BZOJ 1478 Sgu282 Isomorphism 置换

来源：互联网发布：新闻发稿平台源码编辑：程序博客网时间：2024/05/22 01:45

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

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#define N 1010using 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);}

0 0