【BZOJ2438】 [中山市选2011]杀人游戏 tarjan强连通分量+缩点

Description

一位冷血的杀手潜入 Na-wiat，并假装成平民。警察希望能在 N 个人里面，
查出谁是杀手。 
警察能够对每一个人进行查证，假如查证的对象是平民，他会告诉警察，他
认识的人， 谁是杀手， 谁是平民。 假如查证的对象是杀手， 杀手将会把警察干掉。 
现在警察掌握了每一个人认识谁。 
每一个人都有可能是杀手，可看作他们是杀手的概率是相同的。 
问：根据最优的情况，保证警察自身安全并知道谁是杀手的概率最大是多
少？

Input

第一行有两个整数 N,M。 
接下来有 M 行，每行两个整数 x,y，表示 x 认识 y（y 不一定认识 x,例如hjt同志） 。

Output

仅包含一行一个实数，保留小数点后面 6 位，表示最大概率。

Sample Input

5 4
1 2
1 3
1 4
1 5

Sample Output

0.800000

HINT

警察只需要查证 1。假如1是杀手，警察就会被杀。假如 1不是杀手，他会告诉警

察 2,3,4,5 谁是杀手。而 1 是杀手的概率是 0.2,所以能知道谁是杀手但没被杀的概

率是0.8。对于 100%的数据有 1≤N ≤  10 0000,0≤M ≤  30 0000

数据已加强！

tarjan求强连通+缩点后，统计入度为0的点，最后特判一下（并不知道为什么）。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <queue>
 #include <stack>
 #define N 500500
 #define M 500500
 using namespace std;
 struct data{int next,p;}e1[M],e2[M];
 int head1[N],head2[N],belong[N],low[N],dfn[N],mon[N],vscc[N],dis[N],inq[N],vis[N];
 int q1[N],q[N];
 int n,m,S,p,cnt,ans,qt,scc;
 inline int max(int x,int y){return x>y?x:y;}
 inline int min(int x,int y){return x<y?x:y;}
 inline int read()
 {
     char c;
     int anss=;
     while ((c=getchar())==' ' || c=='\n' || c=='\r');
     anss=c-'';
     while (isdigit(c=getchar())) anss=anss*+c-'';
     return anss;
 }
 void se1(int x,int y){cnt++; e1[cnt].next=head1[x]; head1[x]=cnt; e1[cnt].p=y;}
 void se2(int x,int y){cnt++; e2[cnt].next=head2[x]; head2[x]=cnt; e2[cnt].p=y;}
 void tarjan(int x)
 {
     vis[x]=inq[x]=;
     dfn[x]=low[x]=++cnt;
     q1[++qt]=x;
     for (int i=head1[x];i;i=e1[i].next)//i全打成xSB调了好久。
     {
         if (!vis[e1[i].p])
         {
             tarjan(e1[i].p);
             low[x]=min(low[x],low[e1[i].p]);
         }
         else if (inq[e1[i].p])    low[x]=min(low[x],low[e1[i].p]);
     }
     if (low[x]==dfn[x])
     {
         int now=;
         scc++;
         while (now!=x)
         {
             now=q1[qt];qt--;
             belong[now]=scc;
             inq[now]=;
             vscc[scc]+=mon[now];
         }
     }
 }
 void part1_tarjan()
 {
     cnt=;
     for (int i=;i<=n;i++)
         if (!vis[i])    tarjan(i);
 }
 void part2_shr_point()
 {
     cnt=;
     for (int i=;i<=n;i++)
         for (int j=head1[i];j;j=e1[j].next)
             if (belong[i]!=belong[e1[j].p])
                 se2(belong[i],belong[e1[j].p]);
 }
 void part3_spfa()
 {
     memset(inq,,sizeof(inq));
     int t=,w=;
     q[t]=belong[S];
     inq[belong[S]]=;
     dis[belong[S]]=vscc[belong[S]];
     while  (t!=w)
     {
         int now=q[t];t++;
         for (int i=head2[now];i;i=e2[i].next)
             if (dis[e2[i].p]<dis[now]+vscc[e2[i].p])
             {
                 dis[e2[i].p]=dis[now]+vscc[e2[i].p];
                 if (!inq[e2[i].p])
                 {
                     inq[e2[i].p]=;
                     q[w++]=e2[i].p;
                 }
             }
         inq[now]=;
     }
 }
 int main()
 {
     n=read();m=read();
     for (int i=;i<=m;i++)
     {
         int x,y;
         x=read();y=read();
         se1(x,y);
     }
     for (int i=;i<=n;i++)    mon[i]=read();
     part1_tarjan();
     part2_shr_point();
     S=read();p=read();
     part3_spfa();
     for (int i=;i<=p;i++)
     {
         int x=read();
         if (dis[belong[x]]>ans) ans=dis[belong[x]];
     }
     printf("%d\n",ans);
     return ;
 }