bzoj2438: [中山市选2011]杀人游戏(强联通+特判)

2438: [中山市选2011]杀人游戏

简要题意：

　　给出n个点，m条有向边，进行最少的访问并且可以便利（n-1）个点，求这个方案成功的概率

题解：
　　一道非常好的题目！

　　题目要知道最大的存活概率，那么也就是找到直接找到杀手的最小概率

　　那么我们采用强联通缩点：

　　统计每个联通分量的入度，如果入度为0（证明除此联通分量里的点，没有人可以知道连通分量里的信息，那就一定要先选一个人访问），那么sum++（因为依据题意，假如问到连通分量里的任意一个人，只要ta不是杀手，那么一定可以安全的遍历强联通分量里的所有人）

　　接下来就是最关键的特判：

　　对于连通分量里只有一个点的情况，如果它的入度为零，不一定就要访问（这里何前面似乎有些许矛盾）

　　解释：如果有联通分量只有一个家族成员，并且没有入度，因为杀手只有一个，那么我们完全可以在遍历其他的n-1个点时得出答案（排除法嘛）

　　但是它的出度不一定为0，所以还要判断一下在上面的基础上，这个点连出去的边是否能够被其他点访问（入度>1），很好理解吧

　　如果以上的都满足，那么sum--

　　还有一个小槽点：sum--一次就好了，因为如果有两种相同的情况，那么对于两个独立的点，我们还是要选择一个访问才能遍历完成的

　　输出1.0-double(sum)/double(n)

代码：

  1 #include<cstdio>
  2 #include<cstring>
  3 #include<cstdlib>
  4 #include<cmath>
  5 #include<algorithm>
  6 using namespace std;
  7 int n,m;
  8 struct node
  9 {
 10     int x,y,next;
 11 }a[1110000];int len,last[510000];
 12 void ins(int x,int y)
 13 {
 14     len++;
 15     a[len].x=x;a[len].y=y;
 16     a[len].next=last[x];last[x]=len;
 17 }
 18 struct edge
 19 {
 20     int x,y,next;
 21 }e[1110000];int len1,last1[510000];
 22 void inss(int x,int y)
 23 {
 24     len1++;
 25     e[len1].x=x;e[len1].y=y;
 26     e[len1].next=last1[x];last1[x]=len1;
 27 }
 28 int cnt,tp,id;
 29 int belong[510000],dfn[510000],low[510000],sta[510000],size[510000];
 30 bool v[510000];
 31 void dfs(int x)
 32 {
 33     low[x]=dfn[x]=++id;
 34     sta[++tp]=x;v[x]=true;
 35     for(int k=last[x];k;k=a[k].next)
 36     {
 37         int y=a[k].y;
 38         if(dfn[y]==-1)
 39         {
 40             dfs(y);
 41             low[x]=min(low[x],low[y]);
 42         }
 43         else
 44         {
 45             if(v[y]==true)
 46                 low[x]=min(low[x],dfn[y]);
 47         }
 48     }
 49     if(low[x]==dfn[x])
 50     {
 51         int i;cnt++;
 52         do{
 53             i=sta[tp--];
 54             v[i]=false;
 55             belong[i]=cnt;
 56             size[cnt]++;
 57         }while(i!=x);
 58     }
 59 }
 60 int ru[510000],chu[510000],sum1,sum2;
 61 bool check(int x)
 62 {
 63     if(size[x]!=1)return true;
 64     if(last1[x]==0)return false;
 65     for(int k=last1[x];k;k=e[k].next)
 66     {
 67         int y=e[k].y;
 68         if(ru[y]<=1)return true;
 69     }
 70     return false;
 71 }
 72 int main()
 73 {
 74     cnt=tp=id=sum1=sum2=0;
 75     len=0;len1=0;
 76     memset(last,0,sizeof(last));memset(last1,0,sizeof(last1));
 77     scanf("%d%d",&n,&m);
 78     memset(chu,0,sizeof(chu));
 79     memset(ru,0,sizeof(ru));
 80     memset(sta,0,sizeof(sta));
 81     memset(dfn,-1,sizeof(dfn));
 82     memset(low,0,sizeof(low));
 83     memset(v,false,sizeof(v));
 84     for(int i=1;i<=m;i++)
 85     {
 86         int x,y;
 87         scanf("%d%d",&x,&y);
 88         ins(x,y);
 89     }
 90     for(int i=1;i<=n;i++)
 91         if(dfn[i]==-1)
 92             dfs(i);
 93     for(int i=1;i<=m;i++)
 94     {
 95         int st=belong[a[i].x],ed=belong[a[i].y];
 96         if(st!=ed)
 97         {
 98             ru[ed]++;
 99             inss(st,ed);
100         }
101     }
102     for(int i=1;i<=cnt;i++)
103         if(ru[i]==0)sum1++;
104     for(int i=1;i<=cnt;i++)
105         if(ru[i]==0)
106             if(check(i)==false){sum1--;break;}
107     printf("%.6lf
",1.0-double(sum1)/double(n));
108     return 0;
109 }