【BZOJ2438】【中山市选2011】杀人游戏 强连通分量+期望

链接：

#include <stdio.h>
int main()
{
    puts("转载请注明出处[vmurder]谢谢");
    puts("网址：blog.csdn.net/vmurder/article/details/44548895");
}

题解：

SCC：强连通分量。
首先对于一个SCC，你问其中一个点，然后如果它不是杀手，那么你就可以安全地遍历整个SCC——问一个人，是杀手，游戏结束；不是杀手，你可以安全地问他。
然后对于整个图，如果一个SCC有入度，你可以不先问它。
这样我们统计有多少入度为0的SCC就好啦。

但是还需要一个特判，就是其它点都跑完了以后，若都安全，那么它自然就是杀手点，然后被杀的概率就要 - -。

那个，这个逻辑有点诡异，至少我认为它不对。
但是它是AC的……

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 101000
#define M 301000
#define inf 0x3f3f3f3f
using namespace std;
int n,m,ans,ind,top;
struct Eli
{
    int v,next,len;
}e[M],E[M];
int head[N],HEAD[N],d[N],cnt;
inline void add1(int u,int v)
{
    e[++cnt].v=v;
    e[cnt].next=head[u];
    head[u]=cnt;
}
inline void add2(int u,int v)
{
    d[v]++;
    E[++cnt].v=v;
    E[cnt].next=HEAD[u];
    HEAD[u]=cnt;
}
int scc;
int low[N],dfn[N],q[N],id[N],num[N];
bool inq[N],flag[N];
void tarjan(int x)
{
    low[x]=dfn[x]=++ind;
    inq[x]=1;q[++top]=x;
    for(int i=head[x];i;i=e[i].next)
    {
        int v=e[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(inq[v])low[x]=min(low[x],dfn[v]);
    }
    if(low[x]==dfn[x])
    {
        int now=0;scc++;
        while(x!=now)
        {
            now=q[top];top--;
            id[now]=scc;
            inq[now]=0;
            num[scc]++;
        }
    }
}
bool check(int x)
{
    if(d[x]!=0||num[x]!=1)return 0; // 如果一个点没有入度，并且里面人数不为1
    for(int i=HEAD[x];i;i=E[i].next)
        if(d[E[i].v]==1)return 0; // 并且所有出点都存在其它入点
    return 1; // 就ans--
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        add1(u,v);
    }
    for(int i=1;i<=n;i++)if(!dfn[i])tarjan(i);
    cnt=0;
    for(int x=1;x<=n;x++) //重建图（带去重的那种）
    {
        for(int i=head[x];i;i=e[i].next)
        {
            int v=e[i].v;
            if(id[x]!=id[v]&&!flag[id[v]])
                flag[id[v]]=1,add2(id[x],id[v]);
        }
        for(int i=head[x];i;i=e[i].next)
        { 
            int v=e[i].v;
            if(id[x]!=id[v])flag[id[v]]=0;
        }
    }
    for(int i=1;i<=scc;i++)if(!d[i])ans++;
    for(int i=1;i<=scc;i++)if(check(i)){ans--;break;}
    printf("%.6lf",(double)(n-ans)/n);
    return 0;
}