洛谷P1262间谍网络

题目描述

由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果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：

【样例1】
YES
110
【样例2】
NO
3

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define MAXN 3010
#define MAXM 1000010
#define MAX 1000000000
using namespace std;
int n,m,top=1,c=1,d=0,ans=0;
int indegree[MAXN],money[MAXN],deep[MAXN],low[MAXN],cstack[MAXN],minmoney[MAXN],rink[MAXN],head[MAXN];
bool beque[MAXN];
struct node{
       int next,to;
}a[MAXM];
inline int read(){
       int date=0,w=1;char c=0;
       while(c!='-'&&(c<'0'||c>'9'))c=getchar();
       if(c=='-'){w=-1;c=getchar();}
       while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
       return date*w;
}
void add(int x,int y){
     a[c].to=y;
     a[c].next=head[x];
     head[x]=c++;
}
void dfs(int x){
     if(beque[x])return;
     beque[x]=true;
     d++;
     for(int i=head[x];i;i=a[i].next)
     dfs(a[i].to);
}
void work(int x){
     int t;
     beque[x]=true;
     deep[x]=low[x]=d++;
     cstack[top++]=x;
     for(int i=head[x];i;i=a[i].next){
             t=a[i].to;
             if(!deep[t]){
                          work(t);
                          low[x]=min(low[x],low[t]);
                          }
             else if(beque[t])
             low[x]=min(low[x],deep[t]);
             }
     if(low[x]==deep[x]){
                         minmoney[x]=MAX;
                         do{
                             beque[cstack[top-1]]=false;
                             rink[cstack[top-1]]=x;
                             minmoney[x]=min(minmoney[x],money[cstack[top-1]]);
                             }while(cstack[--top]!=x);
                         }
}
int main(){
    int x,y;
    n=read();m=read();
    for(int i=1;i<=n;i++){
            money[i]=MAX;
            beque[i]=false;
            indegree[i]=0;
            }
    for(int i=1;i<=m;i++){
            x=read();y=read();
            money[x]=y;
            }
    m=read();
    for(int i=1;i<=m;i++){
            x=read();y=read();
            add(x,y);
            }
    for(int i=1;i<=n;i++)
    if(money[i]!=MAX)
    dfs(i);
    if(d<=n)
    for(int i=1;i<=n;i++)
    if(!beque[i]){
                  printf("NO\n%d\n",i);
                  return 0;
                  }
    d=1;
    memset(beque,false,sizeof(beque));
    for(int i=1;i<=n;i++)
    if(money[i]!=MAX&&!deep[i])
    work(i);
    for(int i=1;i<=n;i++)
    for(int j=head[i];j;j=a[j].next)
    if(rink[i]!=rink[a[j].to])
    indegree[rink[a[j].to]]++;
    for(int i=1;i<=n;i++)
    if(rink[i]==i&&!indegree[i])
    ans+=minmoney[i];
    printf("YES\n%d\n",ans);
    return 0;
}