codeforce_exercise_r19

高速公路（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

问题分析

解题思路

这个题，打眼一看，emm，传递闭包，floyd5行完事。真简单。然后，看了一眼数据范围，1e4…然后，才反应过来，不对，是强连通分量。然后再仔细一看，实际上是一个裸的强连通分量题。这里的便利城市对就是一个强连通分量中不同的两个点。这样，需要求出强连通分量的节点个数n，然后计算组合数，最后累加即为最后的结果。

参考代码

#include <iostream>
#include <vector>
#include <cstring>

using namespace std;

class edge
{
public:
	int from;
	int to;
	edge(int ff,int tt)
	{
		from=ff;
		to=tt;
	}
};

vector<edge> edges;
vector<int> g1[10010];
vector<int> g2[10010];
vector<int> tp;
int vis[10010];
int m,n;
int ed;
int cnt;

void init()
{
	edges.clear();
	for(int i=0;i<10010;i++)
	{
		g1[i].clear();
		g2[i].clear();
	}
	memset(vis,0,sizeof(vis));
}

void addedge(int u,int v)
{
	edge e1(u,v);
	edges.push_back(e1);
	g1[u].push_back(edges.size()-1);
	edge e2(v,u);
	edges.push_back(e2);
	g2[v].push_back(edges.size()-1);
}

void dfs1(int start_point)
{
    vis[start_point]=1;
	for(int i=0;i<g1[start_point].size();i++)
	{
		edge te=edges[g1[start_point][i]];
		if(vis[te.to]==0) dfs1(te.to);
	}
	tp.push_back(start_point);	
}

int dfs2(int start_point)
{
	cnt++;
	vis[start_point]=1;
	for(int i=0;i<g2[start_point].size();i++)
	{
		edge te=edges[g2[start_point][i]];
		if(vis[te.to]==0) dfs2(te.to);
	}
	return cnt;
}

int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int u,v;
		scanf("%d %d",&u,&v);
		addedge(u,v);
	}
	ed=1;
	for(int i=1;i<=n;i++)
	{
		if(vis[i]==0) dfs1(i);
	}
    int ans=0;
    memset(vis,0,sizeof(vis));
    for(int i=tp.size()-1;i>=0;i--)
    {
    	if(vis[tp[i]]==0)
    	{
    		cnt=0;
    		int num=dfs2(tp[i]);
    		ans+=num*(num-1)/2;
		}
	}
	printf("%d",ans);
	return 0;
}

心得体会

幸亏上学期学了算法，了解了一些强连通分量的性质，要不然这个题是真的做不出来。另外好像上上个学期的程设里面也讲过这个算法，但是由于确实有点复杂，再加上讲的有点快，所以这个算法当时的死记记住的。现在想想，其实感觉也没那么难理解。所以，平时多积累，多掌握一些算法还是挺重要的。