MYF

HDU 3836 Equivalent Sets

题目链接

HDU 3836

题目类型:强联通分量

题目分析

题目大意

n个点,m条边,问,最少再加上多少条边可以使图保持为强联通图。

解析

最少加上多少条边变为强联通图,那么,也就是说我们要消除所有入度为0的点以及出度为0的点,两者之间连线即可消掉一对,所以答案为两者个数的较大者。

如果一开始就是一个强联通分量,需要特判,因为所有点的sccno都是1,最后求出来的in[1]也是0,这样就不对了。

代码

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
#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=20000+5;
int dfs_clock;//时钟
int scc_cnt;//强连通分量总数
vector<int> G[maxn];//G[i]表示i节点指向的所有点
int pre[maxn]; //时间戳
int low[maxn]; //u以及u的子孙能到达的祖先pre值
int sccno[maxn];//sccno[i]==j表示i节点属于j连通分量 sccno[i]的区间为1~scc_cnt
stack<int> S;
void dfs(int u){
pre[u]=low[u]=++dfs_clock;
S.push(u);
for(int i=0;i<G[u].size();i++){
int v=G[u][i];
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 in[maxn],out[maxn];
void work(int n,int m){
for(int i=1;i<=n;i++) G[i].clear();
while(m--){
int u,v;
scanf("%d%d",&u,&v);//index range: 1~n
G[u].push_back(v);
}
find_scc(n);
memset(in, 0, sizeof(in));
memset(out, 0, sizeof(out));
for(int u = 1; u <= n; u++){
for(int i = 0; i < G[u].size(); i++){
int v = G[u][i];
int x = sccno[u], y = sccno[v];
if(x != y){
in[y]++;
out[x]++;
}
}
}
int a = 0, b = 0;
for(int i = 1 ; i <= scc_cnt ; i++ ){
if(in[i] == 0) a++;
if(out[i] == 0) b++;
}
cout<<(scc_cnt == 1? 0 : max(a, b))<<"\n";
}
int main(){
int n,m;
while(scanf("%d%d",&n,&m)!=EOF){
work(n,m);
}
}