HDU - 1827 Summer Holiday

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吧，但是因为刚学了kosaraju就用了一下。还被卡了输入，很蛋疼。

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>


using namespace std;

const int N=100010;

vector<int>mp[N];
vector<int>mpp[N];
int vis[N];
int val[N];
int num[N];
int color[N];
int degree[N];
int n,m,slg;

void dfs(int u)
{
    vis[u]=1;
    for(int i=0;i<mp[u].size();i++)
    {
        int v=mp[u][i];
        if(vis[v]==0)
            dfs(v);
    }
    num[slg++]=u;
}

void dfs2(int u)
{
    vis[u]=1;
    color[u]=slg;
    for(int i=0;i<mpp[u].size();i++)
    {
        int v=mpp[u][i];
        if(vis[v]==0)
            dfs2(v);
    }
}

void kosaraju()
{
    slg=1;
    memset(degree,0,sizeof(degree));
    memset(vis,0,sizeof(degree));
    for(int i=1;i<=n;i++)
    {
        if(vis[i]==0)
            dfs(i);
    }
    slg=0;
    memset(vis,0,sizeof(vis));
    for(int i=n;i>=1;i--)
    {
        if(vis[num[i]]==0)
        {
            slg++;
            dfs2(num[i]);
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<mp[i].size();j++)
        {
            if(color[i]!=color[mp[i][j]])
                degree[color[mp[i][j]]]++;
        }
    }

    int tot=0,output=0;
    for(int i=1;i<=slg;i++)
    {
        if(!degree[i])
        {
            tot++;
            int minn=0x3f3f3f3f;
            for(int j=1;j<=n;j++)
            {
                if(color[j]==i)
                {
                    minn=min(minn,val[j]);
                }
            }
            output+=minn;
        }
    }
    cout<<tot<<' '<<output<<endl;
}
int main()
{
    while(cin>>n>>m)
    {
        for(int i=1;i<=n;i++)
        {
            mp[i].clear();
            mpp[i].clear();
        }
        for(int i=1;i<=n;i++)
            scanf("%d",&val[i]);
        for(int i=0;i<m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            mp[x].push_back(y);
            mpp[y].push_back(x);
        }
        kosaraju();
    }
	return 0;

}