【数学一本通 第一章】Sumdiv [POJ 1845]

题目：

题目链接：Sumdiv
题解：
题目就是在求$A^B$的所有约数和在%mod。

一：整数的唯一分解定理：

$A^B$=（$p_1^{k_1}\ast p_2^{k_2}\ast p_3^{k_3}$……$p_n^{k_n}$）^B

= ($p_1^{k_{1}B}\ast p_2^{k_{2}B}\ast p_3^{k_{3}B}$……$p_n^{k_{n}B}$)

二：约数和定理：

$A^B$=（1+$p_1^1$+$p_1^2$+……+$p_1^{k_{1}B}$) ×（1+$p_2^1$+$p_1^2$+……+$p_2^{k_{2}B}$）× ……×（1+$p_n^1$+$p_1^n$+……+$p_n^{k_{n}B}$）%mod

三、等比数列二分求和 ：

注意上面的这句话：
（1+$p_1^1$+$p_1^2$+……+$p_1^{k_{1}B}$)
就会发现这就是个等比数列，可是如果用等比数列求和公式就死掉了，不用试就炸掉了。
所以可以试一下二分求解：
1 + q + q² + q³ = (1 + q²)+(q + q³) = (1 + q²)(1 + q)
下面就来推一下整个式子：
假设n是奇数：
1+q+……+qⁿ(偶数个项) = (1 + q^n/2+1)+(q + q^n/2+2)……=(1+q^n/2+1)(1 + q …… + q^n/2)
假设n是偶数：
前面和n为奇数的情况是一样的，最后一项暴力得即可，(1+q^n/2)q^n/2
(这样的话就把复杂度直接降成log的了)

#include<bits/stdc++.h>
#define LL long long
#define mod 9901
using namespace std;
const int sea=10001;
int ans=1,k[sea],p[sea],n,a,b;
LL ksm(LL x,LL y)
{
	LL s=1;
	while(y)
	{
		if(y&1) s=(s*x)%mod;
		y>>=1;
		x=x*x%mod;
	}
	return s;
}	
void ycl()//整数的唯一分解定理
{
	n=0;
	for(int i=2;i*i<=a;)
	{
		if(a%i==0)
		{
			p[++n]=i; k[n]=0;
			while(!(a%i)) k[n]++,a/=i;
		}
		if(i==2) i++;
		else i+=2;
	}
	if(a!=1) p[++n]=a,k[n]=1;
}
LL ysh(LL p,LL k)//等比数列二分求和 
{
	if(k==0) return 1;
	if(k%2) return (ysh(p,k/2)*(1+ksm(p,k/2+1)))%mod;
	else return (ysh(p,k/2-1)*(1+ksm(p,k/2+1))+ksm(p,k/2))%mod;
}

int main()
{
	scanf("%d%d",&a,&b);    ycl();
	for(int i=1;i<=n;i++) ans=(ans*(ysh(p[i],k[i]*b)%mod))%mod;)//约数和定理
	printf("%d\n",ans);
	return 0;
}

continue……