题目大意:问一个有向图是否任意两点在两个方向上互相连通。
有向图强连通分量定义:如果一个图中的任意两点在两个方向上都互相连通,则该图为强连通图。极大强连通图为有向图的强连通分量(注意是极大,不是最大。一个图会有多个强连通分量)。感性理解,强连通图就是多个环,或者一个点连接在一起所产生的图。
如何求?定义节点cur->Low,cur的子搜索树节点a中如果存在边(a,b),使得b->DfsN小于cur->DfsN,则cur->Low=min foreach b{b->DfsN},否则为cur->DfsN。显然,cur->Low所对应的节点位于cur所在极大强连通分量中,且是cur所在环中深度最浅的。
新定义一个栈维护一组节点,其使得处理完栈内节点cur以上的节点后,cur以上节点的Low <= cur->Low(它们位于cur子搜索树中的各个枝杈上(而不仅仅存在于一个枝杈中,这是系统栈所做不到的),不存在满足该条件却不在该栈内的节点),即cur及以上节点在一个强连通分量内。如果cur->Low == cur->DfsN,则cur是其所在强连通分量中深度最浅的,再往栈下面走的节点就不属于这个连通分量了,此时便可以输出连通分量,然后将cur及以上节点全部弹出。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
const int MAX_NODE = 10010, MAX_EDGE = 100010;
struct Node;
struct Edge;
struct Node
{
int Id, DfsN, Low;
Edge *Head;
}_nodes[MAX_NODE], *Root;
int _vCount, DfsCnt;
struct Edge
{
Node *From, *To;
Edge *Next;
Edge() {}
}*_edges[MAX_EDGE];
int _eCount;
Edge *NewEdge()
{
_eCount++;
return _edges[_eCount] ? _edges[_eCount] : _edges[_eCount] = new Edge();
}
Edge *AddEdge(Node *from, Node *to)
{
Edge *e = NewEdge();
e->From = from;
e->To = to;
e->Next = from->Head;
from->Head = e;
return e;
}
vector<Node*> Stack;
void Init(int vCount)
{
_vCount = vCount;
_eCount = 0;
DfsCnt = 0;
Root = 1 + _nodes;
Stack.clear();
memset(_nodes, 0, sizeof(_nodes));
}
void Dfs(Node *cur)
{
cur->DfsN = cur->Low = ++DfsCnt;
Stack.push_back(cur);
for (Edge *e = cur->Head; e; e = e->Next)
{
if (!e->To->DfsN)
{
Dfs(e->To);
cur->Low = min(cur->Low, e->To->Low);
}
else
cur->Low = min(cur->Low, e->To->DfsN);
}
if (cur != Root && cur->Low == cur->DfsN)
{
while (Stack.back() != cur)
Stack.pop_back();
Stack.pop_back();
}
}
int main()
{
#ifdef _DEBUG
freopen("c:\\noi\\source\\input.txt", "r", stdin);
freopen("c:\\noi\\source\\output.txt", "w", stdout);
#endif
int totNode, totEdge, uId, vId;
while (scanf("%d%d", &totNode, &totEdge) && (totNode || totEdge))
{
Init(totNode);
for (int i = 1; i <= totEdge; i++)
{
scanf("%d%d", &uId, &vId);
_nodes[uId].Id = uId;
_nodes[vId].Id = vId;
AddEdge(uId + _nodes, vId + _nodes);
}
Dfs(Root);
if (Stack.size()!=totNode)
printf("No\n");
else
printf("Yes\n");
}
return 0;
}