【BZOJ】3925地震后的幻想乡-概率dp

传送门: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));
}