bzoj 1478: Sgu282 Isomorphism && 1815: [Shoi2006]color 有色图【dfs+polya定理】

参考 https://wenku.baidu.com/view/fee9e9b9bceb19e8b8f6ba7a.html?from=search### 的最后一道例题
首先无向完全图是个若干点的置换，但是实际上要染色边，也就是要求边的置换
首先，通过dfs构造一个点的置换，然后再把每个置换分割加起来就是答案（实际上分割方案很少）
那么现在有一个点置换的长度(a1,a2,a3...)，考虑边置换，一条边(pi,pj)，如果pi,pj在不同的置换里，那么显然循环节是lcm(ai,aj)，所以循环个数就是gcd(ao,aj)；
对于pi,pj在同一个置换里，如果a是奇数，那么只有循环节长度为2的循环个数就是(a-1)/2，如果是偶数，除了长度为2的循环节还有长度为a/2的，所以个数是a/2
然后一个拆分的方案数是https://blog.csdn.net/litble/article/details/79116659

#include<iostream>
#include<cstdio>
using namespace std;
const int N=60;
int n,m;
long long fac[N],mod,ans,a[N];
long long gcd(long long a,long long b)
{
	return !b?a:gcd(b,a%b);
}
long long ksm(long long a,long long b)
{
	long long r=1;
	while(b)
	{
		if(b&1)
			r=r*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return r;
}
void dfs(int w,int s,int y)
{
	if(!y)
	{
		int c=0,tot=1;
		long long nw=1;
		for(int i=1;i<w;i++)
			c+=a[i]/2;
		for(int i=1;i<w;i++)
			for(int j=i+1;j<w;j++)
				c+=gcd(a[i],a[j]);
		for(int i=1;i<w;i++)
			nw=nw*a[i]%mod;
		for(int i=2;i<w;i++)
		{
			if(a[i]!=a[i-1])
				nw=nw*fac[tot]%mod,tot=0;
			tot++;
		}
		nw=fac[n]*ksm(nw*fac[tot]%mod,mod-2)%mod;
		ans=(ans+nw*ksm(m,c))%mod;
	}
	if(y<s)
		return;
	for(int i=s;i<=y;i++)
	{
		a[w]=i;
		dfs(w+1,i,y-i);
	}
}
int main()
{
	scanf("%d%d%lld",&n,&m,&mod);
	fac[0]=1;
	for(int i=1;i<=n;i++)
		fac[i]=fac[i-1]*i%mod;
	dfs(1,1,n);
	printf("%lld\n",ans*ksm(fac[n],mod-2)%mod);
	return 0;
}