codeforces 724G - Xor-matic Number of the Graph 线性基+图

题目传送门

题意：给出衣服无向带权图，问有多少对合法的$<u,v,s>$，要求$u$到$v$存在一条路径（不一定是简单路径）权值异或和等于$s$，并且$u<v$。求所有合法三元组的s的和。

思路：

　　参考了一篇大佬的博客。

　　这类题的核心思想就是，两点之间的所有可能的路径，都是由一条简单路径加上若干个环组成的。u，v两点所有路径的异或值的集合，等价于，u，v一条简单路径的异或值，与整个连通图的所有环组成的线性基异或的集合。

　　所以按位考虑每个二进制1给整幅图带来的价值。

　　特别要注意的一点是，在处理线性基的过程中，最多处理到63位，如果处理到64位及以上，当$i>=64$,某些$(x>>i)$，或者$（x<<i）$，就会发出许多诡异的错误，wa了很久。

#pragma GCC optimize (2)
#pragma G++ optimize (2)
#pragma comment(linker, "/STACK:102400000,102400000")
#include<bits/stdc++.h>
#include<cstdio>
#include<vector>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,b,a) for(int i=b;i>=a;i--)
#define clr(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define pii pair<int,int >
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
ll rd()
{
    ll x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
const int maxn=;
const ll mod=1e9+;
int n,m;
struct edge{
    int to;
    ll w;
};
vector<edge >ve[maxn];
int vis[maxn];
ll del[maxn],p[],siz,num[],ans=,tot=;
void init(){
    rep(i,,n){
        ve[i].clear();
        vis[i]=;
    }
    ans=,siz=,clr(p,);
}
bool insert(ll x){
    dep(i,,){
        if((x>>i)&){
            if(!p[i]){
                p[i]=x;
                return true;
            }
            x^=p[i];
        }
    }
    return false;
}
void dfs(int u,ll res){
    del[u]=res;
    for(ll i=;i>=;i--){
        if(((res>>i)&)){
            num[i]++;
        }
    }
    tot++;
    vis[u]=;
    for(auto &st:ve[u]){
        if(!vis[st.to]){
            dfs(st.to,res^st.w);
        }else{
            if(insert(res^st.w^del[st.to])){
                siz++;
            }
        }
    }
}
int main(){
    while(cin>>n>>m){
        init();
        rep(i,,m){
            int u,v;
            ll w;
            u=rd(),v=rd(),w=rd();
            ve[u].pb({v,w});
            ve[v].pb({u,w});
        }
        rep(u,,n){
            if(!vis[u]){
                tot=;
                clr(p,),siz=;
                clr(num,);
                dfs(u,);
                dep(i,,){
                    bool ok=;
                    dep(j,,){
                        if((p[j]>>i)&)ok=;
                    }
                    if(ok){
                        ans=(ans+tot*(tot-)/%mod*((1ll<<(siz-))%mod)%mod*((1ll<<i)%mod)%mod)%mod;
                    }else{
                        ans=(ans+num[i]*(tot-num[i])%mod*((1ll<<siz)%mod)%mod*((1ll<<i)%mod)%mod)%mod;
                    }
                }
            }
        }
        printf("%lld\n",ans);
    }
}