[bzoj 2438--中山市选2011]杀人游戏

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

这道题是一道强联通，是一道比较好的题目，与经典的题目有一点的小改动。首先，设gr为最少查证的人数，那么概率就为(n-gr)/n。之后，这题有个特殊情况，如果有一个点，它的入度为0，并自为一个联通分量，并与它相连的点都可以可以被其他给覆盖的话，那么它就无用，gr-1。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
struct node
{
    int x,y,next;
}a[310000];int len,last[310000];
void ins(int x,int y)
{
    len++;
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
int id,cnt,top;
int low[110000],dfn[110000],sta[110000],belong[110000],s[110000],ss[110000];
int ru[110000];
bool v[110000];
void dfs(int x)
{
    dfn[x]=low[x]=++id;
    sta[++top]=x;v[x]=true;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(dfn[y]==-1)
        {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else 
        {
            if(v[y]==true)low[x]=min(low[x],dfn[y]);
        }
    }
    if(low[x]==dfn[x])
    {
        cnt++;int i;ss[cnt]=x;
        do
        {
            i=sta[top--];
            v[i]=false;
            s[cnt]++;
            belong[i]=cnt;
        }while(i!=x);
    }
}       
bool he(int x)
{
    x=ss[x];
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(ru[belong[y]]==1)return false;
    }
    return true;
}
int main()
{
    len=0;memset(last,0,sizeof(last));
    int n,m,ans=0;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        ins(x,y);
    }
    id=cnt=top=0;
    memset(dfn,-1,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(belong,0,sizeof(belong));
    memset(sta,0,sizeof(sta));
    memset(v,false,sizeof(v));
    for(int i=1;i<=n;i++)
    {
        if(dfn[i]==-1)dfs(i);
    }
    for(int i=1;i<=len;i++)
    {
        if(belong[a[i].x]!=belong[a[i].y])
        {
            ru[belong[a[i].y]]++;
        }
    }
    int gr=0;
    for(int i=1;i<=cnt;i++)
    {
        if(ru[i]==0)gr++;
    }
    for(int i=1;i<=cnt;i++)
    {
        if(s[i]==1 && ru[i]==0 && he(i)==true)
        {
            gr--;break;
        }
    }
    printf("%.6lf\n",(n-gr)*1.0/n*1.0);
    return 0;
}