洛谷1262 间谍网络

题目描述

由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果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缩点，记下每个强连通分量里不能贿赂的最小编号和能贿赂的最小花费，然后从每个入度为零的点开始看能否遍历整张图。

#include<cstdio>#include<cstring>#include<algorithm>#include<vector>using namespace std;const int oo=0x3f3f3f3f;int fir[3010],ne[8010],to[8010],m,n,p,minv[3010],minn[3010],tot,val[3010],dfn[3010],low[3010],sta[3010],bel[3010],top,clo;bool b[3010],ins[3010],need[3010],vis[3010],map[3010][3010];vector<int> to2[3010];void init(){    int i,u,v;    scanf("%d%d",&n,&p);    for (i=1;i<=p;i++)    {        scanf("%d%d",&u,&v);        b[u]=1;        val[u]=v;    }    scanf("%d",&m);    for (i=1;i<=m;i++)    {        scanf("%d%d",&u,&v);        ne[i]=fir[u];        fir[u]=i;        to[i]=v;    }}void dfs(int u){    int i,v,x;    dfn[u]=low[u]=++clo;    ins[u]=1;    sta[++top]=u;    for (i=fir[u];i;i=ne[i])    {        if (!dfn[v=to[i]])        {            dfs(v);            low[u]=min(low[u],low[v]);        }        else        {            if (ins[v]) low[u]=min(low[u],dfn[v]);        }    }    if (low[u]==dfn[u])    {        tot++;        do        {            x=sta[top--];            ins[x]=0;            bel[x]=tot;            if (b[x]) minv[tot]=min(minv[tot],val[x]);            else minn[tot]=min(minn[tot],x);        }        while (x!=u);    }}void tarjan(){    int i;    memset(minv,0x3f,sizeof(minv));    memset(minn,0x3f,sizeof(minn));    for (i=1;i<=n;i++)      if (!dfn[i])        dfs(i);}void dfs2(int u){    int i;    vis[u]=1;    for (i=0;i<to2[u].size();i++)      dfs2(to2[u][i]);}bool solve(int &ret){    int i,u,v;    bool flag;    memset(need,1,sizeof(need));    ret=0;    for (u=1;u<=n;u++)      for (i=fir[u];i;i=ne[i])        if (bel[v=to[i]]!=bel[u]&&!map[bel[u]][bel[v]])        {            map[bel[u]][bel[v]]=1;            to2[bel[u]].push_back(bel[v]);            need[bel[v]]=0;        }    for (i=1;i<=tot;i++)      if (need[i]&&minv[i]<oo)      {        ret+=minv[i];        dfs2(i);      }    flag=1;    for (i=1;i<=tot;i++)      if (!vis[i])      {        if (flag) ret=oo,flag=0;        ret=min(ret,minn[i]);      }    return flag;}int main(){    int x;    init();    tarjan();    if (solve(x)) printf("YES\n");    else printf("NO\n");    printf("%d\n",x);}