bzoj2438 [中山市选2011]杀人游戏(tarjan缩点)

知道谁是杀手相当于知道所有人的身份。因此题目的答案即在有向图中选择最少的点，使得能遍历到至少n-1个点（最后一个点可以推理得到）。设结果为x，则答案为(n-x)/n。
tarjan求强连通缩点，每个强连通我们只需询问一个人即可。而且只询问那些入度为0的即可。特殊情况是有可能存在一个强连通大小为1，入度为0，且它是不必要的（也就是他指向的点的入度都大于1，即使不问他，也可以遍历到那些点），那我们可以不去询问这个人。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
#define M 300010
#define pa pair<int,int>
inline 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;
}
int n,m,h[N],h1[N],num=0,dfn[N],low[N],dfnum=0,bel[N],scc=0,sz[N],in[N];
bool inq[N],f[N];
stack<int>q;
struct edge{
    int to,next;
}data[M],data1[M];
inline void add(int x,int y){
    data1[++num].to=y;data1[num].next=h1[x];h1[x]=num;in[y]++;
}
void tarjan(int x){
    dfn[x]=low[x]=++dfnum;q.push(x);inq[x]=1;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;
        if(!dfn[y]) tarjan(y),low[x]=min(low[x],low[y]);
        else if(inq[y]) low[x]=min(low[x],dfn[y]);
    }
    if(dfn[x]==low[x]){
        ++scc;while(1){
            int y=q.top();q.pop();inq[y]=0;
            bel[y]=scc;sz[scc]++;if(x==y) break;
        }
    }
}
void rebuild(){
    num=0;
    for(int x=1;x<=n;++x){
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;
            if(bel[x]==bel[y]) continue;
            if(!f[bel[y]]) f[bel[y]]=1,add(bel[x],bel[y]);
        }
        for(int i=h1[bel[x]];i;i=data1[i].next) f[data1[i].to]=0;
    }
}
inline bool jud(int x){
    for(int i=h1[x];i;i=data1[i].next)
        if(in[data1[i].to]==1) return 0;
    return 1;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    while(m--){
        int x=read(),y=read();data[++num].to=y;data[num].next=h[x];h[x]=num;
    }
    for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);
    rebuild();int ans=0;//ans--入度为0的点的个数，都得问 
    for(int i=1;i<=scc;++i) if(!in[i]) ans++;
    for(int i=1;i<=scc;++i)
        if(sz[i]==1&&!in[i]&&jud(i)){ans--;break;}
    printf("%.6lf\n",(n-ans)*1.0/n);
    return 0;
}