题目描述
一些学校连入一个电脑网络。那些学校已订立了协议:每个学校都会给其它的一些学校分发软件(称作“接受学校”)。注意即使 �B 在 �A 学校的分发列表中,�A 也不一定在 �B 学校的列表中。
你要写一个程序计算,根据协议,为了让网络中所有的学校都用上新软件,必须接受新软件副本的最少学校数目(子任务 A)。更进一步,我们想要确定通过给任意一个学校发送新软件,这个软件就会分发到网络中的所有学校。为了完成这个任务,我们可能必须扩展接收学校列表,使其加入新成员。计算最少需要增加几个扩展,使得不论我们给哪个学校发送新软件,它都会到达其余所有的学校(子任务 B)。一个扩展就是在一个学校的接收学校列表中引入一个新成员。
输入格式
输入文件的第一行包括一个正整数 �N,表示网络中的学校数目。学校用前 �N 个正整数标识。
接下来 �N 行中每行都表示一个接收学校列表(分发列表),第 �+1i+1 行包括学校 �i 的接收学校的标识符。每个列表用 00 结束,空列表只用一个 00 表示。
输出格式
你的程序应该在输出文件中输出两行。
第一行应该包括一个正整数,表示子任务 A 的解。
第二行应该包括一个非负整数,表示子任务 B 的解。
输入输出样例
输入 #1
5 2 4 3 0 4 5 0 0 0 1 0输出 #1
1 2说明/提示
2≤N≤100。
思路:
对与有向无环图DAG来说,A任务的解显然是入度为0的点的个数,即起点。可以用tarjan缩点。
B任务即最少加多少条边可以将原图变为一个强连通分量。有一个结论,最小加边数为,出度为0的点和入度为0的点较大的哪个,即起点和终点 。证明较为复杂。
Code:
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<unordered_map>
using namespace std;
#define x first
#define y second
#define PII pair<int,int>
#define V vector<int>
#define endl "\n"
#define init(array) memset(array,0,sizeof array)
typedef long long ll;
typedef unsigned long long llu;
const int INF=0x3f3f3f3f;
const int N=1e4+10,M=5e6+10;
int n;
int h[N],e[M],ne[M],idx;
int dfn[N],low[N],timestamp;
int din[N],dout[N];
int id[N],scc;
bool in_stk[N];
stack<int> stk;
void add(int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++timestamp;
stk.push(u);
in_stk[u]=true;
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(!dfn[j])
{
tarjan(j);
low[u]=min(low[u],low[j]);
}
else if(in_stk[j])
low[u]=min(low[u],dfn[j]);
}
if(dfn[u]==low[u])
{
scc++;
int y;
do {
y=stk.top();
stk.pop();
in_stk[y]=false;
id[y]=scc;
} while(y!=u);
}
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
cout.tie(nullptr);
memset(h,-1,sizeof h);
cin>>n;
for(int i=1;i<=n;i++)
{
int x;
while(cin>>x,x) add(i,x);
}
for(int i=1;i<=n;i++)
if(!dfn[i]) tarjan(i);
for(int i=1;i<=n;i++)
{
for(int j=h[i];~j;j=ne[j])
{
int k=e[j];
int a=id[i],b=id[k];
if(a!=b)
{
din[b]++;
dout[a]++;
}
}
}
int in=0,out=0;
for(int i=1;i<=scc;i++)
{
if(!din[i]) in++;
if(!dout[i]) out++;
}
cout<<in<<endl;
if(scc==1) cout<<0; //特判已经是强连通分量
else cout<<max(in,out);
return 0;
}