幻想乡三连B:连在一起的幻想乡
$G[k][x]$表示所有$x$个点的无向图中每一个图的边数的$k$次方之和。
$F[k][x]$就是在$G[k][x]$的基础上加了一个整体连通的性质。
有一个经典的套路就是对于$F$在对应的$G$中刨去枚举$1$号节点所在的连通块大小的答案。
最后一个难点就是对于形如$\sum(x+y)^2$可以转化为$\sum x^2 +2(\sum x)(\sum y)+\sum y^2$。
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> #include<cmath> #define LL long long #define M 2002 using namespace std; int read(){ int nm=0,fh=1; char cw=getchar(); for(;!isdigit(cw);cw=getchar()) if(cw=='-') fh=-fh; for(;isdigit(cw);cw=getchar()) nm=nm*10+(cw-'0'); return nm*fh; } int G[3][M],F[3][M],C[M][M],num; const int n=read(),mod=read(); int mul(int x,int y){return (LL)x*(LL)y%mod;} int add(int x,int y){return (x+y)>=mod?x+y-mod:x+y;} void upd(int &x,int y){x=add(x,y);} int qpow(int x,int sq){ if(sq<0) return 0; int res=1; x%=mod; while(sq){ if(sq&1) res=mul(res,x); x=mul(x,x),sq>>=1; } return res; } int main(){ C[0][0]=1; for(int i=1;i<=n;i++){ num=(i*(i-1)>>1),C[i][0]=C[i][i]=1,G[0][i]=qpow(2,num),G[1][i]=mul(num,qpow(2,num-1)); for(int j=0;j<i;j++)upd(G[2][i],mul(add(add(G[2][i-1],mul(j<<1,G[1][i-1])),mul(j*j,G[0][i-1])),C[i-1][j])); for(int j=1;j<i;j++){ C[i][j]=add(C[i-1][j],C[i-1][j-1]),upd(F[0][i],mul(C[i-1][j-1],mul(F[0][j],G[0][i-j]))); upd(F[1][i],mul(C[i-1][j-1],add(mul(F[1][j],G[0][i-j]),mul(F[0][j],G[1][i-j])))); upd(F[2][i],mul(C[i-1][j-1],add(add(mul(F[2][j],G[0][i-j]),mul(2,mul(F[1][j],G[1][i-j]))),mul(F[0][j],G[2][i-j])))); } for(int j=0;j<3;j++) F[j][i]=mod-F[j][i],upd(F[j][i],G[j][i]); } printf("%d\n",F[2][n]); return 0; }