poj-1845 Sumdiv && nyoj - 928 小M的因子和

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

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

#include<cstdio>
#include<cmath>
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);
    }
}