[BZOJ4036] set - 子集和变换 - 容斥原理/莫比乌斯反演定理

虽然有点长还是建议看这篇博客，大意就是通过子集和变换构造出函数F使得可以快速计算∑[0,+∞)F（当其收敛），然后通过消无穷，最后得出式子就可以辣。题解 

#include"bits/stdc++.h"
using namespace std;
 
const double eps=1e-6;
const int N=1048576+5;
double a[N],ans;int n,m,o[N];
 
#define lowbit(x) ((x)&-(x))
 
int main(){
    scanf("%d",&n);m=1<<n;int i,j;
    for(i=0;i<m;i++)scanf("%lf",&a[i]);
    for(i=0;i<n;i++)for(j=0;j<m;j++)
        if(j&(1<<i))a[j]+=a[j^(1<<i)];
    for(i=0;i<m-1;i++)if(a[i]+eps>1.0)
    {puts("INF");return 0;}o[0]=n&1?1:-1;
    for(i=1;i<m;i++)o[i]=-o[i-lowbit(i)];
    for(i=0;i<m-1;i++)ans+=o[i]*a[i]/(1-a[i]);
    printf("%.10lf\n",ans+1.0);
    return 0;
}