poj2154 Color ——Polya定理

题目：http://poj.org/problem?id=2154

今天学了个高端的东西，Polya定理...

此题就是模板，然而还是写了好久好久...

具体看这个博客吧：https://blog.csdn.net/wsniyufang/article/details/6671122

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int const maxn=1e6+5;
int X,n,p,pri[maxn],cnt,ans;
bool vis[maxn];
int phi(int n)
{int ret=n;for(int i=2;i*i<=n;i++){if(n%i)continue;
//        (ret-=ret/i)%=p;// 有除法，不可随时取模 ret-=ret/i;while(n%i==0)n/=i;}
//    if(n!=1)(ret-=ret/n)%=p;// 有除法，不可随时取模 if(n!=1)ret-=ret/n;return ret%p;
} 
int pw(int a,int b)
{int ret=1; a%=p;//!!!for(;b;b>>=1,(a*=a)%=p)if(b&1)(ret*=a)%=p;return ret;
}
int main()
{scanf("%d",&X);while(X--){scanf("%d%d",&n,&p);ans=0;for(int i=1;i*i<=n;i++){if(n%i)continue;(ans+=phi(n/i)*pw(n,i-1))%=p;if(i*i!=n)(ans+=phi(i)*pw(n,n/i-1))%=p;}printf("%d\n",ans);}
}