poj 1845 Sumdiv

/**
大意: 计算 a^b  的所有因子的和， 和mod 9901
思路； 将a 进行质因子分解，那么所有因子的和为
（2^0+ 2^1 + 2^2 +....+ 2^a1）*(3^0 + 3^1+..+ 3^a2)*.....
注意： 求模n下a的逆，需要 gcd（a，n） =1，，所以此处不适合
用二分求  0-n 次幂的和。。。 又遇到了一次。。。
注意中间过程的取模。，，，不然会溢出。。
**/

#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
using namespace std;
long long ddiv[50000];
long long res[50000];
int cnt;
const int mod = 9901;
void get_dr(long long n){
    memset(res,0,sizeof(res));
    long long m = (long long )sqrt(n+0.5);
    cnt =0;
    for(int i=2;i<=m;i++)if(n%i==0){
        ddiv[cnt] = i;
        while(n%i==0){
            res[cnt]++;
            n = n/i;
        }
        cnt++;
    }
    if(n>1){
        ddiv[cnt] = n;
        res[cnt]++;
        cnt++;
    }
}

long long  pow_mod(long long  a,long long  n){
    if(n==0)
        return 1;
    long long  c =1;
    a =a %mod;
    while(n){
        if(n&1)
            c =c*a%mod;
        a =a*a%mod;
        n>>=1;
    }
    return c;
}

long long sum (long long p,long long n){
        if(n==0)
            return 1; 
        if(n%2){     // 因为是从0开始的。。所以此处才是偶数的情况
            return (sum(p,n/2)*(1+pow_mod(p,n/2+1)))%mod;
        }
        else     // 奇数的情况
            return (sum(p,n/2-1)*(1+pow_mod(p,n/2+1))+pow_mod(p,n/2))%mod;
}

int main()
{
    long long a,b;
    while(cin>>a>>b){
        if(b==0){
            cout<<1<<endl;
        }else if(a==0){
            cout<<0<<endl;
        }else{
            get_dr(a);
            //for(int i=0;i<cnt;i++)
            //    cout<<ddiv[i]<<endl;
            for(int i=0;i<cnt;i++){
                res[i] *= b;
            }
            //for(int i=0;i<cnt;i++)
            //   cout<<res[i]<<endl;
            long long ans =1;
            for(int i=0;i<cnt;i++){
                ans = (ans*(sum(ddiv[i],res[i]))%mod)%mod;
                ans  =ans%mod;
            }
            ans = ans%mod;
            cout<<ans<<endl;
        }
    }
    return 0;
}