[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[],edge2[];
 int head[],num,low[],dfn[],dfscnt,stack[],inStack[];
 int point[],n,m,scnt,top,sccno[],ans,num2,head2[],w[],dist[],q[];
 bool b2[];
 bool b[];
 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]=;
      for (i=head[u];i;i=edge[i].next)
      {
         int v=edge[i].to;
          if (dfn[v]==)
          {
             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+]!=u)
         {
             sccno[stack[top]]=scnt;
             inStack[stack[top--]]=;
         }
      }
 }
 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=;i<=m;i++)
       {
         scanf("%d%d",&u,&v);
         add(u,v);
       }
        for (i=;i<=n;i++)
        scanf("%d",&w[i]);
      cin>>s>>p;
       for (i=;i<=p;i++)
       {
         scanf("%d",&k);
          b[k]=;
       }
       for (i=;i<=n;i++)
        if (dfn[i]==) dfs(i);
        for (i=;i<=n;i++)
        {
          int u=sccno[i];
             point[u]+=w[i];
             if (b[i]) b2[u]=;
           for (j=head[i];j;j=edge[j].next)
            {
             int v=sccno[edge[j].to];
             if (u!=v) add2(u,v);
            }
        }
       q[]=sccno[s];
       dist[sccno[s]]=point[sccno[s]];
        h=;t=;
        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=;i<=scnt;i++)
      if (b2[i])
       ans=max(ans,dist[i]);
       cout<<ans;
 }