洛谷1262间谍网络（tarjan缩点）

tarjan缩点思想：http://m.blog.csdn.net/blog/qq574857122/16361033

题目描述

由于外国间谍的大量渗入，国家安全正处于高度的危机之中。如果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：

【样例1】
3
2
1 10
2 100
2
1 3
2 3
【样例2】
4
2
1 100
4 200
2
1 2
3 4

输出样例#1：

【样例1】
YES
110
【样例2】
NO
3

思想：对于NO的情况，我们先做一遍bfs如果不能全部访问到，输出NO和最小的不能到达的结点。否则，tarjan缩点，将缩完后的点中入度为0的结点值加入ans。

贴代码：

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
const int maxn=3003;
int n,p,m,to[maxn*maxn],nxt[maxn*maxn],fst[maxn],cnt=0,cst[maxn],scc_cnt=0,dfs_clock=0;
int dfn[maxn],low[maxn],ans=0,mon[maxn],du[maxn],nod[maxn];
bool vst[maxn],ins[maxn];
inline int get(){
    char c;while(!isdigit(c=getchar()));
    int v=c-48;while(isdigit(c=getchar()))v=v*10+c-48;
    return v;
}
inline void add(int x,int y){to[++cnt]=y;nxt[cnt]=fst[x];fst[x]=cnt;}
stack<int>s;
inline void tarjan(int x){
	s.push(x);ins[x]=1;
    dfn[x]=low[x]=++dfs_clock;
    for(int i=fst[x];i;i=nxt[i]){
	    if(!dfn[to[i]]){
		    tarjan(to[i]);
		    low[x]=min(low[x],low[to[i]]);
		}
		else if(ins[to[i]])low[x]=min(low[x],low[to[i]]);
	}
	if(dfn[x]==low[x]){
		++scc_cnt;
	    for(;;){
		    int u=s.top();s.pop();
		    nod[u]=scc_cnt;
		    ins[u]=0;
		    if(cst[u])mon[scc_cnt]=min(mon[scc_cnt],cst[u]);
		    if(u==x)break;
		}
	}
}
inline void find_scc(){
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(mon,0x7f,sizeof(mon));
    memset(du,0,sizeof(du));
    memset(ins,0,sizeof(ins));
    for(int i=1;i<=n;++i)if(!dfn[i])tarjan(i);
    for(int i=1;i<=n;++i)
        for(int j=fst[i];j;j=nxt[j])if(nod[i]!=nod[to[j]])++du[nod[to[j]]];
    for(int i=1;i<=scc_cnt;++i)if(!du[i])ans+=mon[i];
    printf("YES\n%d\n",ans);
}
queue<int>q;
inline void bfs(int u){
	vst[u]=1;
    q.push(u);
    bool use[maxn];
    memset(use,0,sizeof(use));
    use[u]=1;
    while(!q.empty()){
		int x=q.front();q.pop();
        for(int i=fst[x];i;i=nxt[i]){
	        int y=to[i];
	        if(!vst[y]){
				vst[y]=1;
	            if(!use[y]){
			        q.push(y);
			        use[y]=1;
			    }
			}
	    }
	    use[x]=0;
    }
}
inline bool check(){
    for(int i=1;i<=n;++i){
	    if(cst[i] && !vst[i])bfs(i);
	}
	for(int i=1;i<=n;++i)if(!vst[i]){
	    printf("NO\n%d\n",i);
	    return 0;
	}
	return 1;
}
int main(){
	int u;
    n=get();p=get();memset(vst,0,sizeof(vst));
    memset(cst,0,sizeof(cst));
    for(int i=1;i<=p;++i)u=get(),cst[u]=get();
    m=get();
    for(int i=1;i<=m;++i){
	    int x=get(),y=get();
	    add(x,y);
	}
	if(check())find_scc();
	return 0;
}