P1412【分类练习5.图论】间谍网络

问题描述

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

样例输入2：
4
2
1 100
4 200
2
1 2
3 4

样例输出

样例输出1：
YES
110

样例输出2：
NO
3

题解

先tajan将强连通分量缩点，在并用一个数组记录一个强连通分量的最小值，再以此判断新生成图中入度为0的点（或带有间谍的点）的值相加即是答案。

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
#include<stack>
using namespace std;
#define maxn1 10005
#define inf 1e9
#define f(x,y) for(i=x;i<=y;i++)
struct node{
    int a,b;
};
node p[maxn1];
int n,Last[maxn1],Next[maxn1],Len[maxn1],End[maxn1];
int be[maxn1],ru[maxn1];
int dfs[maxn1],low[maxn1];
int v[maxn1];
int pp,m,cnt;
int cost[maxn1],mincost[maxn1];
int scc,vis;
bool instack[maxn1],check[maxn1];
stack<int>s;
void tj(int u)
{
    int i,j,k;
    dfs[u]=low[u]=++vis;
    s.push(u);
    instack[u]=true;
    for(i=Last[u];i;i=Next[i])
    {
        int en=End[i];
        if(dfs[en]==0)
        {
            tj(en);
            low[u]=min(low[u],low[en]);
        }
        else if(dfs[en]&&instack[en])
        {
            low[u]=min(low[u],dfs[en]);
                }       
    }
    int v;
    if(dfs[u]==low[u])
    {
      scc++;
      mincost[scc]=inf;
      do{
         v=s.top();
        s.pop();
        be[v]=scc;
        mincost[scc]=min(mincost[scc],cost[v]);
        instack[v]=false;
      }while(u!=v);
     } 
}
void insert(int x,int y)
{
    cnt++;
    Next[cnt]=Last[x];
    Last[x]=cnt;
    End[cnt]=y;
}
int main()
{
    int i,j,k;
    scanf("%d%d",&n,&pp);
    for(i=1;i<=n;i++) cost[i]=inf;
    for(i=1;i<=pp;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        cost[x]=y;
    }
    scanf("%d",&m);
    for(i=1;i<=m;i++) {
        int x,y;
        scanf("%d%d",&x,&y);
        insert(x,y);
        p[i].a=x;p[i].b=y;
        ru[y]++;
    }
    for(i=1;i<=n;i++){
        if(!dfs[i]) tj(i);
    }//tarjan
//    f(1,scc)cout<<mincost[i]<<" ";
    for(i=1;i<=m;i++)
    {
        if(be[p[i].a]!=be[p[i].b])
        {
            check[be[p[i].b]]=true;
        }
    }
    bool flag=false;
    int ans=0;
    for(i=1;i<=n;i++)
    {
        if(check[be[i]]==false&&mincost[be[i]]==inf)
        {
            flag=true;
            break;
        }
        else if(check[be[i]]==false){
            ans+=mincost[be[i]];
            check[be[i]]=true;
        }
    }
    if(flag) cout<<"NO"<<endl<<i;
    else cout<<"YES"<<endl<<ans;
}