hdu 1827 Summer Holiday Tarjan缩点+度数判断 解题报告

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

思路

没错！我又用的是vector存边！感觉挺好用的呢，就是不知道会不会炸空间。邻接表嘛。。。自然是不会抛弃的。
我们首先要Tarjan缩点重新建图，建图的同时统计每个点的入度。我们知道，所需要打电话的人的人数就是新图上面入度为0的点的个数。scc数组在Tarjan内部更新强连通分量的最小代价，最后统计的时候tot+=scc就好了。

我天各位亲朋好友们。千万记得，下次写了函数一定要调用！我的clear()函数就没有调用！！WA了两次才发现。。。。
具体见代码。

代码

话说
为什么用read()读入会挂掉？？

求助

PS：代码是AC的

#include<iostream>
#include<algorithm>
#include<stack>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
using namespace std; 
const int N=1000+5;
const int inf=0x7fffffff;
vector<int> ed[N],edge[N]; 
stack<int> s; 
int n,m,w[N],dfn[N],low[N],bl[N],tim,ans,sum; 
int in[N],num,scc[N],flag[N],instack[N]; 
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void Tarjan(int u)
{ 
    dfn[u]=low[u]=tim++;  
    s.push(u);flag[u]=1; 
    instack[u]=1; 
    for (int i=0;i<ed[u].size();i++) 
    {     
        int v=ed[u][i];       
        if (!flag[v]) 
        {               
            Tarjan(v);       
            low[u]=min(low[u],low[v]);  
        }          
        else if (instack[v]) low[u]=min(low[u],dfn[v]); 
    }    
    if (dfn[u]==low[u]) 
    {             
        scc[num]=inf;
        for (int v;v!=u;)
        {           
            v=s.top();s.pop();              
            instack[v]=0; 
            bl[v]=num;
            scc[num]=min(w[v],scc[num]);
        } 
        num++;  
    }
}
void clear()
{
    memset(flag,0,sizeof(flag)); 
    memset(instack,0,sizeof(instack)); 
    memset(low,0,sizeof(low));
    memset(dfn,0,sizeof(dfn));
    memset(scc,0,sizeof(scc));
    num=0;ans=0;sum=0;
    memset(in,0,sizeof(in)); 
}
int main() 
{
    while(scanf("%d %d",&n,&m)!=EOF) 
    {
        clear();
        for (int i=0;i<=n;i++) 
        {
            ed[i].clear();
            edge[i].clear();
        }
        for (int i=1;i<=n;i++) 
        w[i]=read();
        while(m--) 
        {
            int a,b;
            scanf("%d%d",&a,&b);
            ed[a].push_back(b);
        }
        for (int i=1;i<=n;i++) 
        if (!flag[i]) Tarjan(i);
        for (int i=1;i<=n;i++) 
        for (int j=0;j<ed[i].size();j++) 
        {
            int v=ed[i][j];
            if (bl[i]!=bl[v]) 
            {
                edge[bl[i]].push_back(bl[v]);
                in[bl[v]]++; 
            }
        }
        for (int i=0;i<num;i++) 
        if (in[i]==0) {ans++;sum+=scc[i];}
        printf("%d %d\n",ans,sum); 
    }
    return 0; 
}
/*
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
*/