BZOJ 1478 Sgu282 Isomorphism

Description

给 定一个N 个结点的无向完全图（ 任意两个结点之间有一条边），现在你可以用 M 种颜色对这个图的每条边进行染色，每条边必须染一种颜色。若两个已染色的图，其中一个图可以通过结点重新编号而与另一个图完全相同，就称这两个染色方案相同。现在问你有多少种本质不同的染色方法，输出结果 mod P。P 是一个大于N 的质数。

Input

仅一行包含三个数，N、M、P。

Output

仅一行，为染色方法数 mod P 的结果。

Sample Input

3 4 97

Sample Output

20

HINT

数据范围：1≤N≤53，1≤M≤1000，N<p≤10 。

Source

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

polya~

具体思路见https://wenku.baidu.com/view/6659d7f59e31433239689329.html~

#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long

int t,n,m,modd,a[60],num[60];
ll ans,fac[60];

int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f; 
}

ll gcd(ll u,ll v)
{
	return v ? gcd(v,u%v):u;
}

ll mi(ll u,ll v)
{
	ll ans=1;u%=modd;
	while(v)
	{
		if(v&1) ans=(ans*u)%modd;
		u=(u*u)%modd;v>>=1;
	}
	return ans;
}

void cal(int u)
{
	ll now=1;
	for(int i=1;i<=n;i++) num[i]=0;
	for(int i=1;i<=u;i++) num[a[i]]++;
	for(int i=1;i<=u;i++) now=now*a[i]%modd;
	for(int i=1;i<=n;i++) now=now*fac[num[i]]%modd;
	int c=0;
	for(int i=1;i<=u;i++)
	{
		c+=a[i]/2;
		for(int j=i+1;j<=u;j++) c+=gcd(a[i],a[j]);
	}
	ans=(ans+(mi(m,c)*fac[n]%modd)*mi(now,modd-2)%modd)%modd;
}

void dfs(int u,int v,int k)
{
	if(u==n+1) {cal(v-1);return;}
	for(int i=1;i<=k && u+i<=n+1;i++)
	{
		a[v]=i;dfs(u+i,v+1,i);
	}
}

int main()
{
	n=read();m=read();modd=read();fac[0]=1;
	for(int i=1;i<60;i++) fac[i]=fac[i-1]*i%modd;
	dfs(1,1,n);
	printf("%lld\n",ans*mi(fac[n],modd-2)%modd);
	return 0;
}