POJ 1845 sumdiv(分治、二分递归、快速幂、数论约束和公式)

POJ 1845

解题思路如下：

注意点

1. 小心多组输入
2. 不要太轻易相信板子。我的快速幂重载居然错了。无语无语
3. 记得清零

LL p[maxn],k[maxn],modd=9901;
LL sum(LL p,LL k){
    if(k==0)return 1;
    if(k%2){
        return (sum(p,k/2)*(1+qpow(p,k/2+1)))%modd;
    }else return (sum(p,k/2-1)*(1+qpow(p,k/2+1))+qpow(p,k/2))%modd;
}
LL solve(LL a,LL b){
    LL co=0,ans=1;
    for(LL i=2;i*i<=a;i++){
        if(a%i==0){
            p[++co]=i;
            k[co]=0;
            while(a%i==0){
                a/=i;
                k[co]++;
            }
        }
    }
    if(a!=1){
        p[++co]=a;
        k[co]=1;
    }
    for(int i=1;i<=co;i++){
        ans=(ans*sum(p[i],b*k[i])%modd)%modd;
    }
    return ans;
}
int main(){
    LL a,b,co=0,tt;
    while(scanf("%lld%lld",&a,&b)!=EOF){
        printf("%lld\n",solve(a,b));
    }
    return 0;
}