P3175-[HAOI2015]按位或【min-max容斥,FWT】
正題
題目鏈接:https://www.luogu.com.cn/problem/P3175
題目大意
開始有一個nnn位二進制數(shù)s=0s=0s=0,每次有pip_ipi?概率選取數(shù)字iii讓sss或上這個數(shù)字iii,求期望多少次能夠讓sss的nnn個位都變?yōu)?span id="ze8trgl8bvbq" class="katex--inline">111。
解題思路
因為是或所以我們只關心最后一個選中的數(shù),設第iii位選中的期望次數(shù)為E(i)E(i)E(i)的話答案就是max{E(i)}max\{E(i)\}max{E(i)}。
又是期望又是maxmaxmax所以可以直接上min-max\text{min-max}min-max容斥,答案就是
∑T∈Smin{E(i)}(i∈T)?(?1)∣T∣+1\sum_{T\in S}min\{E(i)\}(i\in T)*(-1)^{|T|+1}T∈S∑?min{E(i)}(i∈T)?(?1)∣T∣+1
算這個東西的話也就是如果我們選中一個與TTT有交集的數(shù)就可以退出了。期望次數(shù)=1/期望概率。所以我們直接算期望概率
也就是我們要算所有∑G∩T≠?pG\sum_{G\cap T\neq \varnothing}p_{G}∑G∩T?=??pG?。GGG和TTT的交集非空就去掉所有交集為空的,交集為空的就是TTT的補集的子集和。
子集和的話就是直接拿ppp出來跑一次ororor的FWT\text{FWT}FWT的結果就是子集和了。
時間復雜度O(n2n)O(n2^n)O(n2n)
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; const int N=1<<21; const double eps=1e-8; int n;double cnt[N],p[N],ans; void FWT_or(double *f,int op){for(int p=2;p<=n;p<<=1)for(int k=0,len=p>>1;k<n;k+=p)for(int i=k;i<k+len;i++)f[i+len]+=f[i]*op;return; } int main() {scanf("%d",&n);cnt[0]=-1;n=1<<n;for(int i=0;i<n;i++)scanf("%lf",&p[i]);FWT_or(p,1);for(int i=0;i<n;i++){if(i)cnt[i]=-cnt[i-(i&-i)];double e=1-p[(n-1)^i];if(fabs(e)<eps)continue;ans+=cnt[i]*(1.0/e);}if(ans<eps)printf("INF");else printf("%.10lf",ans); }總結
以上是生活随笔為你收集整理的P3175-[HAOI2015]按位或【min-max容斥,FWT】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 拼图秀恩爱 带你玩转美图秀秀iPhone
- 下一篇: 变身幽灵海盗的幽灵烈酒如何获得?