[线性基 模板]

网上找到一篇写的非常棒的文章Orz, 看完之后相信你一定也可以学会并理解
点我

个人理解线性基的关键:

1. 线性基可以有多个,但是每个内部一定都是线性无关的,也就是说对于一个可以用线性基表示的x来说, 其组成方案一定是唯一的, 所以一个线性基的任意子集都不会构造出0.
2. 构造线性基的核心就是 : 对于一个基 $a,b$ , 我们可以用 $a$^$b$ 替换任意一个,且不会影响这个基本身. 有了这个我们才可以很巧妙地构造线性基,且可以在这个基础上按照一定规则进行一些操作可以让这个线性基的性质更加的好看

思考

1. 原集合能够亦或出0 ? 如果原集合的大小等于线性基的大小, 那么肯定不能够亦或出0, 因为线性基的任意子集都不会被亦或出0, 如果不等于, 那么原数组中多余的数都可以被线性基表示, 然后这两者亦或一定为0
2. 原集合的亦或第k小 ? 记ai为线性基中最高位的1在第i位的 向量, 可以按照i从大到小的顺序,用ai去异或aj (j>i).
   这样最后得到的重塑后向量组又多了一个优美的性质： 只有ai的第i位是1,其他的第i位都是0. 有了这个性质,就容易证明要求第k小的异或和,只要把k二进制拆分,线性基中第j位是1就异或第j个线性基中的向量。正确性,我们脑补一下, 因为这样 重塑之后, 假设线性基的大小为sz, 那么从最低位到最高位就可以完全控制最后亦或值的大小，那么他可以组合出方案数为(1 << sz) - 1 个，并且按照二进制枚举的方式越来越大．还是无法说明白，ＱＡＱ

HDU-3949 亦或第k小

Code

/***********************************************
Author        :lzs
Created Time  :2018年10月19日 星期五 10时47分05秒
File Name     :xjj.cpp
************************************************/

#include <bits/stdc++.h>
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <string>
#include <math.h>
#include <stack>
#include <vector>
#include <queue>
#include <set>
#include <map>
using namespace std;
#define rep(i, l, r) for(int i = l; i < r; i++)
#define per(i, r, l) for(int i = r; i >= l; i--)
#define dbgln(...) cerr<<"["<<#__VA_ARGS__":"<<(__VA_ARGS__)<<"]"<<"\n"
#define dbg(...) cerr<<"["<<#__VA_ARGS__":"<<(__VA_ARGS__)<<"]"


typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int>pii;

const int N = (int) 1e5 + 11;
const int M = (int) 1e6 + 11;
const int MOD = (int) 1e9 + 7;
const int INF = (int) 0x3f3f3f3f;
const ll INFF = (ll) 0x3f3f3f3f3f3f3f3f;
/*-----------------------------------------------------------*/

int n;
ll xjj[70], xj[70];
bool insert(ll a){ // 如果最后a为0,那么表示 a可以用线性基表示
	for(ll i = 62; i >= 0; i--){
		if(a >> i & 1) {
			if(!xjj[i]) {
				xjj[i] = a;
				break;
			}
			a ^= xjj[i];
		}
	}
	return a == 0;
}
ll getmx(ll x = 0){
	// 不传参数,表示求原集合中的亦或最大值
	for(int i = 62; i >= 0; i--){
		if(x < ( x ^ xjj[i])) x = x ^ xjj[i];
	}
	return x;
}
// 原集合中的亦或最小值 就是线性基中的最低位的值
int sz;
void rebuild(){
	sz = 0;
	for(int i = 62; i >= 0; i--){
		for(int j = i - 1; j >= 0; j--){
			if(xjj[i] >> j & 1) xjj[i] ^= xjj[j];
		}
	}
	for(int i = 0; i <= 62; i++) if(xjj[i]) xj[sz++] = xjj[i];
}
ll query_kth(ll k){
	if(k >= (1ll << sz)) return -1; ll ans = 0;
	for(int i = 62; i >= 0; i--){　// 顺序枚举和倒序枚举都不影响
		if(k >> i & 1) ans ^= xj[i];
	}
	return ans;
} 
int main(){
//	cout << (1ll << 62) <<"\n";	
	int T; scanf("%d", &T);
	for(int cas = 1; cas <= T; cas++){
		scanf("%d" ,&n);
		memset(xjj, 0, sizeof(xjj));
		for(int i = 1; i <= n; i++) {
			static ll t; scanf("%lld", &t);
			insert(t);
		}
		rebuild();
		int q; scanf("%d" ,&q);
		printf("Case #%d:\n", cas);
		while(q--){
			static ll k; scanf("%lld" ,&k);
			if(n != sz) k--; // 如果线性基的大小 等于 原数组大小,则0是不能够被亦或出来的
			printf("%lld\n", query_kth(k));
		}
	}
	return 0;
}