poj_1236 强连通分支

#include<stdio.h>
#include<string.h>
#include<vector>
#include<stack>
#include<set>
#define MAX_NODE 105
#define min(a, b) a < b? a:b
#define max(a, b) a > b? a:b
using namespace std;
vector<vector<int> > gGraph;	//存储图结构
stack<int> gStack;				//用于Tarjan算法

bool gVisited[MAX_NODE];		//判断每个点是否被访问过
bool gInStack[MAX_NODE];		//在Tarjan算法中判断点是否在栈中
int gDfn[MAX_NODE];				//Tarjan算法中标记每个点最开始被访问的时间
int gLow[MAX_NODE];				//Tarjan算法中标记从每个点开始的DFS过程中经过的点所能访问到的点的最小的开始时间

int gClusterOfNode[MAX_NODE];	//标记每个点所在的强连通分支
int gIndex;						//标记点被第一次访问的时间
int gClusterIndex;				//强连通分支的编号

//Tarjan算法求强连通分支
void Tarjan(int u){
	gLow[u] = gDfn[u] = ++gIndex;
	gVisited[u] = true;
	gInStack[u] = true;
	gStack.push(u);

	for (int i = 0; i < gGraph[u].size(); i++){
		int v = gGraph[u][i];
		if (!gVisited[v]){
			Tarjan(v);
			gLow[u] = min(gLow[u], gLow[v]);
		}
		else if (gInStack[v]){
			gLow[u] = min(gLow[u], gDfn[v]);
		}
	}
	if (gLow[u] == gDfn[u]){
		int v;
		do{
			v = gStack.top();
			gStack.pop();
			gInStack[v] = false;
			gClusterOfNode[v] = gClusterIndex;	//标记所属的强连通分支编号
		} while (v != u);
		gClusterIndex++;
	}
}
vector<set<int> > gLinkTo;		//将强连通分支缩成点之后，每个点所连接到的另外的点，用set 使得插入时不重复。用于统计出度
vector<set<int> > gLinkFrom;	//将强连通分支缩成点之后，每个点被其他的点连接，用set 使得插入时不重复。用于统计入度

//将强连通分支缩成点，然后统计每个强连通分支的入度和出度
void ReconstructGraph(int nodes, int clusters){
	gLinkTo.resize(clusters);
	gLinkFrom.resize(clusters);

	for (int u = 1; u <= nodes; u++){
		for (int i = 0; i < gGraph[u].size(); i++){
			int v = gGraph[u][i];
			int uc = gClusterOfNode[u];
			int vc = gClusterOfNode[v];
			if (uc != vc){	//注意！！！自身不能连接到自身！
				gLinkTo[uc].insert(vc);
				gLinkFrom[vc].insert(uc);
			}
		}
	}
}

int main(){
	int n, u, v;
	scanf("%d", &n);
	gGraph.resize(n + 1);
	for (int i = 1; i <= n; i++){
		while (scanf("%d", &v) && v != 0){
			gGraph[i].push_back(v);
		}
	}
	gIndex = gClusterIndex = 0;
	memset(gVisited, false, sizeof(gVisited));
	memset(gInStack, false, sizeof(gInStack));
	for (int u = 1; u <= n; u++){
		if (!gVisited[u]){
			Tarjan(u);
		}
	}
	ReconstructGraph(n, gClusterIndex);
	int zero_outdegree = 0, zero_indegree = 0;
	for (int i = 0; i < gClusterIndex; i++){
		if (gLinkFrom[i].empty()){
			zero_indegree++;
		}
		if (gLinkTo[i].empty()){
			zero_outdegree++;
		}
	}
	if (gClusterIndex == 1){	//如果只有一个强连通分支，则直接输出1 0
		printf("1\n0\n");
	}else
		//最少需要发送的版本个数，等于将强连通分支缩成点之后的图中 入度为0的点的个数，
		//因为，一个有向无环图中，所有入度不为0的点，一定可以从某个入度为0的点出发可达。
		//而这些入度为0的点之间不能可达，故至少需要将所有这些入度为0的点设为起点

		//需要连接多少条有向边才能使得整个图称为一个强连通图
		//入度为0的点设为m个，不能被其他任何点可达，因此，要有 m条边连接到 这m个入度为0的点
		//出度为0的点设为n个，无法到达其他点，因此要有n条边从这n个出度为0的点连接出去
		//可以证明，在一个有向无环图中，所有点的入度和出度均不为0，则该图为强连通的
		//则，至少有 max(m, n)条边。

		printf("%d\n%d\n", zero_indegree, max(zero_outdegree, zero_indegree));
	return 0;
}