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 - 线性基的更多相关文章

  1. 2019牛客暑期多校训练营(第二场) H-Second Large Rectangle(单调栈)

    题意:给出由01组成的矩阵,求求全是1的次大子矩阵. 思路: 单调栈 全是1的最大子矩阵的变形,不能直接把所有的面积存起来然后排序取第二大的,因为次大子矩阵可能在最大子矩阵里面,比如: 1 0 0 1 ...

  2. 2019牛客暑期多校训练营(第九场) D Knapsack Cryptosystem

    题目 题意: 给你n(最大36)个数,让你从这n个数里面找出来一些数,使这些数的和等于s(题目输入),用到的数输出1,没有用到的数输出0 例如:3  4 2 3 4 输出:0 0 1 题解: 认真想一 ...

  3. [题解] 2019牛客暑期多校第三场H题 Magic Line

    题目链接:https://ac.nowcoder.com/acm/contest/883/H 题意:二维平面上有n个不同的点,构造一条直线把平面分成两个点数相同的部分. 题解:对这n个点以x为第一关键 ...

  4. 2019牛客暑期多校训练营(第五场)G - subsequeue 1 (一题我真的不会的题)

    layout: post title: 2019牛客暑期多校训练营(第五场)G - subsequeue 1 (一题我真的不会的题) author: "luowentaoaa" c ...

  5. 2019牛客暑期多校训练营(第一场)A Equivalent Prefixes(单调栈/二分+分治)

    链接:https://ac.nowcoder.com/acm/contest/881/A来源:牛客网 Two arrays u and v each with m distinct elements ...

  6. 2019牛客暑期多校训练营(第一场)A题【单调栈】(补题)

    链接:https://ac.nowcoder.com/acm/contest/881/A来源:牛客网 题目描述 Two arrays u and v each with m distinct elem ...

  7. 2019牛客暑期多校训练营(第一场) B Integration (数学)

    链接:https://ac.nowcoder.com/acm/contest/881/B 来源:牛客网 Integration 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 5242 ...

  8. 2019牛客暑期多校训练营(第一场) A Equivalent Prefixes ( st 表 + 二分+分治)

    链接:https://ac.nowcoder.com/acm/contest/881/A 来源:牛客网 Equivalent Prefixes 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/ ...

  9. 2019牛客暑期多校训练营(第二场)F.Partition problem

    链接:https://ac.nowcoder.com/acm/contest/882/F来源:牛客网 Given 2N people, you need to assign each of them ...

随机推荐

  1. Maven 集成Tomcat插件(引用)

    Maven已经是Java的项目管理标配,如何在JavaEE开发使用Maven调用Web应用,是很多同学关心的问题.本文将介绍,Maven如何介绍Tomcat插件. Maven Tomcat插件现在主要 ...

  2. 总结const、readonly、static三者的区别【收藏、转载】20190614

    总结const.readonly.static三者的区别 const:静态常量,也称编译时常量(compile-time constants),属于类型级,通过类名直接访问,被所有对象共享! a.叫编 ...

  3. 2018-11-19-visualStudio-无法登陆

    title author date CreateTime categories visualStudio 无法登陆 lindexi 2018-11-19 15:24:15 +0800 2018-2-1 ...

  4. Rust(一)介绍 安装

    目录 Rust安装 Rust介绍: Windows 安装步骤: Helle world 创建项目文件夹: 写并执行程序: Rust安装 安装过程简单快捷,直接参照官网即可,Rust安装 Rust介绍: ...

  5. 2018-10-01-weekly

    Algorithm 77. 组合 What 给定两个整数 n 和 k,返回 1 ... n 中所有可能的 k 个数的组合. How 利用递归的思想,当凑够k个数时,就回退回去,remove掉一个数,在 ...

  6. php $_SERVER 中的 QUERY_STRING和REQUEST_URI

    index.php <?php print_r($_GET); parse_str($_SERVER['QUERY_STRING'],$get); print_r($get); print_r( ...

  7. Wannafly挑战赛27 D绿魔法师

    链接Wannafly挑战赛27 D绿魔法师 一个空的可重集合\(S\),\(n\)次操作,每次操作给出\(x,k,p\),要求支持下列操作: 1.在\(S\)中加入\(x\). 2.求\[\sum_{ ...

  8. man uname

    UNAME(1)         User Commands/用户命令        UNAME(1) NAME/名称       uname - print system information/打 ...

  9. 在一个div上增加遮罩

    有一个需求是给一个视频增加遮罩 我研究了下 抽象出来就是给一个div增加遮罩 原理是:最外层的div使用relative定位     然后里面放两个子div  一个是不被遮的  另一个是遮罩(用abs ...

  10. 2019 上海网络赛 F Rhyme scheme (字典树DP)

    题目:https://nanti.jisuanke.com/t/41414 题意:求长度为n的第k个bell number  ,  就是第i位的选取范围在 1-(i-1)位的最大值 +1,第一位固定为 ...