题面：间谍网络

大秦为你打开题目传送门

题目描述

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

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

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

输入格式

第一行只有一个整数 $n$。

第二行是整数 $p$。表示愿意被收买的人数，$1\le p\le n$。

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

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

输出格式

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

样例 #1

样例输入 #1

1
2
3
4
5
6
7
3
2
1 10
2 100
2
1 3
2 3

样例输出 #1

1
2
YES
110

样例 #2

样例输入 #2

1
2
3
4
5
6
7
4
2
1 100
4 200
2
1 2
3 4

样例输出 #2

1
2
NO
3

思路

首先，我们发现一个间谍被发现只有两种可能：

  1. 被贿赂
  2. 被揭发

而一个强连通分量中的所有间谍，只要有一个被贿赂，就会全部被揭发，也就是说我们要计算出每一个强连通分量中最小的贿赂值。但是说所有强连通分量都有必要贿赂吗？显然一个强连通分量还可以被另外的强连通分量揭发，这样我们就少贿赂了一些人。那么必须贿赂的，就是那些入度为 0 的强连通分量，如果存在不可贿赂的入度为 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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
// typedef __int128 int128;
typedef pair<int , int > pii;
typedef unsigned long long ull;

namespace FastIO 
{
    // char buf[1 << 20], *p1, *p2;
    // #define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 20, stdin), p1 == p2) ? EOF : *p1++)
    template<typename T> inline T read(T& x) {
        x = 0;
        int f = 1;
        char ch;
        while (!isdigit(ch = getchar())) if (ch == '-') f = -1;
        while (isdigit(ch)) x = (x << 1) + (x << 3) + (ch ^ 48), ch = getchar();
        x *= f;
        return x;
    }
    template<typename T, typename... Args> inline void read(T& x, Args &...x_) {
        read(x);
        read(x_...);
        return;
    }
    inline ll read() {
        ll x;
        read(x);
        return x;
    }
};
using namespace FastIO;

const int MaxV = 1e5 + 10;
const int MaxE = 5e5 + 10;

template<int N, int M>
class Graph {
    private :
        struct Edge {
            int to, nt, wt;
            Edge() {}
            Edge(int to, int nt, int wt) : to(to), nt(nt), wt(wt) {}
        }e[M];
        int hd[N], cnte;
    public :
        inline void clear() { memset(hd, 0, sizeof(hd)), cnte = 0; }
        inline void AddEdge(int u, int v, int w = 0) {
            e[++cnte] = Edge(v, hd[u], w);
            hd[u] = cnte;
        }
        inline int head(int u) { return hd[u]; }
        inline int nt(int u) { return e[u].nt; }
        inline int to(int u) { return e[u].to; }
        inline int wt(int u) { return e[u].wt; }
};

int n, m, p;
int buy[MaxV];
int uid[MaxE], vid[MaxE];
Graph< MaxV, MaxE >G;

inline void Input() {
    read(n, p);
    for(int i = 1; i <= n; i++) buy[i] = 1e9;
    int u, v;
    for(int i = 1; i <= p; i++) {
        read(u, v);
        buy[u] = min(buy[u], v);
    }
    read(m);
    for(int i = 1; i <= m; i++) {
        read(uid[i], vid[i]);
        G.AddEdge(uid[i], vid[i]);
    }
}

int dfn[MaxV], low[MaxV], cntd;
stack<int >stk; int instk[MaxV];
int cntc, col[MaxV];
vector<int >clb[MaxV];

void Tarjan(int u) {
    dfn[u] = low[u] = ++cntd;
    stk.push(u), instk[u] = 1;
    for(int i = G.head(u); i; i = G.nt(i)) {
        int v = G.to(i);
        if(!dfn[v]) { 
            Tarjan(v);
            low[u] = min(low[u], low[v]);
        }
        else if(instk[v]) { 
            low[u] = min(low[u], dfn[v]);
        }
    }
    if(low[u] == dfn[u]) {
        ++cntc;
        int v = stk.top(); stk.pop();
        while(v != u) {
            col[v] = cntc;
            instk[v] = 0;
            clb[cntc].push_back(v);
            v = stk.top(); stk.pop();
        }
        col[v] = cntc;
        instk[v] = 0;
        clb[cntc].push_back(v);
    }
}

int inD[MaxV], mi[MaxV];

inline void Init() {
    for(int i = 1; i <= m; i++) {
        if(col[uid[i]] == col[vid[i]]) continue;
        inD[col[vid[i]]]++;
    }
    for(int i = 1; i <= cntc; i++) {
        mi[i] = 1e9;
        sort(clb[i].begin(), clb[i].end());
        for(auto &u : clb[i]) {
            mi[i] = min(mi[i], buy[u]);
        }
    }
}

inline void Work() {
    for(int i = 1; i <= n; i++) {
        if(dfn[i] == 0) {
            Tarjan(i);
        }
    }
    Init();
    int flag = 0, ans = 0;
    for(int i = 1; i <= cntc; i++) {
        if(inD[i] != 0) continue;
        if(mi[i] == 1e9) {
            if(flag == 0) ans = n + 1;
            flag = 1;
            ans = min(ans, clb[i][0]);
        }
        else if(flag == 0){
            ans += mi[i];
        }
        // cout << "?! : " << i << ' ' << ans << ' ' << flag << endl;
    }
    if(flag == 1) {
        for(int i = 1; i <= ans; i++) {
            if(buy[i] == 1e9) ans = min(ans, i);
        }
    }
    if(flag == 1) printf("NO\n");
    else printf("YES\n");
    printf("%d", ans);
}

int main() {
    int T = 1;
    while(T--) {
        Input();
        Work();
    }
    return 0;
}