Description
给出一个
n
个点
Input
第一行两个整数
n,m
表示点数和边数,之后
m
行每行三个整数
(1≤n≤105,0≤m≤2⋅105,0≤wi≤1018)
Output
输出所有合法三元组
s
值之和,结果模
Sample Input
4 4
1 2 1
1 3 2
2 3 3
3 4 1
Sample Output
12
Solution
一个连通块内任意两点路径可以表示为这两点的一条简单路径和一个环组合而成,对于
u
点所处连通块的两点
现在考虑每一位对答案的贡献,即有多少点对,其路径权值在第
i
位是
如果存在某个基在第
i
位是
如果没有基在第
i
位是
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=100005;
#define mod 1000000007
int n,m,vis[maxn],tot,head[maxn],cnt[66][2];
ll dis[maxn],base[66];
vector<ll>a;
struct node
{
int v,next;
ll w;
}g[4*maxn];
void add_edge(int u,int v,ll w)
{
g[tot].v=v,g[tot].next=head[u],g[tot].w=w,head[u]=tot++;
}
void dfs(int u,ll d)
{
dis[u]=d,vis[u]=1;
for(int i=0;i<=62;i++)
cnt[i][(dis[u]>>i)&1]++;
for(int i=head[u];~i;i=g[i].next)
{
int v=g[i].v;
ll w=g[i].w;
if(vis[v])a.push_back(dis[u]^dis[v]^w);
else dfs(v,dis[u]^w);
}
}
void add(int &x,int y)
{
x=x+y>=mod?x+y-mod:x+y;
}
int main()
{
scanf("%d%d",&n,&m);
memset(head,-1,sizeof(head));
while(m--)
{
int u,v;
ll w;
scanf("%d%d%I64d",&u,&v,&w);
add_edge(u,v,w),add_edge(v,u,w);
}
int ans=0;
for(int u=1;u<=n;u++)
if(!vis[u])
{
a.clear();
memset(base,0,sizeof(base));
memset(cnt,0,sizeof(cnt));
dfs(u,0);
for(int i=0;i<a.size();i++)
for(int j=62;j>=0;j--)
if((a[i]>>j)&1)
{
if(!base[j])
{
base[j]=a[i];
break;
}
else a[i]^=base[j];
}
int num=0;
for(int i=0;i<=62;i++)
if(base[i])num++;
for(int i=0;i<=62;i++)
{
int res=0,x=cnt[i][0],y=cnt[i][1],flag=0;
for(int j=0;j<=62;j++)
if((base[j]>>i)&1)
{
flag=1;
break;
}
if(flag)
{
int tx=(ll)x*(x-1)/2%mod,ty=(ll)y*(y-1)/2%mod,tz=(1ll<<(num-1))%mod;
add(res,(ll)(tx+ty)%mod*tz%mod);
add(res,(ll)x*y%mod*tz%mod);
}
else
{
int tz=(1ll<<num)%mod;
add(res,(ll)x*y%mod*tz%mod);
}
add(ans,(ll)(1ll<<i)%mod*res%mod);
}
}
printf("%d\n",ans);
return 0;
}