poj2186——双dfs求强连通分量

挑战程序设计原题。
感觉这个方法容易理解。
第一次dfs，先用一个时间戳来标记一下（时间戳就是寻访完这个店结束的时间），不是遍历的时间，
我们可以保证的是 最开始的点 时间戳最大。
然后 根据时间戳，从大到小进行第二次dfs（把图反向）
就可以得到k。
这个方法容易理解。
并且可以练习dfs。。。

#include <cstdio>#include <vector>#include <iostream>#include <cstring>/*双dfs求强联通分量。当然用targin算法也可以，不过我感觉这个也挺好的，具体思想是，先通过dfs来扫一下，再用vs一个数组打一个时间戳，就是扫完点结束的事件（先扎到底，那个就是1，在找枝叶）从时间戳最小的点开始扫，来计算强连通分量。并且强连通分量会得到一个 拓扑排序，这个拓扑排序很好用。先求强连通分量，然后对拓扑序最大的那个，判定一下他是否和所有点都联通。*/using namespace std;const int maxn=10009;vector <int> vs;vector <int>G[maxn];vector <int>G2[maxn];int tp[maxn];bool used[maxn];void add(int u,int v){  G[u].push_back(v);}void add1(int u,int v){  G2[u].push_back(v);}int dfs(int s){   used[s]=true;      for(int i=0;i<G[s].size();i++)      {   if(!used[G[s][i]])           dfs(G[s][i]);      }      vs.push_back(s);//    return 0;}int rdfs(int s,int k)//对反向图进行搜索{   tp[s]=k; used[s]=true;   for(int i=0;i<G2[s].size();i++)   {  if(!used[G2[s][i]])       rdfs(G2[s][i],k);   }   return 0;}int main(){   int m,n;    int a,b;    scanf("%d%d",&m,&n);    for(int i=1;i<=n;i++)        {scanf("%d%d",&a,&b);         add(a,b);         add1(b,a);//购置反向边。        }        //先bfs，对层数为2的在建一个边。     memset(used,false,sizeof(used));     for(int i=1;i<=m;i++)        if(!used[i])           dfs(i);      memset(used,false,sizeof(used));     // for(int i=0;i<vs.size();i++)          //cout<<vs[i]<<endl;      //for(int i=0;i<vs.size();i++)         // cout<<vs[i]<<endl;      int k=1;      for(int i=vs.size()-1;i>=0;i--)      {   if(!used[vs[i]])            {rdfs(vs[i],k);              //cout<<vs[i]<<endl;             k++;            }      }      k--;      int sum=0;      int v;     for(int i=1;i<=m;i++)     {  if(tp[i]==k)          {sum++;           v=i;          }     }     memset(used,false,sizeof(used));     rdfs(v,1);     for(int i=1;i<=m;i++)     {   if(!used[i])          {sum=0;          break;          }     }     cout<<sum<<endl;    return 0;}