POJ1236 tarjan

Network of Schools

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 19613		Accepted: 7725

Description

A number of schools are connected to a computer network. Agreements have been developed among those schools: each school maintains a list of schools to which it distributes software (the “receiving schools”). Note that if B is in the distribution list of school A, then A does not necessarily appear in the list of school B 
You are to write a program that computes the minimal number of schools that must receive a copy of the new software in order for the software to reach all schools in the network according to the agreement (Subtask A). As a further task, we want to ensure that by sending the copy of new software to an arbitrary school, this software will reach all schools in the network. To achieve this goal we may have to extend the lists of receivers by new members. Compute the minimal number of extensions that have to be made so that whatever school we send the new software to, it will reach all other schools (Subtask B). One extension means introducing one new member into the list of receivers of one school. 

Input

The first line contains an integer N: the number of schools in the network (2 <= N <= 100). The schools are identified by the first N positive integers. Each of the next N lines describes a list of receivers. The line i+1 contains the identifiers of the receivers of school i. Each list ends with a 0. An empty list contains a 0 alone in the line.

Output

Your program should write two lines to the standard output. The first line should contain one positive integer: the solution of subtask A. The second line should contain the solution of subtask B.

Sample Input

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

Sample Output

1
2

Source

IOI 1996

思路：求一个有向图从几个点出发可以遍历整个图、以及至少加几条边使整张图强联通。

缩点以后，显然入度为0的点的个数就是第一问的答案。 
然后第二问答案显然是入度为0和出度为0的个数的最大值，即出入度为0的点间连条边就可以了。

代码：

#include"bits/stdc++.h"

#define db double
#define ll long long
#define vl vector<ll>
#define ci(x) scanf("%d",&x)
#define cd(x) scanf("%lf",&x)
#define cl(x) scanf("%lld",&x)
#define pi(x) printf("%d\n",x)
#define pd(x) printf("%f\n",x)
#define pl(x) printf("%lld\n",x)
#define rep(i, a, n) for (int i=a;i<n;i++)
#define per(i, a, n) for (int i=n-1;i>=a;i--)
#define fi first
#define se second
using namespace std;
typedef pair<int, int> pii;
const int N = 1e4 + 5;
const int mod = 1e9 + 7;
const int MOD = 998244353;
const db PI = acos(-1.0);
const db eps = 1e-10;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3fffffffffffffff;
int low[N], dfn[N],head[N],beg[N];
bool ins[N];
int in[N],out[N];
int n;
int cnt, id, num;
stack<int> s;

struct P {
    int to, nxt;
} e[2 * N];

void add(int u, int v) {
    e[cnt].to = v;
    e[cnt].nxt = head[u];
    head[u] = cnt++;
}
void tarjan(int u)
{
    low[u]=dfn[u]=++id;    ins[u]=1;
    s.push(u);//入栈
    for(int i=head[u];~i;i=e[i].nxt){
        int v=e[i].to;
        if(!dfn[v]) tarjan(v),low[u]=min(low[u],low[v]);
        else if(ins[v]) low[u]=min(low[u],dfn[v]);
    }
    if(low[u]==dfn[u]){
        int v;
        do{
            v=s.top();s.pop();
            beg[v]=num;//缩点
            ins[v]=0;
        }while(u!=v);
        num++;
    }
}

void init() {
    memset(ins, 0, sizeof(ins));
    memset(head, -1, sizeof(head));
    memset(dfn, 0, sizeof(dfn));
    memset(low, 0, sizeof(low));
    memset(beg, 0, sizeof(beg));
    memset(in, 0, sizeof(in));
    memset(out, 0, sizeof(out));
    while(!s.empty()) s.pop();
    cnt = id = num = 0;
}
int main() {
    while (scanf("%d",&n) != EOF) {
        init();
        for(int i=1;i<=n;i++){
            int x;
            while(scanf("%d",&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 = head[i]; j != -1; j = e[j].nxt) {
                int v=e[j].to;
                if(beg[i]!=beg[v]) out[beg[i]]++,in[beg[v]]++;
            }
        }
        int I=0,O=0;
        for(int i=0;i<num;i++){
            if(!in[i]) I++;
            if(!out[i]) O++;
        }
        if(num==1) printf("1\n0");
        else pi(I),pi(max(I,O));
    }
    return 0;
}