MYF

HDU 3072 Intelligence System

Posted on | In | | Visitors

题目链接

HDU 3072

题目类型：强联通分量

题目分析

题目大意

给你一个有向网络,网络中的每条有向边都有一个代价,表示从u向v发信息的代价.现在你要从0号点将信息发到所有的其他点去,问你最小代价是多少.其中如果u与v点可以互达(即属于同一个强连通分量),那么他们之间的通信不需要花代价. 输入保证从0号点能到达所有其他点.

解析

DAG，缩点，然后统计出到每个点分量的最短距离，然后从1到scc_cnt统计一下即可。

参考文章：HDU 3072 Intelligence System(强连通分量)

代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <iomanip>
#include <bitset>
#include <cstring>
#include <iostream>
#include <iosfwd>
#include <deque>
#include <algorithm>
#define Memset(a,val) memset(a,val,sizeof(a))
#define PI acos(-1.0)
#define PB push_back
#define MP make_pair
#define rt(n)       (i == n ? '\n' : ' ')
#define hi         printf("Hi----------\n")
#define fre() freopen("data_in.txt","r",stdin); \
			  freopen("data_out.txt","w",stdout);
#define debug(x) cout<<"Debug : ---"<<x<<"---"<<endl;
#define debug2(x,y) cout<<"Debug : ---"<<x<<" , "<<y<<"---"<<endl;
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
typedef pair<int,int> PII;
typedef long long ll;
const int mod=1000000007;
const int INF=0x3f3f3f3f;
const double eps=1e-8;
const int maxn=50000+5;
int dfs_clock;//时钟
int scc_cnt;//强连通分量总数
int pre[maxn]; //时间戳
int low[maxn]; //u以及u的子孙能到达的祖先pre值
int sccno[maxn];//sccno[i]==j表示i节点属于j连通分量 sccno[i]的区间为1~scc_cnt
int w[maxn];
stack<int> S;
struct Edge{
	int v,w,next;
}e[100000+10];
int he[maxn],ecnt= 0;
void adde(int u,int v,int w){
	e[ecnt].v = v;
	e[ecnt].w = w;
	e[ecnt].next = he[u];
	he[u] = ecnt++;
}
void dfs(int u){
	pre[u]=low[u]=++dfs_clock;
	S.push(u);
	for (int i = he[u];i != -1;i = e[i].next) {
		int v = e[i].v;
		if(!pre[v]){
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!sccno[v]){
			low[u]=min(low[u],pre[v]);
		}
	}
	if(low[u] == pre[u]){//u为当前强连通分量的入口
		scc_cnt++;
		while(true){
			int x=S.top(); S.pop();
			sccno[x]=scc_cnt;
			if(x==u) break;
		}
	}
}

//求出有向图所有连通分量
void find_scc(int n){
	scc_cnt=dfs_clock=0;
	memset(sccno,0,sizeof(sccno));
	memset(pre,0,sizeof(pre));
	for(int i=1;i<=n;i++)
		if(!pre[i]) dfs(i);
}
int val[maxn];
void work(int n,int m){
	memset(he, -1, sizeof(he));
	ecnt = 0;
	for (int i = 1;i<=m;i++) {
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);//index range: 1~n
		u++, v++;
		adde(u, v, w);
	}
	find_scc(n);
	memset(val, INF, sizeof(val));
	for(int u = 1; u <= n; u++){
		for(int i = he[u]; i != -1; i = e[i].next){
			int v = e[i].v;
			int w = e[i].w;
			int x = sccno[u] , y = sccno[v];
			if(x != y) val[y] = min(val[y], w);
		}
	}
	ll ans = 0;
	for(int i = 1; i <= scc_cnt; i++){
		ans += val[i];
	}
	ans -= INF;
	cout<<ans<<"\n";
}
int main(){
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF){
		work(n,m);
	}
}