[洛谷1262]间谍网络

传送门洛谷1262

Description

由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果A间谍手中掌握着关于B间谍的犯罪证据，则称A可以揭发B。有些间谍收受贿赂，只要给他们一定数量的美元，他们就愿意交出手中掌握的全部情报。所以，如果我们能够收买一些间谍的话，我们就可能控制间谍网中的每一分子。因为一旦我们逮捕了一个间谍，他手中掌握的情报都将归我们所有，这样就有可能逮捕新的间谍，掌握新的情报。
我们的反间谍机关提供了一份资料，色括所有已知的受贿的间谍，以及他们愿意收受的具体数额。同时我们还知道哪些间谍手中具体掌握了哪些间谍的资料。假设总共有n个间谍(n不超过3000)，每个间谍分别用1到3000的整数来标识。
请根据这份资料，判断我们是否有可能控制全部的间谍，如果可以，求出我们所需要支付的最少资金。否则，输出不能被控制的一个间谍。

Input

一行只有一个整数n。
第二行是整数p。表示愿意被收买的人数，1<=p<=n。
接下来的p行，每行有两个整数，第一个数是一个愿意被收买的间谍的编号，第二个数表示他将会被收买的数额。这个数额不超过20000.
紧跟着一行只有一个整数r，1<=r<=8000。然后r行，每行两个正整数，表示数对（A，B），A间谍掌握B间谍的证据。

Output

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

Sample Input

#1

3
2
1 10
2 100
2
1 3
2 3

#2

4
2
1 100
4 200
2
1 2
3 4

Sample Output

#1

YES
110

#2

NO
3

Solution

• 所有间谍构成有向图，当且仅当有入度为零且不可收买的结点时无解
• 对于一个强连通分量内的间谍，控制一个即可
• 强连通缩点重建图，记录每个scc内最高价，统计入读为零的scc答案

---

Code

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
struct JOKER{
    int b,next;
}e[16007];
int n,m,k,tot,ans,cnt,head[3007],a[3007],init[3007],cost[3007];
int stack[3007],top,scc,Time,dfn[3007],low[3007],id[3007];
bool vis[3007],ins[3007],ct[3007];
void add(int x,int y){
    e[++tot].b=y;
    e[tot].next=head[x];
    head[x]=tot;
}
void Tarjan(int u){
    vis[u]=ins[u]=1;
    dfn[u]=low[u]=++Time;
    stack[++top]=u;
    for (int i=head[u];i;i=e[i].next){
        int v=e[i].b;
        if (!vis[v]){
            Tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if (ins[v]) low[u]=min(low[u],dfn[v]);
    }
    if (dfn[u]==low[u]){
        int v=0; scc++;
        while (v!=u){
            v=stack[top--];
            ins[v]=0;
            id[v]=scc;
            if (a[v]!=-1) cost[scc]=min(cost[scc],a[v]);
        }
    }
}
void dfs(int u){
    vis[u]=1; cnt++;
    for (int i=head[u];i;i=e[i].next){
        int v=e[i].b;
        if (!vis[v]) dfs(v);
    }
}
int main(){
    scanf("%d%d",&n,&k);
    memset(a,-1,sizeof(a));
    for (int i=1;i<=k;i++){
        int x; scanf("%d",&x);
        scanf("%d",&a[x]);
    }
    scanf("%d",&m);
    for (int i=1;i<=m;i++){
        int x,y; scanf("%d%d",&x,&y);
        add(x,y);
    }
    for (int i=1;i<=n;i++)
        if (a[i]!=-1&&!vis[i]) dfs(i);
    if (cnt<n){
        printf("NO\n");
        for (int i=1;i<=n;i++)
            if (!vis[i]){
                printf("%d\n",i);
                return 0;
            }
    }

    for (int i=1;i<=n;i++){vis[i]=0; cost[i]=1e9+7;}    
    for (int i=1;i<=n;i++)
        if (!vis[i]) Tarjan(i);
    for (int u=1;u<=n;u++){
        for (int i=head[u];i;i=e[i].next){
            int v=e[i].b;
            if (id[u]!=id[v]) init[id[v]]++;
        }
    }
    for (int i=1;i<=scc;i++)
        if (!init[i]) ans+=cost[i];
    printf("YES\n");
    printf("%d",ans);
    return 0;
}