nyoj 120 校园网络 【强连通】

校园网络

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 5

描述

南阳理工学院共有M个系，分别编号1~M,其中各个系之间达成有一定的协议，如果某系有新软件可用时，该系将允许一些其它的系复制并使用该软件。但该允许关系是单向的，即：A系允许B系使用A的软件时，B未必一定允许A使用B的软件。

现在，请你写一个程序，根据各个系之间达成的协议情况，计算出最少需要添加多少个两系之间的这种允许关系，才能使任何一个系有软件使用的时候，其它所有系也都有软件可用。

输入

第一行输入一个整数T，表示测试数据的组数(T<10)
每组测试数据的第一行是一个整数M，表示共有M个系(2<=M<=100)。
随后的M行，每行都有一些整数，其中的第i行表示系i允许这几个系复制并使用系i的软件。每行结尾都是一个0，表示本行输入结束。如果某个系不允许其它任何系使用该系软件，则本行只有一个0.

输出

对于每组测试数据，输出最少需要添加的这种允许关系的个数。

样例输入

1
5
2 4 3 0
4 5 0
0
0
1 0

样例输出

2

分析：

强连通问题，代码没有看懂，先保存下来，回去让大神教我一下。

感觉解释的一个比较好的博客地址：

https://www.byvoid.com/blog/scc-tarjan

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<stack> 
using namespace std;

const int maxn =105;
int map[maxn][maxn];
int low[maxn];
int dfn[maxn];
int in[maxn];
int out[maxn];
int instack[maxn];
int t[maxn];
int n,order,res,ans;
stack<int>s;

void init()
{
	memset(map,0,sizeof(map));
	memset(low,0,sizeof(low));
	memset(dfn,0,sizeof(dfn));
	memset(in,0,sizeof(in));
	memset(out,0,sizeof(out));
	memset(instack,0,sizeof(instack));
	memset(t,0,sizeof(t));
	while(!s.empty())
	s.pop();
	res=0;
	order=0;
}

int min(int x,int y)
{
	return x<y?x:y; 
}

void tr(int u)
{
	int v;
	dfn[u]=low[u]=++order;
	instack[u]=1;
	s.push(u);
	for(int i=1;i<=n;i++)
	{
		if(map[u][i])
		{
			if(!dfn[i])
			{
				tr[i];
				low[u]=min(low[u],low[i]);
			}
			else if(instack[i])
			low[u]=min(low[u],dfn[i]);
		}
	}
	if(dfn[u]==low[u])
	{
		++res;
		do{
			v=s.top();
			s.pop();
			instack[v]=0;
			t[v]=res;
		}while(v!=u);
	}
}

void tarjan()
{
	for(int i=1;i<=n;i++)
	{
		if(!dfn[i])
		tr(i);
	}
}

void solve()
{
	for(int i=1;i<=n;i++)
	  {
	  	for(int j=1;j<=n;j++)
	  	   if(map[i][j])
	  	   {
	  	   	++in[t[i]];
	  	   	++out[t[i]];
	  	   }
	  } 
	  int xx,yy;
	  xx=yy=0;
	  for(int i=1;i<=res;i++)
	  {
	  	if(in[i]==0)
	  	xx++;
		  else if(out[i]==0)
		  yy++; 
	  }
	  ans=xx>yy?xx:yy;
}

int main()
{
	int T,x;
	scanf("%d",&T);
	while(T--)
	{
		init();
		scanf("%d",&n);
		for(int i=1;i<=n;i++)
		{
			while(scanf("%d",&x),x)
             map[i][x]=1;
		}
		tarjan();
		solve();
		if(res==1)
		printf("0\n");
		else
		printf("%d\n",ans);
	}
	return 0;
}