UOJ #36 -【清华集训2014】玛里苟斯(线性基+暴搜)
看到 \(k\) 次方的期望可以很自然地想到利用低次方和维护高次方和的套路进行处理,不过。由于这里的 \(k\) 达到 \(5\),直接这么处理一来繁琐,二来会爆 long long
,因此考虑另辟蹊径。注意到答案 \(\le 2^{63}-1\),也就是说当 \(k\) 比较大时值域也不会太大。因此考虑对 \(k\) 分类讨论。
\(k=1\) 时考虑计算每一位的贡献,注意到对于一位 \(i\),如果存在某个 \(a_j\) 满足 \(a_j\) 的 \(2^i\) 位为 \(1\),那么这一位在子集异或和中为 \(1\) 的概率就是 \(\dfrac{1}{2}\)。否则这一位自然不可能异或出 \(0\) 来,对答案的贡献自然就是 \(0\),贡献随便加一下即可。
\(k=2\) 时考虑用低次方和维护高次方和的套路,假设 \(x=2^{b_1}+2^{b_2}+\cdots+2^{b_k}\),那么显然有 \(x^2=(2^{b_1}+2^{b_2}+\cdots+2^{b_k})^2=\sum\limits_{i=1}^k\sum\limits_{j=1}^k2^{b_i}·2^{b_j}\),也就是说我们考虑每两个在异或和中同时为 \(1\) 的位 \(i,j\),答案加上 \(2^i,2^j\) 在 \(x\) 中同时为 \(1\) 的概率乘上 \(2^{i+j}\) 的和。考虑枚举 \(i,j\) 计算 \(2^i,2^j\) 位在 \(x\) 中同时为 \(1\) 的概率,显然满足不存在任何一个 \(a_p\) 的第 \(i\) 位为 \(1\),或者不存在任何一个 \(a_p\) 的第 \(i\) 位为 \(1\),那么 \(2^i,2^j\) 位在 \(x\) 中同时出现的概率为 \(0\)。否则如果对于任意 \(a_p\) 都有 \(a_p\) 的第 \(i,j\) 位同时为 \(0\) 或者同时为 \(1\),那么 \(2^i,2^j\) 在 \(x\) 中同时出现的概率为 \(\dfrac{1}{2}\),否则 \(2^i,2^j\) 在 \(x\) 中同时出现的概率为 \(\dfrac{1}{4}\),贡献累加一下即可。
\(k\ge 3\) 时我们考虑建出原序列的线性基,显然对于不在线性基中的元素,我们删除它们后答案不会发生变化。由于题目保证了答案 \(\le 2^{63}-1\),因此线性基中的个数不会太多。精确地估算一下大概就 \(\sum\limits_{i=0}^{2^p-1}x^3\) 是关于 \(2^p\) 的四次多项式,再除以一个 \(2^p\) 是关于 \(2^p\) 的三次多项式,当 \(p=22\) 时 \(4194304^3\) 就已经超过了 ull
的范围,就算乘上个小常数也超过了 \(2^{63}\),因此暴搜即可。
可以证明此题答案小数点后最多只有一位,证明大概就仿照 \(k=2\) 的思路,考虑枚举 \(k\) 位 \(b_1,b_2,b_3,\cdots,b_k\),答案加上 \(b_1,b_2,\cdots,b_k\) 位同时为 \(1\) 的概率乘上 \(2^{b_1+b_2+\cdots+b_k}\)。我们假设 \(b_1,b_2,b_3,\cdots,b_k\) 中不同数的个数为 \(c\),那么这些位产生的贡献应是某个整数乘上 \(2^{b_1+b_2+\cdots+b_k-c}\),分析一下可知当 \(b_i=0,c=1\) 时 \(2^{b_1+b_2+\cdots+b_k-c}\) 取到最小值 \(0.5\),因此此题小数点后最多只有一位 0.5
。
const int MAXN=1e5;
const int LOG_N=64;
int n,k,flg[LOG_N+2];u64 a[MAXN+5];
void print(u64 x){
printf("%llu",x>>1);
if(x&1) printf(".5");
}
namespace sub{
u64 b[LOG_N+2];
__int128_t res=0;int tot=0;
void insert(u64 x){
for(int i=LOG_N-1;~i;i--) if(x>>i&1){
if(!b[i]) return b[i]=x,void();
x^=b[i];
}
}
void dfs(int x,u64 sm){
if(x==LOG_N+1){
__int128_t mul=1;
for(int i=1;i<=k;i++) mul*=sm;
res+=mul;tot++;return;
} dfs(x+1,sm);
if(b[x]) dfs(x+1,sm^b[x]);
}
void solve(){
for(int i=1;i<=n;i++) insert(a[i]);
dfs(0,0);
if(tot!=1) res/=(tot/2);
else res*=2;
print((u64)res);
}
}
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++) scanf("%llu",&a[i]);
if(k==1){
u64 res=0;
for(int i=0;i<LOG_N;i++){
bool has1=0;
for(int j=1;j<=n;j++) has1|=(a[j]>>i&1);
if(has1) res+=1ull<<i;
} print(res);
} else if(k==2){
u64 res=0;
for(int i=0;i<LOG_N;i++) for(int j=1;j<=n;j++) flg[i]|=(a[j]>>i&1);
for(int i=0;i<LOG_N;i++) for(int j=0;j<LOG_N;j++){
bool sm=1;
if(!flg[i]||!flg[j]) continue;
for(int k=1;k<=n;k++) sm&=(!((a[k]>>i&1)^(a[k]>>j&1)));
if(sm) res+=1ull<<(i+j);
else res+=1ull<<(i+j-1);
} print(res);
} else sub::solve();
return 0;
}
UOJ #36 -【清华集训2014】玛里苟斯(线性基+暴搜)的更多相关文章
- uoj #46[清华集训2014]玄学
uoj 因为询问是关于一段连续区间内的操作的,所以对操作构建线段树,这里每个点维护若干个不交的区间,每个区间\((l,r,a,b)\)表示区间\([l,r]\)内的数要变成\(ax+b\) 每次把新操 ...
- UOJ.41.[清华集训2014]矩阵变换(稳定婚姻)
题目链接 稳定婚姻问题:有n个男生n个女生,每个男/女生对每个女/男生有一个不同的喜爱程度.给每个人选择配偶. 若不存在 x,y未匹配,且x喜欢y胜过喜欢x当前的配偶,y喜欢x也胜过y当前的配偶 的完 ...
- bzoj 3816&&uoj #41. [清华集训2014]矩阵变换
稳定婚姻问题: 有n个男生,n个女生,所有女生在每个男生眼里有个排名,反之一样. 将男生和女生两两配对,保证不会出现婚姻不稳定的问题. 即A-1,B-2 而A更喜欢2,2更喜欢A. 算法流程: 每次男 ...
- uoj 41 【清华集训2014】矩阵变换 婚姻稳定问题
[清华集训2014]矩阵变换 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://uoj.ac/problem/41 Description 给出 ...
- AC日记——【清华集训2014】奇数国 uoj 38
#38. [清华集训2014]奇数国 思路: 题目中的number与product不想冲: 即为number与product互素: 所以,求phi(product)即可: 除一个数等同于在模的意义下乘 ...
- [UOJ#274][清华集训2016]温暖会指引我们前行
[UOJ#274][清华集训2016]温暖会指引我们前行 试题描述 寒冬又一次肆虐了北国大地 无情的北风穿透了人们御寒的衣物 可怜虫们在冬夜中发出无助的哀嚎 “冻死宝宝了!” 这时 远处的天边出现了一 ...
- UOJ#36. 【清华集训2014】玛里苟斯 线性基
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ36.html 题解 按照 $k$ 分类讨论: k=1 : 我们考虑每一位的贡献.若有至少一个数第 $i$ ...
- uoj#36. 【清华集训2014】玛里苟斯(线性基+概率期望)
传送门 为啥在我看来完全不知道为什么的在大佬们看来全都是显然-- 考虑\(k=1\)的情况,如果序列中有某一个\(a_j\)的第\(i\)位为\(1\),那么\(x\)的第\(i\)位为\(1\)的概 ...
- UOJ #36「清华集训2014」玛里苟斯
这怎么想得到啊......... UOJ #36 题意:求随机一个集合的子集的异或和的$k$次方的期望值,保证答案$ \lt 2^{63},1 \leq k \leq 5$ $ Solution:$ ...
随机推荐
- 对epoll机制的学习理解v1
epoll机制 wrk用非阻塞多路复用IO技术创造出大量的连接,从而达到很好的压力测试效果.epoll就是实现IO多路复用的关键. 本节是对epoll的本质的学习总结,进一步的参考资料为: <深 ...
- D:\Software\Keil5\ARM\PACK\Keil\STM32F1xx_DFP\2.1.0\Device\Include\stm32f10x.h(483): error: #5: cannot open source input file "core_cm3.h": No such file or directory
1. 错误提示信息: D:\Software\Keil5\ARM\PACK\Keil\STM32F1xx_DFP\2.1.0\Device\Include\stm32f10x.h(483): erro ...
- NKOJ-2936 城市建设
问题描述: PS国是一个拥有诸多城市的大国,国王Louis为城市的交通建设可谓绞尽脑汁.Louis可以在某些城市之间修建道路,在不同的城市之间修建道路需要不同的花费.Louis希望建造最少的道路使得国 ...
- 《基于SD-SEIR模型的实验室人员不安全行为传播研究》
My Focus:基于SD-SEIR模型的实验室人员不安全行为的传播; 建模与实验仿真 Title: Study on Porpagation of Unsafe Bhavior of Laborat ...
- 洛谷 P6075 [JSOI2015]子集选取
链接:P6075 前言: 虽然其他大佬们的走分界线的方法比我巧妙多了,但还是提供一种思路. 题意: %&¥--@#直接看题面理解罢. 分析过程: 看到这样的题面我脑里第一反应就是DP,但是看到 ...
- 字符串与模式匹配算法(一):BF算法
一.BF算法的基本思想 BF(Brute Force)算法是模式匹配中最简单.最直观的算法.该算法最基本的思想是从主串的第 start 个字符起和模式P(要检索的子串)的第1个字符比较,如果相等,则逐 ...
- AOP源码解析:AspectJAwareAdvisorAutoProxyCreator类的介绍
AspectJAwareAdvisorAutoProxyCreator 的类图 上图中一些 类/接口 的介绍: AspectJAwareAdvisorAutoProxyCreator : 公开了Asp ...
- Java 将Excel转为et和ett格式
以.et结尾的文件格式是属于金山办公软件WPS Office中的电子表格文件,.ett是一种模板文件格式.除了通过WPS软件可以创建该格式的电子表格外,也可以通过格式转换的方法来获得,如将Micros ...
- AXI总线简介、ID分析、DMA、Vivado烧录、系统集成
转载:https://blog.csdn.net/CrazyUncle/article/details/89918030?depth_1-utm_source=distribute.pc_releva ...
- Redis 客户端重试指南
本作品采用知识共享署名-非商业性使用 4.0 国际许可协议进行许可. 在互联网服务中,特别是在云环境下,网络及硬件环境复杂,所有应用程序都可能遇到暂时性故障.暂时性故障包括瞬时的网络抖动,服务暂时不可 ...