poj1845 质因数分解+等比数列求和

题意：求A^B约数之和。
分析：首先可以讲A分解质因数：
A=p1^a1+p2^a2+p3^a3+……+pn^an
则A^B=p1^(a1*B)+p2^(a2*B)……+pn^(an*B)
A^B的所有约数之和sum=[1+p1+p1^2+…+p1^(a1*B)][1+p2+p2^2+…+p2^(a2*B)][1+pn+pn^2+…+pn^(an*B)].
之后就是用二分的方法求解了。
下面是代码（其实getprime（）函数是可以不要的）^_^

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long LL;
const int maxn=10010;
const int mod=9901;
int prime[maxn];
void getprime(){
    int cnt=0;
    for(int i=2;i<=maxn;i++){
        if(!prime[i]) prime[++cnt]=i;
        for(int j=i*i;j<=maxn;j+=i){
            prime[j]=1;
            if(i%prime[j]==0) break;
        }
    }
}
LL pow(LL x,LL y){
    LL ret=1;
    while(y){
        if(y&1) ret=(ret*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return ret;
}
LL zhiyinshu[maxn][2];
int cnt;
void fenjie(LL x){
    cnt=0;
    for(int i=1;prime[i]*prime[i]<=x;i++){
        zhiyinshu[cnt][1]=0;
        if(x%prime[i]==0){
            zhiyinshu[cnt][0]=prime[i];
            while(x%prime[i]==0){
                zhiyinshu[cnt][1]++;
                x/=prime[i];
            }
            cnt++;
        }
    }
    if(x!=1){
        zhiyinshu[cnt][0]=x;
        zhiyinshu[cnt++][1]=1;
    }
}
LL sum(LL x,LL y){
    if(y==0) return 1;
    if(x==0) return 0;
    if(y&1){
        return (((1+pow(x,y/2+1))%mod)*sum(x,y/2)%mod)%mod;
    }else return ((1+pow(x,(y/2+1)))%mod*sum(x,y/2-1)%mod+pow(x,y/2)%mod)%mod;
}
int main(){
    int A,B;
    int i;
    scanf("%d%d",&A,&B);
    getprime();
    fenjie(A);
    LL ans=1;
    for(i=0;i<cnt;i++){
        ans=ans*(sum(zhiyinshu[i][0],zhiyinshu[i][1]*B)%mod)%mod;
    }
    printf("%lld",ans);
    return 0;
}