bzoj 1179 Atm

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1179

题解：

　　一道比较综合的图论题

　　直接讲正解：
　　如果这个图G中存在某个强连通分量，那么这个强连通分量中的所有ATM即可视为都被抢到，所有的酒吧都可视为重点，并且也可以从这个强连通分量的任何结点出发继续向外扩展

　　所以先做一遍Tarjan，找出强连通分量，然后重新构图，把每个强连通分量缩成一个点，此点的权值即为原先强连通分量里所有点权之和，判断此点中有没有酒吧，再将原先所有连接强连通分量的边连接在这个点

　　最后做SPFA，找出权值最大的路径

　　这道题时间限制15s，空间162MB，所以一般不会TLE或MLE

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 #define MAXN 500010
 #define MAXM 500010
 int n,m,time,cnt,head1[MAXN],DFN[MAXN],Low[MAXN],stack[MAXN],top,a1[MAXN],belong[MAXN];
 bool vis[MAXN],bar1[MAXN];
 struct node1
 {
     int v,next;
 }edge1[MAXM];
 //以上变量基本是在进行tarjan（及之前）使用
 int heads[MAXN],mm,a2[MAXN],s,p,d[MAXN],q[MAXN],head,tail,ans;
 bool viss[MAXN];
 struct node2
 {
     int v,next;
 }edge2[MAXM];
 //以上变量是在进行SPFA时使用
 //a存放点权，bar记录是否有酒吧
 bool bar2[MAXN];
 void add(int x,int y)
 {
     edge1[++cnt].next=head1[x];
     head1[x]=cnt;
     edge1[cnt].v=y;
 }
 void buildmap(int k)//重新建图、连边
 {
     for(int i=;i<=n;i++)
     {
         if(belong[i]==k)
         {
             for(int j=head1[i];j;j=edge1[j].next)
             {
                 if(belong[edge1[j].v]==k)continue;
                 edge2[++m].next=heads[k];
                 heads[k]=m;
                 edge2[m].v=belong[edge1[j].v];
             }
         }
     }
 }
 void tarjan(int u)
 {
     DFN[u]=Low[u]=++time;
     vis[u]=true;
     stack[++top]=u;
     for(int i=head1[u];i;i=edge1[i].next)
     {
         int v=edge1[i].v;
         if(!DFN[v])
         {
             tarjan(v);
             Low[u]=min(Low[u],Low[v]);
         }
         else if(vis[v])Low[u]=min(Low[u],DFN[v]);
     }
     if(DFN[u]==Low[u])
     {
         int i;
         cnt++;
         do
         {
             i=stack[top--];
             vis[i]=false;
             belong[i]=cnt;
             a2[cnt]+=a1[i];
             if(bar1[i])bar2[cnt]=true;//先更新点权和是否有酒吧
         }while(u!=i);
     }
 }
 void SPFA()
 {
     head=;tail=;
     q[]=s;
     viss[s]=;
     while(head<tail)
     {
         for(int i=heads[q[head]];i!=;i=edge2[i].next)
         {
             if(d[q[head]]+a2[edge2[i].v]>d[edge2[i].v])//注意，这里不是边权，是点权
             {
                 d[edge2[i].v]=d[q[head]]+a2[edge2[i].v];
                 if(!viss[edge2[i].v])
                 {
                     q[tail++]=edge2[i].v;
                     viss[edge2[i].v]=true;
                 }
             }
         }
         viss[q[head]]=false;
         head++;
     }
 }
 int main()
 {
     scanf("%d%d",&n,&m);
     for(int i=;i<=m;i++)
     {
         int x,y;
         scanf("%d%d",&x,&y);
         add(x,y);
     }
     for(int i=;i<=n;i++)
     {
         scanf("%d",&a1[i]);
     }
     scanf("%d%d",&s,&p);
     for(int i=;i<=p;i++)
     {
         int x;
         scanf("%d",&x);
         bar1[x]=true;
     }
     cnt=;
     m=;
     for(int i=;i<=n;i++)
         if(!DFN[i])tarjan(i);
     for(int i=;i<=cnt;i++)
         buildmap(i);
     s=belong[s];
     for(int i=;i<=cnt;i++)
     {
         d[i]=a2[i];
     }
     SPFA();
     for(int i=;i<=cnt;i++)
     {
         if(bar2[i])ans=max(ans,d[i]);
     }
     printf("%d",ans);
     return ;
 }