题目描述

小w隐藏的心绪已经难以再隐藏下去了。小w有n+ 1(保证n为偶数)个心绪,每个都包含了[1,2n]的一个大小为n的子集。现在他要找到隐藏的任意两个心绪,使得他们的交大于等于n/2。

输入描述

一行一个整数n。接下来每行一个长度为k的字符串,该字符串是一个64进制表示,ASCII码为x的字符代表着x-33,所有字符在33到33+63之间。转为二进制表示有6k位,它的前2n个字符就是读入的集合,第i位为1表示这个集合包含i,为0表示不包含。

输出描述

一行两个不同的整数表示两个集合的编号。如果无解输出“NO Solution”。

样例输入

10
  EVK#
  IH=#
  676"
  R7,#
  74S"
  6V2#
  O3J#
  S-7$
  NU5"
  C[$$
  3N.#

样例输出

1 2

对于20%的数据满足n≤100。
  对于50%的数据满足n≤1×103
  对于100%的数据满足n≤6×103

分析

我觉得应该没有人会想到暴力即是正解

首先先看一眼数据规模,好像用int存字符串会炸内存,于是考虑用bitset

不然就像我一样爆零

然后就可以暴力地O(n3/32)地做啦

我们来验证一下它的正确性

任选出两个子集,考虑每个元素在两个子集中共有的概率,加起来就是期望共有的个数

由于总共有2n个元素,每个子集有n个元素,那么一个元素在子集内的概率就是$\frac 1 2$,

那么一个元素同时在两个子集内的概率就是$\frac 1 4$

所以任选两个子集期望的共有的个数是$\frac 1 4 \times 2n=\frac n 2$

当然,这个期望值是建立在数据纯随机的情况下的

我们接下来可以求一下在出题人控制了数据的情况下我们任选两个子集,子集的共有元素的期望个数

假设第i号元素有Si个子集拥有,那么对于i号元素就有$C_{Si}^2$种情况任选两个子集都包含i号元素,而任选两个子集的情况总数为$C_{n+1}^{2}$

所以对于i号元素,任选两个子集都包含它的概率为

$$\frac {C_{Si}^2}  {C_{n+1}^2}$$

任选两个子集的期望共有元素个数为

$$\sum_{i=1}^{2n} \frac {C_{Si}^2} {C_{n+1}^2}$$

显然出题人可以控制每个Si的大小来卡我们,但他无法控制其它数的大小。

我们只要求出这个东西的最小值就可以知道出题人是否能卡我们

而且我们还知道

$$\sum_{i=1}^{2n} Si=n(n+1)$$

所以我们化简一下这个式子

$$\sum_{i=1}^{2n} \frac {C_{Si}^2} {C_{n+1}^2}$$

$$=\frac { \sum_{i=1}^{2n} {C_{Si}^2} } {C_{n+1}^2}$$

$$=\frac {\sum_{i=1}^{2n} { \frac {Si \times (Si-1) } {2} } } { \frac {n \times (n+1)} {2} }$$

$$=\frac {\sum_{i=1}^{2n} {Si \times (Si-1)}} {n \times (n+1)}$$

$$=\frac {\sum_{i=1}^{2n} {Si^2} - \sum_{i=1}^{2n} {Si}} {n \times (n+1)}$$

$$=\frac {\sum_{i=1}^{2n} {Si^2} - n(n+1)} {n(n+1)}$$

接下来只需要根据$\sum_{i=1}^{2n} Si=n(n+1)$求出$\sum_{i=1}^{2n} {Si^2}$的最小值就好啦

我们可以小的情况推到大的情况,比如知道$a+b=x$,求$min(a^2+b^2)$

可以用均值不等式的思想来证明

因为$(a-b)^2 \geq 0$,所以$a^2+b^2 \geq 2ab$

又因为$(a+b)^2=x^2$,所以$2ab=x^2-a^2-b^2$

代入$a^2+b^2 \geq 2ab$就有$a^2+b^2 \geq {\frac {x^2} {2}}$

同理,因为

$$\sum_{i=1}^{2n}\sum_{j=1}^{2n}(Si-Sj)^2 \geq 0$$

所以

$$\sum_{i=1}^{2n}\sum_{j=1}^{2n}(Si^2+Sj^2-2SiSj) \geq 0$$

$$\sum_{i=1}^{2n}\sum_{j=1}^{2n}Si^2 + \sum_{i=1}^{2n}\sum_{j=1}^{2n}Sj^2 - \sum_{i=1}^{2n}\sum_{j=1}^{2n}2SiSj\geq 0$$

$$\sum_{i=1}^{2n}\sum_{j=1}^{2n}Si^2 + \sum_{i=1}^{2n}\sum_{j=1}^{2n}Sj^2\geq 2\sum_{i=1}^{2n}\sum_{j=1}^{2n}SiSj$$

$$2n\sum_{i=1}^{2n}Si^2 +2n\sum_{j=1}^{2n}Sj^2\geq 2\sum_{i=1}^{2n}\sum_{j=1}^{2n}SiSj$$

$$4n\sum_{i=1}^{2n}Si^2\geq 2\sum_{i=1}^{2n}\sum_{j=1}^{2n}SiSj$$

又因为$\sum_{i=1}^{2n} Si=n(n+1)$

所以

$$\left ( \sum_{i=1}^{2n} Si \right )^2=n^2(n+1)^2$$

$$\sum_{i=1}^{2n}\sum_{j=1}^{2n}SiSj=n^2(n+1)^2$$

代入上面的不等式就有

$$4n\sum_{i=1}^{2n}Si^2\geq 2\sum_{i=1}^{2n}\sum_{j=1}^{2n}SiSj=2n^2(n+1)^2$$

$$\sum_{i=1}^{2n}Si^2\geq \frac {n(n+1)^2} {2}$$

我们终于求出了$\sum_{i=1}^{2n}Si^2$的最小值

把它代入我们之前求出的式子里

$$\frac {\sum_{i=1}^{2n} {Si^2} - n(n+1)} {n(n+1)}\geq \frac {\frac {n(n+1)^2} {2}-n(n+1)} {n(n+1)}$$

$$=\frac {n+1} {2} -1=\frac {n-1} {2}$$

所以不管出题人怎么出数据,任选两个字符串共有元素的期望个数最小都是$\frac {n-1} {2}$,所以直接O(n^3/32)是很正确的

代码?我™手贱重启电脑清空了

【CSP模拟赛】Confess(数学 玄学)的更多相关文章

  1. CSP模拟赛游记

    时间:2019.10.5 考试时间:100分钟(连正式考试时间的一半还没有到)题目:由于某些原因不能公开. 由于第一次接触NOIinux系统所以连怎么建文件夹,调字体,如何编译都不知道,考试的前半小时 ...

  2. 【CSP模拟赛】starway(玄学建边 最小生成树)

    问題描述 小w伤心的走上了 Star way to heaven.   到天堂的道路是一个笛卡尔坐标系上一个n×m的长方形通道(顶点在(0,0))和(n,m)),小w从最左边任意一点进入,从右边任意一 ...

  3. 【CSP模拟赛】坏天平(数学&思维)

    蹭兄弟学校的题目做还不用自己出题的感觉是真的爽 题目描述 nodgd有一架快要坏掉的天平,这架天平右边的支架有问题,如果右边的总重量比左边多太多,天平就彻底坏掉了.现在nodgd手上有n种砝码,质量分 ...

  4. 【CSP模拟赛】方程(数学)

    题目描述 求关于x的方程:x1+x2+……xk=n的非负整数解的个数. 输入格式 仅一行,包含两个正整数n,k. 输出格式 一个整数,表示方程不同解的个数,这个数可能很大,你只需输出mod 20080 ...

  5. 【csp模拟赛5】限制 (restrict.cpp)--数学

    自己看吧: 爆搜代码: //春水初涨-春林初盛-春风十里-不如你 //----hzwer // 这是啥子题,读不懂-- //题意有问题 -- #include<iostream> #inc ...

  6. 20180606模拟赛T4——数学游戏

    数学游戏 题目描述: 小T又发脑残了,没错,她又要求奇怪的东西,这次她想知道[X,Y]之间整数有多少可以表示成K个不同的B的幂的和形势.如\(x,y,k,b=15,20,2,2\),则有: \[17= ...

  7. 【CSP模拟赛】Freda的迷宫(桥)

    题目描述 Freda是一个迷宫爱好者,她利用业余时间建造了许多迷宫.每个迷宫都是由若干房间和走廊构成的,每条走廊都连接着两个不同的房间,两个房间之间最多只有一条走廊直接相连,走廊都是双向通过.  黄昏 ...

  8. CSP模拟赛2游记

    这次由于有课迟到30min,了所以只考了70min. 调linux配置调了5min,只剩下65min了. T1:有点像标题统计,但要比他坑一点,而且我就被坑了,写了一个for(int i=1;i< ...

  9. 【CSP模拟赛】Freda的旗帜

    题目描述  要开运动会了,Freda承担起了制作全校旗帜的工作.旗帜的制作方法是这样的:Freda一共有C种颜色的布条,每种布条都有无数个,你可以认为这些布条的长.宽.厚都相等,只有颜色可能不同.每个 ...

随机推荐

  1. 单IP、网络、别名管道限速的设置

    单IP.网络.别名管道限速的设置 在设备上,控制网络限速最常用的是流量×××中的“限制器”,你可以把它理解为一个管道,20Mit/s的下载管道,那么下载的最大流量就不会超过20Mit/s.流量限制器结 ...

  2. vsftp 安装配置(被动模式)

    vi /etc/vsftpd/vsftpd.conf vsftp配置末尾添加 pasv_enable=YES pasv_min_port=10000 pasv_max_port=10030 防火墙端口 ...

  3. Unicode 字符和UTF编码的理解

    Unicode 编码的由来 我们都知道,计算机的内部全部是由二进制数字0, 1 组成的, 那么计算机就没有办法保存我们的文字, 这怎么行呢? 于是美国人就想了一个办法(计算机是由美国人发明的),也把文 ...

  4. 理解 Cookie,Session,Token 并结合 Redis 的使用

    Http 协议是一个无状态协议, 客户端每次发出请求, 请求之间是没有任何关系的.但是当多个浏览器同时访问同一服务时,服务器怎么区分来访者哪个是哪个呢? cookie.session.token 就是 ...

  5. scrapy随机切换user-agent

    使用github的 scrapy-fake-useragent 不用自己改源码继承自带的userAgent中间件  只需要安装后增加配置即可 https://github.com/alecxe/scr ...

  6. python基础-模块(全是理论,没有代码)

    模块 概念:一系列功能的结合体.相当于模块包着一堆函数与代码.本质上是py文件. 来源: python内置的模块----→ python解释器的模块 第三方的模块 -----→ 其他人编写提供的 自定 ...

  7. MySQL Replication--复制延迟02--exec_time测试

    复制延迟(Seconds_Behind_Master)测试 测试环境: MySQL 5.7.19 测试主从时间差: 检查主从系统时间差,同时在主库和从库执行SELECT NOW()语句: 主库:-- ...

  8. 基于glew,freeglut的imshow

    OpenGL显示图片,这篇博客使用glew + freeglut + gdal来实现imshow. 主要修改: 使用BGR而不是RGB,保持和opencv行为一致 纯C,去掉C++相关的 去掉GDAL ...

  9. Linux系统自动化安装之cobbler实现

    一.cobbler简介 cobbler是快速网络安装linux操作系统的服务,支持众多的Linux发行版本,如redhat|.fedora.centos.debian.ubuntu和suse,也可以支 ...

  10. Linux命令——getfacl、setfacl

    简介 ACL是Access Control List的缩写,传统的Linux权限只能针对一个用户.一个群组及非此群组的其他人设置权限而已,无法针对单一用户或个人来设计权限.ACL可以对权限进行更细致的 ...