[luogu 3175] [HAOI2015]按位或

题面

刚开始你有一个数字0,每一秒钟你会随机选择一个[0,2n-1]的数字,与你手上的数字进行按位或运算。问期望多少秒后,你手上的数字变成2n-1。

分析

前置知识:min-max容斥

记\(\max(S)\)为集合\(S\)中的最大值,\(\min(S)\)为集合\(S\)中的最小值(如果\(S=\emptyset\)

,那\(\max(S)=\min(S)=0\)),那么有

\[\max(S)=\sum _{T\subseteq S}\left( -1\right) ^{|T|-1} \min(T)
\]

这里感性理解一下就好了

前置知识:高维前缀和

戳这里


把二进制数看成一个集合,第i位为1表示集合里有元素i.设全集\(U\)为二进制数\(2^n-1\)对应的集合

设\(\max(S)\)为S中最后一个元素被或为1的期望时间,min就是最先被或为1的元素的期望时间,那么答案就是\(\max(U)\)

根据min-max容斥,我们有

\[\max(S)=\sum_{T \subseteq S} (-1)^{|T|-1} \min(T)
\]

因为是最先被或为1,根据定义我们有

\[\min(T)=\frac{1}{\sum_{X \subseteq U , X \cap T \neq \emptyset }p(x)}
\]

那么

\[\begin{aligned} \min(T) &= \frac{1}{\sum_{X \subseteq U , X \cap T \neq \emptyset }p(x)} \\ &=\frac{1}{1-\sum_{X \subseteq U,X \cap T = \emptyset} p(x)} \\&= \frac{1}{1-\sum_{X \subseteq U-T} p(x)}\end{aligned}
\]

其实就是用了2次补集转化,然后\(\sum_{X \subseteq U-T} p(x)\)显然就是一个高维前缀和,直接套模板就可以了

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 20
#define eps 1e-10
using namespace std;
int n;
double p[(1<<maxn)+5]; int count_1(int x){
int ans=0;
while(x){
if(x&1) ans++;
x>>=1;
}
return ans;
} int main(){
scanf("%d",&n);
for(int i=0;i<(1<<n);i++) scanf("%lf",&p[i]);
for(int i=0;i<n;i++){
for(int j=0;j<(1<<n);j++){
if(j&(1<<i)) p[j]+=p[j^(1<<i)];
}
}
double ans=0;
int all=(1<<n)-1;
for(int i=1;i<=all;i++){
if(fabs(1-p[all^i])<eps) continue;//防止除0错误
ans+=pow(-1,count_1(i)-1)*1/(1-p[all^i]);
}
if(fabs(ans)<eps) printf("INF");
else printf("%.10lf",ans);
}

[luogu 3175] [HAOI2015]按位或(min-max容斥+高维前缀和)的更多相关文章

  1. BZOJ4036:按位或 (min_max容斥&高维前缀和)

    Description 刚开始你有一个数字0,每一秒钟你会随机选择一个[0,2^n-1]的数字,与你手上的数字进行或(c++,c的|,pascal 的or)操作.选择数字i的概率是p[i].保证0&l ...

  2. [HAOI2015]按位或(min-max容斥,FWT,FMT)

    题目链接:洛谷 题目大意:给定正整数 $n$.一开始有一个数字 $0$,然后每一秒,都有 $p_i$ 的概率获得 $i$ 这个数 $(0\le i< 2^n)$.一秒恰好会获得一个数.每获得一个 ...

  3. luoguP3175 [HAOI2015]按位或 min-max容斥 + 高维前缀和

    考虑min-max容斥 \(E[max(S)] = \sum \limits_{T \subset S} min(T)\) \(min(T)\)是可以被表示出来 即所有与\(T\)有交集的数的概率的和 ...

  4. BZOJ4036 [HAOI2015]按位或 【minmax容斥 + 期望 + FWT】

    题目链接 BZOJ4036 题解 好套路的题啊,,, 我们要求的,实际上是一个集合\(n\)个\(1\)中最晚出现的\(1\)的期望时间 显然\(minmax\)容斥 \[E(max\{S\}) = ...

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

    其实也不是FWT--我也不知道刷FWT专题问什么会刷出来这个东西 这是min-max容斥讲解:https://www.zybuluo.com/ysner/note/1248287 总之就是设min(s ...

  6. Codeforces.449D.Jzzhu and Numbers(容斥 高维前缀和)

    题目链接 \(Description\) 给定\(n\)个正整数\(a_i\).求有多少个子序列\(a_{i_1},a_{i_2},...,a_{i_k}\),满足\(a_{i_1},a_{i_2}, ...

  7. 【BZOJ4036】按位或(Min-Max容斥,FWT)

    [BZOJ4036]按位或(Min-Max容斥,FWT) 题面 BZOJ 洛谷 题解 很明显直接套用\(min-max\)容斥. 设\(E(max\{S\})\)表示\(S\)中最晚出现元素出现时间的 ...

  8. [Hdu-6053] TrickGCD[容斥,前缀和]

    Online Judge:Hdu6053 Label:容斥,前缀和 题面: 题目描述 给你一个长度为\(N\)的序列A,现在让你构造一个长度同样为\(N\)的序列B,并满足如下条件,问有多少种方案数? ...

  9. luogu P3175 [HAOI2015]按位或

    传送门 如果每个位置上的数字的意义是这个位置被加进集合的最早时间,那么我们要求的就是集合中最大数的期望,使用Min-Max容斥,\(E(max(S))=\sum_{T\subset S}(-1)^{| ...

随机推荐

  1. jquery delegate()方法 语法

    jquery delegate()方法 语法 作用:delegate() 方法为指定的元素(属于被选元素的子元素)添加一个或多个事件处理程序,并规定当这些事件发生时运行的函数.使用 delegate( ...

  2. [模板] KMP算法/Border

    KMP 算法 KMP (Knuth-Morris-Pratt) 算法是一种在线性时间内匹配文本串和模式串的算法. 称字符串的 Border 集合为 \[ \operatorname {Border} ...

  3. codevs 1077 多源最短路x

                         题目描述 Description 已知n个点(n<=100),给你n*n的方阵,a[i,j]表示从第i个点到第j个点的直接距离. 现在有Q个询问,每个询 ...

  4. (Java多线程系列一)快速入门

    Java多线程快速入门 1.线程和进程的区别 进程是所有线程的集合,每一个线程是进程的一条执行路径. 2.多线程的应用场景 多线程主要体现在提高程序的效率,比如迅雷多线程下载,多线程分批发送短信等. ...

  5. PHP 判断当前日期是否是法定节假日或者休息日 (原)

    //判断日期是不是节假日或者双休日接口 @param time [时间(时间戳或者Y-m-d都可)] public function get_type_by_date(){ $t = $_GET['t ...

  6. Vue之vue中的data为什么是一个函数+vue中路径别名alias设置

    问题描述 为什么在vue组件中,我们的data属性必须是一个函数,new Vue()中的data除外,因为new Vue中只有一个data属性. 原因 因为我们能抽离出来的组件,肯定是具有复用性的,它 ...

  7. Kaggle 网站上的练习

    快速的数据科学教程: https://www.kaggle.com/learn/overview

  8. ES6指北【1】——let、const

    1.如何学习ES6 1.1 js的学习顺序 ES5 -> ES6 -> ES7 -> ES8 以此类推 ES5没学好就别想学好ES6 1.2 边学边用 学了就要用 2.变量声明的方式 ...

  9. .tcc文件

    今天看源码时碰到一个MemoryPool.h文件和MemoryPool.tcc文件,毫不犹豫在vs工程下把.tcc加到了源文件文件夹下, 把.h文件放到了头文件文件夹下.结果闹了笑话: 以下是解释, ...

  10. C++如何限制对象在堆上或栈上生成

    1,限制类的对象只能生成在栈上 将 operator new 各种原型设为私有 #include <iostream> class OnlyOnStack { public: OnlyOn ...