POJ1845 Sumdiv

逆元

题目传送门

题目大意：给你两个正整数A、B，求$A^B$的所有因子之和$mod\ 9901$。

把A分解成质因数形式

$$A^B=p_1^{a_1B}p_2^{a_2B}\dots p_n^{a_nB}$$
一个数的质因子之和
$$\phi(n)=\frac{p_1^{a_1+1}}{p_1-1}\times\frac{p_2^{a_2+1}}{p_2-1}\times\dots\times\frac{p_n^{a_n+1}}{p_n-1}$$
于是
$$A^B=\frac{p_1^{a_1B+1}}{p_1-1}\times\frac{p_2^{a_2B+1}}{p_2-1}\times\dots\times\frac{p_n^{a_nB+1}}{p_n-1}$$
然后因为每个 $p$和$a$的指数都很大 废话，所以我们需要用到快速幂，而快速幂在模的时候需要用到一般求逆元的方法。
即 $ans=\frac{a}{b}mod\ m=a\ mod\ (mb)/b$

然后就差不多了

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 10000
#define WJD 9901//不要在意名字
using namespace std;
typedef long long LL;
LL a,b;
LL prime[MAXN+5],p,t,num[MAXN+5][2];
bool f[MAXN+5];
void Make(){//筛素数
    memset(f,true,sizeof(f));
    for (LL i=2;i<=MAXN;i++){
        if (f[i]) prime[++p]=i;
        for (LL j=1;j<=p&&prime[j]*i<=MAXN;j++){
            f[i*prime[j]]=false;
            if (i%prime[j]==0) break;
        }
    }
}
void fenjie(LL x){//质因数分解
    LL node=1; memset(num,0,sizeof(num));
    while (x>1&&node<=p){
        if (x%prime[node]==0) num[++t][0]=prime[node];
        while (x%prime[node]==0){
            x/=prime[node]; num[t][1]++;
        }
        node++;
    }
    if (x>1) {
        num[++t][0]=x; num[t][1]=1;
    }
}
LL mul(LL a,LL b,LL m){//快速乘
    LL ret=0; a%=m;
    while (b){
        if (b&1)
            ret=(ret+a)%m;
        a=(a+a)%m;
        b/=2;
    }
    return ret;
}
LL qsm(LL a,LL b,LL m){//快速幂
    LL ret=1; a%=m;
    while (b){
        if (b&1)
            ret=mul(ret,a,m);
        a=mul(a,a,m);
        b/=2;
    }
    return ret;
}
int main(){
    Make(); 
    while (scanf("%lld%lld",&a,&b)==2){
        fenjie(a); LL ans=1;
        for (LL i=1;i<=t;i++){
            LL M=(LL)(num[i][0]-1)*WJD;
            ans*=(qsm(num[i][0],num[i][1]*b+1,M)+M-1)/(num[i][0]-1)%WJD;//求逆元把除换成乘
            ans=(ans+WJD)%WJD;
        }
        printf("%lld\n",ans);
    }
}