【洛谷】1593 因子和 组合数学

题目传送门

好像这题并不是组合数学啊……只是一个套公式的题目罢了。

一个正整数$n$有素因子分解$=p_1^{a_1}*p_2^{a_2}*…*p_m^{a_m}$，则有以下两个定理：

$$1.n的因子和=\frac{p_1^{a_1+1}-1}{p_1-1}*\frac{p_2^{a_2+1}-1}{p_2-1}*…*\frac{p_m^{a_m+1}-1}{p_m-1}=\prod_{i=1}^m\frac{p_i^{a_i+1}-1}{p_i-1}$$

$$2.n的因子个数=(a_1+1)*(a_2+1)*…*(a_m+1)=\prod_{i=1}^m(a_i+1)$$

然后这题就变成了直接套第一条定理的水题，直接对于每个素因子快速幂求贡献即可。

然而这题给出的模数9901是一个质数，可能被给出的 n <script type="math/tex" id="MathJax-Element-5">n</script>包含，需要特判一下。

p.s.这题在poj上也有，题目要为1845，只是因为poj最近在维护所以不能交就去洛谷了。

附上AC代码：

#include <cstdio>
using namespace std;

typedef long long ll;
const int mod=9901;
ll a,b,c,ans;

inline ll ksm(ll x,ll y){
    ll sum=1;
    while (y){
        if (y&1) sum=sum*x%mod;
        x=x*x%mod,y>>=1;
    }
    return sum;
}

int main(void){
    scanf("%d%d",&a,&b),ans=1;
    for (int i=2; i*i<=a; ++i)
        if (a%i==0){
            c=0;while (a%i==0) a/=i,++c;
            if (i==9901) continue;
            ans=ans*(ksm(i,c*b+1)-1)%mod*ksm(i-1,mod-2)%mod;
        }
    if (a>1&&a!=9901) ans=ans*(ksm(a,b+1)-1)%mod*ksm(a-1,mod-2)%mod;
    printf("%lld",ans);
    return 0;
}