POJ 1236 Network of Schools 强连通分量

链接：http://poj.org/problem?id=1236

题意：提供一个路径为单向的校园网络，给出两个任务，第一个是至少要多少台电脑接受新软件才能让信息传遍整个校园网络，第二个是任务是问至少要加多少条边让多少点都彼此相通。

思路：第一个是求入度为零的强连通分量，第二个答案是入度为零的强连通分量数和出度为零的强连通分量数更大的那个。（tarjan的思路看看书，模板网上漫天飘！）

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <cstdlib>
#include <cmath>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <vector>
#include <ctype.h>
#include <algorithm>
#include <string>
#define PI acos(-1.0)
#define maxn 105
#define INF 1<<25
#define MAX 0x7fffffff
#define mem(a,b) memset(a,b,sizeof(a))
#define f(i,a,b) for(i=a;i<b;i++)
typedef long long ll;
using namespace std;
struct Edge
{
    int v;
    int next;
} edge[maxn*maxn];
int head[maxn],dfn[maxn],low[maxn],belong[maxn],instack[maxn],in[maxn],ou[maxn];
int tt,top,cnt,v,A,B;
int add_edge(int u,int v)
{
    edge[top].v=v;
    edge[top].next=head[u];
    head[u]=top++;
}
stack <int> S;
int init()
{
    mem(head,-1);
    mem(dfn,0);
    mem(low,0);
    mem(belong,0);
    mem(instack,0);
    mem(in,0);
    mem(ou,0);
    A=B=top=cnt=tt=0;
}
int tarjan(int u)
{
    dfn[u]=low[u]=++cnt;
    S.push(u);
    instack[u]=1;
    for(int i=head[u]; i!=-1; i=edge[i].next)
    {
        int v=edge[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(instack[v]&&dfn[v]<low[u])
            low[u]=dfn[v];
    }
    if(dfn[u]==low[u])
    {
        tt++;
        do
        {
            v=S.top();
            instack[v]=0;
            belong[v]=tt;
            S.pop();
        }
        while(v!=u);
    }
}
int main()
{
    int tot,i;
    scanf("%d",&tot);
    init();
    f(i,1,tot+1)
    {
        int x;
        while(scanf("%d",&x))
        {
            if(!x)
                break;
            else add_edge(i,x);
        }
    }
    f(i,1,tot+1)
    if(!dfn[i])
        tarjan(i);
    f(i,1,tot+1)
    {
        int a=belong[i];
        for(int j=head[i]; j!=-1; j=edge[j].next)
        {
            int b=belong[edge[j].v];
            if(a!=b)
            {
                in[b]++;
                ou[a]++;
            }
        }
    }
    if(tt==1)
    {
        printf("1\n");
        printf("0\n");
    }
    else
    {
        f(i,1,tt+1)
        {
            if(!in[i])
                A++;
            if(!ou[i])
                B++;
        }
        printf("%d\n",A);
        printf("%d\n",max(A,B));
    }
}