洛谷 P1262 间谍网络

题目描述
由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果A间谍手中掌握着关于B间谍的犯罪证据，则称A可以揭发B。有些间谍收受贿赂，只要给他们一定数量的美元，他们就愿意交出手中掌握的全部情报。所以，如果我们能够收买一些间谍的话，我们就可能控制间谍网中的每一分子。因为一旦我们逮捕了一个间谍，他手中掌握的情报都将归我们所有，这样就有可能逮捕新的间谍，掌握新的情报。
我们的反间谍机关提供了一份资料，色括所有已知的受贿的间谍，以及他们愿意收受的具体数额。同时我们还知道哪些间谍手中具体掌握了哪些间谍的资料。假设总共有n个间谍(n不超过3000)，每个间谍分别用1到3000的整数来标识。
请根据这份资料，判断我们是否有可能控制全部的间谍，如果可以，求出我们所需要支付的最少资金。否则，输出不能被控制的一个间谍。

---

【题目分析】
Tarjan。

---

【代码】

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int n,p,x,y,m;
int a[3001];
int h[3001],to[8001],ne[8001],fr[8001],en=0;
int bel[3001],cnt[3001],idx=0,low[3001],dfn[3001],minn[3001],vcnt=0,ins[3001],sta[3001],top=0;
int minnum[3001];
void add(int a,int b)
{
    fr[en]=a;
    to[en]=b;
    ne[en]=h[a];
    h[a]=en++;
}
void tarjan(int k)
{
    low[k]=dfn[k]=++idx;
    ins[k]=1;
    sta[++top]=k;
    for (int i=h[k];i>=0;i=ne[i])
    {
        if (!dfn[to[i]])
        {
            tarjan(to[i]);
            low[k]=min(low[k],low[to[i]]);
        }
        else if (ins[to[i]]) low[k]=min(low[k],dfn[to[i]]);
    }
    if (low[k]==dfn[k])
    {
        int x=-1;
        ++vcnt;
        while (x!=k)
        {
            x=sta[top--];
            bel[x]=vcnt;
            cnt[vcnt]++;
            if (a[x]==0) minnum[vcnt]=min(x,minnum[vcnt]);
            ins[x]=0;
            if (a[x]!=0) minn[vcnt]=min(a[x],minn[vcnt]);
        }
    }
}
int h2[3001],to2[8001],ne2[8001],fr2[8001],en2=0;
int du[3001];
void add2(int a,int b)
{
//    printf("add2 %d %d\n",a,b);
    to2[en2]=b;
    ne2[en2]=h2[a];
    fr2[en2]=a;
    h2[en2]=en2++;
    du[b]++;
}
queue <int> q;
void prt()
{
    int vis[3001]; 
    memset(vis,0,sizeof vis);
//    printf("its no\n");
    int ans=0x3f3f3f3f;
    while (!q.empty())
    {
        int x=q.front(); q.pop(); ans=min(ans,minnum[x]);
//        printf("%d \n",minnum[x]);
        vis[x]=1;
        for (int i=h2[x];i>=0;i=ne2[i]) if (!vis[to2[i]]) q.push(to2[i]);
    }
    printf("NO\n%d\n",ans);
}
int main()
{
    while (!q.empty()) q.pop();
    memset(minnum,0x3f,sizeof minnum);
    memset(h2,-1,sizeof h2);
    memset(minn,0x3f,sizeof minn);
    memset(h,-1,sizeof h);
    scanf("%d%d",&n,&p);
    for (int i=1;i<=p;++i)
    {
        scanf("%d%d",&x,&y);
        a[x]=y;
    }
    scanf("%d",&m);
    for (int i=1;i<=m;++i)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        add(a,b);
    }
    for (int i=1;i<=n;++i) if (!dfn[i]) tarjan(i);
//    for (int i=1;i<=n;++i) printf("%d ",bel[i]); printf("\n");
//    for (int i=1;i<=n;++i) printf("%d ",minnum[i]); printf("\n");
    for (int i=0;i<en;++i)
        if (bel[fr[i]]!=bel[to[i]])
            add2(bel[fr[i]],bel[to[i]]);
    int ans=0,flag=1;
    for (int i=1;i<=vcnt;++i)
        if (!du[i]){
            if (minn[i]==0x3f3f3f3f)
            {
                q.push(i);
                flag=0;
//                prtno();
//                return 0;
            }
            else ans+=minn[i];
        }
    if (flag) printf("YES\n%d\n",ans);
    else prt();
}