poj1236:Network of School（Tarjan）

• Description
  一些学校连入一个电脑网络。那些学校已订立了协议：每个学校都会给其它的一些学校分发软件（称作“接受学校”）。注意如果 B 在 A 学校的分发列表中，那么 A 不必也在 B 学校的列表中。
  你要写一个程序计算，根据协议，为了让网络中所有的学校都用上新软件，必须接受新软件副本的最少学校数目（子任务A）。更进一步，我们想要确定通过给任意一个学校发送新软件，这个软件就会分发到网络中的所有学校。为了完成这个任务，我们可能必须扩展接收学校列表，使其加入新成员。
  请你计算最少需要增加几个扩展，使得不论我们给哪个学校发送新软件，它都会到达其余所有的学校（子任务B）。一个扩展就是在一个学校的接收学校列表中引入一个新成员。

• 题解：

先是几个概念

强连通分量

在有向图G中，如果两个顶点(u,v)间存在一条路径u到v的路径且也存在一条v到u的 路径，则称这两个顶点 (u,v)是强连通的 (strongly connected) 。如果有向图G的每两个顶点都 强连通，称 是一个强连通图。有向非强连通图的极大强连通子图，称为强连通分量(strongly connectedcomponents)。若将有向图中的强连通分量都缩为一个点，则原图会形成一个DAG （有向无环图）。

Tarjan - 一种求强连通分量的算法

（由于字数太多，图片太多，直接放上一张图片。。）
1.Tarjan算法框架

2.代码主体

3.补充
上面的说明仍然有些抽象。可以用一种理解方式，即想象一颗搜索树。

假设A子树中没有前向边可以到达根节点,则：1.A子树中的任何点无法到达B子树中任何点。2.对于B子树中的节点中的横叉边a到A子树中的点肯定已经退栈。那么此时这条横叉边无效。
假设B子树中存在前向边可以到达根节点,则：1.B子树中的点可以到达C子树中点。2.对于C子树中的节点中的横叉边b到B子树中的点没有退栈。那么此时这条横叉边有效。

#include<bits/stdc++.h>
using namespace std;

const int Maxn=1e4+50;

inline int read()
{
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<3)+(i<<1)+ch-'0';ch=getchar();}
    return i*f;
}

int n,dfn[Maxn],low[Maxn],st[Maxn],ins[Maxn],top,ind,cnt,out[Maxn],in[Maxn],bl[Maxn];
vector<int>edge[Maxn];

inline void dfs(int now)
{
    dfn[now]=low[now]=++ind;
    st[++top]=now;
    ins[now]=1;
    for(int e=edge[now].size()-1;e>=0;e--)
    {
        int v=edge[now][e];
        if(!dfn[v])
        {
            dfs(v);
            low[now]=min(low[now],low[v]);
        }
        else if(ins[v])
        {
            low[now]=min(low[now],dfn[v]);
        }
    }
    if(low[now]==dfn[now])
    {
        ++cnt;
        while(st[top]!=now)
        {
            ins[st[top]]=0;
            bl[st[top--]]=cnt;
        }
        ins[st[top]]=0;
        bl[st[top--]]=cnt;
    }
}

int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        int j=read();
        while(j)
        {
            edge[i].push_back(j);
            j=read();
        }
    }
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])dfs(i);
    }
    for(int i=1;i<=n;i++)
    {
        for(int e=edge[i].size()-1;e>=0;e--)
        {
            if(bl[i]!=bl[edge[i][e]])out[bl[i]]++,in[bl[edge[i][e]]]++;
        }
    }
    int a=0,b=0;
    for(int i=1;i<=cnt;i++)
    {
        if(!in[i])++a;
        if(!out[i])++b;
    }
    cout<<a<<endl;
    cnt?(cout<<0<<endl):(cout<<max(a,b)<<endl);
}