强联通 poj 2762

t个样例    (注意清零)

n个点m条边 有向;

任意2点是否能从a->b或者b->a;

Yes  No

 #include<stdio.h>
 #include<algorithm>
 #include<string.h>
 #include<queue>
 #include<math.h>
 #include<stack>

 using namespace std;

 #define MAXN  20000
 #define MAXN1 10000
 int cnt,k,num;
 struct edg
 {
     int next,to,fr;

 }x[MAXN+];
 int head[MAXN1],dfn[MAXN1],low[MAXN1],f[MAXN1];
 bool vis[MAXN1];
 int in[MAXN1];

 void add(int u,int v)
 {
     x[cnt].fr=u;
     x[cnt].next=head[u];
     x[cnt].to=v;
     head[u]=cnt++;
 }
 stack<int>s;

 void dfs(int u)
 {
     low[u]=dfn[u]=k++;
     s.push(u);
     vis[u]=;
     int i;
     for(i=head[u];i!=-;i=x[i].next)
     {
         if(!dfn[x[i].to])
         {
             dfs(x[i].to);
             low[u]=min(low[u],low[x[i].to]);
         }
         else if(vis[x[i].to])
             low[u]=min(low[u],dfn[x[i].to]);
     }
     if(low[u]==dfn[u]) //强联通分量
     {
         num++;
         while(!s.empty())
         {
             int now=s.top();
             s.pop();
             f[now]=num;  //新的点
             vis[now]=;
             if(u==now)break;
         }
     }
 }
 queue<int>q1;

 int judge()
 {
     int i;
     while(!q1.empty())q1.pop();   //很重要
     for(i=;i<=num;i++)
     {
         if(!in[i])
             q1.push(i);
     }
     if(q1.size()>)return ;   //入度为0的点1个

     while(!q1.empty())
     {
         int now=q1.front();
         q1.pop();
         for(i=head[now];i!=-;i=x[i].next)
         {
             in[x[i].to]--;
             if(!in[x[i].to])
                 q1.push(x[i].to);
         }
         if(q1.size()>)return ;
     }
     return ;
 }

 int main()
 {
     int t;
     scanf("%d",&t);

     while(t--)
     {
         int n,m;
         scanf("%d%d",&n,&m);
         int i;
         memset(head,-,sizeof(head));
         cnt=;
         for(i=;i<=m;i++)
         {
             int a,b;
             scanf("%d%d",&a,&b);
             add(a,b);
         }
         k=;
         num=;
         while(!s.empty())s.pop();
         memset(dfn,,sizeof(dfn));
         memset(low,,sizeof(low));
         memset(vis,,sizeof(vis));
         memset(f,,sizeof(f));
         for(i=;i<=n;i++)
         {
             if(!dfn[i])
                 dfs(i);
         }
         int en=cnt;
         memset(head,-,sizeof(head));
         memset(in,,sizeof(in));
         cnt=;
         for(i=;i<en;i++)  //重新建图
         {
             int a,b;
             a=f[x[i].fr];
             b=f[x[i].to];
             if(a!=b)
             {
                 add(a,b);
                 in[b]++;
             }
         }
         if(judge())
             printf("Yes\n");
         else
             printf("No\n");
     }

     return ;
 }