思路:
1.对整数A进行质因数分解,统计质因数的个数即pn。
2.对于本题,如果直接对上面进行等比数列的计算会WA,将其拆分成两种情况:奇数和偶数考虑。则sum(p,c)可以表示为
对偶数:sum(p,c-1)*p+1
对奇数:sum(p,(c-1)/2) * (p(c+1)/2+1)
代码:
#include<iostream>
using namespace std;
const int mod=9901;
//快速幂
int qmi(int a,int b)
{
int res=1%mod;
while(b)
{
if(b&1) res=res*a%mod;
b>>=1;
a=a*a%mod;
}
return res;
}
//使用分治对数列求和
int sum(int p,int c)
{
if(c==0) return 1;
if(c%2==0)
return (p%mod*sum(p,c-1)+1)%mod;//偶数
return (qmi(p,c/2+1)+1)*sum(p,c/2)%mod;//奇数
}
int main()
{
int A,B;
cin>>A>>B;
int res=1;//结果
//计算分解的质因数个数
for(int i=2;i<=A;i++)
{
int count=0;//计算分解的p^c个数
while(A%i==0)
{
count++;
A/=i;
}
if(count) res=res*sum(i,count*B)%mod;//记得乘B
}
if(A==0) res=0;
cout<<res<<endl;
return 0;
}