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;}