[HAOI2015]按位或

一、题目

点此看题

二、解法

设$max(S)$表示组成$S$的最后一个元素的出现时间，$min(S)$表示组成$S$的第一个元素的出现时间，这个是用时间顺序作为大小关系的，所以可以套用$\text{min-max}$容斥，下文的$E$表示期望，$U$表示全集，我们要求的是$E(max(U))$：
$$E(max(U))=\sum _{S\ne \varnothing }(-1{)}^{|s|-1}E(min(S))$$关键在于如果快速求$E(min(S))$了，设$p(min(S)=i)$表示这个出现时间为$i$的概率，那么：
$$E(min(S))=\sum _{i=1}^{+\infty }i\times p(min(S)=i)$$现在考虑如何计算$p(min(S)=i)$，我们前$i-1$个时间选不到，正好在$i$时间选到：
$$p(min(S)=i)=(\sum _{T\cap S=\varnothing }a[T]{)}^{i-1}(1-\sum _{T\cap S=\varnothing }a[T])$$方便推导，我们设$x={\sum }_{T\cap S=\varnothing }a[T]$，把上面的式子带到$E(min(S))$中去，可以发现：
$$E(min(S))=\sum _{i=1}^{+\infty }i(x{)}^{i-1}(1-x)=\sum _{i=1}^{\infty }i{x}^{i-1}-i{x}^i$$仔细观察一下，就会化成这个形式：
$$E(min(S))=\sum _{i=0}^{\infty }{x}^i=\frac{1}{1-x}=\frac{1}{1-{\sum }_{T\cap S=\varnothing }a[T]}$$其实吧，${\sum }_{T\cap S=\varnothing }a[T]={\sum }_{T\subseteq C_{U}S}a[T]$，这显然就是$or$卷积可以求出来的东西，套板子预处理，然后用第一个式子算就行了，时间复杂度$O(n{2}^n)$

#include <cstdio>
#define db double
const int M = 1<<20;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,bit[M];db ans,a[M];
void fwt(db *a,int n,int op)
{
	for(int i=1;i<n;i<<=1)
		for(int p=i<<1,j=0;j<n;j+=p)
			for(int k=0;k<i;k++)
			{
				if(op==1) a[i+j+k]=a[i+j+k]+a[j+k];
				else a[i+j+k]=a[i+j+k]-a[j+k];
			}
}
int main()
{
	n=read();m=1<<n;
	for(int i=0;i<m;i++)
	{
		scanf("%lf",&a[i]);
		bit[i]=bit[i>>1]+(i&1);
	}
	fwt(a,m,1);
	for(int i=1;i<m;i++)
	{
		int Cu=i^(m-1);
		if(1-a[Cu]<1e-8)
		{
			puts("INF");
			return 0;
		}
		db v=1.0/(1.0-a[Cu]);
		ans+=(bit[i]%2?1:-1)*v;
	}
	printf("%.10lf\n",ans);
}