Sumdiv POJ - 1845 A^B的所有约数之和

A^B的所有约数之和为：

     sum = [1+p1+p1^2+...+p1^(a1*B)] * [1+p2+p2^2+...+p2^(a2*B)] *...* [1+pn+pn^2+...+pn^(an*B)].

用递归二分求等比数列1+pi+pi^2+pi^3+...+pi^n：

（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 <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>
#include <cmath>
#include <set>
#include <queue>
using namespace std;

const int INF=1e9+10;
const double EPS = 1e-7;  
typedef long long ll;

ll mypow(ll x,ll k){
	ll sum=1LL;
	while(k){
		if(k&1) sum=1LL*(sum*x)%9901;
		x=1LL*(x*x)%9901;
		k>>=1;
	}
	return sum;
}
ll sum(ll p,ll n)  //递归二分求 (1 + p + p^2 + p^3 +...+ p^n)%mod  
{                          //奇数二分式 (1 + p + p^2 +...+ p^(n/2)) * (1 + p^(n/2+1))  
    if(n==0)               //偶数二分式 (1 + p + p^2 +...+ p^(n/2-1)) * (1+p^(n/2+1)) + p^(n/2)  
        return 1;  
    if(n%2)  //n为奇数,  
        return (1LL*sum(p,n/2)*(1+mypow(p,n/2+1)))%9901;  
    else     //n为偶数  
        return (1LL*sum(p,n/2-1)*(1+mypow(p,n/2+1))+mypow(p,n/2))%9901;  
}
int main(){
	int a,b;
	while(scanf("%d %d",&a,&b)!=EOF){
		ll ans=1LL;
		int m=sqrt(a+0.5);
		for(int i=2;i<=m;){
			if(a%i==0){
				ll cnt=0;
				while(a%i==0) {a/=i;cnt++;}
				ans=(ans*sum(i,1LL*b*cnt)%9901)%9901;
			}
			if(i==2) i++;
			else i+=2;
		}
		if(a>1){
			ans=(ans*sum(a,b)%9901)%9901;
		}
		printf("%lld\n",ans);
	}
	return 0;
}