tarjan求割点 poj 1144

最新推荐文章于 2020-05-31 19:14:59 发布
最新推荐文章于 2020-05-31 19:14:59 发布
阅读量2k 收藏
点赞数
分类专栏： 图论 文章标签： input
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/huyuncong/article/details/6415915
版权

tarjan可真是多才多艺，既求强连通，又求割点割边，还求lca，居然每一种都相差不多。

求割点主要是利用dfs遍历树（李顶龙童鞋一直不理解tarjan，从dfs树理解应该比较直观）。

rel[x]记录x实际在树中深度，low[x]记录x及其子树在树中回边可到达的最浅深度。

如果u是割点，那么在dfs树中

1、u为根，u有一个以上子树（删除u后显然两子树不连通）。

2、u非叶子（易知叶子无关紧要），且low[u的孩子]>=rel[u]（即该子树无法不通过u到达u的祖先）

下面是代码（跑不到0ms）

 

 

附缩点双联通分量

每次点入栈，把回边和父亲边加进去

hdu3686

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
int tail[200000],next[2000000],sora[2000000],id[2000000],po[2000000];
int st[200000],f[50000][18],g[50000][18],rel[200000],low[200000],d[200000],flag[200000];
int top,n,m,ss,w_time,block,b[200000];
void link(int x,int y,int i)
{
	++ss,next[tail[x]]=ss,tail[x]=ss,sora[ss]=y,next[ss]=0,id[ss]=i;
	++ss,next[tail[y]]=ss,tail[y]=ss,sora[ss]=x,next[ss]=0,id[ss]=i;
	po[ss]=ss-1,po[ss-1]=ss;
}
void link2(int x,int y)
{
//	cout<<x<<' '<<y<<endl;
	++ss,next[tail[x]]=ss,tail[x]=ss,sora[ss]=y,next[ss]=0;
	++ss,next[tail[y]]=ss,tail[y]=ss,sora[ss]=x,next[ss]=0;
}
void dfs(int x,int y)
{
	++w_time;
	rel[x]=low[x]=w_time;
	for (int i=x,ne;next[i];) {
		i=next[i],ne=sora[i];
		if (rel[ne]) st[++top]=id[i]+n;
	}
	st[++top]=x;
	for (int i=x,ne;next[i];) {
		i=next[i],ne=sora[i];
		if (po[i]==y) continue;
		if (rel[ne]) low[x]=min(low[x],rel[ne]);
		else {
			int tmp=top;
			dfs(ne,i);
			low[x]=min(low[x],low[ne]);
			if (low[ne]>=rel[x]) {
				flag[x]=flag[x+n+n]=1;
				++block;
				//cout<<x<<' '<<top<<' '<<tmp<<endl;
				for (;top>tmp;top--) {
					b[st[top]]=block;
//					cout<<st[top]<<endl;
					if (st[top]<=n && flag[st[top]]) {
						link2(block,n+n+st[top]);
//						cout<<st[top]<<' '<<block<<endl;
					}
				}
				link2(block,n+n+x);
				top=tmp;
			}
		}
	}
}
void bfs(int s,int fl)
{
	int h,r,ne,na;
	h=r=0;
	int k=(int)(log((double)n+n+n)/log(2.0))+1;
	st[r=1]=s,rel[s]=fl,d[s]=0;
	for (;h<r;) {
		ne=st[++h];
		for (int i=ne;next[i];) {
			i=next[i],na=sora[i];
			if (!rel[na]) {
				rel[na]=fl,d[na]=d[ne]+1;
				st[++r]=na;
				f[na][0]=ne,g[na][0]=flag[na];
			}
		}
	}
	for (int j=1;j<=k;j++)
		for (int i=1;i<=r;i++) 
			f[st[i]][j]=f[f[st[i]][j-1]][j-1],
				g[st[i]][j]=g[st[i]][j-1]+g[f[st[i]][j-1]][j-1];
	///cout<<r<<endl;
	//for (int i=1;i<=r;i++) cout<<st[i]<<' '<<f[st[i]][0]<<endl;cout<<endl;
	//cout<<f[7][0]<<endl;
}
int ask(int x,int y)
{
	int e;
	//cout<<x<<' '<<y<<endl;
	if (d[x]<d[y]) e=x,x=y,y=e;
	e=d[x]-d[y];
	int ans=0;
	for (int i=0;e;e>>=1,i++)
		if (e&1) ans+=g[x][i],x=f[x][i];
	//cout<<x<<' '<<y<<endl;
	if (x==y) return ans+flag[x];
	int k=0;
	for (;;) {
		if (f[x][k]==f[y][k]) {
			if (!k) break;
			k--;
		}
		else {
			ans+=g[x][k]+g[y][k];
			x=f[x][k],y=f[y][k];
		}
	}
	ans+=g[x][k]+g[y][k];
	x=f[x][k],y=f[y][k];
	return ans+flag[x];
}
void origin()
{
	ss=n+n+n;
	for (int i=1;i<=n+n+n;i++) tail[i]=i,next[i]=0;
}
int l[200000],r[200000];
int main()
{
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
	for (;;) {
		scanf("%d%d",&n,&m);
		if (!n && !m) break;
		origin();
		for (int i=1;i<=m;i++) {
			int x,y;
			scanf("%d%d",&x,&y);
			l[i]=x,r[i]=y;
			link(x,y,i);
		}
		for (int i=1;i<=n+n+n;i++) rel[i]=low[i]=flag[i]=d[i]=0;
		for (int i=1;i<=n+m;i++) b[i]=0;
		int k=(int)(log((double)n)/log(2.0))+1;
		for (int i=1;i<=n+n+n;i++)
			for (int j=0;j<=k;j++) f[i][j]=g[i][j]=0;
		top=0,w_time=0,block=n;
		for (int i=1;i<=n;i++)
			if (!rel[i]) dfs(i,0);
	//for (int i=1;i<=n+n+n;i++) cout<<flag[i]<<' ';cout<<endl;
		//for (int i=1;i<=n+m;i++) cout<<b[i]<<' ';cout<<endl;
		for (int i=n+1;i<=n+n;i++) 
			if (!rel[i]) bfs(i,i);
		//for (int i=1;i<=n+n+n;i++) cout<<f[i][0]<<' ';cout<<endl;
		int q;
		scanf("%d",&q);
		for (int i=1;i<=n;i++)
			if (flag[i]) b[i]=i+n+n;
		for (;q;q--) {
			int x,y;
			scanf("%d%d",&x,&y);
			if ((b[x+n]==b[y+n]) || (rel[b[x+n]]!=rel[b[y+n]])) printf("0\n");
			else {
				int ans=ask(b[x+n],b[y+n]);
/*				ans=min(ans,ask(b[l[x]],b[l[y]])-flag[l[x]]);
				ans=min(ans,ask(b[r[x]],b[l[y]])-flag[r[x]]);
				ans=min(ans,ask(b[l[x]],b[r[y]])-flag[l[x]]);
				ans=min(ans,ask(b[r[x]],b[r[y]])-flag[r[x]]);*/
				printf("%d\n",ans);
			}
		}
	}
	return 0;
}


拆边为点，考虑地更少

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cmath>
using namespace std;
int ss,cnt,n,m,w_time,top,K,q;
int st[250000],rel[250000],low[250000],d[250000],tail[250000];
int b[250000],w[250000];
int f[250000][20],g[250000][20];
int next[2000000],sora[2000000];
vector <int> vec[250000];
void link(int x,int y)
{
    ++ss,next[tail[x]]=ss,tail[x]=ss,sora[ss]=y,next[ss]=0;
    ++ss,next[tail[y]]=ss,tail[y]=ss,sora[ss]=x,next[ss]=0;
}
void link2(int x,int y)
{
    vec[x].push_back(y);
    vec[y].push_back(x);
}
void dfs(int x,int y)
{
    ++w_time;
    rel[x]=low[x]=w_time;
    st[++top]=x;
    for (int i=x,ne;next[i];) {
        i=next[i],ne=sora[i];
        if (ne==y) continue;
        if (rel[ne]) low[x]=min(low[x],rel[ne]);
        else {
            int tmp=top;
            dfs(ne,x);
            low[x]=min(low[x],low[ne]);
            if (low[ne]>=rel[x]) {
                b[x]=x;
                if (x<=n) w[x]=1;
                b[++cnt]=cnt;
                link2(x,cnt);
                for (;top!=tmp;top--) {
                    int ne=st[top];
                    if (!b[ne]) b[ne]=cnt;
                    link2(ne,cnt);
                }
            }
        }
    }
}
void dfs2(int x,int s,int dis)
{
    low[x]=s;
    rel[x]=w_time;
    d[x]=dis;
    for (int i=0;i<vec[x].size();i++) {
        int ne=vec[x][i];
        if (rel[ne]!=w_time) {
            f[ne][0]=x,g[ne][0]=w[x];
            dfs2(ne,s,dis+1);
        }
    }
}
void origin()
{
    ss=n+m;
    for (int i=1;i<=n+m;i++) tail[i]=i,next[i]=0;
    for (int i=1;i<=n+m;i++) rel[i]=low[i]=0;
    w_time=0,top=0;
    for (int i=1;i<=cnt;i++) vec[i].clear();
    for (int i=1;i<=cnt;i++) w[i]=0,b[i]=0;
    cnt=n+m;
}
int ask(int x,int y)
{
    if (low[x]!=low[y]) return 0;
    if (d[x]<d[y]) swap(x,y);
    int e=d[x]-d[y];
    int sum=0;
    for (int b=0;e;e>>=1,b++)
        if (e&1) sum+=g[x][b],x=f[x][b];
    if (x==y) return sum;
    for (int i=K;i>=0;i--)
        if (f[x][i]!=f[y][i]) {
            sum+=g[x][i]+g[y][i];
            x=f[x][i],y=f[y][i];
        }
    if (x!=y) sum+=g[x][0];
    return sum;
}
int main()
{
    for (;scanf("%d%d",&n,&m)==2;) {
        if (!n && !m) break;
        origin();
        for (int i=1;i<=m;i++) {
            int x,y;
            scanf("%d%d",&x,&y);
            link(x,i+n);
            link(i+n,y);
        }
        for (int i=1;i<=n+m;i++)
            if (!rel[i]) dfs(i,0);
        ++w_time;
        for (int i=1;i<=cnt;i++) 
            if (rel[i]!=w_time) f[i][0]=0,g[i][0]=0,dfs2(i,i,0);
        K=(int)(log((double)cnt)/log(2.0))+1;
        for (int j=1;j<=K;j++)
            for (int i=1;i<=cnt;i++)
                f[i][j]=f[f[i][j-1]][j-1],
                    g[i][j]=g[i][j-1]+g[f[i][j-1]][j-1];
        scanf("%d",&q);
//        cout<<cnt<<endl;
//        cout<<f[7][0]<<' '<<f[8][0]<<' '<<f[6][0]<<endl;
        for (int i=1;i<=q;i++) {
            int x,y;
            scanf("%d%d",&x,&y);
            x+=n,y+=n;
            int ans=ask(x,y);
            printf("%d\n",ans);
        }
    }
    return 0;
}


huyuncong
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
Tarjan算法求解一个无向图中的割点和桥问题
好记性不如烂笔头的专栏
10-27 5041
基本概念割点：Articulation Point 在无向连通图中，删除一个顶点v及其相连的边后，原图从一个连通分量变成了两个或多个连通分量，则称顶点v为割点，同时也称关节点(Articulation Point)。 双连通的图：一个没有关节点的连通图称为重连通图(biconnected graph)（双连通图）。 连通度：k，若在连通图上至少删去k 个顶点才能破坏图的连通性。算法应用算法应用
ACM快速输入输出
王跃坤的博客
11-17 1392
package 输入输出; import java.io.BufferedReader; import java.io.IOException; import java.io.InputStreamReader; import java.io.OutputStreamWriter; import java.io.PrintWriter; import java.io.StreamTokenize...
Tarjan算法超详细讲解(割点割边强连通)
汪阿少的博客
02-06 6481
本文出自我的掘金博客, 欢迎大家访问传送门 今天我主要介绍Tarjan算法在割点割边以及强连通分量中的应用以及缩点技巧 按照老规矩, 先上两道模板题 【模板】强连通分量 【模板】割点（割顶） 割点割边 一, 离散数学中的定义: 割点: 无向连通图中，去掉一个顶点及和它相邻的所有边，图中的连通分量数增加，则该顶点称为割点。 割边: 无向联通图中，去掉一条边，图中的连通分量数增加，则这条边，称为...
割点（Tarjan算法）
dengshouzi7943的博客
06-07 896
本文可转载，转载请注明出处：www.cnblogs.com/collectionne/p/6847240.html。本文未完，如果不在博客园（cnblogs）发现此文章，请访问以上链接查看最新文章。 前言：之前翻译过一篇英文的关于割点的文章（英文原文、翻译），但是自己还有一些不明白的地方，这里就再次整理了一下。有兴趣可以点我给的两个链接。 割点的概念 在无向连通...
超详细Tarjan算法总结，求强连通分量，割点，割边，有重边的割边
旺旺的博客
06-14 4466
Tarjan是一个人，他一身中发明了很多算法，就这几个算法最为出名。 1、求有向图的强连通分量，那么什么是强连通分量呢，就是一个顶点集合，任意两个顶点间都可以互相到达。如果图中任意两点可以互相到达，则此图强连通。 TARJAN是基于dfs算法的基础上，所以也会得到一颗搜索树。 那么他是怎么运行的呢，首先大家看着两个数组  dfn[i],和low[i] dfn[i]表示第i个节点深
「洛谷3469」「POI2008」BLO-Blockade【Tarjan求割点】
chhokmah的博客
04-21 167
题目链接 【洛谷传送门】 题解 很显然，当这个点不是割点的时候，答案是2∗(n−1)2*(n-1)2∗(n−1) 如果这个点是割点，那么答案就是两两被分开的联通分量之间求组合数。 代码 #include <bits/stdc++.h> #define ll long long using namespace std; const int N = 500005; struct edge...
图论- 图的连通性- Tarjan 求割点与桥.rar
09-16
图论- 图的连通性- Tarjan 求割点与桥.rar
tarjan求割边
07-15
Tarjan算法求割边 Tarjan算法是图论中的一种常见算法，用于寻找有向图中的割点和割边。割边是指删除该边后，图将被分割成两个或多个连通分量的边。 Tarjan算法可以快速地找到所有的割边。 割边的定义与割点的定义...
Tarjan算法求强连通分量
02-03
使用Tarjan算法进行快速计算强连通分量，C++语言实现。
Tarjan算法讲义
07-14
Tarjan 算法是图论中非常实用 / 常用的算法之一，能解决强连通分量，双连通分量，割点和桥，求最近公共祖先（LCA）等问题。 关于 Tarjan 算法，笔者将用一系列文章系统介绍 Tarjan 算法的原理以及其主要解决的问题...
Tarjan 的 LCA 算法
03-19
c++写的Tarjan 的 LCA 算法，最近公共祖先算法，可供算法学习参考
tarjan算法总结 (强连通分量+缩点+割点)，看这一篇就够了~
热门推荐
csyifanZhang的博客
04-08 1万+
文章目录一、tarjan求强连通分量1：算法流程2：模板二、tarjan缩点1：相关定义2：算法流程三、tarjan求割点、桥1、什么是割点2.割点怎么求？割点tarjan模板3.什么是桥4：桥怎么求四、双联通分量（for 无向图） tarjan可以做什么？ 根据 Robert Tarjan 的名字命名的算法Tarjan算法可以在线性时间内求出无向图的割点与桥，再进一步的求出双联通分量，也在数据结...
tarjan 求割点/割边模板
hxc2101的博客
09-08 2122
求割点/割边： //fa为当前结点的父亲结点，dep记录当前深度，初始化为0 //dfn为搜索树中的次序号，初始化为-1 //每个点被访问时都会将dfn置为当前dep，注意是 ++dep，不是dep++，low数组随着搜索更新 vector&lt;int&gt;edg[maxn]; int dfn[maxn],low[maxn]; int dep=0,child=0; void tarjan(...
Tarjan割点
weixin_44959460的博客
04-25 1009
[定义]： 在无向连通图中，如果将其中一个点以及所有连接该点的边去掉，图就不再连通，那么这个点就叫做割点。 在无向连通图中，如果将其中一条边去掉，图就不再连通，那么这条边叫桥。 [思想]： 割点：1）当前节点u为树根的时候，这个节点要有至少两棵子树，否则去掉该节点不会有影响。 2）当前节点u不是树根的时候，要满足low[v]>=dfn[u]，注意u是v的父亲 ...
关于tarjan算法的一些理解（割边和割点）
llin-黎辰
07-19 3279
首先介绍以下tarjan算法： ---------------------------------------------tarjan算法----------------------------------------------------------- tarjan算法思想： tarjan算法在离线求LCA,强连通分量，割边，割点，点双连通分量，边双连通分量很有用 tarjan算法中主
tarjan算法求割点(java)
sidnee的博客
05-31 498
tarjan算法过程 思路很简单，如果某个dfs树里的点，它的所有后代的回边都不能连到这个点或这个点的上面，那么它就是割点。所以用dfn记录dfs序，然后low记录后代中回边里最优秀(连的点最靠上)的那个回边所到达的点。所以一个点是割点当且仅当满足下面条件之一： (1) u为树根，且u有多于一个子树。 (2) u不为树根，且满足存在(u,v)为树枝边(或称父子边， 即u为v在搜索树中的父亲)，使得dfn(u)<=low(v)。 另一种想法 也可以通过每个点必经点集的并集来确定割点，必经点集可以用逐
【模板】最小表示法
zxyoi_dreamer的博客（不定期诈尸）
10-06 588
参考题目：洛谷P1368 解析： 联赛完后统一更所有模板题题解 代码： #include&lt;bits/stdc++.h&gt; using namespace std; #define ll long long #define re register #define gc getchar #define pc putchar #define cs const inline int get...
tarjan算法求割点割边
LanQiLi的博客
12-14 3079
在上一节我们已经知道tarjan算法可以求联通图，在这里我们也运用tarjan的思想求割点与割边，首先我们先来说说割点，那么什么事割点呢，先来看一张图（a）,图片来自网络  在（a）图中，我们将A点以及与A点相连的边全部去除，会发现这个联通图被分成了俩个联通图，一个是节点F，另外一个是余下的所有的节点组成的图，因此我们将A点称为割点，同理我们发现B点也是割点，因此我们可以这样定义割点，在一个...
[BZOJ4651][Noi2016]网格（Tarjan 求割点）
女装的你，如此好看！
04-01 242
Address 洛谷 P1173 BZOJ 4651 UOJ #220 LOJ #2084 Solution 下面把跳蚤视为白点，蛐蛐视为黑点 可以发现答案只可能是 −1-1−1 、 000 、 111 、 222 四种 然后随机输出一个，你可以以 1 / (4 ^ T) 的高概率通过此题 易得， 如果没有黑点，则有三种情况 （1） n×m≤2n\times m\le 2n×m≤2 ： −1...
Tarjan求LCA
最新发布
07-27
Tarjan算法是一种用于解决最近公共祖先（LCA）问题的离线算法。离线算法指的是在读取所有查询之后一次性计算所有查询的答案，而不是每读取一个查询就计算一次。\[1\] 在Tarjan算法中，需要使用并查集来实现。并查集是一种数据结构，用于维护元素之间的集合关系。下面是一个并查集的模板代码： ```cpp int fa\[100000\]; void reset(){ for (int i=1;i<=100000;i++){ fa\[i\]=i; } } int getfa(int x){ return fa\[x\]==x?x:getfa(fa\[x\]); } void merge(int x,int y){ fa\[getfa(y)\]=getfa(x); } ``` 在Tarjan算法的伪代码中，首先标记当前节点为已访问状态。然后遍历当前节点的子节点，递归调用Tarjan函数并合并子节点。接下来，遍历与当前节点有查询关系的节点，如果该节点已经访问过，则输出当前节点和该节点的LCA（通过并查集的查找函数getfa获取）。\[3\] 以上是关于Tarjan算法求解LCA的相关内容。 #### 引用[.reference_title] - *1* [Tarjan 算法解决 LCA 问题](https://blog.csdn.net/chengqiuming/article/details/126878817)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [详解使用 Tarjan 求 LCA 问题（图解）](https://blog.csdn.net/weixin_34315485/article/details/93801193)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值