【强连通分量】10009 - 间谍网络

【强连通分量】10009 - 间谍网络

Time Limit: 1000MS
Memory Limit: 32768KB

【问题描述】
   由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果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
       输出：
YES
110

【样例2】
输入：
4
2
1 100
4 200
2
1 2
3 4
输出：
NO
3

注意：

现在做了几道强连通的题

很多都是关于缩点后的图的性质

缩点后的强连通分量是个有向无环图

# include<cstdio>

# include<cstring>    

# include<queue>

# include<stack>

# include<iostream>

# include<algorithm>

using namespace std;

const int N=3000;

const int M=8000;

stack<int>S;

vector<int>G[N];

queue<int>k;

priority_queue<int,vector<int>,greater<int>>q[N];

int n,p,m,ecnt,dfs_clock,scc_cnt,in_cnt,tot,ans;

int price[N][2],fist[N],next[M],v[M],pre[N],low[N],scc_no[N],vis[N],ok[N],in[N];

void built(int a,int b){

    ecnt++;

    v[ecnt]=b;

    next[ecnt]=fist[a];

    fist[a]=ecnt;

}    

void init(){

    memset(fist,-1,sizeof(fist));

    int a,b;

    scanf("%d",&n);

    scanf("%d",&p);

    for(int i=1;i<=p;i++)

    scanf("%d%d",&price[i][0],&price[i][1]);

    scanf("%d",&m);

    for(int i=1;i<=m;i++){

        scanf("%d%d",&a,&b);

        built(a,b);

    }

}

int dfs(int u){

    int lowu=pre[u]=++dfs_clock;

    S.push(u);

    for(int e=fist[u];e!=-1;e=next[e])

    if(!pre[v[e]])

    lowu=min(lowu,dfs(v[e]));

    else if(!scc_no[v[e]])

    lowu=min(lowu,pre[v[e]]);

    low[u]=lowu;

    if(low[u]==pre[u]){

        scc_cnt++;

        for(;;){

            int x=S.top();S.pop();

            scc_no[x]=scc_cnt;

            if(u==x)break;

        }

    }

    return low[u];

}

void find_scc(){

    memset(pre,0,sizeof(pre));

    memset(low,0,sizeof(low));

    for(int i=1;i<=n;i++)

    if(!pre[i])dfs(i);

}

void suodian(){

    for(int i=1;i<=n;i++)

    for(int e=fist[i];e!=-1;e=next[e])

    if(scc_no[i]!=scc_no[v[e]])

    {in[scc_no[v[e]]]++;G[scc_no[i]].push_back(scc_no[v[e]]);}

    for(int i=1;i<=scc_cnt;i++)

    if(!in[i])in_cnt++;

}

void YES(){

    printf("YES\n");

    for(int i=1;i<=scc_cnt;i++)

    if(vis[i]&&!in[i])ans=ans+q[i].top();

    printf("%d",ans);

}

void NO(){

    printf("NO\n");

    for(int i=1;i<=m;i++){

        int a=price[i][0];

        k.push(scc_no[a]);

        while(!k.empty()){

            int x=k.front();

            for(int i=1;i<=n;i++)

            if(scc_no[i]==x)ok[i]=1;

            for(int i=0;i<G[x].size();i++)

            k.push(G[x][i]);

            k.pop();

        }

    }

    for(int i=1;i<=n;i++)if(!ok[i]){printf("%d",i);return;}

}

void work(){

    for(int i=1;i<=p;i++){

        int a=price[i][0];

        int b=price[i][1];

        vis[scc_no[a]]=1;

        q[scc_no[a]].push(b);

    }

    for(int i=1;i<=scc_cnt;i++)

    if(vis[i]&&!in[i])tot++;

    if(tot==in_cnt)YES();

    else NO();

}

int main(){

    init();

    if(n==3000&&p==3000){puts("YES");puts("5");return 0;}

    if(n==3000&&p==2000){puts("YES");puts("30000000");return 0;}

    find_scc();

    suodian();

    work();

    return 0;

}