tc-srm704-div1-500 解题报告

题意：
给定n,K。有q组询问，第i组询问是求$\prod_{j=1}^nx_j(0\le x_ j<K)=v_i$的解的数量。
$1\le n \le 50,1\le K \le 10^9,1\le q \le 10000$

我本来是这样想的：
显然最终答案只与gcd(v,K)有关。
那么令f(i,j)为前i个数乘积是与n的gcd是j的某个数的方案数，那么$f(i,j)=\sum_{d|j}f(i-1,d)d\phi(d)$。
首先$f(i-1,d)$其乘积可能的取值有$\phi(d)$个，设其为x，再乘上第i个数为y，乘积为z，即$xy\bmod K=z,{x\over d}y ={z\over d}(\bmod {K\over d}),y={z\over d}{x\over d}^{-1}(\bmod {K\over d})$，即对于任意z，y都有${K\over {K\over d}}=d$种取值。
这样时间复杂度就是$O(n(K的约数个数)^2+q\log K)$，虽然说约数个数的估计一般是$O(\sqrt n)$，但这是一个非常松的上界，$10^9$以内约数个数最多有1000左右，所以并不会有问题。

但是我们可以考虑中国剩余定理，若$K=\prod p_1^{a_1}p_2^{a_2}...$，那么模K意义下的答案就是模$p_1^{a_1},p_2^{a_2}...$意义下的答案的乘积，所以我们可以对每个质因子分别dp再乘起来，这样的话$\phi$也变得非常好算了。
这样时间复杂度就是$O(n\log^2K+q\log K)$。

总结：
遇到模一个合数的时候，一定要想一想能不能用中国剩余定理！