按位或「HAOI2015」

题意

刚开始你有一个数字\(0\)，每一秒钟你会随机选择一个\([0,2^n-1]\)的数字，与你手上的数字进行或（c++,c的|,pascal的or）操作。选择数字\(i\)的概率是\(p[i]\)。保证\(0<=p[i]<=1\)，\(Σp[i]=1\)问期望多少秒后，你手上的数字变成\(2^n-1\)。

思路

minmax容斥加上高维前缀和。

minmax容斥就是两个式子：

\[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))\)表示\(S\)中最大的元素的期望，反之亦然。

这个东西的证明需要使用二项式反演，但我并不会，所以只好先背结论。如果要看证明可以去这里。

对于本题，如果我们把所有位被填上的期望时间放到一个集合里，那么求的答案就是这个集合的最大值，发现符合我们式子的左半边。

那么式子右边就是至少有一个数字出现的期望时间。这里需要用到一个结论：\(MIN(T)\)的期望为其概率的倒数。（这个证明不是很复杂，这里不放了）

那么有\(E(MIN(T))=\frac{1}{\sum_{s'\cap T\neq 0}~~P_{s'}}\)（公式好像有点抽搐）

但是求相交不是很容易，所以我们转换为求不相交，这个比较简单，也就是

\[E(MIN(T))=\frac{1}{1-\sum_{s'与T没有相交}~~~~~P_{s'}}\]

分母后面那个东西可以表示为\[\sum_{s'\subseteq (2^n-1-T)}P_{s'}\]

而这个东西可以高维前缀和搞，于是这道题就结束了。

代码

#include <bits/stdc++.h>

using namespace std;

namespace StandardIO{

    template<typename T>inline void read(T &x){
        x=0;T f=1;char c=getchar();
        for(; c<'0'||c>'9'; c=getchar()) if(c=='-') f=-1;
        for(; c>='0'&&c<='9'; c=getchar()) x=x*10+c-'0';
        x*=f;
    }

    template<typename T>inline void write(T x){
        if(x<0) putchar('-'),x*=-1;
        if(x>=10) write(x/10);
        putchar(x%10+'0');
    }

} using namespace StandardIO;

namespace Project{

    const int N=1<<21;

    int n,tot;
    double p[N];
    int num[N];
    double ans;

    inline void MAIN(){
        read(n);
        for(register int i=0; i<(1<<n); ++i){
            scanf("%lf",&p[i]);
            if(p[i]) tot|=i;
        }
        if(tot!=(1<<n)-1) return puts("INF"),void();
        for(register int i=1; i<(1<<n); ++i) num[i]=num[i>>1]+(i&1);
        for(register int i=0; i<n; ++i){
            for(register int j=0; j<(1<<n); ++j){
                if(j&(1<<i)) p[j]+=p[j^(1<<i)];
            }
        }
        for(register int i=1; i<(1<<n); ++i){
            ans+=(num[i]&1?1.0:-1.0)/(1.0-p[((1<<n)-1)^i]);
        }
        printf("%.6lf",ans);
    }

}

int main(){
//    freopen(".in","r",stdin);
//    freopen(".out","w",stdout);
    Project::MAIN();
}

原文地址：https://www.cnblogs.com/ilverene/p/11843381.html