G. Xor-matic Number of the Graph

http://codeforces.com/problemset/problem/724/G

题意:给你一张无向图。定义一个无序三元组(u,v,s)表示u到v的(不一定为简单路径)路径上xor值为s。求出这张无向图所有不重复三元组的s之和。1≤n≤10^5,1≤m≤2*10^5。

想法:

如果做过【Wc2011 xor】这道题目(题解),那么问题变得简单起来了。

①假设我们钦定一个(u,v),设任意一条u->v的路径xor值为X,该连通图所有小环xor值构成的序列为{Ai}。

那么(u,v)的所有路径的xor值可以由X xor {Ai}的子集xor值得到。于是一个(u,v)的 s 之和变成了求X xor{Ai}的子集可以得到多少个不同的数,这些不同的数的和是多少?

如果能知道{Ai}的子集xor值的值域,那么好办了。于是用线性基得到值域{T}。求和的话,按位考虑定义S(i)为{T}中第i为1的个数,为0的个数取个补集就好了。

对于一个(u,v):

②考虑所有的无序点对(u,v)的答案。上面说过任意一条u->v的路径都可以,不如就钦定是DFS遍历得到DFS树的树上路径。

树上两点路径xor值的求法很简单:设dis(i)表示第i个到根节点路径xor值。

dis(a,b)=dis(a) xor dis(lca(a,b)) xor dis(b) xor dis(lca(a,b))=dia(a) xor dis(b)。

根据上面求ans 的式子,ans只与X的第j位是什么有关,所以设cnt(i)表示两点路径xor值第i位为1的个数。cnt(i)可以利用上面dis(a,b)=dia(a) xor dis(b)求。

对于所有(u,v):

于是解决了。

#include<cstdio>
#include<vector>
#define ll long long
const int len(),MP();
struct Node{int nd;ll co;};
std::vector<Node>Edge[len+];
int n,m,u,v,top,ans,much,vis[len+];
ll sum,S[],cnt[],now[];//cnt(i) 统计 i-th =0的个数 now(i):dx^dy i-th==1的个数
ll st[len+],dis[len+],All,t;
struct Base_Linear
{
ll p[];int size;
void ins(ll x)
{
for(int j=;j>=;j--)
if((x>>j)&)
{
if(p[j])x^=p[j];
else {p[j]=x;size++;break;}
}
}
}BL;
template <class T>void read(T &x)
{
x=;int f=;char ch=getchar();
while(ch<''||ch>''){f=(ch=='-');ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
x=f?-x:x;
}
ll power(int a,int b)
{
ll t=,y=a;
for(;b;b>>=)
{
if(b&)t=(t*y)%MP;
y=(y*y)%MP;
}
return t;
}
void add(int a,int b,ll c){Edge[a].push_back((Node){b,c});}
void plus(ll x)
{
for(int j=;j<=;j++)
if(((x>>j)&)==)cnt[j]++;
}
void Dfs(int x)
{
much++; vis[x]=; plus(dis[x]); st[++top]=dis[x];
for(int v=,sz=Edge[x].size();v<sz;v++)
{
Node y=Edge[x][v];
if(vis[y.nd])BL.ins(dis[x]^dis[y.nd]^y.co);
else
{
dis[y.nd]=dis[x]^y.co;
Dfs(y.nd);
}
}
}
void Back()
{
for(int j=;j<=;j++)cnt[j]=now[j]=S[j]=;
for(int j=;j<=;j++)BL.p[j]=; BL.size=;
much=; All=;
}
void Total()
{
for(;top;top--)
{
for(int j=;st[top];j++,st[top]>>=)
if(st[top]&)now[j]=(now[j]+cnt[j])%MP;
}
for(int j=;j<=;j++)All|=BL.p[j];
for(int j=;All;j++,All>>=)
if(All&)S[j]=power(,BL.size-);
All=power(,BL.size); ll C=;
for(int j=;j<=;j++,C<<=,C%=MP)
{
sum=(ll)much*(much-)/;
ll t1=now[j]*(All-S[j])%MP;
ll t2=((sum-now[j])*S[j])%MP;
ans=(ans+C*t1+C*t2)%MP;
}
}
int main()
{
read(n),read(m);
for(int i=;i<=m;i++)
{
read(u),read(v),read(t);
add(u,v,t),add(v,u,t);
}
for(int i=;i<=n;i++)
if(!vis[i])//图可能不连通
{
Back();
Dfs(i);
Total();
}
ans+=ans<?MP:;
printf("%d",ans);
return ;
}

Codeforces 724 G Xor-matic Number of the Graph 线性基+DFS的更多相关文章

  1. codeforces 1101G (Zero XOR Subset)-less 前缀异或+线性基

    题目传送门 题意:给出一个序列,试将其划分为尽可能多的非空子段,满足每一个元素出现且仅出现在其中一个子段中,且在这些子段中任取若干子段,它们包含的所有数的异或和不能为0. 思路:先处理出前缀异或,这样 ...

  2. Intel Code Challenge Final Round (Div. 1 + Div. 2, Combined) G - Xor-matic Number of the Graph 线性基好题

    G - Xor-matic Number of the Graph 上一道题的加强版本,对于每个联通块需要按位算贡献. #include<bits/stdc++.h> #define LL ...

  3. Codeforces.724G.Xor-matic Number of the Graph(线性基)

    题目链接 \(Description\) 给定一张带边权无向图.若存在u->v的一条路径使得经过边的边权异或和为s(边权计算多次),则称(u,v,s)为interesting triple(注意 ...

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

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

  5. 【BZOJ-2115】Xor 线性基 + DFS

    2115: [Wc2011] Xor Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 2142  Solved: 893[Submit][Status] ...

  6. 2115: [Wc2011] Xor (线性基+dfs)

    2115: [Wc2011] Xor Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 5714  Solved: 2420 题目链接:https://w ...

  7. 【BZOJ2115】[Wc2011] Xor 高斯消元求线性基+DFS

    [BZOJ2115][Wc2011] Xor Description Input 第一行包含两个整数N和 M, 表示该无向图中点的数目与边的数目. 接下来M 行描述 M 条边,每行三个整数Si,Ti ...

  8. bzoj 2115: [Wc2011] Xor【线性基+dfs】

    -老是想到最长路上 其实可以这样:把每个环的xor和都存起来,然后任选一条1到n的路径的xor和ans,答案就是这个ans在环的线性基上跑贪心. 为什么是对的--因为可以重边而且是无相连通的,并且对于 ...

  9. 【题解】 bzoj2115: [Wc2011] Xor (线性基+dfs)

    bzoj2115,戳我戳我 Solution: 看得题解(逃,我太菜了,想不出这种做法 那么丢个链接 Attention: 板子别写错了 又写错了这次 \(long long\)是左移63位,多了会溢 ...

随机推荐

  1. linux&nbsp;ip地址自动获取,ip地址…

    linux ip地址自动获取,ip地址手动设置(图文解释) 2011-04-19 16:19:31| 分类: 服务器(appache/n | 标签: |字号大中小 订阅 linux ip地址自动获取( ...

  2. MySQL 之 导出导入数据

    导出数据库(sql脚本)  mysqldump -u 用户名 -p 数据库名 > 导出的文件名mysqldump -u root -p --databases db_name > test ...

  3. WP之样式

    1.定义资源 <Window.Resources> <!--下面用样式--> <Style x:Key="BigFontButtonStyle"> ...

  4. C++开源库(一) ----libConfig详解

    博主天生患有蛋疼疾病,写博不易,转载注明出处http://www.cnblogs.com/liboBlog/,谢谢! 在写程序的时候必不可少的一个部分就是conf文件的解析,但是如果自己解析的话会比较 ...

  5. ms sql server line feed

    多行文本换行: SET ANSI_NULLS ON GO SET QUOTED_IDENTIFIER ON GO -- ======================================== ...

  6. JS中map list 数组的迭代

    后台传给前台一个map 前台如何迭代呢 $.post("getSys.jhtml", function(data){ var temp = ""; $.each ...

  7. es6常用方法

    一.let 和 constlet 声明变量,只在所在的块区有效,不存在变量提升:var 存在变 量提升const 声明常量,只在所在块区有效 二.变量的解构赋值1.数组的解构赋值let [a, b, ...

  8. 51nod1105(二分)

    题目链接:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1105 题意:中文题诶- 思路:直接二分答案,再通过二分找有多少 ...

  9. —Libre#2009. 「SCOI2015」小凸玩密室

    #2009. 「SCOI2015」小凸玩密室 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据   题目描述 ...

  10. 洛咕11月月赛部分题解 By cellur925

    听说是你谷史上最水月赛?我不听我最菜 T1:终于结束的起点 月天歌名好评 给你一个模数 \(M\),请你求出最小的 \(n > 0\),使得\(fib(n)\) \(mod\) \(m=0\), ...