POJ 1236 Network of Schools 强连通分支 .

题目地址：http://poj.org/problem?id=1236

题目意思：

给定一个有向图，求：
1) 至少要选几个顶点，才能做到从这些顶点出
发，可以到达全部顶点
2) 至少要加多少条边，才能使得从任何一个顶
点出发，都能到达全部顶点

解：

1）缩点

2） DAG上面有多少个入度为0的顶点就是1）的答案

3）假定有 n 个入度为0的点，m个出度为0的点，max(m,n)就是第二个问题的解

Korasaju

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100+10;
bool vis[maxn];
int ID[maxn];    //点的颜色编号 
int ncolor;      
vector<int> s;
vector<vector<int> > color(maxn); //同一种颜色的点 
vector<vector<int> > G(maxn);  
vector<vector<int> > GT(maxn);
void dfs(int u){
	vis[u]=true;
	for(int i=0;i<G[u].size();i++){
		int v=G[u][i];
		if(!vis[v]) dfs(v);
	}
	s.push_back(u);
}
void dfs2(int u){
	vis[u]=true;
	ID[u]=ncolor;
	color[ncolor].push_back(u);
	for(int i=0;i<GT[u].size();i++){
		int v=GT[u][i];
		if(!vis[v]) dfs2(v);
	}
}
void Korasaju(int n)
{
	memset(vis,false,sizeof(vis));
	s.clear();
	for(int i=1;i<=n;i++) if(!vis[i]) dfs(i);

	memset(vis,false,sizeof(vis));
	ncolor=0;
	for(int i=s.size()-1;i>=0;i--) {
		int j=s[i];
		if(!vis[j]) ncolor++,dfs2(j);
	}
}
void solve()
{
	int m=0,ans; //m标记入度为0的点有几个 
	for(int i=1;i<=ncolor;i++)
	{
		bool ok=true;
		for(int j=0;j<color[i].size()&&ok;j++)
		{
			int u=color[i][j];
			for(int k=0;k<GT[u].size()&&ok;k++)
			{
				int v=GT[u][k];
				if(ID[u]!=ID[v])  //有出度 
					ok=false; 
			}
		}
		if(ok) m++;
	}
	
	int n=0; //n标记出度为0的点有几个 
	for(int i=1;i<=ncolor;i++)
	{
		bool ok=true;
		for(int j=0;j<color[i].size()&&ok;j++)
		{
			int u=color[i][j];
			for(int k=0;k<G[u].size()&&ok;k++)
			{
				int v=G[u][k];
				if(ID[u]!=ID[v])  //有出度 
					ok=false; 
			}
		}
		if(ok) n++;
	}
	if(ncolor==1) cout<<1<<endl<<0;   //注意，只有一个强连通分支的时候，不需要边 
	else cout<<m<<endl<<max(m,n);
}

int main()
{
	int n;
	cin>>n;
	for(int u=1;u<=n;u++)
	{
		int v;
		while(cin>>v&&v)
		{
			G[u].push_back(v);
			GT[v].push_back(u);
		}
	}
	Korasaju(n);
	solve();
	return 0;
} 

Tarjan

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100+10;
bool vis[maxn];
int ID[maxn];    //µãµÄÑÕÉ«±àºÅ 
int ncolor,index;      
vector<int> st,dfn(maxn),low(maxn);
vector<vector<int> > color(maxn); //ͬһÖÖÑÕÉ«µÄµã 
vector<vector<int> > G(maxn);  
vector<vector<int> > GT(maxn);
void Tarjan(int u)
{
	dfn[u]=low[u]=++index;
	vis[u]=true;
	st.push_back(u);
	for(int i=0;i<G[u].size();i++)
	{
		int v=G[u][i];
		if(!vis[v]){
			Tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else if(find(st.begin(),st.end(),v)!=st.end()) //in stack
			low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u])
	{
		int v; ncolor++;
		do{
			v=st.back(); st.pop_back();
			color[ncolor].push_back(v);
			ID[v]=ncolor;
		}while(v!=u);
	}
}
int getGraph(int n)
{
	index=ncolor=0;
	for(int i=1;i<=n;i++)
		if(!vis[i]) Tarjan(i);
}
void solve()
{
	int m=0,ans; //m±ê¼ÇÈë¶ÈΪ0µÄµãÓм¸¸ö 
	for(int i=1;i<=ncolor;i++)
	{
		bool ok=true;
		for(int j=0;j<color[i].size()&&ok;j++)
		{
			int u=color[i][j];
			for(int k=0;k<GT[u].size()&&ok;k++)
			{
				int v=GT[u][k];
				if(ID[u]!=ID[v])  //Óгö¶È 
					ok=false; 
			}
		}
		if(ok) m++;
	}
	
	int n=0; //n±ê¼Ç³ö¶ÈΪ0µÄµãÓм¸¸ö 
	for(int i=1;i<=ncolor;i++)
	{
		bool ok=true;
		for(int j=0;j<color[i].size()&&ok;j++)
		{
			int u=color[i][j];
			for(int k=0;k<G[u].size()&&ok;k++)
			{
				int v=G[u][k];
				if(ID[u]!=ID[v])  //Óгö¶È 
					ok=false; 
			}
		}
		if(ok) n++;
	}
	if(ncolor==1) cout<<1<<endl<<0;   //×¢Ò⣬ֻÓÐÒ»¸öÇ¿Á¬Í¨·ÖÖ§µÄʱºò£¬²»ÐèÒª±ß 
	else cout<<m<<endl<<max(m,n);
}

int main()
{
	int n;
	cin>>n;
	for(int u=1;u<=n;u++)
	{
		int v;
		while(cin>>v&&v)
		{
			G[u].push_back(v);
			GT[v].push_back(u);
		}
	}
	getGraph(n);
	solve();
	return 0;
}