ZOJ3232--It's not Floyd Algorithm(强连通+缩点+建图+floyd)

题目是给你一个矩阵，1表示u可以到达v，0代表不可到达，问你至少需要多少条边组成的传递闭包符合这个矩阵。

我们可以求出强连通分量，然后在对每个强连通分量进行缩点，每个强连通分量的最少边的数量就是该强连通分量的结点数，再建立新图。对新图中的点用floyd算法，若图中用floyd算法能达到的，且在新图中为1的点，我们将它变为0，则答案就是每个强连通分量内的边数加上新图中为0的点的个数。

#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <string.h>
#include <math.h>
#include <queue>
#include <stack>
#include <vector>
#include <ctype.h>
#define INF 999999999
#define LL long long
#define M 205
using namespace std;

vector<int> G[M];  
int dfn[M],low[M],sccno[M],scc_cnt;
int indx;  
int num[M];
int d[M][M];
int TG[M][M];
stack<int> s;  

void Tarjan(int u)  
{  
	indx++;
	dfn[u]=low[u]=indx; 
	s.push(u);  
	for(int i=0;i<G[u].size();i++)  
	{  
		int v=G[u][i];  
		if(!dfn[v])  
		{  
			Tarjan(v);  
			low[u]=min(low[u],low[v]);  
		}  
		else if(!sccno[v])  
		{  
			low[u]=min(low[u],dfn[v]);  
		}  
	}  
	if(low[u]==dfn[u])    //缩点
	{  
		scc_cnt++;  
		for(;;)  
		{  
			int x=s.top();  
			s.pop();  
			sccno[x]=scc_cnt;  
			num[scc_cnt]++;  //记录每个强连通分量的结点数
			if(x==u)  
				break;  
		}  
	}  
}  

void find_scc(int n)  
{  
	indx=scc_cnt=0;  
	memset(sccno,0,sizeof(sccno));  
	memset(dfn,0,sizeof(dfn));  
	for(int i=0;i<n;i++)  
		if(!dfn[i])  
			Tarjan(i);  
} 

int main()
{
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		for(int i=0;i<n;i++)
			G[i].clear();
		memset(num,0,sizeof(num));
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<n;j++)
			{
				scanf("%d",&d[i][j]);
				if(i!=j && d[i][j]==1)
					G[i].push_back(j);
			}
		}
		find_scc(n);
		int ans=0;
		memset(TG,0,sizeof(TG));
		for(int u=0;u<n;u++)       //建立新图
		{
			for(int i=0;i<G[u].size();i++)
			{
				int v=G[u][i];
				if(sccno[u]!=sccno[v] && d[u][v]==1)
					TG[sccno[u]][sccno[v]]=1;
			}
		}
		for(int k=1;k<=scc_cnt;k++)   //floyd算法
		{
			for(int i=1;i<=scc_cnt;i++)
			{
				for(int j=1;j<=scc_cnt;j++)
				{
					if(TG[i][j]==1 && TG[i][k]==1 && TG[k][j]==1)
						TG[i][j]=0;
				}
			}
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
			{
				if(TG[i][j]==1)
					ans++;
			}
		}
		for(int i=1;i<=scc_cnt;i++)
		{
			if(num[i]>1)
				ans+=num[i];
		}
		printf("%d\n",ans);
	}
	return 0;
}