poj 2154

来源：互联网发布：淘宝订单信息编辑编辑：程序博客网时间：2024/06/13 11:16

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     }