因为在沈阳站的比赛中我们没有做出来H题,后来才知道是高斯消元。于是我先来补一下高斯消元的题目。先做了hdu XOR对线性基有了一定的了解,然后来做一下南阳ccpc的Ba Gua Zhen
分析:这个想法我也很快就想到了。看第二组样例就应该能够看出来了。
在图上,其实只有每个环上的边权的异或值有对答案有贡献。对于图中的任意一条桥边,要么没有经过,要么经过了偶数次数,所以桥边的边权贡献为0。这样的话,问题转化成了给一些数,求从中选取一些异或的最大值。
只要把所有环的异或值先找出来,然后高斯消元求一下线性基,这样最大的值就是线性基全部异或。(我知道这里也许会有童鞋有疑问,看一下XOR那道题,大家就会明白了)
想到这里,原先感觉这到题已经可以解决了,然而死在了怎么找图中的环上了。。。
我一开始试图用点的双联通分量来解环的问题,然后发现了一个大的bug。我想各位菊苣已经在暗中嘲讽我了。。
后来看了别人的博客,发现了一个十分优美的方法。dfs + 时间戳。时间戳保存的是最后一次访问这个节点的时间。然后每经过一条边,就把边权值异或进到达的点里。如果经过了一个环两次,那么说明找到了一个环,对于环上的异或值,就是
val[u]⊗val[v]⊗edge.val
稍微想一想应该就能看出来了。试想从某一个点走到了点
A
,这个时候
补充一点,http://www.cnblogs.com/andyqsmart/p/4968589.html,一张图只有M-N+1个独立回路。但要注意可能在计算的时候会包括由多个小环组成的大环,所以数组还是开得大一点吧。
至此,终于完成了题目。。。
/*****************************************************/
//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
#define offcin ios::sync_with_stdio(false)
#define sigma_size 26
#define lson l,m,v<<1
#define rson m+1,r,v<<1|1
#define slch v<<1
#define srch v<<1|1
#define sgetmid int m = (l+r)>>1
#define ll long long
#define ull unsigned long long
#define mem(x,v) memset(x,v,sizeof(x))
#define lowbit(x) (x&-x)
#define bits(a) __builtin_popcount(a)
#define mk make_pair
#define pb push_back
#define fi first
#define se second
const int INF = 0x3f3f3f3f;
const ll INFF = 1e18;
const double pi = acos(-1.0);
const double inf = 1e18;
const double eps = 1e-9;
const ll mod = 1e9+7;
const int maxmat = 10;
const ull BASE = 133333331;
/*****************************************************/
inline void RI(int &x) {
char c;
while((c=getchar())<'0' || c>'9');
x=c-'0';
while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
}
/*****************************************************/
const int maxn = 5e4 + 5;
const int maxm = 1e5 + 5;
std::vector<pair<int, ll> > G[maxn];
int N, M;
int pre[maxn], dfs_clock;
ll mat[maxn << 2], p[maxn];
int pos;
void Debug() {
for (int i = 1; i <= pos; i ++) {
for (int j = 63; j >= 0; j --)
cout<<((mat[i] >> j) & 1);
cout<<endl;
}
}
void dfs(int u, int fa) {
pre[u] = ++ dfs_clock;
for (unsigned i = 0; i < G[u].size(); i ++) {
int v = G[u][i].first;
if (v == fa) continue;
if (pre[v] != -1){
if (pre[v] < dfs_clock) mat[++ pos] = p[v] ^ p[u] ^ G[u][i].second;
}
else {
p[v] = p[u] ^ G[u][i].second;
dfs(v, u);
}
}
}
void find(int n) {
memset(pre, -1, sizeof(pre));
memset(mat, 0, sizeof(mat));
memset(p, 0, sizeof (p));
dfs_clock = pos = 0;
for (int i = 1; i <= n; i ++)
if (pre[i] == -1) dfs(i, -1);
}
ll work() {
int k = 1;
// Debug();
for (int i = 63; i >= 0; i --) {
int t = 0;
for (int j = k; j <= pos; j ++) {
if ((mat[j] >> i) & 1) {
t = j;
break;
}
}
if (t) {
swap(mat[k], mat[t]);
for (int j = 1; j <= pos; j ++)
if (j != k && (mat[j] >> i) & 1)
mat[j] ^= mat[k];
k ++;
}
}
// cout<<endl;
// Debug();
ll res = 0;
for (int i = 1; i <= pos; i ++)
res ^= mat[i];
return res;
}
int main(int argc, char const *argv[]) {
// freopen("t.txt", "r", stdin);
int T; cin>>T;
for (int kase = 1; kase <= T; kase ++) {
scanf("%d%d", &N, &M);
for (int i = 0; i <= N; i ++) G[i].clear();
for (int i = 0; i < M; i ++) {
int u, v; ll val;
scanf("%d%d%I64d", &u, &v, &val);
G[u].push_back(make_pair(v, val));
G[v].push_back(make_pair(u, val));
}
find(N);
ll res = work();
printf("Case #%d: %I64d\n", kase, res);
}
return 0;
}