洛谷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】321 102 10021 32 3【样例2】421 1004 20021 23 4

输出样例#1：

【样例1】YES110【样例2】NO3

#include<iostream>#include<algorithm>#include<cstdio>#include<cstring>#define MAXN 3010#define MAXM 1000010#define MAX 1000000000using 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;}