题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5544
题目大意:给出一个无向连通图,求一个回路,使得该回路中所有路径权值的亦或最大(某条边经过几次,则亦或几次)
思路:首先找出图中所有的独立环,所谓独立环,就是不能由其它环线性表示出的环。易知,回路是由若干个独立环线性组合而成的。那么我们可以求出每个独立环的亦或值,存到数组中,那么问题就转化为了:给出若干元素,从中选出几个,使它们的总亦或最大。
因为图是连通的,所以从1出发即可找到所有的环。另d[i]表示从1出发到点i的亦或值,当到达点u时,若u的孩子v已经访问过了,说明u和v在同一环中,并且该环的亦或为d[u]^d[v]^dist(u, v)。因为从1到u与从1到v有一段共同的路程,所以亦或后这一段变为0,剩下的是从v到u的路程,所以再加上边(u, v)即组成了一个环。
可以用线性基来求解一个集合中若干元素亦或的最大值。关于线性基的介绍可参照https://blog.sengxian.com/algorithms/linear-basis。简单来说,对于线性基中的每个数,其二进制表示中最多只有1位为1,且线性基中每个数的二进制表示中1的位置都不相同。因此,在本题中,将线性基中所有元素亦或起来就是答案。
#include<cstdio>
#include<cstring>
#include<string>
#include<cctype>
#include<iostream>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<queue>
#include<algorithm>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
typedef long long LL;
using namespace std;
const int maxn = 1e5 + 10;
struct Edge {
int from, to; LL cost;
Edge(int u, int v, LL c):from(u),to(v),cost(c){}
};
LL d[maxn], circle[maxn];
int vis[maxn];
int cnt;
vector<int> G[maxn];
vector<Edge> edges;
void AddEdge(int from, int to, LL cost) {
edges.push_back(Edge(from, to, cost));
edges.push_back(Edge(to, from, cost));
int m = edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
void init(int n) {
memset(d, 0, sizeof d);
memset(vis, 0, sizeof vis);
for(int i = 1; i <= n; i++) G[i].clear();
edges.clear();
cnt = 0;
}
void input_Edge(int m) {
for(int i = 1; i <= m; i++) {
int u, v; LL c;
scanf("%d%d%lld", &u, &v, &c);
AddEdge(u, v, c);
}
}
void dfs(int u, int fa, int t) {
vis[u] = t;
for(int i = 0; i < G[u].size(); i++) {
Edge& e = edges[G[u][i]];
if(e.to == fa) continue;
if(vis[e.to] && vis[e.to] <= t) {
circle[cnt++] = d[u] ^ d[e.to] ^ e.cost;
}
else if(!vis[e.to]) {d[e.to] = d[u] ^ e.cost; dfs(e.to, u, t+1);}
}
}
int Gauss(int n) { //传入一个n个元素的集合,求该集合的线性基,返回线性基的大小
int row = 0;
for(int i = 62; i >= 0; i--) {
int j;
for(j = row; j < n; j++)
if(circle[j] & (1LL << i))
break;
if(j != n) {
swap(circle[row], circle[j]);
for(j = 0; j < n; j++) {
if(j == row) continue;
if(circle[j] & (1LL << i))
circle[j] ^= circle[row];
}
row++;
}
}
return row;
}
int main() {
int T, kase = 0;
scanf("%d", &T);
while(T--) {
int n, m;
scanf("%d%d", &n, &m);
init(n);
input_Edge(m);
dfs(1, -1, 1);
int num = Gauss(cnt);
LL ans = 0;
for(int i = 0; i < num; i++)
ans = max(ans, ans ^ circle[i]);
printf("Case #%d: %lld\n", ++kase, ans);
}
return 0;
}