本文共 954 字，大约阅读时间需要 3 分钟。

题意：求解A^B的因子和 mod 9901

先求解素因子，然后二分求解等比数列

#include
   
    #include
    
     typedef long long LL;const LL mod = 9901 ;LL pow(LL a,LL b){    LL res=1;    while(b)    {        if(b&1) res=(res*a)%mod;        a=(a*a)%mod;        b>>=1;    }    return res;}LL fun(LL a,LL b){    LL s,t;    if(b==0) return 0;    if(b==1) return a%mod;    s=fun(a,b/2)%mod;    if(b&1)    {        t=pow(a,b/2+1)%mod;        return (s*(t+1)+t)%mod;    }    else    {        t=pow(a,b/2)%mod;        return (s*(t+1))%mod;    }}int main(){    LL a,b,c,i,k;    while(~scanf("%lld%lld",&a,&b))    {        if(a == 0) {printf("0\n");continue;}        LL sum=1;        for(i=2;i<=sqrt(a);i++)        {            k=0;            if(a%i==0)            {                while(a%i==0)                {                    a/=i;k++;                }                sum=sum*(fun(i,k*b)+1)%mod;            }        }        if(a>1) sum=sum*(fun(a,b)+1)%mod;        printf("%lld\n",sum);    }}
    
   

转载地址：http://pfsgi.baihongyu.com/