这道题 题目描述及其简单
给出 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;
}