poj1845求n的因子和

此题求a^b的因子和，并且是大数。

一个数n可以写成素数的乘积。

 对A进行素因子分解得

     A = p1^a1 * p2^a2 * p3^a3 *...* pn^an.
  故 A^B = p1^(a1*B) * p2^(a2*B) *...* pn^(an*B);

例如：（15）^3=（3*5）^3=（3^0+3^1+3^2+3^3）*(5^0+5^1+5^2+5^3);

#include<stdio.h>
#include<string.h>
#define maxn 9901
#define tt 10002
int pri[tt],cnt[tt];
__int64 mod(int n,int m) //a^b%maxn.
{
    __int64 sum=1;
    while(m>0)
    {
        if(m&1) sum=(sum%maxn*n%maxn)%maxn;
        m>>=1;
        n=((n%maxn)*(n%maxn))%maxn;
    }
    return sum;
}
__int64 sum(int n,int m) //求和。二分法.
{
    if(m==0) return 1;
    if(m&1) return (((1+mod(n,m/2+1))%maxn)*sum(n,m/2)%maxn)%maxn; //没进行多次求余，wa了好多次。下面没处都加上。
    else return ((1+mod(n,m/2+1))%maxn*sum(n,m/2-1)%maxn+mod(n,m/2)%maxn)%maxn;
}
int main()
{
    int x,y,k,i;
    while(~scanf("%d%d",&x,&y))
    {
        if(x==0){ printf("0\n");continue ;} //加上这两个if，时间即为0ms，
        if(y==0) { printf("1\n");continue;}
        for(k=0,i=2;i*i<=x;i++)  //省去打素数表的部分，认为打表更浪费时间。
          if(x%i==0)
          {
              pri[k]=i;
              while(x%i==0)
              {
                  cnt[k]++;
                  x/=i;
              }
              k++;
          }
          if(x!=1)
          {
              pri[k]=x;
              cnt[k++]=1;
          }
          __int64 ans=1;
          for(i=0;i<k;i++)
            ans=ans*(sum(pri[i],cnt[i]*y)%maxn)%maxn;
        printf("%I64d\n",ans);
    }
    return 0;
}