线性基

线性基是什么？

一个原序列能求出至少一组线性基
线性基是一个集合
第一，原序列中的任意一个数都可以通过线性基中的一些数异或得到
第二，线性基中的个数是一定的（取决于原序列，下面会讲）它不会重复
即线性基中的某些数异或的值不会存在线性基中，或者说线性基中的一些数异或起来不会等于0

线性基又有什么作用呢？为什么要求线性基？

用于减少时间复杂度
线性基的长度大约为log（n）（n为原序列的长度）
原序列和它的线性基在异或操作这个层面上是等价的
所以用线性基来代替它的原序列来处理就可以了

线性基怎么求？

d数组代表原序列的线性基，x代表原序列中的某个元素
d[i]表示x在二进制下的最高位的位置位i+1（i从0开始）
有人会问了
如果有2个数的x和y的最高位的位置都是i+1怎么办？
首先d[i]值能存一个值，那么d[i]先存x，那么可以把y变成x^y， 因为x^ (x^y)=y，再存y，如果这时的y的最高位和z的是最高位的位置一样，且z已经被存进去了怎么办？那把再把y变成 y^z，直到y能存进去为止
那如果把整个d数组都遍历完了都存不进去怎么办？
说明y不能存进去，也就是说如果把y存进去，线性基会出现一些数异或得到0的情况，就会重复，即这个y可以通过已有的线性基中的一些数异或得到
下面是代码

void insert(ll x){
	for(int i=62;i>=0;i--)
	if(x&(1ll<<i))
	{
		if(d[i])x^=d[i];
		else
		{
			d[i]=x;
			return;
		}
	}
	f=1;//标记有没有不能存进去的数
}

下面是关于线性基的一些应用
求在一个序列中，取若干个数，使得它们的异或和最大

ll getmax(){
	ll ans=0;
	for(int i=62;i>=0;i--)
	if(ans^d[i]>ans)ans^=d[i];
	return ans;
}

求最小

ll getmin(){
	if(f)return 0;
	for(int i=0;i<=62;i++)
	if(d[i])return d[i];
}

求第k小的值
从一个序列中取任意个元素进行异或，求能异或出的所有数字中第k小的那个
我们先来看一下
a1=1 0 1 0 0， a2=1 1 1 0 0
----b=1 1 1 1, -------b=1 1 1 1
a1异或b会得到比a1大的数，而a2异或b会得到比a2小的数
只要a的第4位（b的最高位的位置为4）为0那么a异或b的值一定大于a
那么我们需要重组一个线性基使他们变成这样

那么得到的d数组中的任意2个数异或都会变大
d[0]<d[1]<d[2]<d[3]<…
最小的数为d[0];
第二小的数为d[1];
第三小的数为d[0]^d[1];
第四小的数为d[2];
第五小的数为d[0]^d[2]
以此类推能得到如果k=10110；
那么第k小的数为d[1]^d[2] ^d[4];
代码

void rebuild(){
	for(int i=62;i>=1;i--)
	for(int j=i-1;j>=0;j--)
	if(d[i]&(1ll<<j))d[i]^=d[j];
	for(int i=0;i<=62;i++)
	if(d[i])p[cnt++]=d[i];
}
ll getkth(ll k){
	if(f)--k;
	if(!k)return 0;
	if(k>=(1ll<<cnt))return -1;
	ll ans=0;
	for(int i=0;i<cnt;i++)
	if(k&(1ll<<i))ans^=p[i];
	return ans;
}

P4570 [BJWC2011]元素

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll d[100],sum;
int t;
struct node{
	ll n,m;
}a[2000];
bool cmp(node x,node y)
{
	return x.m>y.m;
}
bool insert(ll x)
{
	for(int i=62;i>=0;i--)
	if(x&(1ll<<i))
	{
		if(d[i])x^=d[i];
		else{d[i]=x;return 1;}
	}
	return 0;
}
int main()
{
	scanf("%d",&t);
	for(int i=0;i<t;i++)scanf("%lld%lld",&a[i].n,&a[i].m);
	sort(a,a+t,cmp);
	for(int i=0;i<t;i++)if(insert(a[i].n))sum+=a[i].m;
	cout<<sum<<endl;
	return 0;
}

P3857 [TJOI2008]彩灯

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=2008;
ll d[200];
int n,m,cnt;
string s;
void insert(ll x)
{
	for(int i=62;i>=0;i--)
	{
		if(x&(1ll<<i))
		{
			if(d[i])x^=d[i];
			else{d[i]=x;return;}
		}
	}
}
void rebuild()
{
	for(int i=62;i>=1;i--)
	for(int j=i-1;j>=0;j--)
	if(d[i]&(1ll<<j))d[i]^=d[j];
	for(int i=0;i<=62;i++)
	if(d[i])cnt++;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=0;i<m;i++)
	{
		cin>>s;
		ll a=0;
		for(int j=0;j<n;j++)
		{
			a<<=1;
			if(s[j]=='O')a++;
		}
		insert(a);
	}
	rebuild();
	cout<<(1ll<<cnt)%mod<<endl;
	return 0;
}

nowcoder牛牛数数
用个二分求这个数是不是大于k，找到第一个大于k的位置
再减一下就能得到答案

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+10;
ll d[100],p[100],k,a;
int n,cnt;
bool f=0;
void insert(ll x){
	for(int i=62;i>=0;i--)
	if(x&(1ll<<i))
	{
		if(d[i])x^=d[i];
		else{d[i]=x;return;}
	}
	f=1;
}
void rebuild(){
	for(int i=62;i>=1;i--)
	for(int j=i-1;j>=0;j--)
	if(d[i]&(1ll<<j))d[i]^=d[j];
	for(int i=0;i<=62;i++)
	if(d[i])p[cnt++]=d[i];
}
ll getkth(ll k){
	if(f)--k;
	if(!k)return 0;
	if(k>=(1ll<<cnt))return -1;
	ll ans=0;
	for(int i=0;i<cnt;i++)
	if(k&(1ll<<i))ans^=p[i];
	return ans;
}
int main()
{
	scanf("%d%lld",&n,&k);
	for(int i=1;i<=n;i++)scanf("%lld",&a),insert(a);
	rebuild();
	ll l=1,r=(1ll<<cnt)+f-1;
	while(l<r)
	{
		ll m=l+r>>1;
		if(getkth(m)>k)r=m;
		else l=m+1;
	}
	cout<<(1ll<<cnt)-l+f<<endl;//由于线性基得不到0所以f用来表示0存不存在
	return 0;
}