传送门:bzoj3925
题解
粘一波题解
代码
#include<bits/stdc++.h>
#define db double
using namespace std;
typedef long long ll;
const int N=(1<<11);
ll f[N][50],g[N+10][50],c[50][50];
db ans;
int n,m,e[50],cnt[N],sz[N];
int main(){
scanf("%d%d",&n,&m);
for(int x,y,i=0;i<m;++i){
scanf("%d%d",&x,&y);x--;y--;
e[x]|=(1<<y);e[y]|=(1<<x);
}c[0][0]=1;
for(int j,i=1;i<=m;++i){
c[i][0]=c[i][i]=1;
for(j=1;j<i;++j)
c[i][j]=c[i-1][j-1]+c[i-1][j];
}
for(int t,s=0;s<(1<<n);++s){
sz[s]=sz[s>>1]+(s&1);
if(sz[s]==1){g[s][0]=1;continue;}
for(int i=0;i<n;++i)
if(s & (1<<i))
cnt[s]+=sz[e[i]&s];
cnt[s]>>=1;t=(s&(-s));
for(int i,j,k=(s-1)&s;k;k=(k-1)&s)
if(k & t)
for(i=0;i<=cnt[k];++i)
for(j=0;j<=cnt[s^k];++j)
f[s][i+j]+=g[k][i]*c[cnt[s^k]][j];
for(int i=0;i<=cnt[s];++i)
g[s][i]=c[cnt[s]][i]-f[s][i];
}
n=(1<<n)-1;
for(int i=0;i<=m;i++) ans+=(db)f[n][i]/(db)c[cnt[n]][i];
printf("%.6lf\n",ans/(db)(m+1));
}