Sumdiv poj1845

题目链接
解：
1、首先考虑质因数分解
2、每个质因数可以取从0到该质因数对应的幂次
3、所以答案就是若干个等比数列求和的结果。
4、因为考虑mod，所以除以一个数等于乘以这个数的逆元，9901是质数，费马小定理求逆元
5、就是关于mod运算的性质，要注意该模的时候模，不该模的时候不要乱模
一篇关于模运算的文章
代码如下：

#include<iostream>
#include<cstdio>
using namespace std;

typedef long long ll;
#define mod 9901

ll a,b,ans;ll pri[50000010],aa[50000010],tot[50000010];

ll quick_pow(ll ml,ll nl){
	ll s=1;
	while(nl>0){
		if(nl%2==1){
			s=(s%mod)*(ml%mod)%mod;
		}
		ml=ml*ml%mod;
		nl=nl>>1;
	}
	return s%mod;
}

int Decomposition(ll x){
	ll cnt=0,Cnt=0;
	for(int i=2;i*i<=x;++i){
		for(;x%i==0;x/=i){
			aa[++cnt]=i;
		}
	}
	if(x>1)aa[++cnt]=x;
	for(int i=1;i<=cnt;++i,tot[Cnt]++){
		if(aa[i]!=aa[i-1])pri[++Cnt]=aa[i];
	}
	return Cnt;
}

int main(){
	cin>>a>>b;ans=1;
	if(a==0){
	    printf("0\n");
	    return 0;  
	}
	ll mm=Decomposition(a);
	for(int i=1;i<=mm;++i){
		if((pri[i]-1)%mod==0){
			ans=ans*(1+tot[i]*b)%mod;
		}
		else 
		ans=ans%mod*quick_pow((pri[i]-1)%mod,mod-2)*(quick_pow(pri[i]%mod,(tot[i]*b+1))-1)%mod;
	}
	printf("%lld\n",((ans%mod)+mod)%mod);
	return 0;
}