洛谷 P1262 间谍网络 Label: Kosarajn强联通

题目描述

由于外国间谍的大量渗入,国家安全正处于高度的危机之中。如果A间谍手中掌握着关于B间谍的犯罪证据,则称A可以揭发B。有些间谍收受贿赂,只要给他们一定数量的美元,他们就愿意交出手中掌握的全部情报。所以,如果我们能够收买一些间谍的话,我们就可能控制间谍网中的每一分子。因为一旦我们逮捕了一个间谍,他手中掌握的情报都将归我们所有,这样就有可能逮捕新的间谍,掌握新的情报。

我们的反间谍机关提供了一份资料,色括所有已知的受贿的间谍,以及他们愿意收受的具体数额。同时我们还知道哪些间谍手中具体掌握了哪些间谍的资料。假设总共有n个间谍(n不超过3000),每个间谍分别用1到3000的整数来标识。

请根据这份资料,判断我们是否有可能控制全部的间谍,如果可以,求出我们所需要支付的最少资金。否则,输出不能被控制的一个间谍。

输入输出格式

输入格式:

第一行只有一个整数n。

第二行是整数p。表示愿意被收买的人数,1≤p≤n。

接下来的p行,每行有两个整数,第一个数是一个愿意被收买的间谍的编号,第二个数表示他将会被收买的数额。这个数额不超过20000。

紧跟着一行只有一个整数r,1≤r≤8000。然后r行,每行两个正整数,表示数对(A, B),A间谍掌握B间谍的证据。

输出格式:

如果可以控制所有间谍,第一行输出YES,并在第二行输出所需要支付的贿金最小值。否则输出NO,并在第二行输出不能控制的间谍中,编号最小的间谍编号。

输入输出样例

输入样例#1:
【样例1】
3
2
1 10
2 100
2
1 3
2 3
【样例2】
4
2
1 100
4 200
2
1 2
3 4

代码

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4 #include<algorithm>
 5 #include<vector>
 6 #define MAX 100005
 7 #define INF 0x3f3f3f3f
 8 using namespace std;
 9 vector<int> G[MAX],rG[MAX],vs;//G存边,rG存反向边
10 int used[MAX],cmp[MAX],ans[MAX],c[MAX];//cmp映射所属联通块 ,c记录贿赂值,ans记录联通块最小值
11 int V,m,cant=INF,tot=0;
12 int k=1;
13 
14 void add(int from,int to){//反向的反向边 
15     G[from].push_back(to);
16     rG[to].push_back(from);
17 }
18 
19 void dfs(int v){
20     used[v]=1;
21     for(int i=0;i<G[v].size();i++)
22         if(!used[G[v][i]]) dfs(G[v][i]);
23     vs.push_back(v);
24 }
25 
26 void rdfs(int v,int k){
27     used[v]=1;
28     cmp[v]=k;
29     int tmp_cost=c[v];
30     for(int i=0;i<rG[v].size();i++)
31         if(!used[rG[v][i]]) rdfs(rG[v][i],k);
32     //该函数以上为Kosarajn核心代码 
33     for(int i=0;i<rG[v].size();i++){
34         if(cmp[rG[v][i]]==k)
35             tmp_cost=min(c[rG[v][i]],tmp_cost);
36         else{
37             c[rG[v][i]]=0;//反正该联通块会影响到此节点,直接标记
38             ans[cmp[rG[v][i]]]=0;//反正该联通块会影响到此节点,直接标记
39         }
40     }
41     ans[k]=min(ans[k],tmp_cost);
42 }
43 
44 int scc(){
45     memset(used,0,sizeof(used));vs.clear();
46     for(int v=1;v<=V;v++)
47         if(!used[v]) dfs(v);
48     memset(used,0,sizeof(used));
49     for(int i=vs.size()-1;i>=0;i--)
50         if(!used[vs[i]]) rdfs(vs[i],k++);
51     //该函数以上为Kosarajn核心代码 
52     for(int i=vs.size()-1;i>=0;i--)//查找无法到达的顶点 
53         if(ans[cmp[vs[i]]]==INF) cant=min(cant,vs[i]);
54 }
55 
56 int main(){
57     memset(ans,0x3f,sizeof(ans));
58     memset(c,0x3f,sizeof(c));
59     scanf("%d%d",&V,&m);
60     
61     for(int i=1;i<=m;i++){
62         int cost,to;
63         scanf("%d%d",&to,&cost);
64         c[to]=cost;
65     }
66     scanf("%d",&m);
67     for(int i=1;i<=m;i++){
68         int from,to;
69         scanf("%d%d",&from,&to);
70         add(to,from);//反向添加边,可以保证搜索的时候是正向的 
71     }
72     scc();
73     
74     if(cant!=INF){puts("NO");printf("%d
",cant);return 0;}
75     
76     puts("YES");
77     for(int i=1;i<k;i++)//计算总花费 
78         tot+=ans[i];
79     
80     printf("%d
",tot);
81     
82     return 0;
83 }

详见注释啦

版权声明: 本博客所有文章除特别声明外,均采用 CC BY-NC-SA 3.0 许可协议。转载请注明出处!
原文地址:https://www.cnblogs.com/radiumlrb/p/5874859.html