组合数取模

组合数取模

求

$$C_m^n\mod mo$$
mo=p1^k1p2^k2….
n,m<=10^9,pi^ki<=10001

首先，我们把$mo$分解质因数，写成$p_1^{k_1}p_2^{k_2}...p_r^{k_r}$的形式。
然后，我们只需要算出$C_m^n\mod p_i^{k_i}$的值，然后用中国剩余定理组合起来就行了。

中国剩余定理

同余方程组$X≡a_i\mod (m_i)$在模$M=\prod_{i=1}^{n}m_i$的意义下有唯一解。
设$M_i=M/m_i,t_i=M_i^{-1}\mod m_i$
则

$$X=\sum_{i=1}^{n}a_iM_it_i\mod M$$
（不要问为什么）

继续

$$C_m^n={m!\over{n!(m-n)!}}$$
如果我们能把$x!$表示成$a*p^b$次方的形式的话，那么就可以直接用分式上下相除而得到，因为a和p互质。
那么我们看一看，把1~x中与p的倍数提出来，这样的数有$\lfloor {x\over p}\rfloor$个，也就得到了$p^{\lfloor {x\over p}\rfloor}$，剩下的数就变成了1~$\lfloor {x\over p}\rfloor$，递归计算就可以了。
剩下的东西，我们把它分组为1~$p^k-1$，$p^k+1$~$2*p^k-1$….可以证明每一组中的数的乘积都相等（在模$p^k$的意义下），由于$p^k$不会太大，预处理出然后快速幂就好了。剩下一段也是预处理的。
然后就可以解决了。