POJ 1236 Network of Schools（强连通分量）

题目大意：

给定一个有向图，求：

1、至少要选多少个顶点，才能做到从这些顶点出发，能够到达全部的顶点

2、至少要加多少条边，才能使得从任何一个顶点出发，能都到达全部顶点

补充：

有向无环图中所有入度不为0的点，一定可以由某个入度为0的点出发可达。（由于无环，所以从任何入度不为0的点往回走，必然终止于一个入度为0的点）

解题思路：

1、求出所有的强连通分量

2、每个强连通分量缩成一点，则形成一个有向无环图DAG

3、DAG上面有多少个入度为0的顶点，问题1的答案就是多少

4、在DAG上要加几条边，才能使得DAG变成前连通的，问题2的答案就是多少

加边的方法：

要为每个入度为0的点添加入边，为每个出度为0的点添加出边

假定有n个入度为0的点，m个出度为0的点，max(m, n)就是问题2的答案

代码：

#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#define M 105
using namespace std;

struct School{
    int starT;
    int endT;
    int mark;
}sch[M];

struct Order{   //记录在按照endT从大到小排序后的元素所对应的排序前的序号
    int th;
    int endT;
}order[M];

int n, t, lable;
int adj[M][M], adjT[M][M];
int vis[M];

int cmp(Order o1, Order o2) {
    return o1.endT > o2.endT;
}

void dfs1(int v) {
    vis[v] = 1;
    for(int i=0; i<n; i++) {
        if(adj[v][i] == 1 && vis[i] == 0) {
            sch[i].starT = t++;
            dfs1(i);
            order[i].endT = t;
            sch[i].endT = t++;
        }
    }
    return ;
}

int max(int a, int b) {
    return a > b ? a : b;
}

void dfs2(int v) {
    vis[v] = 1;
    for(int j=0; j<n; j++) {
        if(adjT[v][j] == 1 && vis[j] == 0) {
            sch[j].mark = lable;
            dfs2(j);
        }
    }
    return ;
}

int main() {
    int v;
    int i, j, k;
    while(~scanf("%d", &n)) {
        memset(adj, 0, sizeof(adj));
        memset(adjT, 0, sizeof(adjT));
        for(i=0; i<n; i++) {
            order[i].th = i;
            while(scanf("%d", &v), v) {
                adj[i][v-1] = 1;
                adjT[v-1][i] = 1;
            }
        }
        t = 1;
        memset(vis, 0, sizeof(vis));
        for(i=0; i<n; i++) {
            if(vis[i] == 0) {
                sch[i].starT = t++;
                dfs1(i);
                order[i].endT = t;
                sch[i].endT = t++;
            }
        }
        sort(order, order+n, cmp);  //排序后点的顺序出现絮乱，一定要和之前的点的顺序对应正确！
        lable = 1;
        memset(vis, 0, sizeof(vis));
        for(i=0; i<n; i++) {
            if(vis[order[i].th] == 0) {     //注意vis数组下标的值
                sch[order[i].th].mark = lable;
                dfs2(order[i].th);
                lable ++;
            }
        }
        int inN, inA[M];    //缩点后的DAG上的入度为0的点
        int outN, outA[M];  //缩点后的DAG上的出度为0的点
        inN = 0; outN = 0;
        memset(inA, 0, sizeof(inA));
        memset(outA, 0, sizeof(outA));
        for(i=0; i<n; i++) {
            for(j=0; j<n; j++) {
                if(adj[i][j] == 1 && sch[i].mark != sch[j].mark) {
                    outA[sch[i].mark] = 1;
                    inA[sch[j].mark] = 1;
                }
            }
        }
        if(lable == 2) {    //当原图为强连通图时，要特殊处理！！
            printf("1\n0\n");
            continue;
        }
        for(i=1; i<lable; i++) {
            if(inA[i] == 0) inN ++;
            if(outA[i] == 0) outN ++;
        }
        printf("%d\n%d\n", inN, max(inN, outN));
    }
    return 0;
}