[hdu2460]network(依次连边并询问图中割边数量) tarjan边双联通分量+lca
题意:
给定一个n个点m条边的无向图,q个操作,每个操作给(x,y)连边并询问此时图中的割边有多少条。(连上的边会一直存在)
n<=1e5,m<=2*10^5,q<=1e3,多组数据。
题解:
用tarjan求边双连通分量并缩点,缩点后组成一棵树,记录此时割边共有sum条。
连接(x,y),设c[i]为缩点后i所在的新点(边双连通分量),则c[x]-->lca-->c[y]形成一个环,环上的所有边都不再是割边,走一遍并标记,如果遇到没标记过的就sum--。
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std; const int N=,M=,K=;
int n,m,al,bl,sum,num,sl,cnt;
int afirst[N],bfirst[N],fa[N][K],dfn[N],low[N],is[N],c[N],s[N],dep[N];
struct node{
int x,y,tmp,next;
}a[*M],b[*M]; int minn(int x,int y){return x<y ? x:y;} void ins_a(int x,int y)
{
a[++al].x=x;a[al].y=y;a[al].tmp=;
a[al].next=afirst[x];afirst[x]=al;
} void ins_b(int x,int y)
{
b[++bl].x=x;b[bl].y=y;
b[bl].next=bfirst[x];bfirst[x]=bl;
} void tarjan(int x)
{
dfn[x]=low[x]=++num;
s[++sl]=x;
for(int i=afirst[x];i;i=a[i].next)
{
if(a[i].tmp) continue;
a[i].tmp=;
a[i%== ? i-:i+].tmp=;
int y=a[i].y;
if(!dfn[y])
{
tarjan(y);
low[x]=minn(low[x],low[y]);
if(dfn[x]<low[y])
{
sum++;//printf("%d -- > %d\n",x,y);
cnt++;
int z;
while()
{
z=s[sl--];
c[z]=cnt;
if(z==y) break;
}
}
}
else low[x]=minn(low[x],dfn[y]);
}
} void dfs(int x)
{
for(int i=bfirst[x];i;i=b[i].next)
{
int y=b[i].y;
if(y==fa[x][]) continue;
fa[y][]=x;dep[y]=dep[x]+;
dfs(y);
}
} void prepare_lca()
{
memset(fa,,sizeof(fa));
memset(dep,,sizeof(dep));
memset(is,,sizeof(is));
dep[]=;
dfs();
for(int j=;j<=;j++)
for(int i=;i<=n;i++)
fa[i][j]=fa[fa[i][j-]][j-];
} int get_lca(int x,int y)
{
if(dep[x]<dep[y]) swap(x,y);
for(int i=;i>=;i--)
if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
if(x==y) return x;
for(int i=;i>=;i--)
{
if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
}
return fa[x][];
} void del(int x,int y,int lca)
{
while(x!=lca)
{
if(!is[x]) sum--;
is[x]=;
x=fa[x][];
}
while(y!=lca)
{
if(!is[y]) sum--;
is[y]=;
y=fa[y][];
}
} int main()
{
freopen("a.in","r",stdin);
int q,x,y,T=;
while()
{
scanf("%d%d",&n,&m);
if(!n && !m) return ;
printf("Case %d:\n",++T);
al=;bl=;sum=;
memset(afirst,,sizeof(afirst));
memset(bfirst,,sizeof(bfirst));
num=;sl=;cnt=;
memset(dfn,,sizeof(dfn));
memset(c,,sizeof(c));
for(int i=;i<=m;i++)
{
scanf("%d%d",&x,&y);
ins_a(x,y);ins_a(y,x);
}
for(int i=;i<=n;i++)
{
if(!dfn[i])
{
tarjan(i);
if(sl)
{
cnt++;
for(int i=;i<=sl;i++) c[s[i]]=cnt;
sl=;
}
}
} for(int i=;i<=*m;i++)
{
x=a[i].x,y=a[i].y;
if(c[x]!=c[y]) ins_b(c[x],c[y]);
}
// for(int i=1;i<=bl;i++)
// printf("%d -- > %d\n",b[i].x,b[i].y);
// for(int i=1;i<=n;i++) printf("c [ %d ] = %d\n",i,c[i]);
prepare_lca();
scanf("%d",&q);
for(int i=;i<=q;i++)
{
scanf("%d%d",&x,&y);
x=c[x];y=c[y];
int lca=get_lca(x,y);
del(x,y,lca);
printf("%d\n",sum);
}
printf("\n");
}
return ;
}
[hdu2460]network(依次连边并询问图中割边数量) tarjan边双联通分量+lca的更多相关文章
- [J]computer network tarjan边双联通分量+树的直径
https://odzkskevi.qnssl.com/b660f16d70db1969261cd8b11235ec99?v=1537580031 [2012-2013 ACM Central Reg ...
- POJ3694 Network —— 边双联通分量 + 缩点 + LCA + 并查集
题目链接:https://vjudge.net/problem/POJ-3694 A network administrator manages a large network. The networ ...
- 图连通性【tarjan点双连通分量、边双联通分量】【无向图】
根据 李煜东大牛:图连通性若干拓展问题探讨 ppt学习. 有割点不一定有割边,有割边不一定有割点. 理解low[u]的定义很重要. 1.无向图求割点.点双联通分量: 如果对一条边(x,y),如果low ...
- poj 3177 Redundant Paths 求最少添加几条边成为双联通图: tarjan O(E)
/** problem: http://poj.org/problem?id=3177 tarjan blog: https://blog.csdn.net/reverie_mjp/article/d ...
- UML类图中的关系和表示方法
类图是用来描述程序中的类以及它们之间的关系的,使用类图可以帮助我们简化对系统的理解.在UML类图中比较常见的关系有六种,它们分别是:依赖.关联.聚合.组合.泛化.实现,这六种关系中类之间的紧密程度是依 ...
- Ex 3_25 图中每个顶点有一个相关价格..._十一次作业
(a)首先对有向无环图进行拓扑排序,再按拓扑排序的逆序依次计算每个顶点的cost值,每个顶点的cost值为自身的price值与相邻顶点间的cost值得最小值 (b)求出图中的每一个强连通分量,并把所有 ...
- Dijkstra 算法,用于对有权图进行搜索,找出图中两点的最短距离
Dijkstra 算法,用于对有权图进行搜索,找出图中两点的最短距离,既不是DFS搜索,也不是BFS搜索. 把Dijkstra 算法应用于无权图,或者所有边的权都相等的图,Dijkstra 算法等同于 ...
- UML类图中的各种箭头代表的含义(转自:http://www.cnblogs.com/damsoft/archive/2016/10/24/5993602.html)
1.UML简介Unified Modeling Language (UML)又称统一建模语言或标准建模语言. 简单说就是以图形方式表现模型,根据不同模型进行分类,在UML 2.0中有13种图,以下是他 ...
- UML类图中的关系表示
UML类图中的关系和表示方法 类图是用来描述程序中的类以及它们之间的关系的,使用类图可以帮助我们简化对系统的理解.在UML类图中比较常见的关系有六种,它们分别是:依赖.关联.聚合.组合.泛化.实现,这 ...
随机推荐
- Internet History
Alan Turing and Bletchley Park Top secret breaking effort(二战破译希特勒密码) 10,000 people at the peak(team ...
- UVALive - 6868 Facility Locations 想法题
题目链接: http://acm.hust.edu.cn/vjudge/problem/88634 Facility Locations Time Limit: 3000MS 题意 给你一个m*n的矩 ...
- lintcode-167-链表求和
167-链表求和 你有两个用链表代表的整数,其中每个节点包含一个数字.数字存储按照在原来整数中相反的顺序,使得第一个数字位于链表的开头.写出一个函数将两个整数相加,用链表形式返回和. 样例 给出两个链 ...
- 安装配置erlang_db_driver
erlang-db-driver是北京融易通公司开源的一个erlang支持众多数据库的一个驱动类库,据其wiki介绍,其支持MySQL, Oracle, Sybase, DB2 and Informi ...
- phpcms开启在线编辑模版 方法
目录:\caches\configs\system.php 将:第20行 'tpl_edit'=> 0 修改为 'tpl_edit'=> 1 (0:默认的,不开启: 1: ...
- 统计VS2013中有效行数
将鼠标放在解决方案处,按下ctrl+shift+F b*[^:b#/]+.*$(带前面的using)^b*[^:b#/]+.*$
- 禁止移动端input弹出软键盘
在做三级联动,或者一些时间插件的时候总是弹出软键盘,用下面的方法就可以禁用掉,废话不多说直接上代码. HTML代码 <div class=""> <div> ...
- 2011 Multi-University Training Contest 6 - Host by JLU
打了4hours,做出一道题...太菜了.rank:45/107 开场看B,题目看不懂...3hours半才发现i<=N-1,不是i<=x-1.然而还是不会. 看到J有人过了,发现是个简单 ...
- BZOJ4737 组合数问题(卢卡斯定理+数位dp)
不妨不管j<=i的限制.由卢卡斯定理,C(i,j) mod k=0相当于k进制下存在某位上j大于i.容易想到数位dp,即设f[x][0/1][0/1][0/1]为到第x位时是否有某位上j> ...
- POJ2891:Strange Way to Express Integers——题解
http://poj.org/problem?id=2891 题目大意: k个不同的正整数a1,a2,...,ak.对于一些非负m,满足除以每个ai(1≤i≤k)得到余数ri.求出最小的m. 输入和输 ...