BZOJ 2438 [中山市选2011] 杀人游戏 （图的联通，tarjan）

题目链接：BZOJ 2438

首先用tarjan缩点，然后重新建图，找到入度为0的强联通分量，ans为其大小之和。要特判一种情况，即缩点之后的这个点的入度为0，但是它不存在出度，或者是以它引出去的点可以被其他点判断，那么ans--。注意这种情况只能减去一次。代码后面写的有点烦= =。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<stack>
using namespace std;

#define maxn (100000+10)
int N,M,idx=0,tot=0,ans=0,W=0;
int head[maxn],belong[maxn],dfn[maxn],low[maxn],in[maxn],sai[maxn];
bool vis[maxn];
struct node{int v,next;}e[3*maxn];

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 k=0;
void adde(int u,int v){
	e[k].v=v; e[k].next=head[u]; head[u]=k++;
}

stack<int>q;
void tarjan(int x){
	dfn[x]=low[x]=++idx; vis[x]=1;
	q.push(x);
	for(int i=head[x];i!=-1;i=e[i].next){
		int v=e[i].v;
		if(!dfn[v]){
			tarjan(v); low[x]=min(low[x],low[v]);
		}
		else if(vis[v])low[x]=min(low[x],dfn[v]);
	}
	if(dfn[x]==low[x]){
		tot++; int now=0,t=0;
		while(now!=x){
			now=q.top(); q.pop(); t++;
			vis[now]=0; belong[now]=tot;
		}
		if(t==1)sai[++W]=x;
	}
}

void rebuilt(){
	for(int i=1;i<=N;i++){
		for(int j=head[i];j!=-1;j=e[j].next){
			int v=e[j].v;
			if(belong[i]!=belong[v])in[belong[v]]++;
		}
	}
}

int main(){
	freopen("sai.in","r",stdin);
	freopen("sai.out","w",stdout);
	
	memset(head,-1,sizeof(head));
	N=read(); M=read();
	for(int i=1;i<=M;i++){
		int x=read(),y=read();
		adde(x,y);
	}
	for(int i=1;i<=N;i++)if(!dfn[i])tarjan(i);
	rebuilt();
	for(int i=1;i<=tot;i++){
		if(!in[i])ans++;
	}
	for(int i=1;i<=W;i++){
		if(in[belong[sai[i]]]!=0)continue;
		int col=1;
		for(int j=head[sai[i]];j!=-1;j=e[j].next){
			int v=e[j].v;
			if(in[belong[v]]<=1){col=0;break;}
		}
		if(col){ans--;break;}
	}
	printf("%.6lf",(double)(N-ans)/N);
	return 0;
}