POJ 1845 约数和

这道题 题目描述及其简单

给出 a ,b 求 a^b 所有约数 % 9901 的数是多少。

首先，a 可能不是个质数，可能有 f1 f2 .. fn 个因数

所以要先分解a, 然后 使用约数和公式

公式及其证明 http://blog.csdn.net/u012935901/article/details/20144899

最后 在约束和公式中要用到等比数列求和，因为等比数列和涉及到除法，而本体又有mod运算，所以，不能直接用公式。

而因为等比数列良好的结构，可以分治一下。就是先求左边 n/2 项，然后 乘个 q^(n/2) 次方 就是后面了，当然，奇偶不同，应当处理一下这些细节。

最后，a =0 的情况要特殊处理。

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
#include <map>
#include <set>
#include <cmath>
using namespace std;
#define MAX 1100000 
#define INF 0x3f3f3f3f
#define MS(x) memset(x,0,sizeof(x))
#define ll long long
#define P pair<int,int>
#define fst first
#define sec second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define MOD 9901
vector<int> prim;
bool noprim[2500000];
map<int,int> fac;
void getPrim()
{	
	for(int i=2;i<250000;i++)
	{
		if(!noprim[i])
		{
			prim.push_back(i);
			for(int j=2;i*j<25000;j++)
				noprim[i*j]=1;
		}
	}
}
void getFac(ll n)
{
	if(n==1)
		return ;
	int ma=sqrt(n);
	for(int i=0;prim[i]<=ma;i++)
	{
		if(n%prim[i]==0)
		{
			fac[prim[i]]++;
			getFac(n/prim[i]);
			return ;
		}
	}
	fac[n]++;
}

ll pow(ll a,ll b)
{
	ll res=1;
	while(b>0)
	{
		if(1&b)
			res=(a*res)%MOD;
		a=a*a%MOD;
		b>>=1;
	}
	return res;
}

ll arrSum(ll a,ll n)//0 ~ n
{
	if(n==0)
		return 1;
	ll m=n>>1;
	ll tmp;
	if(n%2)
	{
		tmp=arrSum(a,m);
		tmp=(pow(a,m+1)*tmp+tmp)%MOD;
	}else
	{
		tmp=arrSum(a,m-1);
		tmp=(pow(a,m+1)*tmp+pow(a,m)+tmp)%MOD;
	}
	return tmp;
}

int main()
{
	getPrim();
	ll a,b;
	while(cin>>a>>b)
	{
		if(a==0)
		{
			cout<<0<<endl;
			continue;
		}
		fac.clear();
		getFac(a);
		ll ans=1;
		for(map<int,int>::iterator it=fac.begin();it!=fac.end();it++)
		{
			//cout<<it->fst<<endl;
			ans*=arrSum(it->fst,(it->sec)*b);
			ans%=MOD;
		}
		cout<<ans<<endl;
	}	
	return 0;
}