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;
}