UOJ 题面传送门

看到 \(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】玛里苟斯(线性基+暴搜)的更多相关文章

  1. uoj #46[清华集训2014]玄学

    uoj 因为询问是关于一段连续区间内的操作的,所以对操作构建线段树,这里每个点维护若干个不交的区间,每个区间\((l,r,a,b)\)表示区间\([l,r]\)内的数要变成\(ax+b\) 每次把新操 ...

  2. UOJ.41.[清华集训2014]矩阵变换(稳定婚姻)

    题目链接 稳定婚姻问题:有n个男生n个女生,每个男/女生对每个女/男生有一个不同的喜爱程度.给每个人选择配偶. 若不存在 x,y未匹配,且x喜欢y胜过喜欢x当前的配偶,y喜欢x也胜过y当前的配偶 的完 ...

  3. bzoj 3816&&uoj #41. [清华集训2014]矩阵变换

    稳定婚姻问题: 有n个男生,n个女生,所有女生在每个男生眼里有个排名,反之一样. 将男生和女生两两配对,保证不会出现婚姻不稳定的问题. 即A-1,B-2 而A更喜欢2,2更喜欢A. 算法流程: 每次男 ...

  4. uoj 41 【清华集训2014】矩阵变换 婚姻稳定问题

    [清华集训2014]矩阵变换 Time Limit: 20 Sec  Memory Limit: 256 MB 题目连接 http://uoj.ac/problem/41 Description 给出 ...

  5. AC日记——【清华集训2014】奇数国 uoj 38

    #38. [清华集训2014]奇数国 思路: 题目中的number与product不想冲: 即为number与product互素: 所以,求phi(product)即可: 除一个数等同于在模的意义下乘 ...

  6. [UOJ#274][清华集训2016]温暖会指引我们前行

    [UOJ#274][清华集训2016]温暖会指引我们前行 试题描述 寒冬又一次肆虐了北国大地 无情的北风穿透了人们御寒的衣物 可怜虫们在冬夜中发出无助的哀嚎 “冻死宝宝了!” 这时 远处的天边出现了一 ...

  7. UOJ#36. 【清华集训2014】玛里苟斯 线性基

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ36.html 题解 按照 $k$ 分类讨论: k=1 : 我们考虑每一位的贡献.若有至少一个数第 $i$ ...

  8. uoj#36. 【清华集训2014】玛里苟斯(线性基+概率期望)

    传送门 为啥在我看来完全不知道为什么的在大佬们看来全都是显然-- 考虑\(k=1\)的情况,如果序列中有某一个\(a_j\)的第\(i\)位为\(1\),那么\(x\)的第\(i\)位为\(1\)的概 ...

  9. UOJ #36「清华集训2014」玛里苟斯

    这怎么想得到啊......... UOJ #36 题意:求随机一个集合的子集的异或和的$k$次方的期望值,保证答案$ \lt 2^{63},1 \leq k \leq 5$ $ Solution:$ ...

随机推荐

  1. spring cloud hystrix的隔离策略和dashboard

    随着服务的拆分,各个服务有着明确的职责,服务之间通过轻量级的协议进行通讯.但有时候我们完成一个功能需要同时调用多个微服务,比如完成订单的创建,那么获取用户信息需要调用用户微服务,获取商品信息需要调用商 ...

  2. postman_参数关联

    注:postman中的参数引用符号是{{变量名}},两个大括号 1.创建登录接口请求,在Tests中,添加如下代码,并send: // 获取响应的json数据 var jsonData = pm.re ...

  3. Zabbix 5.0:监控阿里云RDS

    Blog:博客园 个人 由于近期压测,需要频繁登录阿里云查看RDS监控,每次登录查看监控步骤较为繁琐,故将监控接入到zabbix. 概述 由于阿里云已做了RDS的监控,我们只需要通过阿里云SDK把这些 ...

  4. [luogu2973]driving out the piggies 驱逐猪猡【高斯消元+概率DP】

    看到题面的那一刻,我是绝望的ORZ 图论加概率期望加好像不沾边的高斯消元???我人直接傻掉 还没学过概率期望的我果断向题解屈服了(然后还是傻掉了两节课来找线性方程.. Description 奶牛们建 ...

  5. Linux文件IO操作

    来源:微信公众号「编程学习基地」 目录 文件操作 Linux文件类型 Linux文件权限 修改文件权限 Linux error 获取系统调用时的错误描述 打印错误信息 系统IO函数 open/clos ...

  6. Vivado Synth/Place Faild但是没有给出error信息

    最近遇到一个现象,以前可以编译通过的工程,修改之后发现Synthesis编译报错,而且没有给出error信息,以前也出现过无故place 失败但是没有给出error信息的现象,查看错误日志输出文件,出 ...

  7. notepad++ 替换回车换行

    以" | "为分隔符,换行 结果如下图:

  8. oracle 归档日志:db_recovery_file_dest、log_archive_dest和log_archive_dest_n的区别和使用

    概念: db_recovery_file_dest:默认的指定闪回恢复区路径 log_archive_dest:指定归档文件存放的路径,所有归档路径必须是本地的,默认为''.log_archive_d ...

  9. LOTO实践【干货】电压比较器的快速应用

    LOTO实践[干货]电压比较器的快速应用 话不多说先上图: 最近出差去客户那边做设备调试,现场如上图,其中我负责的技术部分包含了一个掉落物监测的功能.硬件上的原理比较简单,发射板上的红外发射头阵列,以 ...

  10. Piakchu之RCE漏洞

    一.Ping(远程系统命令执行) 首先正常输入一个ip,查看页面的返回值.发现有乱码,但是能看出执行了ping命令. 查看源代码,可以看到只是对操作系统进行了判断,而对输入内容是否为ip地址并没有判断 ...