jzoj2473. 【NOI2011中山联考7.18】杀人游戏（强连通分量）

2473. 【NOI2011中山联考7.18】杀人游戏

Description
　　一位冷血的杀手潜入Na-wiat，并假装成平民。警察希望能在N个人里面，查出谁是杀手。
　　警察能够对每一个人进行查证，假如查证的对象是平民，他会告诉警察，他认识的人，谁是杀手，谁是平民。假如查证的对象是杀手，杀手将会把警察干掉。
　　现在警察掌握了每一个人认识谁。
　　每一个人都有可能是杀手，可看作他们是杀手的概率是相同的。
　　问：根据最优的情况，保证警察自身安全并知道谁是杀手的概率最大是多少？

Input

输入文件killer.in，第一行有两个整数N,M。
　　接下来有M行，每行两个整数x,y，表示x认识y（y不一定认识x,例如胡锦涛同志）。

Output
　　输出文件killer.out仅包含一行一个实数，保留小数点后面6位，表示最大概率。

Sample Input
5 4
1 2
1 3
1 4
1 5

Sample Output
0.800000

Hint
【样例解释】
　　警察只需要查证1。假如1是杀手，警察就会被杀。假如1不是杀手，他会告诉警察2,3,4,5谁是杀手。而1是杀手的概率是0.2,所以能知道谁是杀手但没被杀的概率是0.8。
【数据规模】
　　对于30%的数据有1≤N ≤ 10,0≤M ≤10
　　对于100%的数据有1≤N ≤ 10 0000,0≤M ≤ 30 0000

分析：显然，这是一个有向有环图，对于每个环只要知道其中一个人的身份即可知道所有人的身份，也就是说只需要问其中一个人，因此我们可以把环缩成一个点，得到一张有向无环图，那么容易发现，只要把所有入度为0的点询问一次就可以知道所有点的身份，统计一下入度为0的点的个数ans，1-ans/n即为答案。

代码

#include <cstdio> 
#include <cstring>
#include <string>
#include <stack>
#define N 500000
using namespace std;

struct arr
{
	int to,nxt;
}a[N],c[N];
int n,m,l,ls[N],lc,lsc[N],d[N],size;
int dfn[N],low[N],b[N],siz[N],tot,cnt;
bool v[N],vis[N],fl;
stack<int> s;

double fmax(double x, double y){return x>y?x:y;}
int fmin(int x, int y){return x<y?x:y;}


int read()
{
	int x = 0, f = 1;
	char ch = getchar();
	while (ch < '0' || ch > '9') {if (ch == '-') f = -1; ch = getchar();}
	while (ch >= '0' && ch <= '9') x = x * 10 + ch - '0', ch = getchar();
	return x * f;
}

void add(int x, int y)
{
	a[++l].to = y;
	a[l].nxt = ls[x];
	ls[x] = l;
}
void addc(int x, int y)
{
	c[++lc].to = y;
	c[lc].nxt = lsc[x];
	lsc[x] = lc;
}

void tarjan(int x)
{
	dfn[x] = low[x] = ++tot;
	s.push(x);
	v[x] = true;
	for (int i = ls[x]; i; i = a[i].nxt)
		if (!dfn[a[i].to])
		{
			tarjan(a[i].to);
			low[x] = fmin(low[x], low[a[i].to]);
		}
		else if (v[a[i].to]) low[x] = fmin(low[x], dfn[a[i].to]);
	if (low[x] == dfn[x])
	{
		++cnt;
		int j;
		do
		{
			j = s.top();
			s.pop();
			siz[cnt]++;
			b[j] = cnt;
			v[j] = false;
		}while (x != j);
	}
}

void dfs(int x)
{
	vis[x] = true;
	size += siz[x];
	for (int i = lsc[x]; i; i = c[i].nxt)
		if (!vis[c[i].to]) dfs(c[i].to);
}

int main()
{
	freopen("investigation.in","r",stdin);
	freopen("investigation.out","w",stdout);
	scanf("%d%d", &n, &m);
	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++)
		for (int j = ls[i]; j; j = a[j].nxt)
			if (b[i] != b[a[j].to]) addc(b[i],b[a[j].to]),d[b[a[j].to]]++;
	int ans = 0;
	for (int i = 1; i <= cnt; i++)
		if (!d[i])
		{
			if (siz[i] == 1 && !fl) 
			{
				for (int j = lsc[i]; j; j = c[j].nxt)
					if (d[c[j].to] != 1) fl = true;
			}
			ans++;
		}
	if (fl || ans == n) ans--;
	double xx = 1.0 * ans / (1.0 * n);
	printf("%.6lf", 1.0 - xx);
	fclose(stdin);
	fclose(stdout);
}