「WC2011」Xor - 线性基

给一个无向连通图,求一条从 的路径(可以不是简单路径),使经过的边权的异或和最大。

链接

BZOJ 2115

题解

首先,一个显然的结论是,图中所有简单环的异或和,都是可以直接获得的 —— 因为我们可以从起点走到一个环,由 进入这个环,走一圈回到 然后原路返回起点,这样起点和 之间的边都经过了两次,相当于没有经过这些边,而环上的边权被留在了答案中。

我们可以 BFS 求出所有简单环的边权异或和,并求出它们的线性基。任意找一条从 的路径,然后从大到小考虑线性基中的每个元素加入到答案中会不会使答案更大即可。

**证明:**设当前答案为 ,考虑到的线性基中的当前元素为 控制着线性基中的第 位,那么分为三种情况:

  1. 中第 位为 中第 位为 ,此时不选
  2. 中第 位为 中第 位为 ,此时选
  3. 中第 位为 ,此时一定有 ,不会影响。

根据归纳法,最终答案是正确的。

代码

#include <cstdio>
#include <vector>
#include <queue>

const int MAXN = 50000;
const int MAXL = 60;

struct Node {
    std::vector<struct Edge> e;
    Node *fa;
    long long dist;
    int dep;
} N[MAXN + 1];

struct Edge {
    Node *s, *t;
    long long w;

    Edge(Node *s, Node *t, long long w) : s(s), t(t), w(w) {}
};

inline void addEdge(int s, int t, long long w) {
    N[s].e.push_back(Edge(&N[s], &N[t], w));
    N[t].e.push_back(Edge(&N[t], &N[s], w));
}

inline std::vector<long long> bfs() {
    std::queue<Node *> q;
    q.push(&N[1]);
    N[1].dep = 1;

    std::vector<long long> circles;
    while (!q.empty()) {
        Node *v = q.front();
        q.pop();

        for (Edge *e = &v->e.front(); e <= &v->e.back(); e++) {
            if (!e->t->dep) {
                e->t->dep = v->dep + 1;
                e->t->dist = v->dist ^ e->w;
                e->t->fa = v;
                q.push(e->t);
            } else if (e->t != v->fa) {
                circles.push_back(e->t->dist ^ v->dist ^ e->w);
            }
        }
    }

    return circles;
}

struct LinearBasis {
    std::vector<long long> v;

    void build(std::vector<long long> x) {
        std::vector<long long> a(MAXL + 1);

        for (int i = 0; i < int(x.size()); i++) {
            long long t = x[i];

            for (int j = MAXL; j >= 0; j--) {
                if (!(t & (1ll << j))) continue;

                if (a[j]) t ^= a[j];
                else {
                    for (int k = 0; k < j; k++) if (t & (1ll << k)) t ^= a[k];
                    for (int k = j + 1; k <= MAXL; k++) if (a[k] & (1ll << j)) a[k] ^= t;
                    a[j] = t;
                    break;
                }
            }
        }

        v.clear();
        for (int i = 0; i <= MAXL; i++) if (a[i]) v.push_back(a[i]);
    }

    long long queryMaxWith(long long x) {
        long long res = x;
        for (int i = 0; i < int(v.size()); i++) if ((res ^ v[i]) > res) res ^= v[i];
        return res;
    }
} lb;

int main() {
    int n, m;
    scanf("%d %d", &n, &m);

    while (m--) {
        int u, v;
        long long w;
        scanf("%d %d %lld", &u, &v, &w);
        addEdge(u, v, w);
    }

    std::vector<long long> circles = bfs();
    lb.build(circles);

    printf("%lld\n", lb.queryMaxWith(N[n].dist));

    return 0;
}