【POJ1845】Sumdiv【数论】

$link$

分析：

考虑唯一分解定理 $A$ 可以表示为 ${\prod }_{i=1}^n{p}_i^{c_i}$
约数之和即 $(1+p_1+p_1^2+...+p_1^{c_1})\times (1+p_2+p_2^2+...+p_2^{c_2})\times ...\times (1+p_n+p_n^2+...+p_n^{c_n})$

那么 $A^B$可表示为 ${\prod }_{i=1}^n{p}_i^{c_i\times B}$
此时约数之和即 $(1+p_1+p_1^2+...+p_1^{c_1\times B})\times (1+p_2+p_2^2+...+p_2^{c_2\times B})\times ...\times (1+p_n+p_n^2+...+p_n^{c_n\times B})$
每部分可以等比数列求和 即 ${\prod }_{i=1}^n\frac{p_i^{c_i+1}-1}{p_i-1}$ 所以快速幂乘 $p_i-1$ 的逆元即可

注意到 $9901$ 是质数 若 $p_i-1$ 是 $9901$ 的倍数 逆元会不存在 此时 $p_i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}9901)$
所以等比数列之和为 $(1+1+1^2+...+1^{c_n\times B+1})=c_n\times B+1$ 答案即 ${\prod }_{i=1}^n{c}_i\times B+1$

CODE：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define reg register
using namespace std;
typedef long long ll;
const ll Mod=9901,N=1e5+5;
ll A,B,p[N],c[N],tot,ans=1;
void calc(ll x)
{
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{
			p[++tot]=i;
			c[tot]=1;
			x/=i;
			while(x%i==0)
				x/=i,c[tot]++;
		}
	}
	if(x>1) p[++tot]=x,c[tot]=1;
}
ll ksm(ll a,ll k)
{
	ll res=1;
	while(k)
	{
		if(k&1) (res*=a)%=Mod;
		(a*=a)%=Mod;
		k>>=1;
	}
	return res;
}
int main(){
	scanf("%lld%lld",&A,&B);
	calc(A);
	for(int i=1;i<=tot;i++)
	{
		int a=p[i],k=B*c[i];
		((a-1)%Mod==0)?ans=ans%Mod*(k+1)%Mod:ans=ans%Mod*(ksm(a,k+1)%Mod-1+Mod)%Mod*ksm(a-1,Mod-2)%Mod;
	}
	printf("%lld",ans);
	return 0;
}