省选专练【WC2011】最大XOR和路径（异或问题学习笔记）

我不会线性基啊！！！

我还搞个屁

但是也不是不可以做。

考虑两个线性基的模板：给定n个非负整数选任意多个求最大异或值。

朴素解法（不是暴力）
二分答案（所谓二分答案啊，不是说二分最大异或值，而是按位异或）

我先进来：枚举最高位1可否？

怎么枚举：问题转化：给定n个非负整数任选多个求异或值是否可以达到S

解决方案：这就是开关灯问题啊：给定操作策略然后求能否关完灯。

解线性异或方程组。

X1^X2^X3^0^x4^0=1

X1表示原本这个对应操作要关上灯对吧。

然后按位拆数。

一个数二进制表示下可以得到是否关这位的灯的含义。

于是高斯消元，按位查就好了。

复杂度：O（N*logn^3）

但是考虑类似增广的思想

不需要每次重新构图，继续在残图增广就好了。

复杂度：O（N*logn^2）;

这依旧卡常数

在矩阵上做高斯消元

当前列为自由元：无法改变这一位

否则可以改变消去其余数的这一位

剩余一个三角矩阵，向上累乘。

然后如果是求这n个数和一个给定数S的异或值最大？

还是这样搞就行了。

只是累乘的时候判断S对应位。

那么这个题怎么搞？

一条路径转化为另一条路径的代价是原代价抑或上一个环。

异或的哲学性质啊！！！！一个数异或自己为0，0异或任何数为其本身。

然后预处理环

转变为：

已知一个S可以异或一些数然后最大。

完，就是如上解法。

#include<bits/stdc++.h>
using namespace std;
typedef int INT;
#define int long long
const int N=5e4+10000;
const int maxn=1e5+100;
inline void read(long long &x){
	x=0;
	int f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-'){
			f=-1;
		}
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		x=x*10+ch-'0';
		ch=getchar();
	}
	x*=f;
}
struct Front_star{
	long long u,v,w,nxt;
}e[maxn*4];
long long cnt=0;
long long first[N]={};
void add(long long u,long long v,long long w){
	cnt++;
	e[cnt].u=u;
	e[cnt].v=v;
	e[cnt].w=w;
	e[cnt].nxt=first[u];
	first[u]=cnt;
}
long long n,m;
long long _xor[N]={0};//_xor just record the way to the end
long long vis[N]={};
long long S[maxn*4]={};
long long top=0;
void dfs(long long u){
	vis[u]=1;
	for(long long i=first[u];i;i=e[i].nxt){
		long long v=e[i].v;
		if(vis[v]){//if visit v already means there must be a Circle so push the sum into the Stack in order to judge the max xor  
			S[++top]=e[i].w^_xor[u]^_xor[v];
		}
		else{
			_xor[v]=e[i].w^_xor[u];
			dfs(v);	
		} 
	}
}
void Gaussian_Elimination(){
	long long i,k=0;
	for(long long j=((long long)1<<60);j;j>>=1){
		for(i=k+1;i<=top;i++){//不需要层层增广,记录上一次的位置。 
			if(S[i]&j)break;  
		}
		if(i==top+1)continue;
		swap(S[++k],S[i]);
		for(i=1;i<=top;i++){
			if(i!=k){
				if(S[i]&j){
					S[i]^=S[k];
				}
			}
		}
	}
}
INT main(){
	memset(_xor,0,sizeof(_xor));
	read(n);
	read(m);
	for(long long i=1;i<=m;i++){
		long long u,v,w;
		read(u);
		read(v);
		read(w);
		add(u,v,w);
		add(v,u,w);
	}
	dfs(1);
	Gaussian_Elimination();
	long long ans=_xor[n];
	for(long long i=1;i<=top;i++){
		if((ans^(S[i]))>ans){
			ans^=S[i];
		}
	}
	cout<<ans;
}