[kosaraju] Csp201509-4高速公路

题意

某国有n个城市，为了使得城市间的交通更便利，该国国王打算在城市之间修一些高速公路，由于经费限制，国王打算第一阶段先在部分城市之间修一些单向的高速公路。
　　现在，大臣们帮国王拟了一个修高速公路的计划。看了计划后，国王发现，有些城市之间可以通过高速公路直接（不经过其他城市）或间接（经过一个或多个其他城市）到达，而有的却不能。如果城市A可以通过高速公路到达城市B，而且城市B也可以通过高速公路到达城市A，则这两个城市被称为便利城市对。
　　国王想知道，在大臣们给他的计划中，有多少个便利城市对。

样例

样例输入：

输入的第一行包含两个整数n, m，分别表示城市和单向高速公路的数量。
接下来m行，每行两个整数a, b，表示城市a有一条单向的高速公路连向城市b。

5 5
1 2
2 3
3 4
4 2
3 5

样例输出:

输出一行，包含一个整数，表示便利城市对的数量。

3

---

思路

1.一道裸的求强连通分量的题
强连通分量:原图中的一个子图,该子图内任意两点互相可达(强连通)
2.由于题目给定,求互相可达的城市对,而互相可达的城市对只可能出现在同一个强连通分量中
故我们找出每个成熟所属的强连通分量,然后计算每个SCC中有多少对城市组合即可

---

总结

强连通分量的模板题练习
现在原图中 dfs求出后序序列
再在反图中 按照逆后序序列的顺序进行染色

---

代码

#include<iostream>
#include<string.h>
#include<string>
#include<algorithm>
#include<math.h>
#include<stdio.h>
#include<queue>
#define MAXN 10010
#define INF 0x7ffffff
#define ll long long
using namespace std;
int n,m,c[MAXN],scc_num[MAXN],dfn[MAXN],vis[MAXN],dcnt,scnt;
vector<int> G1[MAXN],G2[MAXN];
void dfs1(int x)
{
    vis[x]=1;
    for(auto y:G1[x])   if(!vis[y]) dfs1(y);
    dfn[++dcnt]=x;
}

void dfs2(int x)
{
    c[x]=scnt;scc_num[scnt]++;
    for(auto y:G2[x])   if(!c[y])   dfs2(y);
}
void kosaraju()
{
    for(int i=1;i<=n;i++)   if(!vis[i]) dfs1(i);//原图中后序遍历 
    for(int i=n;i>=1;i--)   if(!c[dfn[i]])  ++scnt,dfs2(dfn[i]);//反图中通过逆后序序列 染色
    ll ans=0;
    for(int i=1;i<=scnt;i++)    ans+=(ll)scc_num[i]*(scc_num[i]-1)/2;
    cout<<ans<<endl;
}
int main()
{
    cin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        int u,v;
        cin>>u>>v;
        G1[u].push_back(v);//原图
        G2[v].push_back(u);//反图
    }
    kosaraju();
    system("pause");
    return 0;
}