题目背景
感谢@hehe_54321 提供的3组hack数据
题目描述
输入两个正整数aa和bb，求a^ba
b
的因子和。由于结果太大，只要输出它对9901的余数。
输入格式
仅一行，为两个正整数aa和bb (0leq a,b leq 500000000≤a,b≤50000000)。
质因数分解
(1+p1^1+p1^2+...+p1^c1)(1+p2^1+p2^2+...+p2^c2)...*(1+pn^1+pn^2+...+pn^cn)
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int N=3e6+10,M=5e7+1,mod=9901;
#define int long long
int a,b,tmp,ans=1;
inline int ksm(int x,int y){
int ans=1;
while(y){
if(y&1)ans=(ans*x)%mod;
x=(x*x)%mod;
y=y>>1;
}
return ans;
}
inline int sum(int x,int y){
if(x==0)return 1;
if(x%2==1)return sum(x/2,y)*(1+ksm(y,x/2+1))%mod;
return (sum(x/2-1,y)*(1+ksm(y,x/2+1))%mod + ksm(y,x/2))%mod;
}
signed main(){
cin>>a>>b;
int lim=sqrt(a);
for(int i=2;i<=lim;i++)
if(!(a%i)){
tmp=0;
while(!(a%i)){a/=i;tmp++;};
ans=ans*sum(tmp*b,i)%mod;
}
if(a!=1)ans=ans*sum(b,a)%mod;
cout<<ans<<endl;
}