Summer Holiday 【有向图中连通最少的点来使其整个图 连通的】+【tarjan求SCC +缩点】

Summer Holiday
Time Limit: 10000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2088 Accepted Submission(s): 966

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

学习了 —– 缩点和scc的入度处理
思路： tarjan 扫一遍，然后缩点 ；对于每个scc ，如果其入度为0.则只能从内部 打通一个点，否则入度不为零，就一定 会有入度为0的点来 使其连通。。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<set>
#define CLR(a,b) memset((a),(b),sizeof(a))
#define inf 0x3f3f3f3f
#define mod 100009
#define LL long long
#define M  2000+100
#define ll o<<1
#define rr o<<1|1
#define lson o<<1,l,mid
#define rson o<<1|1,mid+1,r
using namespace std;
void read(int &x){
    x=0;char c;
    while((c=getchar())<'0');
    do x=x*10+c-'0';while((c=getchar())>='0');
}

struct Edge {
    int from,to,next;
}edge[M];
int head[M],top;
vector<int >G[M]; // 缩点 后的新图 
vector<int>scc[M]; //  存储每个SCC中的点
int low[M],dfn[M];
int dfs_clock;
int Instack[M];
stack<int>S;
int sccno[M],scc_cnt; // i 点是属于哪个scc ，scc的计数器 
int in[M]; // 每个scc的入度 
int cost[M];//每个点权值
int n,m;

void init()
{
    top=0;
    memset(head,-1,sizeof(head));
    while(!S.empty()) S.pop();  
} 
void addedge(int a,int b)
{
    Edge e={a,b,head[a]};
    edge[top]=e;head[a]=top++;
}
void getmap()
{
    int i,j,a,b;
    for(i=1;i<=n;i++)
    read(cost[i]);
    while(m--)
    {
        read(a);read(b);
        addedge(a,b);
    }
}
void tarjan(int now)
{
    int nexts;
    dfn[now]=low[now]=++dfs_clock;
    S.push(now);Instack[now]=1;
    for(int i=head[now];i!=-1;i=edge[i].next)
    {
        Edge e=edge[i];
        if(!dfn[e.to])  
        {
            tarjan(e.to);
            low[now]=min(low[e.to],low[now]);
        }
        else if(Instack[e.to])
        low[now]=min(low[now],dfn[e.to]);
    }
    if(low[now]==dfn[now])
    {
        scc_cnt++;
        scc[scc_cnt].clear();   //  这个只要在这里清空就可以了 
        for(;;)
        {
            nexts=S.top();S.pop();
            sccno[nexts]=scc_cnt;
            Instack[nexts]=0;
            scc[scc_cnt].push_back(nexts);  // 存储每个scc里的点 
            if(nexts==now) break;
        }
    }
}
void find_cut(int le,int ri)
{
    int i,j;
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(sccno,0,sizeof(sccno));
    memset(Instack,0,sizeof(Instack));
    scc_cnt=dfs_clock=0;
    for(i=le;i<=ri;i++)
    if(!dfn[i]) tarjan(i);
 } 
void suodian()
{
    for(int i=1;i<=scc_cnt;i++) G[i].clear(),in[i]=0;

    for(int i=0;i<top;i++)  // 一个scc表示一个点来进行缩点 
    {
        int now=sccno[edge[i].from];
        int nexts=sccno[edge[i].to];
        if(now!=nexts) 
        {
            G[now].push_back(nexts); 
            in[nexts]++; 
        }
    }
}
void solve()
{
    if(scc_cnt==1) // 则这个图就是一个强联通图，找到一个人就可以传遍所有人，当然找最便宜的 
    {
        sort(cost+1,cost+1+n);
        printf("%d %d\n",1,cost[1]); 
     } 
    else 
    {
        int ans=0;
        int mincost=0;
        for(int i=1;i<=scc_cnt;i++)         
        {
            if(in[i]) continue; // 对于入度为为1的 scc 那么一定可以 入度为0的scc 传消息进来，就不用花钱了  
            ans++;
            int money=inf;// 对于入度为0的来说 就必须要找一个人来打电话，当然找最便宜的 了 
            for(int j=0;j<scc[i].size();j++)  // 遍历 第 i scc，找到最便宜 
            money=min(money,cost[scc[i][j]]);
            mincost+=money;     
        }
        printf("%d %d\n",ans,mincost);
    }
}
int main()
{
    while(~scanf("%d%d",&n,&m)) 
    {
        init();
        getmap();
        find_cut(1,n);
        suodian();
        solve();
     }
    return 0;
}