[POJ2154] Color [Polya定理][莫比乌斯反演]

题意:

从N种颜色的珠子中选出N个串成一圈，问有多少种方案。（不考虑翻转）

读入X(X <= 3500)
对于X组数据，每一组给出N(1 <= N <= 1000000000)和P(1 <= P <= 30000)，要求输出Ans%P

---

环是可以旋转的（本题不考虑翻转）

显然是置换。

有N种置换（包括不转动）；

第i种置换（也就是旋转了i格）有$r=gcd(N,i)$个循环

根据Polya定理有$Ans=\frac{\sum N^r}{N}=\frac{\sum N^{(gcd(N,i))}}{N}=\sum N^{(gcd(N,i)-1)}$

每次Euclid求gcd复杂度$O(\lg i)$

快速幂复杂度$O({\log }_2\sqrt{N})$

总复杂度$O(N\ast (\lg i+{\log }_2\sqrt{N})$

好像十分完美

可是我们想想这道题还有个X

总复杂度应该是$O(X\ast N\ast (\lg i+{\log }_2\sqrt{N})$

凉凉

---

我们前面得到了$Ans=\sum N^{(gcd(N,i)-1)}$

有sum有gcd，思考利用莫比乌斯反演进行优化

如果$F(n)={\sum }_{d|N}f(d)$，有$f(n)={\sum }_{d|N}\mu (d)F(\frac{N}{d})$

由于N为正整数，有${\sum }_{d|N}\frac{\mu (d)}{d}=\frac{\varphi (N)}{N}$

显然 A n s = ∑ ( d ∣ N ) {   N d ∗ ∑ k = 0 N [   d = g c d ( K , N )   ]   } Ans=\sum_{(d|N)}\{\ N^{d}*\sum_{k=0}^N[\ d=gcd(K,N)\ ]\ \} Ans=∑(d∣N)​{ Nd∗∑k=0N​[ d=gcd(K,N) ] }

$F(n)={\sum }_{(d|N)}f(n),f(n)=N^d$

$N^d={\sum }_{(d|N)}\frac{\varphi (N)}{N}\ast {\sum }_{(d|N)}\frac{N^d}{d}$

---

$Ans={\sum }_{(d|N)}{N}^{d-1}\ast {\sum }_{k=0}^{(N-1)}(\frac{gcd(K,N)}{d}=1)$

$={\sum }_{(d|N)}{N}^{d-1}\ast {\sum }_{k=0}^{(N-1)}(gcd(\frac{K}{d},\frac{N}{d})=1)$

$={\sum }_{(d|N)}{N}^{d-1}\ast {\sum }_{k=0}^{(\frac{N}{d}-1)}(gcd(K,\frac{N}{d})=1)$

注意到右边等于$\varphi (\frac{N}{d})$

$Ans={\sum }_{(d|N)}{N}^{d-1}\ast \varphi (\frac{N}{d})$

求欧拉函数因为要mod p就没办法线性筛了有点遗憾（

不过我们可以线性筛素数，用分解质因数的方法枚举素数求出$\Phi$

最后mod不是问题，直接%就好，因为前面已经把除N搞掉了

和p也不一定互质，没办法用什么逆元（

---

感觉做了一道假的polya

Code:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<cctype>
using namespace std;
#define LL long long
const int Limit=sqrt(1e9);
int Primetot=0;
int X,N,P;
int Ans;
LL Phi[32767]={};
int Prime[32767]={};
bool notPrime[32767]={};
int Quick_Pow(int Num,int Pow)
{
	register int Base=Num%P;
	int Res=1;
	while(Pow)
	{
		if(Pow&1)Res=Res*Base%P;
		Base=Base*Base%P;
		Pow>>=1;
	}
	return Res;
}
int phi(int Num)
{
	register int Res=Num;
	for(int i=1;Prime[i]<=sqrt(Num);++i)
	{
		if(Num%Prime[i])continue;
		Res-=Res/Prime[i];
		while(Num%Prime[i]==0)Num/=Prime[i];
	}
	if(Num!=1)Res-=Res/Num;
	return Res%P;
}
int main()
{
	scanf("%d",&X);
	for(int i=2;i<=Limit;++i)
	{
		if(notPrime[i])continue;
		for(int j=2;i*Prime[j]<=Limit&&j<=Primetot;++j)
		{
			notPrime[i*Prime[j]]=1;
			if(i%Prime[j]==0)
			{
				notPrime[i]=1;
				break;
			}
		}
		if(!notPrime[i])Prime[++Primetot]=i;
	}
	while(X--)
	{
		register int i;
		Ans=0;
		scanf("%d%d",&N,&P);
		for(i=1;i*i<=N;++i)
		{
			if(N%i==0)
			{
				Ans+=Quick_Pow(N,i-1)*phi(N/i);
				if(i*i!=N)Ans+=Quick_Pow(N,N/i-1)*phi(i);
				Ans%=P;
			}
		}
		printf("%d\n",Ans);
	}
	return 0;
}

---

上面是一两年前写的
然后呢 里面的莫反就是个套路
$$\sum _{k=0}^{n-1}{n}^{(k,n)-1}$$
$$\sum _{d|n}{n}^{d-1}\sum _{k=0}^{n-1}[d=(k,n)]$$
$$\sum _{d|n}{n}^{d-1}\sum _{k=0}^{\frac{n}{d}-1}[(k,\frac{n}{d})=1]$$
$$\sum _{d|n}\varphi (\frac{n}{d})n^{d-1}$$
吸收一下思想？考虑到暴力计算的重复计算部分然后就大概有莫反的思路了
假如做得熟练，看出莫反就是一眼的事情了。。
然后本题的另一部分（环置换）也很经典啊。所以这道题是一道一眼大水题