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);
}