【HDOJ】【4336】Card Collector

概率DP/数学期望/状压DP/容斥原理

　　kuangbin总结中的第14题

　　好神奇的做法……题解看kuangbin的代码好了……

//HDOJ 4336
#include<cstdio>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
const int N=22;

double p[N],f[1<<N];
int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
        double tt=0.0;
        rep(i,n){
            scanf("%lf",&p[i]);
            tt+=p[i];
        }
        tt=1-tt;
        f[(1<<n)-1]=0;
        D(i,(1<<n)-2,0){
            double x=0,sum=1.0;
            rep(j,n)
                if(i & (1<<j))x+=p[j];
                else sum+=p[j]*f[i|(1<<j)];
            f[i]=sum/(1-tt-x);
        }
        printf("%.5lf
",f[0]);
    }
    return 0;
}

//HDOJ 4336
#include<cstdio>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)

double p[22];
int main(){
    int n;
    while(scanf("%d",&n)!=EOF){
        rep(i,n) scanf("%lf",&p[i]);
        double ans=0;
        for(int i=1;i<(1<<n);++i){
            int cnt=0; double sum=0;
            rep(j,n) if(i&(1<<j)){
                sum+=p[j];
                cnt++;
            }
            if (cnt&1) ans+=1.0/sum;
            else ans-=1.0/sum;
        }
        printf("%.5lf
",ans);
    }
    return 0;
}