组合数取模

组合数取模

求

$$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$不会太大，预处理出然后快速幂就好了。剩下一段也是预处理的。
然后就可以解决了。