BZOJ2438[中山市选2011] 杀人游戏

原题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=2438

杀人游戏

Description

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

Input

第一行有两个整数 N,M。

接下来有 M 行，每行两个整数 x,y，表示 x 认识 y（y 不一定认识 x,例如某个伟人） 。

Output

仅包含一行一个实数，保留小数点后面 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。

对于 100%的数据有 1≤N ≤ 10 0000,0≤M ≤ 30 0000

数据已加强！

题解

对于一个强连通分量，只需要问一个人就能掌握整个$\mathcal{S}\mathcal{C}\mathcal{C}$的身份，而通过一个$\mathcal{S}\mathcal{C}\mathcal{C}$又能到达其他$\mathcal{S}\mathcal{C}\mathcal{C}$，所以我们实际需要询问的就是入度为$0$的$\mathcal{S}\mathcal{C}\mathcal{C}$，所以我们缩一波点，统计一下入度为$0$的$\mathcal{S}\mathcal{C}\mathcal{C}$的个数即可。

呵呵，血$\mathcal{W}\mathcal{A}$。

当存在一个大小为$1$，入度为$0$的$\mathcal{S}\mathcal{C}\mathcal{C}$，它没有出边或者通向的$\mathcal{S}\mathcal{C}\mathcal{C}$可以被其他点到达（即入度$\ge 2$），我们在询问完其他点以后就不必询问这个点，此时询问数$-1$.

这下是真的$\mathcal{A}$了。

代码

#include<bits/stdc++.h>
using namespace std;
const int M=1e5+5;
int n,m,tot,cot,dfn[M],low[M],sta[M],siz[M],top,col[M],ru[M],ans,f,t;
bool vis[M];
vector<int>mmp[M],scc[M];
void in()
{
	int x,y;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;++i)scanf("%d%d",&x,&y),mmp[x].push_back(y);
}
void dfs(int v)
{
	int to;dfn[v]=low[v]=++tot;vis[v]=1;sta[++top]=v;
	for(int i=mmp[v].size()-1;i>=0;--i)
	{
		to=mmp[v][i];
		if(!dfn[to])dfs(to),low[v]=min(low[v],low[to]);
		else if(vis[to])low[v]=min(low[v],dfn[to]);
	}
	if(dfn[v]==low[v])
	{
		col[v]=++cot,vis[v]=0;
		for(;sta[top]!=v;col[sta[top]]=cot,vis[sta[top--]]=0);
		--top;
	}
}
void ac()
{
	for(int i=1;i<=n;++i)if(!dfn[i])dfs(i);
	for(int i=1;i<=n;++i)++siz[col[i]],scc[col[i]].push_back(i);
	for(int i=1;i<=cot;++i)
	{
		vis[i]=1;
		for(int j=scc[i].size()-1;j>=0;--j)
		{
			f=scc[i][j];
			for(int k=mmp[f].size()-1;k>=0;--k)
			{
				t=mmp[f][k];if(vis[col[t]])continue;
				++ru[col[t]],vis[col[t]]=1;
			}
		}
		for(int j=scc[i].size()-1;j>=0;--j)for(int k=mmp[scc[i][j]].size()-1;k>=0;--k)vis[col[mmp[scc[i][j]][k]]]=0;
		vis[i]=0;
	}
	for(int i=1;i<=cot;++i)if(!ru[i])++ans;
	bool flag;
	for(int i=1;i<=cot;++i)
	{
		if(siz[i]>1||ru[i])continue;
		flag=1;
		for(int j=mmp[scc[i][0]].size()-1;j>=0;--j)
		if(ru[col[mmp[scc[i][0]][j]]]<2){flag=0;break;}
		if(flag){--ans;break;}
	}
	printf("%.6lf",1.0*(n-ans)/n);
}
int main(){in();ac();}