[Apio2009]Atm

题目描述

输入

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号

输出

输出一个整数，表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。

样例输入

6 7 1 2 2 3 3 5 2 4 4 1 2 6 6 5 10 12 8 16 1 5 1 4 4 3 5 6

样例输出

47

提示

50%的输入保证N, M<=3000。所有的输入保证N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

抢劫是不对的，要弘扬社会主义核心价值观。

缩点裸题

Tarjan求联通分量缩点，再spfa

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct Messi
{
    int next,to;
}edge[500001],edge2[500001];
int head[500001],num,low[500001],dfn[500001],dfscnt,stack[500001],inStack[500001];
int point[500001],n,m,scnt,top,sccno[500001],ans,num2,head2[500001],w[500001],dist[500001],q[5000001];
bool b2[500001];
bool b[500001];
void add(int u,int v)
{
    num++;
    edge[num].next=head[u];
     edge[num].to=v;
     head[u]=num;
}
void add2(int u,int v)
{
    num2++;
    edge2[num2].next=head2[u];
     edge2[num2].to=v;
     head2[u]=num2;
}
void dfs(int u)
{int i,j;
    low[u]=dfn[u]=++dfscnt;
    stack[++top]=u;
     inStack[u]=1;
     for (i=head[u];i;i=edge[i].next)
     {
        int v=edge[i].to;
         if (dfn[v]==0)
         {
            dfs(v);
            low[u]=min(low[u],low[v]);
         }else if (inStack[v]) low[u]=min(low[u],dfn[v]);
     }
     if (dfn[u]==low[u])
     {
        ++scnt;
        while (top&&stack[top+1]!=u)
        {
            sccno[stack[top]]=scnt;
            inStack[stack[top--]]=0;
        }
     }
}
int main()
{int i,j,u,v,s,p,k,h,t;
//freopen("2.in","r",stdin);
//freopen("2.out","w",stdout);
    cin>>n>>m;
     for (i=1;i<=m;i++)
      {
        scanf("%d%d",&u,&v);
        add(u,v);
      }
       for (i=1;i<=n;i++)
       scanf("%d",&w[i]);
     cin>>s>>p;
      for (i=1;i<=p;i++)
      {
        scanf("%d",&k);
         b[k]=1;
      }
      for (i=1;i<=n;i++)
       if (dfn[i]==0) dfs(i);
       for (i=1;i<=n;i++)
       {
         int u=sccno[i];
            point[u]+=w[i];
            if (b[i]) b2[u]=1;
          for (j=head[i];j;j=edge[j].next)
           {
            int v=sccno[edge[j].to];
            if (u!=v) add2(u,v);
           }
       }
      q[1]=sccno[s];
      dist[sccno[s]]=point[sccno[s]];
       h=0;t=1;
       while (h<t)
       {
        h++;
        u=q[h];
         for (i=head2[u];i;i=edge2[i].next)
         {
            v=edge2[i].to;
             if (dist[u]+point[v]>dist[v])
              {
                  dist[v]=dist[u]+point[v];
              t++;
              q[t]=v;
             }
         }
       }  
    for (i=1;i<=scnt;i++)
     if (b2[i])
      ans=max(ans,dist[i]);
      cout<<ans;
}