【BZOJ3811】玛里苟斯(线性基)

题面

BZOJ

题解

\(K=1\)很容易吧,拆位考虑贡献,所有存在的位出现的概率都是\(0.5\),所以答案就是所有数或起来的结果除二。

\(K=2\)的情况,我们直接拆开平方式,平方项的贡献直接算,出现的概率还是\(0.5\),然后\(2ab\)这样子的东西出现的概率是\(0.5*0.5=0.25\),然而注意到有一些位直接两两之间存在联系,即选择了第\(i\)位的时候必定会同时选择第\(j\)位,那么此时两位之间的概率就是\(0.5\),这一部分要特殊判断一下。

接下来考虑\(K\ge 3\)的情况。因为答案不会超过\(2^{63}\),所以任何一个数字都不会从超过\(2^{22}\)。我们知道如果一个数可以被线性基中的其他所有数给表示出来,那么这个数出现的概率一定是\(1/{2^{|G|}}\),其中\(|G|\)是线性基内的元素个数。既然数字大小不大,那么线性基中的元素个数也不多,所以可以构建线性基之后直接\(2^{|G|}\)枚举所有可能出现的数字,\(K\)次方之后求和即可。

注意一下,最终答案除掉总的数字个数之后在\(2^{63}\)以内,那么在除之前是可能爆掉的,所以这里可以拆位什么的存一下就好了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
#define ll long long
#define ull unsigned long long
inline ull read()
{
ull x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,K;ull a[100100];
namespace Task1
{
void Solve()
{
ull ans=0;
for(int i=1;i<=n;++i)ans|=a[i];
if(ans&1)printf("%llu.5\n",ans>>1);
else printf("%llu\n",ans>>1);
}
}
namespace Task2
{
bool vis[50];
bool check(int x,int y)
{
if(!vis[x]||!vis[y])return false;
for(int i=1;i<=n;++i)
{
if((a[i]&(1ll<<x))&&!(a[i]&(1ll<<y)))return false;
if((a[i]&(1ll<<y))&&!(a[i]&(1ll<<x)))return false;
}
return true;
}
void Solve()
{
ull ans=0;
for(int i=0;i<33;++i)
for(int j=1;j<=n;++j)
if(a[j]&(1ll<<i)){vis[i]=true;break;}
for(int i=0;i<33;++i)if(vis[i])ans+=1ll<<(i+i);
for(int i=0;i<33;++i)
for(int j=i+1;j<33;++j)
if(vis[i]&&vis[j])ans+=1ll<<(i+j);
for(int i=0;i<33;++i)
for(int j=i+1;j<33;++j)
if(check(i,j))ans+=1ll<<(i+j);
if(ans&1)printf("%llu.5\n",ans>>1);
else printf("%llu\n",ans>>1);
}
}
namespace Task3
{
int p[50];
void insert(int x)
{
for(int i=30;~i;--i)
if(x&(1<<i))
{
if(!p[i]){p[i]=x;return;}
x^=p[i];
}
}
int S[55],top;
ull ans=0,r=0;
void Calc(int val)
{
int MOD=1<<top;
ull D=0,R=1;
for(int i=1;i<=K;++i)
{
D*=val;R*=val;
D+=R/MOD;R%=MOD;
}
ans+=D;r+=R;
ans+=r/MOD;r%=MOD;
}
void dfs(int x,int val)
{
if(x==top+1){Calc(val);return;}
dfs(x+1,val);dfs(x+1,val^S[x]);
}
void Solve()
{
for(int i=1;i<=n;++i)insert(a[i]);
for(int i=0;i<=30;++i)if(p[i])S[++top]=p[i];
dfs(1,0);
if(r)printf("%llu.5\n",ans);
else printf("%llu\n",ans);
return;
}
}
int main()
{
n=read(),K=read();
for(int i=1;i<=n;++i)a[i]=read();
if(K==1){Task1::Solve();return 0;}
if(K==2){Task2::Solve();return 0;}
if(K>=3){Task3::Solve();return 0;}
return 0;
}

【BZOJ3811】玛里苟斯(线性基)的更多相关文章

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

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

  2. BZOJ.3811.玛里苟斯(线性基)

    BZOJ UOJ 感觉网上大部分题解对我这种数学基础差的人来说十分不友好...(虽然理解后也觉得没有那么难) 结合两篇写的比较好的详细写一写.如果有错要指出啊QAQ https://blog.csdn ...

  3. [清华集训2015 Day1]玛里苟斯-[线性基]

    Description Solution 考虑k=1的情况.假设所有数中,第i位为1的数的个数为x,则最后所有的子集异或结果中,第i位为1的个数为$(C_{k}^{1}+C_{k}^{3}+...)$ ...

  4. BZOJ3811 玛里苟斯(线性基+概率期望)

    k=1的话非常好做,每个有1的位都有一半可能性提供贡献.由组合数的一些性质非常容易证明. k=2的话,平方的式子展开可以发现要计算的是每一对位提供的贡献,于是需要计算每一对位被同时选中的概率.找出所有 ...

  5. 【BZOJ 3811】玛里苟斯 大力观察+期望概率dp+线性基

    大力观察:I.从输出精准位数的约束来观察,一定会有猫腻,然后仔细想一想,就会发现输出的时候小数点后面不是.5就是没有 II.从最后答案小于2^63可以看出当k大于等于3的时候就可以直接搜索了 期望概率 ...

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

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

  7. bzoj 3811: 玛里苟斯【线性基+期望dp】

    这个输出可是有点恶心啊--WA*inf,最后抄了别人的输出方法orz 还有注意会爆long long,要开unsigned long long 对于k==1,单独考虑每一位i,如果这一位为1则有0.5 ...

  8. UOJ #36 -【清华集训2014】玛里苟斯(线性基+暴搜)

    UOJ 题面传送门 看到 \(k\) 次方的期望可以很自然地想到利用低次方和维护高次方和的套路进行处理,不过.由于这里的 \(k\) 达到 \(5\),直接这么处理一来繁琐,二来会爆 long lon ...

  9. 洛谷P3389 高斯消元 / 高斯消元+线性基学习笔记

    高斯消元 其实开始只是想搞下线性基,,,后来发现线性基和高斯消元的关系挺密切就一块儿在这儿写了好了QwQ 先港高斯消元趴? 这个算法并不难理解啊?就会矩阵运算就过去了鸭,,, 算了都专门为此写个题解还 ...

随机推荐

  1. Caffe源码中caffe.proto文件分析

    Caffe源码(caffe version:09868ac , date: 2015.08.15)中有一些重要文件,这里介绍下caffe.proto文件. 在src/caffe/proto目录下有一个 ...

  2. Windows 10 中 VMware 要求禁用 Device Guard 问题

    今天在打开虚拟机的时候,突然出现下面这个错误.网上给了很多教程,基本上都是禁用 Device Guard 和关闭 Hyper-v,博主按照其方法操作,依旧出现下面错误.后来经过不懈努力,终于找到解决办 ...

  3. react/React Native 在 import 导入时,有的带花括号{},有的不带原理解析

    在使用import引用模块时,如何正确使用{} 例如:有两个文件,home.js.user.js 一:不使用{}: 当需要在home.js中引入user.js的时候 //home.js 文件中impo ...

  4. 作业20171102 alpha-review 成绩

    申诉 对成绩有疑问或不同意见的同学,请在群里[@杨贵福]. 申诉时间截止2017年12月12日 17:00. 成绩 review NABCD-评论 SPEC-评论 例行报告 附加分数 合计 本周归一化 ...

  5. 猫咪记单词Beta版使用说明

    猫咪记单词Beta版使用说明 一.项目背景 英语四级考试.六级考试.托福.雅思等英语方面的考试是现在大学生必须面对的问题.同时因为学生对手机的使用越来越频繁,而且仅仅通过书本背诵单词又比较无聊坚持的时 ...

  6. IT行业的创新的读后感

    一.什么是创新 创新是以新思维.新发明和新描述为特征的一种概念化过程.它原意有三层含义,第一,更新:第二,创造新的东西:第三,改变.创新是人类特有的认识能力和实践能力,是人类主观能动性的高级表现形式, ...

  7. leetcode 730 Count Different Palindromic Subsequences

    题目链接: https://leetcode.com/problems/count-different-palindromic-subsequences/description/ 730.Count ...

  8. 四则运算安卓版ver.mk3

    在原有的基础上做了些许改动以及添加了一点小功能,以下是代码: package com.example.add; import java.io.File; import com.example.add. ...

  9. node基础 npm、module、exports、require

    module 模块.包:可以认为是一个代码包,package,提供特定的功能(暴露给外界接口,让外界调用) exports 输出.导出:导出模块中的各种类型的变量,以及各种方法,导出之后,才可以被外界 ...

  10. IIS7开启gZip动态压缩

    1.安装动态压缩模块: 安装过程可能的报错:This application has requested the Runtime to terminate it in an unusual way. ...