Hdu-2421 Deciphering Password

题目地址：

http://acm.hdu.edu.cn/showproblem.php?pid=2421

大意：

给你2个数A,B，范围为1到100万。N是A的B次方。设a1，a2,……ai为N的因子，设k1,k2……,ki为a1，a2,……ai的因子个数。求k1,k2……,ki的立方和。

思路：

算数基本定理推论：

任一个大于1的整数a都能够唯一的写成

a=p1^a1*p2^a2……pk^ak>0.

其中p1<p2<……<pk.且pk为素数。

代码如下：

#include<iostream> #include<cstdio> #include<cstring> #include<string> #include<cmath> #include<algorithm> using namespace std; #define MAXN 1010 bool is_prime[MAXN]; int prime[170]; void pre() //素数打表 { int i, j, k = 0; memset(is_prime, 1, sizeof(is_prime)); for(i = 2; i < 35; ++i) { if(is_prime[i]) for(j = i; j * i < MAXN; ++j) is_prime[i * j] = 0; } for(i = 2; i < MAXN; ++i) if(is_prime[i]) prime[k++] = i; } int main() { int m, n, i,t = 1,count; int sum; pre(); while(scanf("%d%d", &m, &n) != EOF) { sum = 1; for(i = 0; i < 169; ++i) //算数基本定理推论 { count = 0; while(m % prime[i] == 0 && m != 1) { m /= prime[i]; count++; } count =( count * n + 1) % 10007; count=count * (count + 1) / 2 % 10007; sum = sum * (count * count % 10007) % 10007; if(m == 1) break; } if(m != 1) { count=(n + 1) % 10007; count=count * (count + 1) / 2 % 10007; sum = sum * (count * count % 10007) % 10007; } printf("Case %d: %d\n", t++, sum); } return 0; }