HDU -1827 Summer Holiday（Tarjan，缩点）

Summer Holiday

Time Limit: 10000/1000 MS (Java/Others)	Memory Limit: 32768/32768 K (Java/Others)

Problem Description
To see a World in a Grain of Sand
And a Heaven in a Wild Flower,
Hold Infinity in the palm of your hand
And Eternity in an hour.
—— William Blake

听说lcy帮大家预定了新马泰7日游，Wiskey真是高兴的夜不能寐啊，他想着得快点把这消息告诉大家，虽然他手上有所有人的联系方式，但是一个一个联系过去实在太耗时间和电话费了。他知道其他人也有一些别人的联系方式，这样他可以通知其他人，再让其他人帮忙通知一下别人。你能帮Wiskey计算出至少要通知多少人，至少得花多少电话费就能让所有人都被通知到吗？

Input
多组测试数组，以EOF结束。
第一行两个整数N和M（1<=N<=1000, 1<=M<=2000），表示人数和联系对数。
接下一行有N个整数，表示Wiskey联系第i个人的电话费用。
接着有M行，每行有两个整数X，Y，表示X能联系到Y，但是不表示Y也能联系X。

Output
输出最小联系人数和最小花费。
每个CASE输出答案一行。

Sample Input
12 16
2 2 2 2 2 2 2 2 2 2 2 2
1 3
3 2
2 1
3 4
2 4
3 5
5 4
4 6
6 4
7 4
7 12
7 8
8 7
8 9
10 9
11 10

Sample Output
3 6

直接套模板呗，相较于一般的Tarjan，多了一个缩点，个人理解，每个强连通分量都是一个整体，那么一张图可以看作由很多个“整体”组成，将这些“整体”看作点并编号，处理点与点就方便多了。
“整体”的编号：cat
每个“整体”所含元素：point_list[cat]
每个“整体”的出度：out[cat]

#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <vector>
#include <stack>
using namespace std;
#define inf 99999999
#define N 1005

int Low[N];//最小子节点编号
int Bfn[N];//访问顺序，时间戳
int m, n;
int cnt, cat;//时间戳、强连通分量的数量
int vit[N];//标记
vector<int>Map[N];//每个点连接的点
stack<int>S;
int point_setnum[N];//每个点所在强连通分量的编号，别的程序中的scc[]
vector<int> point_list[N];//每个强连通分量所包含的元素
int out[N];//出度
int price[N];//价钱

void tarjan(int x)
{
    Low[x]=Bfn[x]=++cnt;
    vit[x]=1;
    S.push(x);
    for(int i=0; i<Map[x].size(); i++)
    {
        int t=Map[x][i];
        if(!Bfn[t])
        {
            tarjan(t);
            Low[x]=min(Low[x], Low[t]);
        }
        else if(vit[t])
            Low[x]=min(Low[x], Bfn[t]);
    }
    if(Low[x]==Bfn[x])
    {
        cat++;
        while(!S.empty())
        {
            int t=S.top();
            S.pop();
            vit[t]=0;
            point_list[cat].push_back(t);//将元素压入所在强连通分量
            point_setnum[t]=cat;//标记点所在强连通分量的编号
            if(x==t)
                break;
        }
    }
}

void solve()
{
    for(int i=1; i<=n; i++)
        if(!Bfn[i])
            tarjan(i);
}

void init()
{
    memset(Low, 0, sizeof(Low));
    memset(Bfn, 0, sizeof(Bfn));
    memset(vit, 0, sizeof(vit));
    memset(out, 0, sizeof(out));
    memset(point_setnum, 0, sizeof(point_setnum));
    cnt=cat=0;
    while(!S.empty())
        S.pop();
    for(int i=0; i<=n; i++)
        Map[i].clear(), point_list[i].clear();
}

int Mapmin(int i)//查找每个强连通分量中价钱最小的，联系他
{
    int ans=inf;
    for(int j=0; j<point_list[i].size(); j++)
        ans=min(ans, price[point_list[i][j]]);//第i个强连通分量中的第j个元素的价钱
    return ans;
}

int main()
{
    while(scanf("%d%d", &n, &m)!=EOF)
    {
        init();
        for(int i=1; i<=n; i++)
            scanf("%d", &price[i]);
        for(int i=0; i<m; i++)
        {
            int a, b;
            scanf("%d%d", &a, &b);
            Map[a].push_back(b);
        }
        solve();
        for(int i=1; i<=n; i++)//将每个强连通分量看作一个整体，统计每个整体的出度
        {
            for(int j=0; j<Map[i].size(); j++)
            {
                int t=Map[i][j];
                if(point_setnum[i]!=point_setnum[t])
                    out[point_setnum[t]]++;
            }
        }
        int pricesum=0, kidsum=0;;
        for(int i=1; i<=cat; i++)
            if(out[i]==0)
                pricesum+=Mapmin(i), kidsum++;

        printf("%d %d\n", kidsum, pricesum);
    }
    return 0;
}