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
"); else printf("No
");
    }
    return 0;
}