CCF认证 2015-09 高速公路

要找可以两两连通的点个数，可以先找整个图的强连通分量，每个强连通分量中的所有点两两之间都可以互相到达

先用Tarjan算法找到所有强连通分量，假如一个强连通分量中有n个点，那么就有((n-1)+1)*n/2对互通点（等差求和）

找了一个模板，套用进来写下主函数就OK了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
// Tarjan算法   复杂度O(N+M)  
const int MAXN = 10010;//点数 
const int MAXM = 100010;//边数 
struct Edge
{
	int to,next;
}edge[MAXM];
int head[MAXN],tot; 
int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN];//Belong数组的值是1~scc 
int Index,top; 
int scc;//强连通分量的个数
bool Instack[MAXN]; 
int num[MAXN];//各个强连通分量包含点的个数，数组编号1~scc 
//num数组不一定需要，结合实际情况 
void addedge(int u,int v) 
{  
	edge[tot].to = v;edge[tot].next = head[u];head[u] = tot++; 
}
void Tarjan(int u) 
{
	int v; 
	Low[u] = DFN[u] = ++Index; 
	Stack[top++] = u; 
	Instack[u] = true; 
	for(int i = head[u];i != -1;i = edge[i].next) 
	{
		v = edge[i].to; 
		if( !DFN[v] )
		{
			Tarjan(v); 
			if( Low[u] > Low[v] )Low[u] = Low[v]; 
		}
		else if(Instack[v] && Low[u] > DFN[v]) 
		Low[u] = DFN[v]; 
	}
	if(Low[u] == DFN[u]) 
	{
		scc++; 
		do 
		{
			v = Stack[--top];
			Instack[v] = false; 
			Belong[v] = scc; 
			num[scc]++; 
		}while( v != u);
	}
}
void solve(int N) 
{
	memset(DFN,0,sizeof(DFN)); 
	memset(Instack,false,sizeof(Instack)); 
	memset(num,0,sizeof(num)); 
	Index = scc = top = 0;  
	for(int i = 1;i <= N;i++) 
	if(!DFN[i]) Tarjan(i); 
}
void init() 
{  
	tot = 0;  
	memset(head,-1,sizeof(head)); 
} 
int main()
{
	init();
	int n,m,a,b;
	scanf("%d%d",&n,&m);
	for(int i=0;i<m;i++)
	{
		scanf("%d%d",&a,&b);
		addedge(a,b);
	}
	solve(n); //求出强连通分量，个数为scc,每个分量点的个数保存在num[i]中 
	int ans=0;
	for(int i=1;i<=scc;i++)
	ans+=(num[i]-1)*num[i]/2; //等差数列求和 
	printf("%d\n",ans);
	return 0;
}