首页 > 技术文章 > 1517:间谍网络

lyc-lb-blogs 2021-07-05 16:04 原文

1517:间谍网络


时间限制: 1000 ms         内存限制: 562144 KB
提交数: 668     通过数: 333

【题目描述】

由于外国间谍的大量渗入,国家安全正处于高度危机之中。如果 $A$ 间谍手中掌握着关于 $B$ 间谍的犯罪证据,则称 $A$ 可以揭发 $B$。有些间谍接受贿赂,只要给他们一定数量的美元,他们就愿意交出手中掌握的全部情报。所以,如果我们能够收买一些间谍的话,我们就可能控制间谍网中的每一分子。因为一旦我们逮捕了一个间谍,他手中掌握的情报都将归我们所有,这样就有可能逮捕新的间谍,掌握新的情报。

我们的反间谍机关提供了一份资料,包括所有已知的受贿的间谍,以及他们愿意收受的具体数额。同时我们还知道哪些间谍手中具体掌握了哪些间谍的资料。假设总共有 $n$ 个间谍,每个间谍分别用 $1$ 到 $n$ 的整数来标识。

请根据这份资料,判断我们是否可能控制全部的间谍,如果可以,求出我们所需要支付的最少资金。否则,输出不能被控制的一个间谍。

【输入】

第一行只有一个整数 $n$。第二行是整数 $p$。表示愿意被收买的人数。

接下来的 $p$ 行,每行有两个整数,第一个数是一个愿意被收买的间谍的编号,第二个数表示他将会被收买的数额。

紧跟着一行只有一个整数 $r$。然后 $r$ 行,每行两个正整数,表示数对 $(A,B)$,$A$ 间谍掌握 $B$ 间谍的证据。

【输出】

如果可以控制所有间谍,第一行输出$YES$,并在第二行输出所需要支付的贿金最小值。否则输出$NO$,并在第二行输出不能控制的间谍中,编号最小的间谍编号。

【输入样例】

2 
1 
2 512 
2 
1 2 
2 1

【输出样例】

YES
512

【提示】

$1≤n≤3000,1≤p≤n,1≤r≤8000$, 每个收买的费用为非负数且不超过 $20000$。

缩点。对于入度为0的点就必需要收买。

初始化 cost[i] = w[i] = 0x3f3f3f3f,表示控制第i个联通分量 和 收买 第 i 个 间谍的花费。

从 1 枚举到 n ,如果能收买该间谍就用 tarjan算法 求强连通分量

最后统计所有点,如果存在孤立点不被访问到,则无解

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int N = 1e5 + 10;
int h[N], e[N], ne[N], idx;
int hs[N];
int w[N];
int dfn[N],low[N],timestamp;
int stk[N],in_stk[N],top;
int id[N],scc_cnt;
int cost[N];
int n, p, r;
int din[N];
void add(int a,int b)
{
    e[idx] = b; ne[idx] = h[a]; h[a] = idx++;
}

void add2(int a,int b,int c)
{
    e[idx] = b; w[idx] = c; ne[idx] = hs[a]; hs[a] = idx++;
}

void tarjan(int u)
{
    dfn[u] = low[u] = ++timestamp;
    stk[++top] = u; in_stk[u] = 1;
    for(int i = h[u]; ~i ; i = ne[i])
    {
        int j = e[i];
        if(!dfn[j])
        {
            tarjan(j);
            low[u] = min(low[u], low[j]);
        }
        else if(in_stk[j])
        {
            low[u] = min(low[u], dfn[j]);
        }
    }

    if(dfn[u] == low[u])
    {
        ++scc_cnt;
        int y;
        do
        {
           y = stk[top--];
           in_stk[y] = 0;
           id[y] = scc_cnt;
           cost[scc_cnt] = min(cost[scc_cnt], w[u]);
        } while (y != u);
    }

}

int main()
{
    scanf("%d %d",&n,&p);
    memset(h, -1, sizeof(h));
    memset(w, 0x3f,sizeof(w));
    memset(cost,0x3f,sizeof(cost));
    for(int i = 1; i <= p; i++)
    {
        int x;
        scanf("%d",&x);
        scanf("%d",&w[x]);
    }
    scanf("%d",&r);
    for(int i = 1; i <= r; i++)
    {
        int a, b;
        scanf("%d %d",&a,&b);
        add(a, b);
    }
    for(int i = 1; i <= n; i++)
    {
        if(!dfn[i] && w[i] < 0x3f3f3f3f) tarjan(i);
    }
    for(int i = 1 ;i <= n; i++)
    {
        if(!dfn[i])
        {
            printf("NO\n%d",i);
            return 0;
        }
    }
    for(int i = 1; i <= n; i++)
    {
        for(int j = h[i]; ~j; j = ne[j])
        {
            int k = e[j];
            int a = id[i] , b = id[k];
            if(a != b)
            {
                din[b]++;
            }
        }
    }
    int ans = 0;
    for(int i = 1; i <= scc_cnt; i++)
    {
        if(!din[i]) ans += cost[i];
    }
    printf("YES\n%d",ans);
    return 0;
}

推荐阅读