Sumdiv POJ - 1845(分解质因数+小学奥数题+快速幂+等比数列求和递归公式推导)

题意：传送门
题解：要算$A^B$的所有约数之和$\%9901$，首先可以对$A$进行质因数分解$$A=p_1^{k_1}\ast p_2^{k_2}\ast \cdots \ast p_n^{k_n}$$约数之和的取法有多少种呢？$$num=(k_1+1)\ast (k_2+1)\ast \cdots \ast (k_n+1)$$也就是$A$的所有约数之和可以表示为$$sum=(1+p_1^1+p_1^2+\cdots +p_1^{k_1})\ast (1+p_2^1+p_2^2+\cdots +p_2^{k_2})\ast \cdots \ast (1+p_n^1+p_n^2+\cdots +p_n^{k_n})$$这其实也就是个小学奥数题，然后有个幂$B$，那么质因数分解完相对应的每个都乘上$B$，变成：$$A^B=p_1^{k_1\ast B}\ast p_2^{k_2\ast B}\ast \cdots \ast p_n^{k_n\ast B}$$然后对应的约数之和取法就有：$$num=(k_1\ast B+1)\ast (k_2\ast B+1)\ast \cdots \ast (k_n\ast B+1)$$对应的$A^B$的所有约数之和就是：$$sum=(1+p_1^1+p_1^2+\cdots +p_1^{k_1\ast B})\ast (1+p_2^1+p_2^2+\cdots +p_2^{k_2\ast B})\ast \cdots \ast (1+p_n^1+p_n^2+\cdots +p_n^{k_n\ast B})$$知道了这个以后还得了解一个等比数列求和的公式推导：$$sum(p,c)=1+p+p^2+\cdots +p^c$$如果c为奇数$$sum(p,c)=(1+p+\cdots +p^{\frac{c-1}{2}})+(p^{\frac{c+1}{2}}+\cdots +p^c)$$$$sum(p,c)=(1+p+\cdots +p^{\frac{c-1}{2}})+p^{\frac{c+1}{2}}\ast (1+p+\cdots +p^{\frac{c-1}{2}})$$$$sum(p,c)=(1+p^{\frac{c+1}{2}})\ast sum(p,\frac{c-1}{2})$$如果c为偶数的话，那么也可以推导公式，可以差分为$$sum(p,c)=p\ast sum(p,c-1)+1$$或者推导为公式：$$sum(p,c)=(1+p^{\frac{c}{2}})\ast sum(p,\frac{c}{2}-1)+p^c$$这个可自行推导。
附上代码：

#include<iostream>
using namespace std;
const int mod=9901;
int a,b;
int qpow(int a,int b)
{
    a%=mod;
    int res=1;
    for(int i=b;i;i>>=1,a=a*a%mod)if(i&1)res=res*a%mod;
    return res;
}
int sum(int p,int k)
{
    if(k==0)return 1;
    else if(k%2==0)return ((1+qpow(p,k/2))%mod*sum(p,k/2-1)%mod+qpow(p,k))%mod;
    else return ((1+qpow(p,(k+1)/2))%mod*sum(p,(k-1)/2)%mod)%mod;
}
int main()
{
    cin>>a>>b;
    int res=1;
    for(int i=2;i<=a;i++){
        int s=0;
        while(a%i==0){
            s++;
            a/=i;
        }
        if(s)res=res*sum(i,s*b)%mod;
    }
    if(!a)res=0;
    cout<<res<<endl;
    return 0;
}