【洛谷P1262】 间谍网络

题目描述

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

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

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

一个环内的点当做一个点，进行有向图缩点，缩点之后的点权就是该点包括的点的点权最小值。

之后重新建图，对于每个点如果他的点权不是无限大，那么就可以从他开始扩展，dfn就可以解决。不过有一点问题，就是如果可以遍历所有点，那么我们是需要输出最小代价的。

那么尝试hack一下现在的做法。

可以发现，确实会有地方多计算了点权。

比如下面这个图。

按照当前做法，1号点的点权10也会被计算到答案里，但是我们只需要5号点的点权20就可以了。

所以就有了一个优化，就是统计重新建图之后每个点的入度，然后从入度为零的点开始遍历。

之后再去从入度不为零的点遍历。

code：

#include <iostream>
#include <cstdio>
#include <cstring>using namespace std;const int wx=50017;inline int read(){int sum=0,f=1; char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}while(ch>='0'&&ch<='9'){sum=(sum<<1)+(sum<<3)+ch-'0'; ch=getchar();}return sum*f;
}int tot,st[wx],top,col,n,p,m;
int head[wx],h[wx],num,Num;
int dfn[wx],low[wx],belong[wx],size[wx],a[wx],v[wx],vis[wx];
int in[wx];struct e{int nxt,to;
}edge[wx*2];void add(int from,int to){edge[++num].nxt=head[from];edge[num].to=to;head[from]=num;
}void Tarjan(int u){dfn[u]=low[u]=++tot;st[++top]=u;for(int i=head[u];i;i=edge[i].nxt){int v=edge[i].to;if(!dfn[v]){Tarjan(v);low[u]=min(low[u],low[v]);}else if(!belong[v]){low[u]=min(low[u],dfn[v]);}}if(low[u]==dfn[u]){belong[u]=++col;size[col]++;while(st[top]!=u){belong[st[top]]=col;size[col]++;top--;}top--;}
}struct node{int nxt,to;
}e[wx*2];void Add(int from,int to){e[++Num].nxt=h[from];e[Num].to=to;h[from]=Num;
}void dfs(int u){vis[u]=1;for(int i=h[u];i;i=e[i].nxt){int v=e[i].to;if(vis[v])continue;dfs(v);}
}int main(){n=read(); p=read();memset(a,0x3f,sizeof a);memset(v,0x3f,sizeof v);for(int i=1;i<=p;i++){int x; x=read(); a[x]=read();}m=read();for(int i=1;i<=m;i++){int x,y;x=read(); y=read();add(x,y);}for(int i=1;i<=n;i++)if(!dfn[i])Tarjan(i);for(int i=1;i<=n;i++)v[belong[i]]=min(v[belong[i]],a[i]);for(int u=1;u<=n;u++){for(int i=head[u];i;i=edge[i].nxt){int v=edge[i].to;if(belong[v]!=belong[u])in[belong[v]]++,Add(belong[u],belong[v]);}}int ans=0;for(int i=1;i<=col;i++){if(v[i]!=0x3f3f3f3f&&vis[i]==0&&in[i]==0){if(i==1)vis[i]=1;ans+=v[i];dfs(i);}}for(int i=1;i<=col;i++){if(v[i]!=0x3f3f3f3f&&vis[i]==0){if(i==1)vis[i]=1;ans+=v[i];dfs(i);}}for(int i=1;i<=n;i++){if(!vis[belong[i]]){puts("NO");printf("%d\n",i); return 0;}}puts("YES");printf("%d\n",ans);return 0;
}