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

题目传送门

题面描述：

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

Input
第一行有两个整数 N,M。
接下来有 M 行，每行两个整数 x,y，表示 x 认识 y（y 不一定认识 x）。

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

题意转化：有n个点，每个点都有一个身份。其中有 $m$ 个关系，每个关系给出 $a$ 和 $b$ ，表示 $a$ 认识 $b$ （ $b$ 不一定认识 $a$ ），求最少要询问多少个人，才能知道所有人的身份（把那个杀手找出来）。

分析：

我们要询问最少的人有什么用呢？
假设有 $n$ 个人，我们询问了 $k$ 个点知道了所有人的身份。但是，杀手有可能在这k个人之中，也就是说，我们在询问这 $k$ 个人的时候就被杀掉了。对于每个点，他是杀手的概率是 $\frac{1}{n}$ ,我们询问了 $k$ 个人,概率就是 $\frac{k}{n}$ 。所以，我们被杀的概率是 $\frac{k}{n}$ ,能够存活的概率就是 $1-\frac{k}{n}$ 。

题意理解了之后，我们再思考做法。

考虑样例：
在这里插入图片描述
我们只要访问 $1$ 号节点，就能知道全部的身份，1号节点是杀手的概率是 $\frac{1}{5}$ ，所以安全的概率就是 $\frac{4}{5}$ 。

那么样例太水，尝试来一个复杂一点的样例。
在这里插入图片描述
对于这个样例，我们只要访问1号节点知道2,3,4号点的身份和访问6号点知道5号点的身份就可以了。概率是 $\frac{2}{3}$ 。

现在我们已经大致理清楚了题意，但是题目上的图不免有些复杂了，我们能否尝试将图变得简单一点呢？
对于一个集合而言，如果访问了集合内的一个点就知道了集合内所有人的身份，那么这个集合中的点都是等价的，相当于一个点，那么我们只要把他们当成一个点即可。又联想到有向图，我们想到了什么？强连通分量!

没错，这道题就是一道强联通分量缩点的例题。

在这里插入图片描述

图就变得简洁多了，我们只要访问1号点和六号点即可。

第一个猜想：缩点后入度为0的点的个数就是所求的最小点数
恭喜你，那么整整21分！

来组反例：
在这里插入图片描述

那么对于这幅图，我们实际上只要访问一号点知道1,2,3,4号点的身份，通过排除法就可以得到5号点的身份。即我们只需要访问一个人就可以知道所有人的身份。

我们在猜想：对于缩点后的一个点，如果他的大小为1，入读为0，出度>1，那么这个点就可以不用访问，不过我们最多只能不访问这么一个节点（用排除法的性质就可以想）

我们在屡一下思路：

强联通分量缩点
找出入度为0的点数k
找出是否存在上述的特属于的点，记为布尔f
答案为 $1-\frac{k-f}{n}$

那么具体代码如下：

#include<bits/stdc++.h>
using namespace std;
struct node{
    int x,y,Next; 
}e[600010],e1[600010];
int linkk[600010];
int in[600010];
int linkk1[600010];
int dfn[600010];
int low[600010];
int stackk[600010];
int fam[600010];
int size[600010];
bool vis[600010];
int len=0,cnt=0,num=0,top=0,ans=0;
int n,m;
void insert(int x,int y){
    e[++len].Next=linkk[x];
    linkk[x]=len;
    e[len].y=y;
    e[len].x=x;
}
void insert2(int x,int y){
    e1[++len].Next=linkk1[x];
    linkk1[x]=len;
    e1[len].y=y;
}
void insert_new(){
    for (int i=1;i<=m;i++)
      if (fam[e[i].x]!=fam[e[i].y])
        insert2(fam[e[i].x],fam[e[i].y]),in[fam[e[i].y]]++;
}
void tarjan(int x){
    dfn[x]=low[x]=++cnt;
    stackk[++top]=x;vis[x]=1;
    for (int i=linkk[x];i;i=e[i].Next)
      if (!dfn[e[i].y])
        tarjan(e[i].y),low[x]=min(low[x],low[e[i].y]);
      else if (vis[e[i].y]) low[x]=min(low[x],dfn[e[i].y]);
    if (dfn[x]==low[x]){
	    ++num;
	    int t;
	    do{
		    t=stackk[top--];
		    vis[t]=0;
		    size[num]++;
		    fam[t]=num;
		}while (x!=t);
	}
}
int main(){
    scanf("%d %d",&n,&m);
    for (int i=1,x,y;i<=m;i++)
      scanf("%d %d",&x,&y),insert(x,y);
    for (int i=1;i<=n;i++)
      if (!dfn[i]) tarjan(i);
    len=0;
    insert_new();
    bool f=0;
    for (int i=1;i<=num;i++){
    	bool ff=0;
	    if (!f&&!in[i]&&size[i]==1){
	        ff=1;
		    for (int j=linkk1[i];j;j=e1[j].Next) if (in[e1[j].y]==1){ff=0;break;}
	   }
		if (ff) f=1;
		if (!in[i]) ans++;
	}
    printf("%.6lf",1.0-double(double(ans-f)/double(n)));
    fclose(stdin);
    fclose(stdout);
    return 0;
}