线性基 Codeforces724G Xor-matic Number of the Graph

最新推荐文章于 2022-09-06 17:54:51 发布

逍遥丶綦

最新推荐文章于 2022-09-06 17:54:51 发布

阅读量1.5k

点赞数 1

分类专栏： ACM_线性基

本文链接：https://blog.csdn.net/qwb492859377/article/details/52944803

版权

ACM_线性基专栏收录该内容

2 篇文章 0 订阅

订阅专栏

该博客探讨了Codeforces 724G问题，涉及一个含有1e5节点的无向图，目标是找到所有三元组(u, v, s)，其中u到v的路径上边的异或和为s。文章介绍了如何处理图中的环以求得所有三元环的异或和总和，并提出了利用线性基解决含有环的连通块的方法。博主通过分析树形结构和DFS搜索树，详细解释了如何计算答案。" 124261052,7614210,checkra1n越狱问题解决方案：修复失败与AFC2失效,"['ios开发', '越狱工具', 'mobile开发', 'macos应用']

摘要由CSDN通过智能技术生成

传送门：点击打开链接

题意：1e5个点的无向图，求三元环(u,v,s),u<v，满足u经过路径到达v，把路径上经过的边的值求异或和等于s。三元环的价值为s。现在求所有三元环的总价值之和。

思路：可能有很多个连通块，对每个连通块考虑。

首先，我们考虑，如果这个连通块只是一棵树的话。我们对每一位去考虑，直接通过计数就能得到答案。

但是现在图上还有很多环，我们可以发现，环的答案可以直接异或到链上去，我只要走到环上，再绕环一圈，再走回来，就能恰好只异或了环上的值。

所以我们求出所有的环，对环求一遍线性基。

之后的做法，我们前面说的树相当于现在的DFS搜索树，直接对这个树去讨论计数，只要认真讨论一下，答案非常简单就能推导出来的.

#include <map>
#include <set>
#include <cmath>
#include <ctime>
#include <stack>
#include <queue>
#include <cstdio>
#include <cctype>
#include <bitset>
#include <string>
#include <vector>
#include <cstring>
#include <iostream>
#include <iomanip>
#include <algorithm>
#include <functional>
#define fuck(x) cout<<"["<<x<<"]";
#define FIN freopen("input.txt","r",stdin);
#define FOUT freopen("output.txt","w+",stdout);
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;

const int MX = 5e5 + 5;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;

struct Edge {
    LL val;
    int v, nxt;
} E[MX];
int Head[MX], erear;
void edge_init() {
    erear = 0;
    memset(Head, -1, sizeof(Head));
}
void edge_add(int u, int v, LL val) {
    E[erear].v = v;
    E[erear].val = val;
    E[erear].nxt = Head[u];
    Head[u] = erear++;
}

int n, m, sz, r;
LL A[MX], P[62], dis[MX];
vector<int> path;

LL power(LL a, LL b) {
    LL ret = 1;
    while(b) {
        if(b & 1) ret = ret * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ret;
}

void DFS(int u, LL s) {
    if(dis[u] == -1) {
        dis[u] = s;
    } else {
        A[++sz] = dis[u] ^ s;
        return;
    }
    path.push_back(u);
    for(int i = Head[u]; ~i; i = E[i].nxt) {
        int v = E[i].v;
        DFS(v, s ^ E[i].val);
    }
}

void Guass_base() {
    memset(P, 0, sizeof(P));
    for(int i = 1; i <= sz; i++) {
        for(int j = 62; j >= 0; j--) {
            if(!(A[i] >> j & 1)) continue;
            if(!P[j]) {
                P[j] = A[i]; break;
            }
            A[i] ^= P[j];
        }
    }
    r = 0;
    for(int i = 62; i >= 0; i--) {
        if(P[i]) r++;
    }
}

LL solve() {
    Guass_base();
    LL ans = 0;
    for(int i = 62; i >= 0; i--) {
        int cnt[2] = {0}, sign = 0;
        for(int j = 62; j >= 0; j--) {
            if(P[j] >> i & 1) sign = 1;
        }
        for(int j = 0; j < path.size(); j++) {
            int u = path[j];
            if(sign) ans += power(2, i + r - 1) * j % mod;
            else ans += power(2, i + r) * cnt[!(dis[u] >> i & 1)] % mod;
            ans %= mod;
            cnt[dis[u] >> i & 1]++;
        }
    }
    return ans;
}

int main() {
    // FIN;
    edge_init();
    scanf("%d%d", &n, &m);
    memset(dis, -1, sizeof(dis));
    for(int i = 1; i <= m; i++) {
        int u, v; LL val;
        scanf("%d%d%I64d", &u, &v, &val);
        edge_add(u, v, val);
        edge_add(v, u, val);
    }

    LL ans = 0;
    for(int i = 1; i <= n; i++) {
        if(dis[i] == -1) {
            path.clear(); sz = 0;
            DFS(i, 0);
            ans += solve();
            ans %= mod;
        }
    }
    printf("%I64d\n", ans);
    return 0;
}