『追捕盗贼 Tarjan算法』

<更新提示>

<第一次更新>

---

<正文>

追捕盗贼(COCI2007)

Description

为了帮助警察抓住在逃的罪犯，你发明了一个新的计算机系统。警察控制的区域有N个城市，城市之间有E条双向边连接，城市编号为1到N。

警察经常想在罪犯从一个城市逃亡另一个城市的过程中抓住他。侦查员在仔细研究地图，以决定在哪个城市设置障碍，或者断掉某条路。你的计算机系统必须回答以下两种问题：

1、 如果连接城市G1和G2的路被封掉，罪犯能否从城市A逃到城市B？

2、 如果城市C被封掉，罪犯能否从城市A逃到城市B？

编程实现上述计算机系统。

Input Format

输入第一行包含两个整数N和E(2<=N<=100000,1<=E<=500000)，表示城市和边的数量。

接下来E行，每行包含两个不同的整数Ai和Bi，表示城市Ai和Bi之间有一条边直接相连，任意两个城市之间最多只有一条边相连。

接下来一行包含一个整数Q(1<=Q<=300000)，表示询问次数。

接下来Q行，每行包含4个或5个整数，第一个数为1或2，表示询问问题的种类。

如果问题种类是1，后面跟着4个整数A,B,G1,G2，如上描述表示询问如果G1和G2之间的路封掉罪犯能否从城市A逃到城市B，保证A和B不同，G1和G2之间一定存在路。

如果问题种类是2，后面跟着三个整数A,B,C，三个数互不相同，表示询问如果封掉城市C，罪犯能否从城市A逃到城市B。

输入数据保证一开始任意两个城市都可以相互到达。

Output Format

每个询问输出一行“yes”或“no”。

Sample Input

13 15
1 2
2 3
3 5
2 4
4 6
2 6
1 4
1 7
7 8
7 9
7 10
8 11
8 12
9 12
12 13
5
1 5 13 1 2
1 6 2 1 4
1 13 6 7 8
2 13 6 7
2 13 6 8

Sample Output

yes
yes
yes
no
yes

解析

这是一道\(Tarjan\)算法的综合运用题，考察了了有关割点割边等内容。

由于询问的数量较多，我们考虑离线的\(Tarjan\)算法，将原图化为一棵搜索树。

先看一下第一问，和割边的格式有些相似。我们先考虑封掉的边是不是割边。如果不是割边，则对连通性没有影响。如果是割边，我们可以利用如下方法判定。

设封掉的边连接的两个点为\(G1\)，\(G2\)。(深度较深的为\(G1\))

1.点\(A\)在以\(G1\)为根的子树中，但点\(B\)不在以\(G1\)为根的子树中。

2.点\(B\)在以\(G1\)为根的子树中，但点\(A\)不在以\(G1\)为根的子树中。

由于一条割边只能将原图分为两个不同的联通块，所以当两个条件满足一个时，割边必然将\(A\)，\(B\)两点分在了不同的联通块中。

割边我们可以在做\(Tarjan\)时求出，考虑如何判断某个点\(P1\)是否在以点\(P2\)为根的子树中。

我们重新记录一个值\(Maxdfn_x\)代表以x为根的子树中的最大\(dfn\)值。由于一个节点的\(dfn\)值必然比他的祖先节点的\(dfn\)值小，我们就可以利用\(Maxdfn\)数组得到\(subtree\)函数。

inline bool subtree(int p1,int p2)//代表p1是否在以p2为根的子树中
{
	return dfn[p1]>=dfn[p2]&&dfn[p1]<=Maxdfn[p2];
}

一遍\(Tarjan\)就可以解决第一问，考虑如何解决第二问。

我们仍然要求\(C\)点为割点，由于一个割点会将图分为若干个连通分量，我们需要做特殊的记录。对于每一个\(2\)号询问，设\(yA\)代表点\(C\)到点\(A\)路径上的第一个点的编号，\(yB\)代表点\(C\)到点\(B\)路径上的第一个点的编号。

如果\(A\)，\(B\)不在同一条路径上，即\(yA≠yB\)，并且\(yA\)没有能够向上跳超过\(C\)点的返祖边，或者\(yB\)没有能够向上跳超过\(C\)点的返祖边\(B\)，那么\(A\)，\(B\)必然不能互达。

只要求出\(yA\)，\(yB\)就可以了，可以用如下方法：

设\(askA_x\)代表对于某个\(2\)号询问中\(A\)点为\(x\)的询问编号，\(askB_x\)代表对于某个\(2\)号询问中\(B\)点为\(x\)的询问编号，\(askCa_x\)代表某个\(2\)号询问中\(A\)点要求封掉的点\(C\)编号为\(x\)的询问编号，\(askCb_x\)代表某个\(2\)号询问中\(B\)点要求封掉的点\(C\)编号为\(x\)的询问编号。

对于\(Tarjan\)访问到每一个点x时，我们将\(askA_x,askB_x\)一一取出，并以\(C\)节点作为下标将询问编号统计加入\(askCa,askCb\)数组中，在回溯到父节点时，我们就可以方便地统计得到每一个询问的\(yA\)，\(yB\)值。

这样问题就得到了解决。

\(Code:\)

#include<bits/stdc++.h>
using namespace std;
const int N=100000+200,E=500000+200,Q=300000+200;
int n,m,q,Last[E*2],t,dfn[N],low[N],Maxdfn[N],cnt=0,yA[Q],yB[Q];
struct edge{int ver,next;}e[E*2];
struct question{int a,b,c,x,y,index;}ques[Q];
vector < int > askA[N],askB[N],askCa[N],askCb[N];
inline void insert(int x,int y)
{
	e[++t].ver=y;e[t].next=Last[x];Last[x]=t;
}
inline void input(void)
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int u,v;
		scanf("%d%d",&u,&v);
		insert(u,v);
		insert(v,u);
	}
	scanf("%d",&q);
	for(int i=1;i<=q;i++)
	{
		int index;
		scanf("%d",&ques[i].index);
		if(ques[i].index==1)
			scanf("%d%d%d%d",&ques[i].a,&ques[i].b,&ques[i].x,&ques[i].y);
		if(ques[i].index==2)
		{
			scanf("%d%d%d",&ques[i].a,&ques[i].b,&ques[i].c);
			askA[ques[i].a].push_back(i);
			askB[ques[i].b].push_back(i);
		}
	}
}
inline void Tarjan(int x,int fa)
{
	dfn[x]=low[x]=++cnt;
	while(askA[x].size())
	{
		int num=askA[x].back();
		if(dfn[ques[num].c])
			askCa[ques[num].c].push_back(num);
		askA[x].pop_back();
	}
	while(askB[x].size())
	{
		int num=askB[x].back();
		if(dfn[ques[num].c])
			askCb[ques[num].c].push_back(num);
		askB[x].pop_back();
	}
	for(int i=Last[x];i;i=e[i].next)
	{
		int y=e[i].ver;
		if(y==fa)continue;
		if(!dfn[y])
		{
			Tarjan(y,x);
			low[x]=min(low[x],low[y]);
			while(askCa[x].size())
			{
				yA[askCa[x].back()]=y;
				askCa[x].pop_back();
			}
			while(askCb[x].size())
			{
				yB[askCb[x].back()]=y;
				askCb[x].pop_back();
			}
		}
		else low[x]=min(low[x],dfn[y]);
	}
	Maxdfn[x]=cnt;
}
inline bool subtree(int p1,int p2)
{
	return dfn[p1]>=dfn[p2]&&dfn[p1]<=Maxdfn[p2];
}
inline void answer(void)
{
	for(int i=1;i<=q;i++)
	{
		if(ques[i].index==1)
		{
			if(ques[i].x==ques[i].y)printf("yes\n");
			else
			{
				if(dfn[ques[i].x]<dfn[ques[i].y])swap(ques[i].x,ques[i].y);
				if(  (low[ques[i].x]==dfn[ques[i].x])  &&  (subtree(ques[i].a,ques[i].x)^subtree(ques[i].b,ques[i].x))  )
					printf("no\n");
				else printf("yes\n");
			}
		}
		if(ques[i].index==2)
		{
			if(yA[i]==yB[i])
				printf("yes\n");
			else
			if( (yA[i]&&low[yA[i]]>=dfn[ques[i].c]) || (yB[i]&&low[yB[i]]>=dfn[ques[i].c]) )
				printf("no\n");
			else printf("yes\n");
		}
	}
}
int main(void)
{
	freopen("data.in","r",stdin);
	freopen("data.out","w",stdout);
	input();
	for(int i=1;i<=n;i++)
		if(!dfn[i])Tarjan(i,0);
	answer();
	return 0;
}

---

<后记>