BZOJ4036 HAOI2015按位或（概率期望+容斥原理）

　　考虑min-max容斥，改为求位集合内第一次有位变成1的期望时间。求出一次操作选择了S中的任意1的概率P[S]，期望时间即为1/P[S]。

　　考虑怎么求P[S]。P[S]=∑p[s] (s&S>0)=1-∑p[s] (s&S==0)。做一个高维前缀和即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define N (1<<20)
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
const double eps=1E-;
int n,t,s[N];
double p[N],ans;
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4036.in","r",stdin);
    freopen("bzoj4036.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=;i<(<<n);i++)
    {
        scanf("%lf",&p[i]);
        if (p[i]>eps) t|=i;
    }
    if (t!=(<<n)-) {cout<<"INF";return ;}
    for (int i=;i<n;i++)
        for (int j=;j<(<<n);j++)
        if (j&(<<i)) p[j]+=p[j^(<<i)];
    for (int i=;i<(<<n);i++)
    {
        s[i]=s[i^(i&-i)]+;
        if (s[i]&) ans+=/(-p[(<<n)-^i]);
        else ans-=/(-p[(<<n)-^i]);
    }
    printf("%.6f",ans);
}