题意就是输入 a,b 输出 a的b次方的因子求和并对9902 取模

这题可以对因子化简

由于唯一分解定理 可以把a表示成 p1^q1*p2^q2...*pn^qn

也就是说 a^b = p1^(q1*b)*p2^(q2*b)...*pn^(qn*b)

由于约数和公式可得a^b得约数求和为(1+p1+p1^2+p1^3+...p1^(q1*b))(1+p2+p2^2+p2^3+...p2^(q2*b))...(1+pn+pn^2+pn^3+...pn^(qn*b))

故问题转化成了如何求等比数列得问题

我们再看 1+p1+p1^2+p1^3+...p1^(q1*b) 可以对其因式分解 从而简化计算

若q*b是偶数 那么一共为奇数项 最终可以化简为 （1+p^(n/2+1) )*（1+p+p^2+...p^(n/2-1)）+p^(n/2)

若q*b是奇数 那么一共为偶数项 化简为--> （1+p+p^2+...p^(n/2)*(1+p^(n/2+1))

其中都可以利用一个二分的思想 就是把数列拆开来算

code：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;const int maxn = 50000000;
const int mod = 9901;
vector<int>p;
int size,a,b;
bool vis[20005];
void prime()
{for(int i=2;i*i<=maxn;i++){if(!vis[i]){p.push_back(i);for(int j=i+i;j*j<=maxn;j+=i)vis[j]=1;}}
}
ll fun(ll aa,ll bb)
{ll ans = 1;aa%=mod;while(bb){if(bb&1)ans = ans*aa %mod;bb>>=1;aa = aa*aa%mod;}return ans%=mod;
}
ll bi(ll x,ll n)
{if(n==0)return 1;else{if(n%2)return (bi(x,n/2)*( 1 + fun(x,n/2+1 ) ))%mod;else return ( bi(x,n/2-1)*( 1 + fun(x,n/2+1 ) )+fun(x,n/2))%mod;}
}
ll solve()
{ll ans=1; ll t = a;for(int i=0;p[i]*p[i]<=a;i++){ll tem =0; //cout<<p[i]<<endl;//system("pause");while(t%p[i]==0)tem++,t/=p[i];tem*=b;ans*=bi(p[i],tem),ans%=mod;	}if(t!=1)ans*= bi(t,b)%mod,ans%=mod;return ans%mod;
}int main()
{ios::sync_with_stdio(0);prime();cin>>a>>b;if(b==0||a<=1)cout<<1<<endl;else cout<<solve()<<endl;return 0;
}

 

#include<iostream>
#include<cstdio>
using namespace std;
const int maxn = 50000000;
const int lim = 10005;
const int mod = 9901;
typedef long long ll;
ll a,b,p[lim],sze;
bool vis[50005];
void prime()
{for(int i=2;i<50005;i++){if(!vis[i]){p[sze++] = i;for(int j=i+i;j<50005;j+=i)vis[j]=1;}}
}
ll fun(ll aa,ll bb,ll m)
{ll ans = 1;aa%=m;while(bb){if(bb&1)ans = ans*aa%m;bb>>=1;aa = aa*aa%m;}return ans%m;
}
ll solve()
{ll ans=1;ll c =a;for(int i=0;p[i]*p[i]<=c;i++){ll tem =0;while(c%p[i]==0)tem++,c/=p[i];if(tem){
//			cout<<p[i]<<" "<<tem<<" "<<c<<endl;ll m = (mod*(p[i]-1));ans = ans*(((fun(p[i],(tem)*b+1,m)-1+m))/(p[i]-1)),ans%=mod;}}
//	cout<<c<<endl;if(c>1){ll m = (mod*(c-1));//cout<<c<<endl;ans = ans*(((fun(c,b+1,m)-1+m))/(c-1)),ans%=mod;}return ans;
}
int main()
{prime();cin>>a>>b;cout<<solve()<<endl;		return 0;
}