BZOJ 1815: [Shoi2006]color 有色图

这题时限不知道怎么搞的，我的代码跑了4000+MS A掉了，可是时限是4S

然后做法什么的不是很懂，好像很玄学的样子，具体看08年的论文。

大概是利用划分对置换进行分类，每一类的循环节个数是一样的，这样统计出每一类的置换个数，然后就可做了。

代码好短……

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll p;
ll qmul(ll a,ll b){
	ll ans=1;
	for(;b;b>>=1,a=a*a%p)if(b&1)ans=ans*a%p;
	return ans;
}
ll inv(ll x){
	return qmul(x,p-2);
}
ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
ll ans,l[60],tot,m,n;
ll f(){
	ll ans=0;
	for(int i=1;i<=tot;i++){
		ans+=l[i]/2;
		for(int j=i+1;j<=tot;j++)
		ans+=gcd(l[i],l[j]);
	}
	return ans%(p-1);
}
void dfs(int last,int k,int n,ll s){
	l[++tot]=n;
	int t=s*inv(n)%p;
	if(n==last)t=t*inv(k+1)%p;
	else t=t*inv(1)%p;
	ans=(ans+t*qmul(m,f())%p)%p;
	tot--;
	for(int j=last;j<=n-j;j++){
		l[++tot]=j;
		int t=s*inv(j)%p;
		if(j==last)dfs(j,k+1,n-j,t*inv(k+1)%p);
		else dfs(j,1,n-j,t*inv(1)%p);
		tot--;
	}
}
int main(){
	//freopen("a.in","r",stdin);
	scanf("%lld%lld%lld",&n,&m,&p);
	ll s=1;
	for(int i=1;i<=n;i++)s=s*i%p;
	dfs(1,0,n,s);
	for(int i=1;i<=n;i++)ans=ans*inv(i)%p;
	printf("%d\n",ans);
	return 0;
}