记录编号	583469	评测结果	AAAAAAAAAAAAA
题目名称	[BZOJ 2438]杀人游戏	最终得分	100
用户昵称	┭┮﹏┭┮	是否通过	通过
代码语言	C++	运行时间	0.330 s
提交时间	2023-10-15 15:21:48	内存使用	6.99 MiB

#include <bits/stdc++.h> 
using namespace std;
//tarjan求强连通分量 + 概率知识！ 
const int N = 1e5+10,M = 3e5+10;
int n,m;
struct made{
    int ver,nx;
}e[M];
int hd[N],tot,cnt,top,num;
int low[N],dfn[N],st[N],color[N],size[N];
vector<int>scc[N];// 
int v[N];
map<pair<int,int>,bool>mp;
void add(int x,int y){
    tot++;
    e[tot].ver = y,e[tot].nx = hd[x],hd[x] = tot;
}
void tarjan(int x){
    low[x] = dfn[x] = ++cnt;
    st[++top] = x,v[x] = 1;
    for(int i = hd[x];i;i = e[i].nx){
        int y = e[i].ver;
        if(!dfn[y])tarjan(y),low[x] = min(low[x],low[y]);
        else if(v[y])low[x] = min(low[x],dfn[y]);
    }
    if(low[x] == dfn[x]){
        num++;int y = 0;
        do{
            y = st[top--];
            v[y] = 0,color[y] = num,size[num]++;
            scc[num].push_back(y);
        }while(x != y);
    }
}//正常tarjan 
void push(int x){
    for(int i = hd[x];i;i = e[i].nx){
        int y = e[i].ver;
        if(color[x] == color[y])continue;
        if(mp[make_pair(color[x],color[y])])continue;//map判重边 
        v[color[y]]++;
        mp[make_pair(color[x],color[y])] = 1;//必须要判重边，不然会多记录入度，影响结果 
    }
}//求缩点后强连通分量入度 
bool check(int x){
    for(int i = hd[x];i;i = e[i].nx){
        int y = e[i].ver;
        if(v[color[y]] == 1)return 0;
    }
    return 1;
}
int main(){
    freopen("kill_game.in","r",stdin);
    freopen("kill_game.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i = 1;i <= m;i++){
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);
    }
    for(int i = 1;i <= n;i++)
        if(!color[i])tarjan(i);
    memset(v,0,sizeof(v));// 
    for(int i = 1;i <= n;i++)
        push(i);
    double ans = 0;// 
    for(int i = 1;i <= num;i++)
        if(!v[i])ans++;//求入度为零的强连通分量 
    for(int i = 1;i <= num;i++){
        if(v[i] || size[i] != 1)continue;
        int x = scc[i][0];
        if(check(x)){
            ans--;
            break;
        }
    }// 
    double fn = n;//double类型 
    printf("%.6lf\n",1.0-1.0*ans/n);// 
    
    
    return 0;
    
}