【模板】Min-Max容斥-bzoj4036按位或

最新推荐文章于 2021-03-05 16:41:49 发布

ccosi

最新推荐文章于 2021-03-05 16:41:49 发布

阅读量248

点赞数

分类专栏： -------MATH------- FMT 子集卷积容斥原理

本文链接：https://blog.csdn.net/corsica6/article/details/86660299

版权

-------MATH------- 同时被 3 个专栏收录

80 篇文章 1 订阅

订阅专栏

容斥原理

8 篇文章 0 订阅

订阅专栏

FMT

2 篇文章 0 订阅

订阅专栏

传送门:bzoj4036

Min-Max容斥

设 $\max(S)$ 为集合 $S$ 中的最大值， $\min(S)$ 为集合 $S$ 中的最小值， $∣ S ∣$ 为集合 $S$ 的大小，则存在：
$\max(S)=\sum_{T\subseteq S}(-1)^{|T|+1}\min(T) \\ \min(S)=\sum_{T\subseteq S}(-1)^{|T|+1}\max(T)$

同样地， $M i n - M a x$ 容斥在期望下成立，即：
$E(\max(S))=\sum_{T\subseteq S}(-1)^{|T|+1}E(\min(T)) \\ E(\min(S))=\sum_{T\subseteq S}(-1)^{|T|+1}E(\max(T))$

(证明略)

题解

套用 $E(\max(S))=\sum_{T\subseteq S}(-1)^{|T|+1}E(\min(T))$

其中 $E(\min(T))=\dfrac{1}{\sum\limits_{G\cap T\neq \phi}p[G]}$

求出 $P[S]=\sum\limits_{T\subseteq S}p[T]$ ，则 $E(\min(T))=\dfrac{1}{1-P[(2^n-1)\oplus T]}$

求 $P$ 做个 $F W T$ 即可。

代码

#include<bits/stdc++.h>
using namespace std;
typedef double db;
const int N=(1<<20)+5;
const db eps=1e-10;

int n,s,sz[N];
db p[N],ans;

int main(){
    int i,j,k;
    scanf("%d",&n);s=1<<n;
    for(i=0;i<s;++i) scanf("%lf",&p[i]),sz[i]=sz[i>>1]+(i&1);
    for(i=1;i<s;i<<=1)
        for(j=0;j<s;j+=(i<<1))
            for(k=0;k<i;++k)
              p[i+j+k]+=p[j+k];
    for(i=0,s--;i<=s;++i) if((1-p[s^i])>eps) ans+=((sz[i]&1)?1:(-1))/(1-p[s^i]);
    if(ans<eps) printf("INF");else printf("%.10lf",ans);
    return 0;
}