題意：求a^b的所有約數的和對9901取余。

分析：我們轉化為a^b的所有質因子的等比數列的成積，例如100^1，轉化為(1+2+4) * (1 + 5 ＋ 25)。由于a^b的質因子與a的質因子相同，只是每個的數量是a的質因子的b倍。具體做法是先求所有素數，求a的所有質因子，對于每個質因子求num[i]*b＋1項的等比數列。并求乘積。

View Code #include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
usingnamespace std;

#define maxn 10000
#define w 9901

int a, b;

boolis[maxn];
int prm[maxn];
int fac[maxn], num[maxn];

int getprm(int n)
{
int i, j, k =0;
int s, e = (int) (sqrt(0.0+ n) +1);
memset(is, 1, sizeof(is));
prm[k++] =2;
is[0] =is[1] =0;
for (i =4; i < n; i +=2)
is[i] =0;
for (i =3; i < e; i +=2)
if (is[i])
{
prm[k++] = i;
for (s = i *2, j = i * i; j < n; j += s)
is[j] =0;
}
for (; i < n; i +=2)
if (is[i])
prm[k++] = i;
return k;
}

int power(int a, int n)
{
int ret =1;
int m = a % w;
while (n)
{
if (1& n)
ret *= m;
ret %= w;
n >>=1;
m *= m;
m %= w;
}
return ret;
}

int cal(int d, int num)
{
d %= w;
if (num ==0)
return1;
if (num ==1)
{
return (d +1) % w;
}
int mid;
if (num &1)
{
mid = power(d, num /2+1);
return (mid +1) * cal(d, num /2) % w;
}
int x = power(d, num /2);
mid = power(d, num /2+1);
return ((mid +1) * cal(d, num /2-1) + x) % w;
}

int main()
{
//freopen("t.txt", "r", stdin);
scanf("%d%d", &a, &b);
int n = getprm(int(sqrt(a)) +1);
int temp =0;
for (int i =0; i < n; i++)
if (a % prm[i] ==0)
{
fac[temp] = prm[i];
while (a % prm[i] ==0)
{
a /= prm[i];
num[temp]++;
}
temp++;
}
n = temp;
int ans =1;
for (int i =0; i < n; i++)
ans = ans * cal(fac[i], num[i] * b) % w;
if (a !=1)
ans = ans * cal(a, b) % w;
printf("%d\n", ans);
return0;
}

總結

以上是生活随笔為你收集整理的poj1845的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。