poj2154Color polya定理+欧拉函数优化

没想到贱贱的数据居然是错的。。搞得我调了一中午+晚上一小时（哦不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;//循环节为gcdif(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;
}

转载于:https://www.cnblogs.com/AKCqhzdy/p/7598571.html

poj2154Color polya定理+欧拉函数优化相关推荐

poj2154(Polya+欧拉函数优化模版）
#include <cstdio> #include <cstring> #include<iostream> using namespace std; const ...
poj2154-color-polyan次二面体+欧拉函数优化
N<=1e9,O(nlogn)的做法会超时.从枚举置换转变为枚举轮换长度,然后可以利用欧拉函数,把复杂度变为O(√n * logn) 1 /*-------------------------- ...
【数学知识】三种方法求 [1,n] 中所有数欧拉函数（线性筛欧拉函数优化至 O(n) ）
整理的算法模板合集: ACM模板 ①直接求小于或等于n,且与n互质的数个数(求[1,n]中所有数的欧拉函数时间复杂度:O(nn)O(n\sqrt{n})O(nn)) ②求[1,n]之间每个数的质因数 ...
POJ2154(Pólya定理与欧拉函数优化)
题目:Color 题意:将正n边形的n个顶点用n种颜色染色,问有多少种方案(答案mod p,且可由旋转互相得到的算一种) 先说说Pólya定理设Q是n个对象的一个置换群,用m种颜色涂染这n个对象,一 ...
BZOJ-1951-古代猪文-SDOI2010-费马小定理+欧拉函数+lucas定理+中国剩余定理
描述 =>G∑(ni),i|nmodP =>G^{\sum {{n\choose i}\text{,i|n}}} modP 分析 k=∑Cin,i|n(modP) k=\sum {C_n^ ...
nyoj-Color the necklace（Ploya定理 + 欧拉函数 + 扩展欧几里得（求逆元））
题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=688 此题题解不太懂,因为对这些概念,定理太模糊,理解起来比较困难,不过想想还是应该把代 ...
数论一之定理证明——裴蜀/威尔逊/费马/扩展欧几里得/[扩展]欧拉/[扩展]中国剩余定理，欧拉函数，逆元，剩余系，筛法
打死没想到会在H老师处学懂数论同余,整除模运算埃式筛法欧拉筛法最大公约数和最小公倍数辗转相除法更相减损术裴蜀定理威尔逊定理费马定理同余等价类.剩余系.缩系欧拉函数欧拉定理扩 ...
初等数论--同余--欧拉函数、欧拉定理、费马小定理
初等数论--同余--欧拉函数.欧拉定理.费马小定理概念同余类,既约同余类欧拉函数完全剩余系,既约剩余系关于完全剩余系.既约剩余系一些比较简单的定理欧拉定理.费马小定理博主是初学初等数论( ...
素数的有关性质（二）欧拉函数的一些定理证明与计算
文章目录写在前面内容回顾模mmm剩余类环定理模ppp剩余类域定义欧拉函数的定义欧拉函数的性质命题1:欧拉函数等于与mmm互素整数个数命题2:取值为素数ppp的欧拉函数等于p−1p- ...

poj2154Color polya定理+欧拉函数优化

poj2154Color polya定理+欧拉函数优化相关推荐

最新文章

热门文章