【JZOJ5239】【GDOI2018模拟8.7】图的异或

Description

Data Constraint

Solution

这道题要用线性基来做，看到一个博客讲线性基感觉挺好[学习笔记]线性基
线性基的目的就是用log个数的任意异或起来的集合表示一堆数任意异或起来的集合。
对于这道题，我们跑一个s出发的dfs树，然后对于每一条返租边或自环，我们将他们放进一个集合，他们任意异或后的值异或上s到t的异或值的和就是答案。
我们求出线性基和s到t的异或值dt后，假如dt第i位为1，若线性基中所有数的第i为均为0，那么线性基的所有异或出来的数这一位都肯定是0，所以均有贡献。否则只有一半的数这一位都是0有贡献。dt第i位为0的情况分析同理。

Code

#include<iostream> 
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const ll maxn=2e5+5,mo=1e9+7;
ll first[maxn],last[maxn],next[maxn],value[maxn],bz[maxn],bz1[maxn],a[maxn],b[maxn];
ll n,m,i,t,j,k,l,x,y,z,s,st,er[maxn],num,ans;
void insert(ll x){
    ll i,j;
    for (i=60;i>=0;i--)
        if (er[i]&x){
            if (!b[i]){
                b[i]=x,num++;break;
            }else x^=b[i];
        }
}
void lian(ll x,ll y,ll z){
    last[++num]=y;next[num]=first[x];first[x]=num;value[num]=z;bz[num]=i;
}
void dg(ll x,ll y){
    int t;bz1[x]=1;
    for (t=first[x];t;t=next[t])
        if (!bz1[last[t]]) a[last[t]]=(a[x]^value[t]),dg(last[t],bz[t]);
        else if (bz[t]!=y)insert(a[x]^value[t]^a[last[t]]);
}
int main(){
    freopen("xor.in","r",stdin);freopen("xor.out","w",stdout);
    scanf("%lld%lld",&n,&m);
    for (i=1;i<=m;i++)
        scanf("%lld%lld%lld",&x,&y,&z),lian(x,y,z),lian(y,x,z);
    scanf("%lld%lld",&s,&st);er[0]=1;
    for (i=1;i<=60;i++)er[i]=er[i-1]*2;num=0;
    dg(s,0);
    for (i=0;i<=60;i++){
        x=er[i]%mo;
        if (er[i]&a[st]){       
            if (!b[i]){
                t=0;
                for (j=0;j<=60;j++)
                    if (b[j]&er[i]) t++;
                if (!t) y=er[num]%mo,ans=(ans+x*y%mo)%mo;
                else y=er[num-1]%mo,ans=(ans+x*y%mo)%mo;
            }
            else y=er[num-1]%mo,ans=(ans+x*y%mo)%mo;
        }else{
            t=0;
            for (j=i;j<=60;j++)
                if (b[j]&er[i]) t++;
            if (t) y=er[num-1]%mo,ans=(ans+x*y%mo)%mo;
        }
    }
    printf("%lld\n",ans);
}