Network of Schools（POJ 1236）(强连通+缩点)（Tarjan算法）

http://acm.hust.edu.cn/vjudge/problem/17001

题意：N个学校之间有单向的网络，每个学校得到一套软件后，可以通过单向网络向周边的学校传输，问题1：初始至少需要向多少个学校发放软件，使得网络内所有的学校最终都能得到软件。2，至少需要添加 几条传输线路(边)，使任意向一个学校发放软件后，经过若干次传送，网络内所有的学校最终都能得到软件。 （来自北京大学暑期集训课件）

题解：
由于有向无环图中所有入度不为0的点，一定可以由某个入度为0的点出发可达。 所以对于问题1，只用求有向无环图中一共有多少入度为0的点，具体思路与方法可以看一下我的上一篇写缩点的那篇博客，里面有计算出度为0的详细解释，入度为0类似，答案就是num(in == 0)。对于问题2，就需要把所有点的出度和入度都变为非零值，由于问的是至少，所以先连min（num（in == 0）, num(out == 0)）条线，把出度和入度为0的连一起一部分，如果出入度为0点数目不同，则还需要把多的那些点和出（入）度不为0的点连起来，这样的话整个图就变成了一个连通图，问题2答案就是max（num（in == 0）, num(out == 0)）

思考：
1、有向无环图中所有入度不为0的点，一定可以由某个入度为0的点出发可达。

#include <iostream>
#include <cstdio>
#include <string>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <map>
#include <set>
#include <queue>

using namespace std;

const int maxn = 200;
const int maxm = 10010;

struct Edge{
    int to, next;
}edge[maxm];

int head[maxn], tot;
int low[maxn], dfn[maxn], Stack[maxn], Belong[maxn];
int Index, top;
int scc;
bool Instack[maxn];
int num[maxn];
int n;

void addedge(int u, int v)
{
    edge[tot].to = v;
    edge[tot].next = head[u];
    head[u] = tot++;
}

void Tarjan(int u)
{
    int v;
    low[u] = dfn[u] = ++Index;
    Stack[top++] = u;
    Instack[u] = true;
    for (int i = head[u]; i != -1; i = edge[i].next) {
        v = edge[i].to;
        if (!dfn[v]) {
            Tarjan (v);
            if (low[u] > low[v]) low[u] = low[v];
        } else if (Instack[v] && low[u] > dfn[v]) low[u] = dfn[v];
    }
    if (low[u] == dfn[u]) {
        scc++;
        do{
            v = Stack[--top];
            Instack[v] = false;
            Belong[v] = scc;
            num[scc]++;
        } while (v != u);
    }
}

void solve(int n)
{
    memset (dfn, 0, sizeof(dfn));
    memset (Instack, false ,sizeof(Instack));
    memset (num, 0, sizeof(num));
    Index  = scc = top = 0;
    for (int i = 1; i <= n; i++) {
        if (!dfn[i]) Tarjan (i);
    }
    if (scc == 1) {
        printf ("1\n0\n");
    } else {
        int in[maxn], out[maxn];
        memset (in, 0, sizeof(in));
        memset (out, 0, sizeof(out));

        //Á´Ê½Ç°ÏòÐÇ
        for (int u = 1; u <= n; u++) {
            for (int i = head[u]; i != -1; i = edge[i].next) {
                int v = edge[i].to;
                if (Belong[u] != Belong[v]) {
                    in[Belong[v]]++;
                    out[Belong[u]]++;
                    //printf ("x%d->x%d\n", u, v);
                }
            }
        }
        int ans1 = 0, ans2= 0;
        for (int i = 1; i <= scc; i++) {
            if (!in[i]) ans1++;
            if (!out[i]) ans2++;
        }
        printf ("%d\n%d\n", ans1, max (ans1, ans2));
    }
}

void init()
{
    tot = 0;
    memset (head, -1, sizeof(head));
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen ("in.txt", "r", stdin);
    #endif // ONLINE_JUDGE
    scanf ("%d", &n);
    init ();
    for (int i = 1; i <= n; i++) {
        int Re;
        while (scanf ("%d", &Re) != EOF) {
            if (Re == 0) break;
            addedge (i, Re);
        }
    }
    solve (n);
    return 0;
}