【HAOI2015】按位或【Min-Max容斥】【FWT】

传送门

题意:开始时你有一个数000,每次选出[0,2n−1][0,2^n-1][0,2n1]中的一个数进行按位或,每个数选中的概率给定。求得到2n−12^n-12n1的期望操作次数。

1≤n≤201\leq n\leq 201n20

神仙题

首先发现每一位都是独立的,可以分开考虑。

对于一个集合SSS,记F(S)F(S)F(S)表示SSS的每一个元素出现的期望操作次数组成的可重集,记Max(S)=max(F(S)),Min(S)=min(F(S))Max(S)=max(F(S)),Min(S)=min(F(S))Max(S)=max(F(S)),Min(S)=min(F(S))

S=2n−1S=2^n-1S=2n1,我们要求的就是E(Max(S))E(Max(S))E(Max(S))

这个不好求,但可以Min-Max容斥一下?

E(Max(S))=∑T⊆S(−1)∣T∣+1E(Min(T))E(Max(S))=\sum_{T\subseteq S}(-1)^{|T|+1}E(Min(T))E(Max(S))=TS(1)T+1E(Min(T))

可以直接枚举子集,现在考虑怎么求E(Min(T))E(Min(T))E(Min(T))

这玩意的意义是或到和TTT有交集的期望次数

设每次选出一个和TTT有交集的数的概率是ppp

E(Min(T))=∑i=1∞i(1−p)i−1pE(Min(T))=\sum_{i=1}^{\infin}i(1-p)^{i-1}pE(Min(T))=i=1i(1p)i1p

=p∑i=1∞i(1−p)i−1=p\sum_{i=1}^{\infin}i(1-p)^{i-1}=pi=1i(1p)i1

s=∑i=1∞i(1−p)i−1s=\sum_{i=1}^{\infin}i(1-p)^{i-1}s=i=1i(1p)i1

(1−p)s=∑i=2∞(i−1)(1−p)i−1(1-p)s=\sum_{i=2}^{\infin}(i-1)(1-p)^{i-1}(1p)s=i=2(i1)(1p)i1

=∑i=1∞(i−1)(1−p)i−1=\sum_{i=1}^{\infin}(i-1)(1-p)^{i-1}=i=1(i1)(1p)i1

ps=∑i=1∞(1−p)i−1ps=\sum_{i=1}^{\infin}(1-p)^{i-1}ps=i=1(1p)i1

=∑i=0∞(1−p)i=\sum_{i=0}^{\infin}(1-p)^i=i=0(1p)i

(1−p)ps=∑i=1∞(1−p)i(1-p)ps=\sum_{i=1}^{\infin}(1-p)^i(1p)ps=i=1(1p)i

p2s=1p^2s=1p2s=1

s=1p2s=\frac{1}{p^2}s=p21

E(Min(T))=ps=1pE(Min(T))=ps=\frac{1}{p}E(Min(T))=ps=p1

ppp仍然不好求。

正难则反,我们求和TTT没有交集的概率

TTT的补集的子集的概率之和

第一次写的子集标记,发现会被算多次

然后发现不是个FWTFWTFWT板子吗

然后没了

注意不枚举空集 如果遇到无穷大直接输出

复杂度O(n2n)O(n2^n)O(n2n)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
using namespace std;
double p[1<<20];
const int d[]={0,1,1,2,1,2,2,3};
inline int count(int x)
{int ans=0;while (x){ans+=d[x&7];x>>=3;}return ans;
}
int main()
{int n;scanf("%d",&n);for (int i=0;i<(1<<n);i++) scanf("%lf",&p[i]);for (int mid=1;mid<(1<<n);mid<<=1)for (int s=0;s<(1<<n);s+=(mid<<1))for (int k=0;k<mid;k++)p[s+mid+k]+=p[s+k];double sum=0;for (int i=1;i<(1<<n);++i){double t=1-p[(~i)&((1<<n)-1)];if (t<1e-10){puts("INF");return 0;}t=1/t;sum+=((count(i)&1)? t:-t);}printf("%.10f\n",sum);return 0;
}

wtcl

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315227.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Visual Basic 兴衰记

这是一篇暴露年龄的文章。为什么这么说&#xff1f;因为现在年轻的程序员可能没有接触过当年红极一时的 Visual Basic。28 年前的 1991 年 4 月&#xff0c;Microsoft 发布了 Visual Basic 1.0 for Windows&#xff0c;次年 9 月&#xff0c;发布了 Visual Basic 1.0 for DOS。…

Ozon Tech Challenge 2020 (Div.1 + Div.2, Rated) D. Kuroni and the Celebration 交互 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一颗树&#xff0c;每次可以询问两个点的lcalcalca&#xff0c;询问次数不能超过⌊n2⌋\left \lfloor \frac{n}{2} \right \rfloor⌊2n​⌋&#xff0c;求它的根。 n≤1e3n\le1e3n≤1e3。 思路&#xf…

【CF1047D】Little C Loves 3 II【构造】【赛瓦维斯特定理】

传送门 题意&#xff1a;给一个NMN \times MNM的空棋盘&#xff0c;每次选取两个曼哈顿距离为3的空格子放上棋子&#xff0c;问最多能放多少个。 1≤N,M≤1e91 \leq N,M \leq 1e91≤N,M≤1e9 暴力讨论 假装N≤MN \leq MN≤M ①N1N1N1 容易得到&#xff0c;详见代码 ②N2N…

函数式编程里的Materialization应该翻译成什么?

Materialization是函数式编程里的一个专业术语, 用于特指函数式编程中查询被实际执行并生成结果的这一过程.首先, 搜了一下中文资料, 暂时没有对该词的中文翻译, CSDN\博客园\阿里云上所有关于它的文档都没有做中文翻译, 直接沿用这个英文单词.难道轮到我来创造这个词的中文翻…

Ozon Tech Challenge 2020 (Div.1 + Div.2) E.Kuroni and the Score Distribution 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 不难想到&#xff0c;长度为nnn的数组最多的满足条件的三元组序列是1,2,3....,n1,2,3....,n1,2,3....,n&#xff0c;对于每一个位置贡献为i−12\frac{i-1}{2}2i−1​&#xff0c;那么如果m…

【CF1209E】Rotate Columns【状压DP】【位运算】【贪心】

题意&#xff1a;给一个NMN \times MNM的矩阵&#xff0c;可以进行任意多次操作将一列轮换&#xff0c;求每一行的最大值之和的最大值。多组数据。 Easy VersionN≤4N \leq 4N≤4,M≤100M \leq100M≤100 Hard VersionN≤12N \leq 12N≤12,M≤2000M \leq2000M≤2000 看这数据…

Exceptionless(二) - 使用进阶

作者&#xff1a;markjiang7m2原文地址&#xff1a;https://www.cnblogs.com/markjiang7m2/p/11100563.html官网地址&#xff1a;http://letyouknow.net在上一篇文章Exceptionless - .Net Core开源日志框架中就说到如何对Exceptionless进行本地化部署&#xff0c;不过我也跟大家…

Ozon Tech Challenge 2020 (Div.1 + Div.2) F. Kuroni and the Punishment 随机化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个数&#xff0c;每次操作可以选择将某个数1,−11,-11,−1&#xff0c;求最少进行多少次操作使得所有数都为正数且gcd>1gcd>1gcd>1。 思路&#xff1a; 考虑gcd2gcd2gcd2的情况&#xff0…

【洛谷P4841】城市规划【指数型生成函数】【麦克劳林级数】【多项式对数】

传送门 题意&#xff1a;求NNN个点的带标号无向连通简单图的个数。 N≤130000N \leq 130000N≤130000 这个问题的主要矛盾在于连通 这个并不好表示&#xff0c;但可以用这个表示出不要求连通的方案数 由于带标号&#xff0c;先构造答案的EGF f(x)∑i0∞aii!xif(x)\sum_{i0}…

基于 Kong 和 Kubernetes 的 WebApi 多版本解决方案

前言大家好&#xff0c;很久没有写博客了&#xff0c;最近半年也是比较的忙&#xff0c;所以给关注我的粉丝们道个歉。去年和朱永光大哥聊的时候提了一下我们的这个方案&#xff0c;他说让我有空写篇博客讲一下&#xff0c;之前是非常的忙&#xff0c;所以这次趁着有些时间就写…

CodeCraft-20 (Div. 2) D. Nash Matrix 构造 + dfs

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给定一个n∗nn*nn∗n的矩阵&#xff0c;每个点上面都有一个xi,yix_i,y_ixi​,yi​&#xff0c;表示这个点到xi,yix_i,y_ixi​,yi​这个点停下&#xff0c;当xi−1,yi−1x_i-1,y_i-1xi​−1,yi​−1的时候代表…

【HAOI2018】染色【反向二项式反演】【NTT卷积】

传送门 题意&#xff1a;NNN个位置染MMM种颜色&#xff0c;恰好出现SSS次的颜色数量恰好为kkk时的愉悦度为wkw_kwk​,求所有方案的愉悦度之和。对100453580910045358091004535809取模。 N≤1e7N \leq 1e7N≤1e7,M≤1e5M \leq 1e5M≤1e5,S≤150S \leq 150S≤150 本题的恶心之处…

从严治码-别人在项目中下毒,我该怎么治?

01 从软考说起从4月份开始&#xff0c;由于备考《系统集成项目管理工程师》的原因&#xff0c;博客没有持续更新&#xff0c;在上半年考试结束之后&#xff0c;又对项目进行了一些收尾的工作。下面就这段时间的学习作一个记录和总结吧。在学习的过程中&#xff0c;提炼了一些自…

CodeCraft-20 (Div. 2) C. Primitive Primes 思维 + 数论

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度分别为n,mn,mn,m的多项式&#xff0c;将他们乘起来&#xff0c;问系数modp0\bmod p 0modp0的项的指数是多少&#xff0c;两个多项式所有项的系数gcd1gcd1gcd1。 n,m<1e6n,m<1e6n,m<1e6…

【LOJ166】拉格朗日插值2【拉格朗日插值】【NTT卷积】

传送门 题意&#xff1a;给定n,m,f(0),f(1),......,f(n)n,m,f(0),f(1),... ...,f(n)n,m,f(0),f(1),......,f(n)&#xff0c;求f(m),f(m1),......f(mn)f(m),f(m1),... ...f(mn)f(m),f(m1),......f(mn) 模998244353998244353998244353 n≤100000,m≤1e8,n<mn \leq 100000,m \…

学习MVVM设计模式后第一次用于生产

WPF的MVVM设计模式从winform转变到WPF的过程&#xff0c;难点主要还是在MVVM的设计模式。当然&#xff0c;如果依然采用winform的涉及方式&#xff0c;在每个控件背后绑定事件的方式运用在wpf中&#xff0c;依然可行&#xff0c;但是假如GUI改版&#xff0c;其背后绑定的特别为…

剑指 Offer 14- II. 剪绳子 II

给你一根长度为 n 的绳子&#xff0c;请把绳子剪成整数长度的 m 段&#xff08;m、n都是整数&#xff0c;n>1并且m>1&#xff09;&#xff0c;每段绳子的长度记为 k[0],k[1]...k[m - 1] 。请问 k[0]*k[1]*...*k[m - 1] 可能的最大乘积是多少&#xff1f;例如&#xff0c;…

Educational Codeforces Round 108 (Rated for Div. 2) D. Maximum Sum of Products 思维 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度为nnn的数组a,ba,ba,b&#xff0c;你可以至多反转一段连续区间&#xff0c;求∑i1nai∗bi\sum _{i1}^n a_i*b_i∑i1n​ai​∗bi​最大是多少。 n<5e3n<5e3n<5e3 思路&#xff1a; 首…

【CF1215E】Marbles【状压DP】

传送门 题意&#xff1a;给一个长为NNN的序列aaa&#xff0c;每次操作交换两个相邻位置&#xff0c;求最少操作次数使得所有相同的值连成一片。 N≤400000N \leq 400000N≤400000,ai≤20a_i \leq20ai​≤20 我们发现aia_iai​很小&#xff0c;盲猜单独考虑 我们重新确认一个…

netcore mvc快速开发系统(菜单,角色,权限[精确到按钮])开源

基于netcore2.0 mvc 开发的 快速搭建具有如下特色的后台管理系统用户管理菜单管理角色管理权限管理[精确到按钮]&#xff09;代码生成器代码克隆到本地 用vs2017或以上版本 打开工程。项目结构如下&#xff1a;找到DbModel下面的初始化db脚本里面包含4张表的schema和初始化数据…