【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=;

 double p[N],f[<<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=-tt;
         f[(<<n)-]=;
         D(i,(<<n)-,){
             double x=,sum=1.0;
             rep(j,n)
                 if(i & (<<j))x+=p[j];
                 else sum+=p[j]*f[i|(<<j)];
             f[i]=sum/(-tt-x);
         }
         printf("%.5lf\n",f[]);
     }
     return ;
 }

（概率DP）

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

（容斥原理）