CodeForces 724 G.Xor-matic Number of the Graph(组合数学+dfs+线性基)

125 篇文章 0 订阅
33 篇文章 0 订阅

Description

给出一个 n 个点m条边的无向图,每条边有边权,定义三元组 (u,v,s) 合法当且仅当 1u<vn 且存在一条从 u v的路径使得路径上边的权值异或和为 s ,问所有合法三元组的s值之和

Input

第一行两个整数 n,m 表示点数和边数,之后 m 行每行三个整数ui,vi,wi表示 ui,vi 之间有一条权值为 wi 的边

(1n105,0m2105,0wi1018)

Output

输出所有合法三元组 s 值之和,结果模109+7

Sample Input

4 4
1 2 1
1 3 2
2 3 3
3 4 1

Sample Output

12

Solution

一个连通块内任意两点路径可以表示为这两点的一条简单路径和一个环组合而成,对于 u 点所处连通块的两点v,w,如果求出 v,w u 的简单路径权值dis[v] dis[w] ,那么可知 v w的一条路径权值为 dis[v] ^ dis[w] ,然后异或上一个环的权值可以得到所有 v w的路径权值,而任意环的权值可以通过求出简单环的权值线性基表示出

现在考虑每一位对答案的贡献,即有多少点对,其路径权值在第 i 位是1,假设对于 u 点所处连通块,一遍dfs即可求出任一点 v u的简单路径权值 dis[v] ,且知道 dis[v] 中在第 i 位为0 1 的个数num[i][0] num[i][1] ,且已经得到简单环的权值,求出其线性基为 a[1],...,a[m]

如果存在某个基在第 i 位是1,那么有 2m1 种情况形成权值在第 i 位为0的环,同样的也有 2m1 种情况形成权值在第 i 位为1的环,这样一来,如果 dis[v] dis[w] 在第 i 位相同,则异或完在该位就是0,加上一个权值在该位是 1 的环即可得到一条从v w 的权值在第i位为 1 的路径,如果dis[v] dis[w] 在第 i 位不同,则异或完在该位就是1,加上一个权值在该位是 0 的环同样可以得到一条从v w 的权值在第i位为 1 的路径,故对答案的贡献就是2i2m1(C2num[i][0]+C2num[i][1]+num[i][0]num[i][1])

如果没有基在第 i 位是1,那么这些基所能表示出的 2m 个环的权值在第 i 位都是0,故若要得到一条从 v w的权值在第 i 位为1的路径, dis[v] dis[w] 在第 i 位必然不同,对答案贡献就是2i2mnum[i][0]num[i][1]

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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值