2021 ICPC沈阳 B.Bitwise Exclusive-OR Sequence（位运算+图论）

题目描述

题目链接

题目大意

有n个点和m个约束，每个约束包含三个数u，v，w，表示a[u]^a[v]=w.
求满足m个约束 并且 和最小的序列的和。

题目分析

$这很明显是一个图论问题，对于每个约束u-v-w，我们可以在u和v之间连一条权值为w的边。$

$因为题目对于约束并没有太多的限制。因此这个图是可能存在环和不连通的。$

$首先，对于图的不连通，我们可以分别处理每个连通块。并将每一块的权值进行累加。$

$当一个连通块中存在环的时候，我们可以算出环上所有边的权值异或和，如果环上的异或和不为0，则无合法解，$$输出-1。因为环上所有边的异或和是把环上所有点都用了两遍，因此环上的异或和应该是0。$

$除去这两种情况，剩下的就只有链了。最后我们来看一下如何处理链。$
$假设一个链上的点权为（a1,a2,a3,\dots ,an）$

$a_1\wedge a_2=w_1$
$a_2\wedge a_3=w_2$
$\dots \dots$
$a_{n-1}\wedge a_n=w_{n-1}$

$我们对其做一个转换：$

$a_1\wedge a_2=w_1$
$a_1\wedge a_3=w_1\wedge w_2$
$a_1\wedge a_4=w_1\wedge w_2\wedge w_3$
$\dots \dots$
$a_1\wedge a_n=w_1\wedge w_2\wedge \dots \wedge w_{n-1}$

$设s_i为w_{1-i}的异或前缀和，我们再移一下项$

$a_2=s_1\wedge a_1$
$a_3=s_2\wedge a_1$
$\dots \dots$
$a_n=s_{n-1}\wedge a_1$

$因此，我们只需要确定第一个数a_1，即可算出所有的其它的点的点权。$
$如果我们枚举a_1的所有可能值，那一定是会超时的。所以我们需要寻找其它方法。$

$我们可以按位枚举每一个s_i，并记录每一位1的数量cnt[j]。$
$对于每一位j，1的个数是cnt[j]，一条链上的点数为n。$

$如果cnt[j]>n/2，说明有一半以上的s_i在第j位上都是1。因为a_{i+1}=s_i\wedge a_1，所以只要让a_1在第j位$$也为1，就能让一半以上的点权减少(1<<j)，而剩下的点权增加(1<<j)，对于总和来说是减小了。$

$用这种方式就能在O(nlogn)的复杂度内，求出一条链的和最小值。$

代码如下

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <set>
#include <map>
#include <queue>
#include <vector>
#include <algorithm>
#include <iomanip>
#define LL long long
#define ULL unsigned long long
#define PII pair<int,int>
#define PLL pair<LL,LL>
#define PDD pair<double,double>
#define x first
#define y second
using namespace std;
const int N=1e5+5,mod=1e9+7;
vector<PII> h[N];
bool vis[N],flag=1;			//flag=0表示存在非法环，无解
LL f[N],cnt[35];			//f[v]表示从根节点到v点这条链上的异或和
vector<int> a;				//记录这条链上的点（除第一个节点外）
void dfs(int u)
{
	vis[u]=true;
	for(auto t:h[u])
	{
		int v=t.x,w=t.y;
		if(vis[v])			//如果v被访问过，说明存在环
		{
			if((f[u]^w)!=f[v]) flag=0;		//如果环上的异或和不为0，用flag记录
		}
		else {
			f[v]=f[u]^w;			//更新f[v]的值
			a.push_back(v);			//将v加入链中
			dfs(v);					//继续递归
		}
	}
}
LL get()
{
	memset(cnt,0,sizeof cnt);		//先清空cnt[]
	int x=0;					//表示链上的第一个节点的点权
	for(int i=0;i<30;i++)
	{
		for(int u:a)		//记录链上所有数第i位是否为1（f[u]其实和讲解里的s[]等价）
			if((f[u]>>i)&1) cnt[i]++;
		if(cnt[i]>a.size()/2) x+=(1<<i);
	}
	LL ans=x;
	for(int u:a) ans+=f[u]^x;		//记录链上的点权和
	return ans;
}
int main()
{
	cin.tie(0);
	ios::sync_with_stdio(false);
	int n,m;
	cin>>n>>m;
	while(m--)				//建图
	{
		int u,v,w;
		cin>>u>>v>>w;
		h[u].push_back({v,w});
		h[v].push_back({u,w}); 
	}
	LL ans=0;
	for(int i=1;i<=n;i++)			//遍历所有点
	{
		if(!vis[i]&&flag)			//如果点i没被访问过，且不存在非法环
		{
			a.clear(); 				//清空数组
			dfs(i);					//dfs求f[]
			ans+=get();				//计算当前链的最小和
		}
	}
	if(flag) cout<<ans<<endl;
	else cout<<-1<<endl;
	return 0; 
}