标签：块大小 printf efi inf bsp str sum cpp ret
$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;
}
标签：块大小 printf efi inf bsp str sum cpp ret
原文地址：https://www.cnblogs.com/OYJason/p/9651783.html
