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);
}