hdu 4336 Card Colletor 动态规划-概率dp

题意：有n个卡片，你现在买一包方便面，每包方便面出现其中一个卡片的概率为 p[i] ，问你集齐一套卡片需要的张数的数学期望。每包里至多只有一张卡片。

设E[i]为到达目标状态需要的步数，i为2进制数，2的20次方表示20张卡片是否得到的状态

E[i]=Po(E[i]+1)+    P[j](E[i|(1<<j)]+1)(j从0到n-1的求和）

Po表示拿到的包里没有卡片的概率，因为多拿了一包，所以步数++

P[J]表示拿到的包里是卡片j的概率，因为多拿了一包，所以步数++

可能拿到的是原有的卡片，即还是为i，也有可能拿到了新的卡片

#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
#define maxn 100010
#define rep(i,n) for(int i=1;i<=n;i++)
using namespace std;

double dp[1<<20],p[22];

int main()
{
    int n,tot;
    double tmp;
    while(scanf("%d",&n)!=EOF)
    {
        for(int k=0;k<n;k++)
            scanf("%lf",&p[k]);
        tot=(1<<n)-1;
        dp[tot]=0.0;//即从末状态到目标状态的步数为0
        for(int i=tot-1;i>=0;i--)
        {
            dp[i]=1;//从目标状态的前一步到目标状态的步数为1
            tmp=0.0;
            for(int j=0;j<n;j++)
            {
                if(i&(1<<j)) continue;
                dp[i]+=dp[i|(1<<j)]*p[j];
                tmp+=p[j];
            }
            dp[i]/=tmp;//可以到达i这种状态的状态都找到了
   //在循环里累加的是期望值 要除概率和
        }
        printf("%lf\n",dp[0]);//即从0到目标的步数
    }
    return 0;
}