hdu 5544 Ba Gua Zhen

link

因为在沈阳站的比赛中我们没有做出来H题，后来才知道是高斯消元。于是我先来补一下高斯消元的题目。先做了hdu XOR对线性基有了一定的了解，然后来做一下南阳ccpc的Ba Gua Zhen

分析：这个想法我也很快就想到了。看第二组样例就应该能够看出来了。
在图上，其实只有每个环上的边权的异或值有对答案有贡献。对于图中的任意一条桥边，要么没有经过，要么经过了偶数次数，所以桥边的边权贡献为0。这样的话，问题转化成了给一些数，求从中选取一些异或的最大值。
只要把所有环的异或值先找出来，然后高斯消元求一下线性基，这样最大的值就是线性基全部异或。（我知道这里也许会有童鞋有疑问，看一下XOR那道题，大家就会明白了）
想到这里，原先感觉这到题已经可以解决了，然而死在了怎么找图中的环上了。。。
我一开始试图用点的双联通分量来解环的问题，然后发现了一个大的bug。我想各位菊苣已经在暗中嘲讽我了。。
后来看了别人的博客，发现了一个十分优美的方法。dfs + 时间戳。时间戳保存的是最后一次访问这个节点的时间。然后每经过一条边，就把边权值异或进到达的点里。如果经过了一个环两次，那么说明找到了一个环，对于环上的异或值，就是$val[u] \otimes val[v] \otimes edge.val$
稍微想一想应该就能看出来了。试想从某一个点走到了点$A$，这个时候$val[a] = t$，那么当再次从某一点$B$走到点$A$的时候，$val[b]$保存了从$A$到$B$的一路走来的异或值，那么$val[b] ^ edge.val$就是从某一点先经过$A$，然后到$B$再回到$A$的权值，异或一下$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;
}