codeforces 724G

异或运算求权值和

求任意两点a,b 从a到b路径上边权xor值的总和

算出每一位的贡献求和

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
const int N=1e6+10;
const int M=2e6+10;
const int mod=1e9+7;
typedef __int64 LL;
struct Edge {
    int v,next; LL w;
    Edge(){}
    Edge(int v,int next,LL w):v(v),next(next),w(w){}
}e[M*2];
int head[N],total;
void init(){
    memset(head,-1,sizeof(head));total=0;
}
void adde(int u,int v,LL w){
    e[total]=Edge(v,head[u],w);head[u]=total++;
}

LL dis[N];
LL zz[N];int cnt;
int sta[N],top;
int dep[N];
void dfs(int u,int f){
    sta[top++]=u;
    for(int i=head[u];i!=-1;i=e[i].next){
        int v=e[i].v;
        if(v==f)continue;
        if(dis[v]==-1){
            dis[v]=dis[u]^e[i].w;
            dep[v]=dep[u]+1;
            dfs(v,u);
        }
        else if(dep[v]<dep[u]-1){
            zz[cnt++]=dis[v]^dis[u]^e[i].w;
        }
    }
}
LL mi[66];
int main(){
#ifdef DouBi
freopen("in.cpp","r",stdin);
#endif // gh546
    int n,m;
    mi[0]=1;
    for(int i=1;i<=60;i++){
        mi[i]=mi[i-1]*2%mod;
    }
    while(scanf("%d%d",&n,&m)!=EOF){
        init();
        for(int i=0;i<m;i++){
            int a,b;LL c;scanf("%d%d%I64d",&a,&b,&c);a--;b--;
            adde(a,b,c);
            adde(b,a,c);
        }
        memset(dis,-1,sizeof(dis));
        LL ans=0;
        for(int i=0;i<n;i++)if(dis[i]==-1){
            cnt=0,top=0;dis[i]=0;dep[i]=0;
            dfs(i,-1);
            int r=0;
            if(cnt){
                for(int j=59;j>=0;j--){
                    int id=r;
                    for(int k=r+1;k<cnt;k++)if(zz[k]>zz[id]){
                        id=k;
                    }
                    swap(zz[r],zz[id]);
                    if(!zz[r])break;
                    for(int k=r+1;k<cnt;k++)if((zz[k]^zz[r])<zz[k])
                        zz[k]=zz[k]^zz[r];
                    r++;
                }
            }
            LL x=0;
            for(int j=0;j<r;j++)x|=zz[j];
            for(int j=59;j>=0;j--){
                if(x>>j&1){
                    ans+=1LL*top*(top-1)/2%mod*mi[r-1]%mod*mi[j]%mod;
                    ans%=mod;
                }
                else {
                    int cnt1=0,cnt0=0;
                    for(int k=0;k<top;k++){
                        int u=sta[k];
                        if(dis[u]>>j&1)cnt1++;
                        else cnt0++;
                    }
                    ans+=1LL*cnt1*cnt0%mod*mi[r]%mod*mi[j]%mod;
                    ans%=mod;
                }
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}