BZOJ 1478 Sgu282 Isomorphism 置换

题意:
一个无向完全图，m种颜色对所有边进行染色，求任意更改点的顺序不同构的图的个数。
解析:
做完1488再看这题简直sb题。
1488是边连或者不连，其实就是两种颜色在染色。
所以这题就是1488的强化…
然并卵直接水就行。
代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 1010
using namespace std;
typedef long long ll;
ll n,cnt,ans,m,mod;
ll powtwo[N],factor[N],invfac[N],num[N],val[N],gcd[N][N],inv[N];
ll get_gcd(ll x,ll y)
{
 while(y)
 {
 ll t=y;
 y=x%y;
 x=t;
 }
 return x;
}
ll get_inv(ll x,ll y)
{
 ll ret=1;
 while(y)
 {
 if(y&1)ret=(ret*x)%mod;
 x=(x*x)%mod;
 y>>=1;
 }
 return ret;
}
void init()
{
 powtwo[0]=factor[0]=invfac[0]=1;
 for(ll i=1;i<=1000;i++)
 {
 factor[i]=factor[i-1]*i%mod;
 }
}
void dfs(ll now_num,ll left)
{
 if(left==0)
 {
 ll retnow=0;
 ll bot=1;
 for(ll i=1;i<=cnt;i++)
 {
 retnow+=num[i]*(num[i]-1)/2*val[i]+val[i]/2*num[i];
 for(ll j=i+1;j<=cnt;j++)
 retnow+=num[i]*num[j]*get_gcd(val[i],val[j]);
 }
 for(ll i=1;i<=cnt;i++)
 {
 bot=(bot*get_inv(val[i],num[i])%mod*factor[num[i]])%mod; 
 } 
 bot=get_inv(bot,mod-2)*factor[n]%mod; 
 ans=(ans+get_inv(m,retnow)*bot%mod)%mod; 
 }
 if(now_num>left)return;
 dfs(now_num+1,left);
 for(ll i=1;i*now_num<=left;i++)
 {
 val[++cnt]=now_num,num[cnt]=i;
 dfs(now_num+1,left-i*now_num);
 cnt--;
 }
}
int main()
{
 scanf("%lld%lld%lld",&n,&m,&mod);
 init();
 dfs(1,n);
 ans=ans*get_inv(factor[n],mod-2)%mod;
 printf("%lld\n",ans);
}