线性基
1.算法分析
线性基的本质是用一个数x的尽可能高位的1来代表x
线性基是向量空间的一组基,通常可以解决有关异或的一些题目。
通俗一点的讲法就是由一个集合构造出来的另一个集合,它有以下几个性质:
- 线性基的元素能相互异或得到原集合的元素的所有相互异或得到的值。
- 线性基是满足性质 1 的最小的集合。
- 线性基没有异或和为 0 的子集。
- 线性基中每个元素的异或方案唯一,也就是说,线性基中不同的异或组合异或出的数都是不一样的。
- 线性基中每个元素的二进制最高位互不相同
2.模板
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
struct Linear_Basis{
LL d[61],p[61]; // d记录原先的,p记录重构后的
int cnt, flag; // cnt记录重构后有多少个,flag记录是否出现过0
Linear_Basis() {
memset(d, 0 ,sizeof(d));
memset(p, 0, sizeof(p));
cnt = 0; // 当前线性基内元素个数
flag = 0; // 不存在0
}
// 构造线性基:逐个元素插入
bool insert(LL val) {
for (int i = 60;i >= 0; i--) {
if (val & (1ll << i)) {
// 判断当前i位能否插入
if (!d[i]) {
// 如果当前i位没插入过元素
d[i] = val; // 插入
return true; // 插入成功
}
val ^= d[i];
}
}
flag = 1; // 存在0
return false; // 插入失败
}
// 查询第最大异或和
LL query_max() {
LL ret = 0;
for (int i = 60; i >= 0; i--)
if ((ret ^ d[i]) > ret)
ret ^= d[i];
return ret;
}
// 查询最小异或和
LL query_min() {
for (int i = 0;i <= 60; i++)
if (d[i])
return d[i];
return 0;
}
// 重构,使之形成对角矩阵
void rebuild() {
for (int i = 60;i >= 0; i--) // 当前i位时,把后面的j位情况都和i位情况做异或
for (int j = i - 1;j >= 0; j--)
if (d[i] & (1ll << j))
d[i] ^= d[j];
for (int i = 0; i <= 60; i++) // 记录重构结果
if (d[i])
p[cnt++] = d[i];
}
// 查询第k小
LL kthquery(LL k){
LL ret = 0;
if (flag) {
// 存在0
k--;
if (!k) return (LL)0;
}
if (k >= (1ll << cnt)) // 只能由2^cnt - 1个
return -1;
for (int i = 60; i >= 0; i--)
if (k & (1LL << i)) // 每位的贡献为2^i
ret ^= p[i];
return ret;
}
};
// 合并线性基
Linear_Basis merge(const Linear_Basis &n1, const Linear_Basis &n2) {
Linear_Basis ret = n1;
for (int i = 60; i >= 0; i--)
if (n2.d[i])
ret.insert(n2.d[i]);
return ret;
}
int main(){
int n;
scanf("%d",&n);
Linear_Basis lb;
for(int i=0;i<n;i++){
LL a;
scanf("%lld",&a);
lb.insert(a);
}
printf("%lld\n",lb.query_max());
return 0;
}
3.例题
P4570 [BJWC2011]元素
题意: 给n个整数,每个整数带一个权值v,求一个权值和最大的线性基。
题解: 按权值v从大->小排序,依次插入线性基。要是能够插入,那么累加到答案
代码:
#include<bits/stdc++.h>
using namespace std;