poj 2762 Going from u to v or from v to u?

题目描述：
为了让他们的儿子变得更勇敢些，Jiajia和Wind将他们带到一个大洞穴中。洞穴中有n个房
间，有一些单向的通道连接某些房间。每次，Wind选择两个房间x和y，要求他们的一个儿子从
一个房间走到另一个房间，这个儿子可以从x走到y，也可以从y走到x。Wind保证她布置的任
务是可以完成的，但她确实不知道如何判断一个任务是否可以完成。为了使Wind下达任务更容
易些，Jiajia决定找这样的一个洞穴，每对房间（设为x和y）都是相通（可以从x走到y，或者
可以从y走到x）的。给定一个洞穴，你能告诉Jiajia，Wind是否可以任意选择两个房间而不用
担心这两个房间可能不相通吗？
//
求解的是单连通性，但首先要转换成强连通分量的求解。这是因为，强连通分量中
的顶点间存在双向的路径，因此可以将每个强连通分量收缩成一个新的顶点。在有向图的处理中
经常需要将强连通分量收缩成一个顶点。

强连通分量收缩后，再求其拓扑排序。假设求得的拓扑序存储在topo[MAX]中，topo[i]与
topo[i+1]存在边连通（i 到i+1 或i+1 到i），则定有i 到i+1 的边。而如果每个topo[i]与topo[i+1]
都存在边连通（即有i 到i+1 的边）时，topo[i]到任意topo[j]便都有边连通。
// topsort时偷懒 用了邻接矩阵 复杂度成了 n^2  写成邻接表应该要快好多
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#include <math.h>
#include <stdio.h>
#include <string.h>
using namespace std;
#define MOD 1000000007
#define maxn 6100
#define maxm 1010
struct Edge{
int to;
int next;
Edge(){};
Edge(int u,int v){to=u;next=v;}
}E[maxn];
stack<int> S;
int V[maxm],num;
int belong[maxm];
int pre[maxm];
int dfst,scc;
int ans;
bool G[maxm][maxm];
int in[maxm];
void init(int n){
dfst=scc=0;
num=0;
ans=0;
while(!S.empty())
S.pop();
for(int i=1;i<=n;i++){
V[i]=-1;
pre[i]=0;
belong[i]=0;
}
}
void add(int u,int v){
E[num].to=v;
E[num].next=V[u];
V[u]=num++;
}
int tarjan(int u){
int lowu=pre[u]=++dfst;
int v,e;
S.push(u);
for(e=V[u];e!=-1;e=E[e].next){
v=E[e].to;
if(!pre[v]){
int lowv=tarjan(v);
lowu=min(lowu,lowv);
}
else if(!belong[v]) lowu=min(lowu,pre[v]);
}
if(lowu==pre[u]){
scc++;
for(;;){
int x=S.top();S.pop();
belong[x]=scc;
if(x==u) break;
}
}
return lowu;
}
int top[maxm],tn;
void topsort(){
int i,j,k;
tn=0;
bool vi[maxm]={0};
while(1){
for(i=1;i<=scc;i++)
if(!in[i]&&!vi[i]) break;
vi[i]=true;
//  printf("%d ",i);
if(i>scc) break;
top[tn++]=i;
for(j=1;j<=scc;j++)
if(G[i][j]) in[j]--;
}
}
int main()
{
int n,m,T;
int u,v;
int i,j=1;
scanf("%d",&T);
while(T--){
scanf("%d %d",&n,&m);
init(n);
for(i=1;i<=m;i++){
scanf("%d %d",&u,&v);
add(u,v);
}
for(i=1;i<=n;i++)
if(!pre[i]) tarjan(i);
// for(i=1;i<=n;i++) printf("%d ",belong[i]);
for(i=1;i<=scc;in[i]=0,i++)
for(j=1;j<=scc;j++)
G[i][j]=0;
int e,u,v;

for(i=1;i<=n;i++)
{
for(e=V[i];e!=-1;e=E[e].next){
u=belong[i];
v=belong[E[e].to];
if(u!=v){
G[u][v]=1;
in[v]++;
//  printf("%d ",v);
}
}
}

topsort();//printf("\\");
int flag=1;
for(i=1;i<tn;i++)
if(!G[top[i-1]][top[i]]){
flag=0;break;
}
if(flag) printf("Yes\n"); else printf("No\n");
}
return 0;
}