BZOJ 2438 中山市选2011 杀人游戏

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

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2815  Solved: 808
[Submit][Status][Discuss]

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

数据已加强！ 

Source

容斥一下，最优的的方式就是去询问那些入度为0的点

那种情况是我们按照这种方式完不成任务的呢，就是入度为0的店是杀手

那么我们先对图进行Tarjan缩一下点，有入度的点，显然我们可以不用直接访问，那么我们访问每个入度为0的点

不过这里有个特殊情况，如果存在一个被搁置的点，他最后是不用访问的比如：3个人ABC，A认识B，那么访问A后，A,B和C的身份都能得知

这样就可以少询问一个，但是注意，这种情况的条件是：

入度为0，且只包含1个点，且这个点指向的SCC的入度>=2（缩点前）【并不仅仅是出入度为0】<-特别容易出错

比如：3个人ABC，A认识B，C认识B，那么访问A或C后都可以得到所有人身份；

#include <bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
    int x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int MAXN=1e6+10;
struct node{
    int y,next;
}e[MAXN],E[MAXN];
int linkk[MAXN],Linkk[MAXN],len,stark[MAXN],n,m,cnt[MAXN],sum[MAXN],dfs_clock=0,top,mark[MAXN],ine[MAXN],tot,vis[MAXN],dfn[MAXN],low[MAXN];
inline void insert(int xx,int yy){
    e[++len].y=yy;e[len].next=linkk[xx];linkk[xx]=len;
}
inline void insertt(int xx,int yy){
    E[++len].y=yy;E[len].next=Linkk[xx];Linkk[xx]=len;
}
inline void tarjin(int st){
    dfn[st]=low[st]=++dfs_clock;
    vis[st]=1;stark[++top]=st;
    for(int i=linkk[st];i;i=e[i].next){
        if(!dfn[e[i].y]){
            tarjin(e[i].y);
            low[st]=min(low[e[i].y],low[st]);
        }
        else if(vis[e[i].y]) low[st]=min(low[st],dfn[e[i].y]);
    }
    if(dfn[st]==low[st]){
        int k;tot++;
        do{
            k=stark[top--];
            vis[k]=0;
            ine[k]=tot;
            sum[tot]++;
        }while(k!=st);
    }
}
void init(){
    n=read();m=read();
    for(int i=1;i<=m;i++){
        int xx=read();int yy=read();
        insert(xx,yy);
    }
}
void rebuild(){
    len=0;
    for(int i=1;i<=n;i++){
        for(int j=linkk[i];j;j=e[j].next){
            if(ine[i]!=ine[e[j].y]&&!mark[ine[e[j].y]]){
               cnt[ine[e[j].y]]++;mark[ine[e[j].y]]=1;
               //cout<<ine[i]<<' '<<ine[e[j].y]<<endl;
               insertt(ine[i],ine[e[j].y]);
            }
        }
        for(int j=linkk[i];j;j=e[j].next){
            mark[ine[e[j].y]]=0;
        }
    }
}
inline bool dp(int st){
    if(cnt[st]!=0||sum[st]!=1) return 0;
    for(int i=Linkk[st];i;i=E[i].next){
        if(cnt[E[i].y]==1) return 0;
    }
    return 1;
}
void solve(){
    for(int i=1;i<=n;i++){
        if(!dfn[i]) tarjin(i);
    }
    rebuild();
    int ans=0;
    for(int i=1;i<=tot;i++){
        if(!cnt[i]){
            ans++;
        }
    }
    for(int i=1;i<=tot;i++){
        if(dp(i)){
            ans--;break;
        }
    }
    printf("%lf\n",1.0*(n-ans)/n);
}
int main(){
    init();
    solve();
    return 0;
}