洛谷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，并在第二行输出不能控制的间谍中，编号最小的间谍编号。

我们可以把这个图抽象成几个连通分量连在一起。对于每一个单个的连通分量，搞定它需要花的代价是其中最小的。而如果这些连通分量连在一起，那些没有入边的连通分量一定是要搞定的。

而搞定这些联通分量以后，后面的就不用考虑了。所以我们先Tarjan缩点，对于每一个入度为0的联通分量，ans+它的值。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cstring>
#include <stack>
#include <queue> 
#define in(a) a=read()
#define MAXN 10010
#define REP(i,k,n) for(int i=k;i<=n;i++)
using namespace std;
inline int read(){
    int x=0,f=1;
    char ch=getchar();
    for(;!isdigit(ch);ch=getchar())
        if(ch=='-')
            f=-1;
    for(;isdigit(ch);ch=getchar())
        x=x*10+ch-'0';
    return x*f;
}
int n,p,r,a,b,num,ans;
int val[MAXN],indu[MAXN];
int total,head[MAXN],to[MAXN],nxt[MAXN];
int dfn[MAXN],low[MAXN],cnt=0,vis[MAXN],bel[MAXN];
int minn[MAXN],book[MAXN];
inline void adl(int a,int b){
    total++;
    to[total]=b;
    nxt[total]=head[a];
    head[a]=total;
    return ;
}
stack <int> S;
inline void tarjan(int u){
    low[u]=dfn[u]=++cnt;
    S.push(u);
    vis[u]=1;
    for(int e=head[u];e;e=nxt[e]){
        if(!dfn[to[e]]){
            tarjan(to[e]);
            low[u]=min(low[u],low[to[e]]);
        }
        else  if(vis[to[e]])  low[u]=min(low[u],dfn[to[e]]);
    }
    if(dfn[u]==low[u]){
        num++;
        while(!S.empty() && S.top()!=u)  bel[S.top()]=num,minn[num]=min(minn[num],val[S.top()]),vis[S.top()]=0,S.pop();
        if(!S.empty())  bel[S.top()]=num,minn[num]=min(minn[num],val[S.top()]),vis[S.top()]=0,S.pop();
    }
    return ;
}
int main(){
    in(n),in(p);
    REP(i,1,n)  minn[i]=val[i]=2147483647;
    REP(i,1,p)  in(val[in(a)]);
    in(r);
    REP(i,1,r)  in(a),in(b),adl(a,b);
    REP(i,1,n)  if(!dfn[i] && val[i]!=2147483647)  tarjan(i);
    REP(i,1,n)  if(!dfn[i]){
        cout<<"NO"<<endl<<i;
        return 0;
    }
    cout<<"YES"<<endl;
    REP(u,1,n)
        for(int e=head[u];e;e=nxt[e])
            if(bel[u]!=bel[to[e]])
                indu[bel[to[e]]]++;
    REP(i,1,num)
        if(!indu[i])  ans+=minn[i];
    cout<<ans;
    return 0;
}