【模板】线性基

TIPS

1ll警告

线性基插入

void insert(int rt,ll x)
{
    int i;
    fd(i,63,0)
    if (x & (1ll<<i))
        if (basis[rt][i] == 0)
        {
            basis[rt][i] = x;
            return;
        }   else x = x ^ basis[rt][i];
}

线性基查询

bool check(int rt,ll x)
{
    int i;
    fd(i,63,0)
    if (x & (1ll << i))
        if (basis[rt][i])
            x = x ^ basis[rt][i];
        else return false;
    return true;
}

线性基求交

假设有两个线性基$A,B$，设一个新的线性基$X$，一开始$X=A$，然后不断插入$B$的元素。对于$B_i$，如果它能被$X$线性表示，那么将“$A$中用来组成$B_i$的元素”加入答案线性基，否则将这个线性基加入$X$。（也就是说$X$中的元素不但要记录值，还要记录由$A$贡献的部分）

void cross(int a,int b,int c)
{
    int i,j; ll u,k;
    ll x[100],y[100];
    fo(i,0,63) x[i] = y[i] = basis[a][i];
    fo(i,0,63) c[i] = 0;
    fo(i,0,63)
    if (basis[b][i])
    {
        u = basis[b][i]; k = 0;
        fd(j,63,0)
        if (u & (1ll<<j))
            if (x[j])
            {
                u = u ^ x[j];
                k = k ^ y[j];
            }   else
            {
                x[j] = u;
                y[j] = k;
                break;
            }
        if (u == 0) insert(c,k);
    }
}

线性基求并

暴力将线性基$B$中的每一个数插入到线性基$A$中
实际上就是调用若干次insert操作

void merge(int a,int b,int c)
{
	int i;
	fo(i,0,63) c[i] = a[i];
	fd(i,63,0) if (basis[b][i]) insert(c,basis[b][i]);
}

查询一个数是否存在

和插入操作类似，最后异或结果为0表示存在

bool exist(int rt,ll x)
{
    int i;
    fd(i,63,0)
    if (x & (1ll<<i))
        if (basis[rt][i] == 0)
        {
        	return false;
		}	else
		{
			x = x ^ basis[rt][i];
		}
	if (x == 0) return true; else return false;
}

查询最大值

从高位往低位组合，如果能使值更大就就异或

ll maximum(int rt)
{
    int i; ll x = 0;
    fd(i,63,0)
    if ((x ^ basis[rt][i]) > x) x = x ^ basis[rt][i];
    return x;
}

查询最小值

直接取出线性基内最小的元素就好了

ll minimum(int rt)
{
    int i;
    fo(i,0,63) if (basis[rt][i]) return basis[rt][i];
}

查询k小值

首先通过类似高斯消元的方法使得每列只有一个1（也就是把所有可以消掉的1消掉）
然后把k分解二进制，从线性基中从小到大按位取出即可
注意判断0是否可行（不取产生0/组合产生0）

int trans(int rt)
{
	int i,j;
	fo(i,0,63)
		fd(j,i-1,0)
			if (basis[rt][i] & (1<<j)) basis[rt][i] = basis[rt][i] ^ basis[rt][j];
	int num = 0;
	fo(i,0,63) if (basis[rt][i]) num++;
	return num;
}
ll get_kth(int rt,ll k)
{
	int i;
	ll x = 0;
	fo(i,0,63)
	if (basis[rt][i])
	{
		if (k&1) x ^= basis[rt][i];
		k = k / 2;
	}
	if (k > 0) return -1; else return x;
}
/*.........*/
main
{
	num = trans(1);
	if (num != n) k--; //num==n时表示无法组合出0
	ans = get_kth(1,k);
}