hdoj1827Summer Holiday【scc+缩点】

Summer Holiday Time Limit: 10000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2371    Accepted Submission(s): 1111 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   Author 威士忌   Source HDOJ 2007 Summer Exercise（3）- Hold by Wiskey

求出scc后缩点求出入度为零的点的个数即为最少，求出到入度为零的scc的最小花费和即为最小花费

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<stack>
#include<vector>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100010;
int MIN(int a,int b){
	return a<b?a:b;
}
int value[maxn];
int scc_cnt,dfs_clock;
int sccno[maxn];
int dfn[maxn];
int low[maxn];
int head[maxn];
int num[maxn];
int in[maxn];
int out[maxn];
bool instack[maxn];
stack<int>S;
vector<int>scc[maxn];
vector<int>G[maxn];
struct Node{
	int from,to,next;
}A[maxn];
void init(){
	memset(head,-1,sizeof(head));
	memset(low,0,sizeof(low));
	memset(dfn,0,sizeof(dfn));
	memset(instack,false,sizeof(instack));
	memset(sccno,0,sizeof(sccno));
	memset(in,0,sizeof(in));
	memset(out,0,sizeof(out));
	memset(value,0x3f,sizeof(value));
	scc_cnt=dfs_clock=0;
}
void tarjan(int u,int pre){
	low[u]=dfn[u]=++dfs_clock;
	S.push(u);instack[u]=true;
	int v;
	for(int k=head[u];k!=-1;k=A[k].next){
		v=A[k].to;
		if(!dfn[v]){
			tarjan(v,u);
			low[u]=MIN(low[u],low[v]);
		}
		else if(instack[v]){
			low[u]=MIN(low[u],dfn[v]);
		}
	}
	if(low[u]==dfn[u]){
		scc_cnt++;scc[scc_cnt].clear();
		while(1){
			v=S.top();S.pop();
			instack[v]=false;
			sccno[v]=scc_cnt;
			G[scc_cnt].clear();
			scc[scc_cnt].push_back(v);
			value[scc_cnt]=MIN(value[scc_cnt],num[v]);
			if(v==u)break;
		
		}
	}
}
void suodian(int m){
	for(int i=0;i<m;++i){
		int u=sccno[A[i].from];
		int v=sccno[A[i].to];
		if(u!=v){
			G[u].push_back(v);
			in[v]++;out[u]++;
		}
	}
}
int main()
{
	int n,m,i,j,k;
	while(scanf("%d%d",&n,&m)!=EOF){
		init();
		for(i=1;i<=n;++i){
			scanf("%d",&num[i]);
		}
		int a,b;
		for(i=0;i<m;++i){
			scanf("%d%d",&a,&b);
			A[i].from=a;
			A[i].to=b;
			A[i].next=head[a];
			head[a]=i;
		}
		for(i=1;i<=n;++i){
			if(!dfn[i])tarjan(i,-1);
		}
		suodian(m);
		int ans1=0,ans2=0;
		for(i=1;i<=scc_cnt;++i){
			if(in[i]==0){
				ans1++;
				ans2+=value[i];
			}
		}
		printf("%d %d\n",ans1,ans2);
	}
	return 0;
}