当我们拆分完数据以后，
A^B的所有约数之和为：
sum = [1+p1+p1^2+...+p1^(a1*B)] * [1+p2+p2^2+...+p2^(a2*B)] *...*[1+pn+pn^2+...+pn^(an*B)].
当时面对等比数列的时候，想到了求和公式，因为直接算超时了，但是带膜除法不能直接除，所以又想到了乘法逆元，但是逆元的使用条件是除数和mod互质的时候，题目给我们的膜不够大，然后我就方了，不知道该怎么去处理了，后来看到网上，才学会了等比快速求和的方法。
它的思想是二分法+递归，规律如下：
（1）若n为奇数,一共有偶数项，则：
1 + p + p^2 + p^3 +...+ p^n
= (1+p^(n/2+1)) + p * (1+p^(n/2+1)) +...+ p^(n/2) * (1+p^(n/2+1))
= (1 + p + p^2 +...+ p^(n/2)) * (1 + p^(n/2+1))
（2）若n为偶数,一共有奇数项,则:
1 + p + p^2 + p^3 +...+ p^n
= (1+p^(n/2+1)) + p * (1+p^(n/2+1)) +...+ p^(n/2-1) * (1+p^(n/2+1)) + p^(n/2)
= (1 + p + p^2 +...+ p^(n/2-1)) * (1+p^(n/2+1)) + p^(n/2);
至于幂的求法，可以用快速幂去求。代码如下：
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
///sqrt(50000000) = 7071.07;///数据足够
/// num = a1(q^n - 1)/ (q-1);///方法难以使用
const long long Mod = ;
#define maxn 8000
#define LL long long
LL a,b,p[maxn],e[maxn],tot;
void split()
{
int d = sqrt(a*1.0);///素数因子在它的根号之下
tot = ;
memset(e,,sizeof(e));
for(int i = ; i <= d; i+=)
{
if(a == ) break;
if(a%i == )
{
p[tot] = i;
while(a % i == )
{
a /= i;
e[tot]++;
}
tot++;
}
if(i == ) i--;///这种方法求素数很高效
}
if(a != )
{
p[tot] = a;
e[tot]++;
tot++;
}
for(int i = ; i < tot; i++)
e[i] *= b;
/*for(int i = 0;i < tot;i++){
printf("p[%d] = %lld e[%d] = %lld\n",i,p[i],i,e[i]);
}*/
}
LL mypow(LL a,LL b)
{
if(b == ) return ;
if(b == ) return a % Mod;
if(b % == ) return mypow(((a%Mod)*(a%Mod))%Mod,b/)%Mod;
else return ((a%Mod) * mypow(a%Mod,b-)) % Mod;
}
LL get_sum(LL a,LL b)
{
if(b==) return ;
if(b % ) return ((get_sum(a,b/)%Mod)*(+mypow(a,b/+))%Mod) % Mod;
else return ((get_sum(a,b/-)%Mod * (+mypow(a,b/+)%Mod))%Mod + mypow(a,b/)%Mod) % Mod; }
int main()
{
while(~scanf("%I64d %I64d",&a,&b))
{
split();
LL ans = ;
for(int i = ;i < tot;i++)
{
ans = ans * get_sum(p[i],e[i])%Mod;///这里不可以省略
}
printf("%I64d\n",ans);
}
return ;
}
注意：这里有一个很难发现的错误，在取膜的时候不可以使用 ans ×= 的形式，优先级的不同会让他溢出。