poj 2154

polya定理+欧拉函数优化，1641ms水过，具体参考的解题报告:http://www.cnblogs.com/staginner/archive/2012/03/08/2385052.html

 1    #include <iostream>
 2     #include <cstdio>
 3     using namespace std;
 4     typedef long long LL;
 5     int n,p;
 6     int get_Eu(int x)
 7     {
 8         int i;
 9         int ans=x;
10         for(i=2;i*i<=x;i++)
11         {
12             if(x%i==0)
13             {
14                 while(x%i==0) x/=i;
15                 ans=ans/i*(i-1);
16             }
17         }
18         if(x>1) ans=ans/x*(x-1);
19         return ans;
20     }
21     int pow(LL x,LL n)
22     {
23         if(n==0) return 1;
24         if(n==1) return x%p;
25         LL tem=pow(x,n/2);
26         tem=(tem*tem)%p;
27         if(n%2) tem=(tem*(x%p))%p;
28         return (int)tem;
29     }
30     int main()
31     {
32         int cas;
33         cin>>cas;
34         while(cas--)
35         {
36             cin>>n>>p;
37             int i,j,k,tem,sum=0;
38             for(i=1;i*i<n;i++)
39             {
40                 if(n%i==0)
41                 {
42                     tem=get_Eu(n/i);
43                     if(tem) sum=(sum+((tem%p)*pow((LL)n,(LL)(i-1)))%p)%p;
44                     tem=get_Eu(i);
45                     if(tem) sum=(sum+((tem%p)*pow((LL)n,(LL)(n/i-1)))%p)%p;
46                 }
47             }
48             if(i*i==n)
49             {
50                 tem=get_Eu(n/i);
51                 if(tem) sum=(sum+((tem%p)*pow((LL)n,(LL)(i-1)))%p)%p;
52             }
53             printf("%d\n",sum);
54         }
55         return 0;
56     }