HDU 4336 Card Collector 状态压缩+概率DP

state是压缩的状态，0代表已取，1代表未取, dp[state]代表在state这个状态取完全部卡片的期望步数，然后枚举在这个状态下取的卡片，就可以实现转移，这个dp过程就结束了，答案为dp[(1<<n)-1]。

下面我们思考另外一种dp思路，dp[state]代表达到state状态所需的期望步数,然后枚举能够到达state状态的状态，这样理论上也能实现转移，答案是dp[0]。但是我个人数学能力有限，算不出这种方法中转移的概率，只能放弃这种思路...求指点...

下面是方法一的代码:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define eps 1e-4
double P[25];
double Self[1<<20];
double DP[1<<20];
void Pre(int n)
{
    int i,j,lim=(1<<n);
    double tmp;
    for(i=0;i<lim;++i)
    {
        tmp=P[0];
        for(j=0;j<n;++j)
        {   
            if((i&(1<<j))==0)
            {
                tmp+=P[j+1];
            }
        }
        Self[i]=tmp;
    }
}
double Work(int n)
{
    int i,j,lim=(1<<n);
    double tmp;
    DP[0]=0;
    for(i=1;i<lim;++i)
    {
        tmp=1.0;
        for(j=0;j<n;++j)
        {
            if(i&(1<<j))
            {
                tmp+=P[j+1]*DP[i^(1<<j)];
            }
        }
        DP[i]=tmp/(1-Self[i]);
    }
    return DP[lim-1];
}
int main()
{
    int n,i;
    double tmp;
    while(scanf("%d",&n)!=EOF)
    {
        tmp=0.0;
        for(i=1;i<=n;++i)
        {
            scanf("%lf",&P[i]);
            tmp+=P[i];
        }
        P[0]=1.0-tmp;
        Pre(n);
        printf("%.5f\n",Work(n));
    }
    return 0;
}