BZOJ 压力 tarjan 点双联通分量+树上差分+圆方树
题意
如今,路由器和交换机构建起了互联网的骨架。处在互联网的骨干位置的核心路由器典型的要处理100Gbit/s的网络流量。
他们每天都生活在巨大的压力之下。小强建立了一个模型。这世界上有N个网络设备,他们之间有M个双向的链接。这个世界是连通的。
在一段时间里,有Q个数据包要从一个网络设备发送到另一个网络设备。
一个网络设备承受的压力有多大呢?很显然,这取决于Q个数据包各自走的路径。
不过,某些数据包无论走什么路径都不可避免的要通过某些网络设备。
你要计算:对每个网络设备,必须通过(包括起点、终点)他的数据包有多少个?
对于40%的数据,N,M,Q≤2000
对于60%的数据,N,M,Q≤40000
对于100%的数据,N≤100000,M,Q≤200000
分析
这道题很有意思,如果它要是问可能通过的点有多少个,那问题就麻烦多了,所以它问必须经过的点,对于一个点来说,如果它不是割点,那么它一定不是必经的,除非它是端点,所以我们想到了什么?对点双联通分量,但是如果全缩了还是有问题,必须经过的点怎么办?如果变换一下这个图,让它两个点之间有唯一的路径,是不是就很好办了,对于唯一路径,我们很容易想到树,而又不能全缩了成为树,那就势必要用到别的树。
所以引入一个新的树——圆方树
这个树其实很好理解,就是原图的每个点都是圆点,方点是什么呢?一个方点对应一个点双联通分量,每个点双上的点连到这个这个方点上,这样就形成了一棵树,然后怎么办呢?如果对于每次修改都dfs一遍,显然效率是不高的,而每次发送一回数据包,就是将这一个经过的路径区间val+1,涉及区间加减,不就是差分数组嘛,所以直接u++,v++,lca--,falca--(点差分是falca--因为lca的值也要更新,边差分是lca-=2因为lca上边连着父亲节点的边是不改变的),最后直接dfs跑一遍统计答案就ok。
这里特别强调一个坑,点双连树的时候不要直接pop到割点,而是要pop到E.to,看起来是差不多的,但是在割点到E.to之间,还可能会有点,这里说的有点不是图上边,而是栈里边有没pop的值,点没有pop的原因是没有判断到割点,而那些点是显然不是这个点双里边的,连进来就是,WA,所以这里要特别注意一下。
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=4e5+;
struct Edge{
int to,next;
}e[N],E[N];
int Head[N],len;
void Ins(int a,int b){
e[++len].to=b;e[len].next=Head[a];Head[a]=len;
}
int H[N],l;
void I(int a,int b){
E[++l].to=b;E[l].next=H[a];H[a]=l;
}
int low[N],dfn[N],stk[N],num,top,fang;
void tarjan(int u){
dfn[u]=low[u]=++num;
stk[++top]=u;
for(int x=Head[u];x;x=e[x].next){
int v=e[x].to;
if(!dfn[v]){
tarjan(v);
low[u]=min(low[u],low[v]);
if(low[v]>=dfn[u]){
fang++;int temp;
do{
temp=stk[top--];
I(temp,fang);I(fang,temp);
}while(temp!=v);//是这里不能pop到u
I(fang,u);I(u,fang);
}
}else low[u]=min(low[u],dfn[v]);
}
}
int p[N][],dep[N];
//跑lca
void dfs(int x){
for(int i=;p[x][i];i++)
p[x][i+]=p[p[x][i]][i];
for(int i=H[x];i;i=E[i].next){
int v=E[i].to;
if(v!=p[x][]){
p[v][]=x;
dep[v]=dep[x]+;
dfs(v);
}
}
}
int Lca(int u,int v){
if(dep[u]<dep[v])swap(u,v);
int d=dep[u]-dep[v];
for(int i=;d;i++,d>>=)
if(d&)u=p[u][i];
if(u==v)return u;
for(int i=;i>=;i--)//这里又写错了刚开始,i>=0不是i,i->i>0少了个0
if(p[u][i]!=p[v][i])
u=p[u][i],v=p[v][i];
return p[u][];
}
//lca结束
int val[N];
int calc(int u){
for(int i=H[u];i;i=E[i].next){
int v=E[i].to;
if(v!=p[u][])
val[u]+=calc(v);
}
return val[u];
}
int main(){
// freopen("a.txt","r",stdin);
int m,n,q;
scanf("%d%d%d",&n,&m,&q);
fang=n;
for(int i=;i<=m;i++){
int a,b;
scanf("%d%d",&a,&b);
Ins(a,b);Ins(b,a);
}
tarjan();//图是联通的tarjan一遍就好
dfs();
while(q--){
int a,b;
scanf("%d%d",&a,&b);
int lca=Lca(a,b);
val[a]++;val[b]++;
val[lca]--;val[p[lca][]]--;//树上差分
}
calc();
for(int i=;i<=n;i++)
printf("%d\n",val[i]);
}
BZOJ 压力 tarjan 点双联通分量+树上差分+圆方树的更多相关文章
- POJ2942 Knights of the Round Table【Tarjan点双联通分量】【二分图染色】【补图】
LINK 题目大意 有一群人,其中有一些人之间有矛盾,现在要求选出一些人形成一个环,这个环要满足如下条件: 1.人数大于1 2.总人数是奇数 3.有矛盾的人不能相邻 问有多少人不能和任何人形成任何的环 ...
- [J]computer network tarjan边双联通分量+树的直径
https://odzkskevi.qnssl.com/b660f16d70db1969261cd8b11235ec99?v=1537580031 [2012-2013 ACM Central Reg ...
- POJ 3694Network(Tarjan边双联通分量 + 缩点 + LCA并查集维护)
[题意]: 有N个结点M条边的图,有Q次操作,每次操作在点x, y之间加一条边,加完E(x, y)后还有几个桥(割边),每次操作会累积,影响下一次操作. [思路]: 先用Tarjan求出一开始总的桥的 ...
- BZOJ 2959: 长跑 lct 双联通分量 并查集 splay
http://www.lydsy.com/JudgeOnline/problem.php?id=2959 用两个并查集维护双联通分量的编号和合并. #include<iostream> # ...
- POJ3177 Redundant Paths【tarjan边双联通分量】
LINK 题目大意 给你一个有重边的无向图图,问你最少连接多少条边可以使得整个图双联通 思路 就是个边双的模板 注意判重边的时候只对父亲节点需要考虑 你就dfs的时候记录一下出现了多少条连向父亲的边就 ...
- tarjan求双联通分量(割点,割边)
之前一直对tarjan算法的这几种不同应用比较混淆...我太弱啦! 被BLO暴虐滚过来 用tarjan求点双,很多神犇都给出了比较详细的解释和证明,在这里就不讲了(其实是这只蒟蒻根本不会orz) 这里 ...
- [hdu2460]network(依次连边并询问图中割边数量) tarjan边双联通分量+lca
题意: 给定一个n个点m条边的无向图,q个操作,每个操作给(x,y)连边并询问此时图中的割边有多少条.(连上的边会一直存在) n<=1e5,m<=2*10^5,q<=1e3,多组数据 ...
- 【洛谷 SP2878】Knights of the Round Table(双联通分量)
先放这吧,没时间写,明天再补 "明天到了" 题目链接 题意:求不在任何奇环内的点的数量. Tarjan求点双联通分量,然后再染色判断是不是二分图就好了. 只是不懂为什么Tarjan ...
- BZOJ3331 [BeiJing2013]压力[圆方树+树上差分]
圆方树新技能get.具体笔记见图连通性问题学习笔记. 这题求无向图的必经点,这个是一个固定套路:首先,一张连通的无向图中,每对点双和点双之间是以一个且仅一个割点连接起来的(如果超过一个就不能是割点了) ...
随机推荐
- 11--PHP中的类和对象
PHP类和对象 类是面向对象程序设计的基本概念,通俗的理解类就是对现实中某一个种类的东西的抽象, 比如汽车可以抽象为一个类,汽车拥有名字.轮胎.速度.重量等属性,可以有换挡.前进.后退等操作方法. 通 ...
- GDB调试系列之了解GDB
想要熟练利用GDB进行程序调试,首先要了解什么是GDB. 1. 什么是GDB GDB (the GNU Project Debugger) 是一个可以运行在大多数常见的UNIX架构.Windows.M ...
- linux同步当前网络时间
[root@root ~]# yum install -y ntpdate 执行:ntpdate[root@root ~]# ntpdate 120.24.81.91或[root@root ~]# n ...
- Android Base64图片无法长按保存 问题解决
踩了一个巨坑. 目前微信ios/android 均能长按保存src=base64的图片 (微信android x5 专门解决了这个问题); 但是android其他App没有针对解决这个系统问题(姑且 ...
- 【MySQL 原理分析】之 Explain & Trace 深入分析全模糊查询走索引的原理
一.背景 今天,交流群有一位同学提出了一个问题.看下图: 之后,这位同学确实也发了一个全模糊查询走索引的例子: 到这我们可以发现,这两个sql最大的区别是:一个是查询全字段(select *),而一个 ...
- 如何在Flutter中使用flutter_markdown
很多博客,论坛都支持markdown语法,flutter也有支持markdown语法的插件flutter_markdown 安装依赖 dependencies: flutter: sdk: flutt ...
- CJSON的封装API
为了更方便使用C的JSON库,对其进行了一层封装. H文件: #ifndef __JSONHELPER__ #define __JSONHELPER__ #ifdef __cplusplus exte ...
- 什么是Servlet?Servlet的周期和方法
1.什么是Servlet? Servlet是运行在web服务器或应用服务器的程序,它是作为来自web浏览器或其他http客户端的请求和HTTP服务器上的数据库或应用程序之间的中间层! 2.Servl ...
- 微信小程序开发(一)开发准备
微信小程序开发(一)开发准备 微信小程序,近几年,越来越火,应用场景越来越多,好大学习起来,比较容易,学习曲线平缓,你要是会前端开发,那简直可以用易如反掌来形容. 小程序,开发工具官方也提供了,他兼容 ...
- Java 并发同步工具(转)
转自:https://www.jianshu.com/p/e80043ac4115 在 java 1.5 中,提供了一些非常有用的辅助类来帮助我们进行并发编程,比如 CountDownLatch,Cy ...