bzoj 4036: [HAOI2015]按位或【min-max容斥+FWT】

其实也不是FWT……我也不知道刷FWT专题问什么会刷出来这个东西
这是min-max容斥讲解:https://www.zybuluo.com/ysner/note/1248287
总之就是设min(s),max(s)分别表示集合s里最早和最晚出现的元素,显然E(amx(全集))就是答案
然后有这样的式子:

[E(max(s))=sum_{s'in s}E(min(s'))*(-1)^{|s'|+1} ]

[E(min(s))=frac{1}{sum_{s'cap s!=phi }p(s')} ]

当然不会证
然后就是先求min再求max,求min的时候注意到正安则反,取补集就是orFWT的正变换,O(nlogn)求出后枚举子集求max(全集)即可

#include<iostream>
#include<cstdio>
using namespace std;
const int N=2000005;
const double eps=1e-8;
int n,m,b[N];
double a[N],ans;
int main()
{
	scanf("%d",&n);
	m=(1<<n);
	for(int i=0;i<m;i++)
		scanf("%lf",&a[i]),b[i]=b[i>>1]+(i&1);
	for(int i=1;i<m;i<<=1)
		for(int j=0;j<m;j+=(i<<1))
			for(int k=0;k<i;k++)
				a[i+j+k]+=a[j+k];
	for(int i=1;i<m;i++)
		if(1-a[(m-1)^i]>eps)
			ans+=((b[i]&1)?1:-1)/(1-a[(m-1)^i]);
	if(ans<eps)
		puts("INF");
	else
		printf("%.8f
",ans);
	return 0;
}	
原文地址:https://www.cnblogs.com/lokiii/p/10052737.html