标签：poj font 公式 sig ret ems a* 做了 target
没想到贱贱的数据居然是错的。。搞得我调了一中午+晚上一小时（哦不d飞LJH掉RP毕竟他是BUFF）结果重判就对了五次。。
回归正题，这题傻子都看得出是polya定理（如果你不是傻子就看这里），还没有翻转，就一个旋转，结果我就欢快的打完交上去了。傻子都知道会TLE，n<=1e9啊，O(n)都原地爆炸，那怎么办。。。一脸懵逼（然后就膜题解了）
可以发现，这题公式就是sigma(gcd(k,n))(k=1~n)，然后该怎么优化呢，我(??)发现gcd(k,n)里面肯定有一些k和n的gcd是相同的，那我们设n=i*gcd，k=j*gcd，那i肯定和j互质并且1<=j<=i，而且可以发现，gcd(i*gcd,j*gcd)=gcd,只要知道j有多少个，就让power(n,n/i)乘上这个个数，那gcd=n/i的所有情况就都解决了，那具体j有多少个呢？显而易见(??)就是欧拉函数值(然而我不会)了，那我们O(sqrt(n))枚举i，然后就可以得出gcd，然后就可以求出欧拉函数值，那就是phi(i)*power(n,n/i)
#include<cstdio> #include<cstring> #include<cmath> using namespace std; bool v[1100000]; int pr,prime[1100000]; void linear_prime() { memset(v,true,sizeof(v)); for(int i=2;i<=1000000;i++) { if(v[i]==true)prime[++pr]=i; for(int j=1;j<=pr&&i*prime[j]<=1000000;j++) { v[i*prime[j]]=false; if(i%prime[j]==0)break; } } } int n,mod; int power(int A,int k) { int ans=1;A%=mod; while(k!=0) { if(k%2==1)ans=(ans*A)%mod; A=(A*A)%mod;k/=2; } return ans; } int phi(int x)//求欧拉函数值，即j的种数 { int ans=x; for(int i=1;prime[i]*prime[i]<=x;i++) { if(x%prime[i]==0) { ans=ans-ans/prime[i]; while(x%prime[i]==0)x/=prime[i]; } } if(x!=1)ans=ans-ans/x; return ans%mod; } int main() { linear_prime(); int T; scanf("%d",&T); while(T--) { scanf("%d%d",&n,&mod); int ans=0; //设n=i*gcd 而 k=j*gcd，用欧拉函数将所有j的值求出 for(int i=1;i*i<=n;i++) { if(n%i==0) { ans=(ans+phi(i)*power(n,n/i-1))%mod;//循环节为gcd if(i*i!=n)ans=(ans+phi(n/i)*power(n,i-1))%mod; //这里两个power为什么要-1？由于要%mod，所以求值的时提早将/G（G=n）给做了 } } printf("%d\n",ans); } return 0; }
标签：poj font 公式 sig ret ems a* 做了 target
原文地址：http://www.cnblogs.com/AKCqhzdy/p/7598571.html
