bzoj4036 [HAOI2015]按位或（概率与期望+fmt）

vfk的论文题。具体可以参见国家集训队2015论文集 吕凯风的《集合幂级数的性质与应用及其快速算法》
这题就是求集合并卷积，用快速莫比乌斯变换即可。
答案就是$\sum\limits_{k=1}^{oo}k*(p^k[U]-p^{k-1}[U])$
莫比乌斯变换之后，化简一下，就得到了一个等比数列，因为$\sum_p=1$，所以除了最后一项为0外，其余项均为$-\frac{1}{1-p}$
再快速莫比乌斯反演回来即可。答案就是U的系数。
复杂度$O(n2^n)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 25
#define eps 1e-7
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,bin[N];
double a[1100000];
int main(){
//  freopen("a.in","r",stdin);
    n=read();bin[0]=1;
    for(int i=1;i<=n;++i) bin[i]=bin[i-1]<<1;
    for(int i=0;i<bin[n];++i) scanf("%lf",&a[i]);
    for(int i=0;i<n;++i)
        for(int s=0;s<bin[n];++s)
            if(s&bin[i]) a[s]+=a[s^bin[i]];
    for(int i=0;i<bin[n];++i){
        if(a[i]>1-eps){
            if(i!=bin[n]-1){puts("INF");return 0;}
            a[i]=0;continue;
        }a[i]=1.0/(a[i]-1);
    }for(int i=0;i<n;++i)
        for(int s=0;s<bin[n];++s)
            if(s&bin[i]) a[s]-=a[s^bin[i]];
    printf("%.7lf\n",a[bin[n]-1]);
    return 0;
}