2019牛客暑期多校训练营(第一场) - H - XOR - 线性基
https://ac.nowcoder.com/acm/contest/881/H
题意:
给定n个整数,求其中异或和为 \(0\) 的子集的大小的和。
题解思路:
首先转化为每个可以通过异或表示 \(0\) 的数贡献它参与的子集数。
思考的过程分两步。一开始不管三七二十一先对 \(n\) 个整数求一次线性基 \(B_1\) ,记其秩为 \(r\) 。
第一步:
先考虑线性基 \(B_1\) 外的数(假如有的话)产生的贡献。枚举每一个数,记这个数为 \(x\) ,除去这个数,线性基 \(B_1\) 外还有 \(n-r-1\) 个数,他们可以自由组合出 \(2^{n-r-1}\) 种子集,加上 \(x\) 本身之后,必定能被线性基 \(B_1\) 表示,故这些子集都会使得 \(x\) 发生贡献。这样的 \(x\) 共有 \(n-r\) 个。
第二步:
再考虑线性基 \(B_1\) 中的数产生的贡献。枚举每一个数,记这个数为 \(x\) ,这个数产生的贡献来源于他在第一步中被使用的次数。那怎么判断这个数被使用了多少次呢?把除去这个数 \(x\) 的 \(n-1\) 个数求出线性基 \(B_3\) 。假如 \(x\) 能被线性基 \(B_3\) 表示,此时 \(B_3\) 的秩必定为 \(r\) 。选中\(x\)之后,再从其他的 \(n-r-1\) 个数自由组合出子集,再从 \(B_3\)中选出对应的数配成 \(0\) 。否则这个 \(x\) 没有用。
第三步:
最后是快速得到 \(B_3\) 的办法,注意到 \(B_3\) 必定先通过除去 \(B_1\) 的 \(n-r\) 个数构造,可以把 \(n-r\) 先构造出 \(B_2\) ,然后再插入 \(B_1\) 内除x的 \(r-1\) 个数构成。
注意MN不能取64否则会爆掉ll。
qls说是hdu5544 http://acm.hdu.edu.cn/showproblem.php?pid=5544
题意:给定n个整数,求其中异或和为0的子集的大小的和。
题解思路:
首先转化为每个可以通过异或表示0的数贡献它参与的子集数。
对n个数进行一次消元,得到一组大小为\(r\)的线性基B1。
那么剩下的\(n-r\)个满足以下的规律:
设现在要计算的元素为X1,其他\(n-r-1\)个元素任意组合,总能得到一个整数,再加入X1也是一个整数,这个整数能被线性基B1表示,则这个X1可以对他们贡献总计\(2^{n-r-1}\)个次。
接下来算线性基B1里的数的贡献。
再对这\(n-r\)个数进行一次消元, 得到另一组线性基B2。枚举B1的一个线性基X2,对其他的r-1+B2个数消元,得到某个线性基B3。
csl:你必须钦定一个数必须选,如果这个数在线性基外面那没事,在里面你就得搞事情了。
意思是:这个数在线性基外面说明没有集合需要它的参与就可以自己组成异或和为0的。要是在线性基里面就说明存在这个数可以对某些集合产生贡献。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MN=63;
const int mod=1000000007;
ll num[100005];
struct LinearBase {
ll base[MN];
bool flag;//该线性基能否表示0
int cnt;
void Copy(LinearBase b) {
cnt=b.cnt;
flag=b.flag;
memcpy(base,b.base,sizeof(base));
}
void Clear() {
cnt=0;
flag=false;
memset(base,0,sizeof( base));
}
//尝试向线性基中插入一个值
void Insert(ll x) {
for(int i=MN; ~i; i--)
if(x&(1ll<<i))
if(!base[i]) {
base[i]=x;
cnt++;
return;
} else
x^=base[i];
flag=true;
}
//判断该线性基能否表示x
bool Check(ll x) {
for(int i=MN; ~i; i--)
if(x&(1ll<<i)) {
if(!base[i])
return false;
else
x^=base[i];
}
return true;
}
} B1,B2,B3;
ll qpow(ll x,int n) {
ll res=1;
while(n) {
if(n&1)
res=res*x%mod;
x=x*x%mod;
n>>=1;
}
return res;
}
ll ans,p2;
vector<int>B1ID;
int main() {
#ifdef Yinku
freopen("Yinku.in", "r", stdin);
//freopen("Yinku.out", "w", stdout);
#endif // Yinku
int n;
while(~scanf("%d", &n)) {
for(int i=1; i<=n; i++) {
scanf("%lld",&num[i]);
}
B1.Clear();
B2.Clear();
B1ID.clear();
for(int i=1; i<=n; i++) {
if(B1.Check(num[i])) {
B2.Insert(num[i]);
} else {
B1.Insert(num[i]);
B1ID.push_back(i);
}
}
ans=0;
if(n!=B1.cnt) {
p2=qpow(2,n-B1.cnt-1);
ans+=p2*(n-B1.cnt)%mod;
}
for(ll i:B1ID) {
B3.Copy(B2);
for(ll j:B1ID) {
if(i!=j) {
B3.Insert(num[j]);
}
}
if(B3.Check(num[i])) {
//num[i]能被其他数表示
ans=(ans+p2)%mod;
}
}
printf("%lld\n", ans);
}
}
2019牛客暑期多校训练营(第一场) - H - XOR - 线性基的更多相关文章
- 2019牛客暑期多校训练营(第二场) H-Second Large Rectangle(单调栈)
题意:给出由01组成的矩阵,求求全是1的次大子矩阵. 思路: 单调栈 全是1的最大子矩阵的变形,不能直接把所有的面积存起来然后排序取第二大的,因为次大子矩阵可能在最大子矩阵里面,比如: 1 0 0 1 ...
- 2019牛客暑期多校训练营(第九场) D Knapsack Cryptosystem
题目 题意: 给你n(最大36)个数,让你从这n个数里面找出来一些数,使这些数的和等于s(题目输入),用到的数输出1,没有用到的数输出0 例如:3 4 2 3 4 输出:0 0 1 题解: 认真想一 ...
- [题解] 2019牛客暑期多校第三场H题 Magic Line
题目链接:https://ac.nowcoder.com/acm/contest/883/H 题意:二维平面上有n个不同的点,构造一条直线把平面分成两个点数相同的部分. 题解:对这n个点以x为第一关键 ...
- 2019牛客暑期多校训练营(第五场)G - subsequeue 1 (一题我真的不会的题)
layout: post title: 2019牛客暑期多校训练营(第五场)G - subsequeue 1 (一题我真的不会的题) author: "luowentaoaa" c ...
- 2019牛客暑期多校训练营(第一场)A Equivalent Prefixes(单调栈/二分+分治)
链接:https://ac.nowcoder.com/acm/contest/881/A来源:牛客网 Two arrays u and v each with m distinct elements ...
- 2019牛客暑期多校训练营(第一场)A题【单调栈】(补题)
链接:https://ac.nowcoder.com/acm/contest/881/A来源:牛客网 题目描述 Two arrays u and v each with m distinct elem ...
- 2019牛客暑期多校训练营(第一场) B Integration (数学)
链接:https://ac.nowcoder.com/acm/contest/881/B 来源:牛客网 Integration 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 5242 ...
- 2019牛客暑期多校训练营(第一场) A Equivalent Prefixes ( st 表 + 二分+分治)
链接:https://ac.nowcoder.com/acm/contest/881/A 来源:牛客网 Equivalent Prefixes 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/ ...
- 2019牛客暑期多校训练营(第二场)F.Partition problem
链接:https://ac.nowcoder.com/acm/contest/882/F来源:牛客网 Given 2N people, you need to assign each of them ...
随机推荐
- head 显示文件头部内容
1. 命令功能 head 默认显示文件前10行内容. 2.语法格式 head option file 参数说明 参数 参数说明 -n 指定显示行数 -c 指定显示的字节数 -v 总是显示文件名的文件头 ...
- Tunnel connection failed: 407 Proxy Authentication Required
报错信息 : Retrying (Retry(total=4, connect=None, read=None, redirect=None, status=None)) after connecti ...
- [python 学习] IO操作之读写文件
一.读取全部文件: # -*- coding: utf-8 -*- f = open('qq_url.txt','r'); print f.read(); f.close(); 二.读取规定长度文件 ...
- Java中最基本的集合接口:初识Collection
Collection接口 Collection是最基本的集合接口,一个Collection代表一组Object,即Collection的元素(Elements). 一些 Collection允许相同的 ...
- python基本数据预处理语法函数(1)
numpy包: ####数组###########from numpy import * shape #获取维度 size #获取长度 arange(0,5,1) #生成数组函数,从0到5以1为间隔 ...
- 如何从word文档复制内容到富文本编辑器
在之前在工作中遇到在富文本编辑器中粘贴图片不能展示的问题,于是各种网上扒拉,终于找到解决方案,在这里感谢一下知乎中众大神以及TheViper. 通过知乎提供的思路找到粘贴的原理,通过TheViper找 ...
- POJ 1502 MPI MaeIstrom ( 裸最短路 || atoi系统函数 )
题意 : 给出 N 个点,各个点之间的路径长度用给出的下三角矩阵表示,上上角矩阵和下三角矩阵是一样的,主对角线的元素都是 0 代表自己到达自己不用花费,现在问你从 1 到 N 的最短路,矩阵的 x 代 ...
- Codeforces 846F - Random Query
原题链接:http://codeforces.com/contest/846/problem/F 题意:给一个数列,任意取区间[l, r],问区间内不同数字的个数的期望是多少. 思路: 对于第i个数a ...
- 彻底关闭Postprocess
即使场景中没有postprocess volume,场景中也会有默认的postprocess volume效果,如果需要彻底关闭postprocess, 可以使用命令: sg.PostProcessQ ...
- 学习笔记11 EF查询相当于sql 中的 where in
两种写法 1. int[] Ids={1,2,3} DBContainer db=new DBContainer(); var list=db.表明.where(a=>Ids.Contains( ...