【NOIP提高】礼物

Description

夏川的生日就要到了。作为夏川形式上的男朋友，季堂打算给夏川买一些生日礼物。
商店里一共有种礼物。夏川每得到一种礼物，就会获得相应喜悦值Wi（每种礼物的喜悦值不能重复获得）。
每次，店员会按照一定的概率Pi（或者不拿出礼物），将第i种礼物拿出来。季堂每次都会将店员拿出来的礼物买下来。
众所周知，白毛切开都是黑的。所以季堂希望最后夏川的喜悦值尽可能地高。
求夏川最后最大的喜悦值是多少，并求出使夏川得到这个喜悦值，季堂的期望购买次数。

Solution

比赛的时候怎么都看不懂，果断弃疗。
原来题目也不过如此。
题意就是说，每次可以按照概率取出每个数或者不去，最终每个物品都取到的期望步数。

状压DP

由于正着不方便统计，那么试着倒着做。
设f[i]为当前取的数状态为i的期望步数。
那么找到所有的位置j满足i or （1<<(j-1)）>0，就是说j这个位置可以把i的一个零的位置填成1，那么很显然f[i]+=(f[i|(1<<(j-1))]+1)*p[j]。
但是也会有选了已经选过的情况，那么f[i]+=f[i]*p[j]，这样的情况要把f[i]转到左边去，所以要统计出所有f[i]的系数。
还有不选的情况，那么也跟上面的一样，不选的状态也是f[i]，只用把系数减掉的不选的情况就可以了。
最后等式的右边在除以f[i]的系数就求出f[i]的值了。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=21;
int i,j,k,l,t,n,m;
long long ans;
double p[maxn],chu,f[1<<21],o,zuo,you;
long long a[maxn],er[maxn];
int main(){
    scanf("%d",&n);
    fo(i,1,n){
        scanf("%lf%d",&p[i],&a[i]);
        ans+=a[i];chu+=p[i];
    }
    chu=1-chu;
    fo(i,1,n)er[i]=(1<<(i-1));
    fod(i,(1<<n)-2,0){
        zuo=1,f[i]=1;
        fo(j,1,n){
            if((i|er[j])==i)zuo-=p[j];
            else f[i]+=p[j]*f[i|er[j]];
        }
        zuo-=chu;
        f[i]/=zuo;
    }
    printf("%lld\n",ans);
    printf("%.3lf\n",f[0]);
}