【每日一题】洛谷P4819 [中山市选]杀人游戏

eastcloud6

已于 2022-12-27 18:53:42 修改

阅读量129

点赞数

分类专栏：算法每日一题文章标签：算法 c++ 图论

于 2022-12-27 18:43:58 首次发布

本文链接：https://blog.csdn.net/eastcloud6/article/details/128459983

版权

算法每日一题专栏收录该内容

1 篇文章 0 订阅

订阅专栏

[中山市选]杀人游戏

题目描述

一位冷血的杀手潜入Na-wiat，并假装成平民。警察希望能在 $N$ 个人里面，查出谁是杀手。警察能够对每一个人进行查证，假如查证的对象是平民，他会告诉警察，他认识的人，谁是杀手，谁是平民。假如查证的对象是杀手，杀手将会把警察干掉。现在警察掌握了每一个人认识谁。每一个人都有可能是杀手，可看作他们是杀手的概率是相同的。

问：根据最优的情况，保证警察自身安全并知道谁是杀手的概率最大是多少？

输入格式

第一行有两个整数 $N, M$ 。
接下来有 $M$ 行，每行两个整数 $x, y$ ，表示 $x$ 认识 $y$ （ $y$ 不一定认识 $x$ ,例如President同志）。

注：原文zz敏感内容已替换

输出格式

仅包含一行一个实数，保留小数点后面 $6$ 位，表示最大概率。

样例 #1

样例输入 #1

样例输出 #1

0.800000

提示

题解

警察只需要查证 $1$ 。假如 $1$ 是杀手，警察就会被杀。假如 $1$ 不是杀手，他会告诉警察 $2, 3, 4, 5$ 谁是杀手。而 $1$ 是杀手的概率是 $0.2$ ，所以能知道谁是杀手但没被杀的概率是 $0.8$ 。

对于 $100\%$ 的数据有 $1 \leq N \leq 100000, 0 \leq M \leq 300000$ 。

题意简述

给出一个有向图，表示每个人的认识关系，所有人中有一个杀手，询问正常人时他会说出自己认识的人的身份，询问杀手则会被杀死，问活下来并知道杀手的最大概率是多少。

题目分析

在此题中，由于每个人是杀手的概率均等，假设抽查了 $p$ 个人，活下来的概率则为 $\frac{p}{n}$ ，若要让该概率最大，我们可以转化为求最小的 $p$ 。

由于每个人的认识关系为一张无向图，因此只要对一个人发出询问，就能知道他所连边的人的身份，进而推出他能到达的人的身份，因此我们可以想到一个贪心策略：不选有入度点，否则选取能到达他的点一定更优。

但图中可能有环的存在，加上图不连通，就可能出现必须选一个有入度点的情况，于是我们考虑对图进行操作，注意到如果两点在一个强连通分量中，那么选取这两个点是等效的，因此我们可以将图缩点，进而统计无入度点的个数，因为此处将环也缩成了一个节点，因此不会有上述情况且更方便操作。

但我们还需要注意一种特殊情况，当一个强连通分量的大小为 1 且能直接到达的点都起码有两个入度时，这个点不用选，因为其他点检查过后只剩他一个没有检查，可以直接推测出这个点的身份，但是因为我们使用了推测法，这样的点只能存在一个。

总上，我们建图后求 scc，统计无入度点和特殊情况即可。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<map>
using namespace std;
map<pair<int,int> ,int > t;
struct Node{
	int u,v;
}edg[300001];
vector<int> l[300001];
vector<int> tr[300001];
int dfn[300001],low[300001];
int c[300001],q[300001],ins[300001];
int in[300001],ou[300001];
int siz[300001];
int tail,tot,cnt;
void tarjan(int x){
	dfn[x]=low[x]=++tot;
	q[++tail]=x;
	ins[x]=1;
	for(int i=0;i<l[x].size();i++){
		int v=l[x][i];
		if(!dfn[v]){
			tarjan(v);
			low[x]=min(low[x],low[v]);
		}
		else if(ins[v])low[x]=min(low[x],dfn[v]);
	}
	if(low[x]==dfn[x]){
		int y;
		cnt++;
		do{
			y=q[tail--];
			c[y]=cnt;
			ins[y]--;
			siz[cnt]++;
		}while(y!=x);
	}
}
int main(){
	int n,m,x,y;
	cin>>n>>m;
	for(int i=1;i<=m;i++){
		cin>>x>>y;
		l[x].push_back(y);
		edg[i].u=x;edg[i].v=y;
	}
	for(int i=1;i<=n;i++) if(!dfn[i])tarjan(i);
	for(int i=1;i<=m;i++){
		if(c[edg[i].u]==c[edg[i].v])continue;
		pair<int,int> tmp;
		tmp.first=c[edg[i].u];tmp.second=c[edg[i].v];
		if(t[tmp])continue;
		else t[tmp]=1;
		in[c[edg[i].v]]++;
		ou[c[edg[i].u]]++;
		tr[c[edg[i].u]].push_back(c[edg[i].v]);
	}
	int ans=0;
	for(int i=1;i<=cnt;i++) if(in[i]==0) ans++;
	for(int i=1;i<=cnt;i++){
		if(in[i] || siz[i]>1) continue;
		int flag=0;
		for(int j=0;j<tr[i].size();j++){
			if(in[tr[i][j]]<2){
				flag=1;
				break;
			}
		}
		if (!flag){
			ans--;
			break;
		}
	}
	double p=(double)(n-ans)*1.0/n;
	printf("%.6lf",p);
}