BZOJ3925: [Zjoi2015]地震后的幻想乡【状压DP】

3925: [Zjoi2015]地震后的幻想乡

我们定义

$f[i][j]$表示当前点集状态为$i$，处理了前$j$条边，不连通的方案数

$g[i][j]$表示当前点集状态为$i$，处理了前$j$条边，连通的方案数

显然$f[i][j]+g[i][j]={(}_{cnt[i]}^j)$,$cnt[i]$表示点集为$i$的边的数量

考虑如何用$g$来计算$f$，$f[S][j]={\sum }_{T\subset S}g[T][k]\ast {(}_{cnt[S-T]}^{j-k})$

但是我们会发现这样子会有重复，所以我们可以随便选一个点，强制$S$和$T$都包含它，这样可以做到不重不漏。

最后答案就是$\frac{1}{m+1}{\sum }_{i=0}^m\frac{f[U][i]}{{(}_m^i)}$

#include<cstdio>
using namespace std;
typedef long long LL;
int n,m,Siz[1<<11];bool mp[14][14];LL f[1<<11][50],g[1<<11][50],C[50][50];double Ans;
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1,x,y;i<=m;i++) scanf("%d%d",&x,&y),x--,y--,mp[x][y]=mp[y][x]=1;
	for(int i=1;i<(1<<n);i++)
	for(int j=0;j<n;j++)
	for(int k=0;k<j;k++)
	Siz[i]+=(i&(1<<j))&&(i&(1<<k))&&mp[j][k];
	C[0][0]=C[1][0]=C[1][1]=1;
	for(int i=2;i<=m;i++){
		C[i][0]=1;
		for(int j=1;j<=i;j++) C[i][j]=C[i-1][j]+C[i-1][j-1]; 
	}
	for(int i=0;i<(1<<n);i++) f[i][0]=1;
	for(int i=0;i<n;i++) g[1<<i][0]=1,f[1<<i][0]=0;
	for(int j=1;j<=m;j++)
	for(int i=1;i<(1<<n);i++){
		int p=0;for(;((i>>p)&1)==0;p++);
		for(int s=(i-1)&i;s;s=(s-1)&i)
		if(s&(1<<p))
			for(int k=0;k<=j;k++) f[i][j]+=g[s][k]*C[Siz[i^s]][j-k];
		g[i][j]=C[Siz[i]][j]-f[i][j];
	}
	for(int i=0;i<=m;i++) Ans+=1.0*f[(1<<n)-1][i]/C[m][i];
	printf("%.6lf\n",Ans/(m+1));
	return 0;
}