P2255【L1 SOLO 第五场 APIO2009】抢掠计划

问题描述

Siruseri城中的道路都是单向的。不同的道路由路口连接。按照法律的规定，在每个路口都设立了一个Siruseri银行的ATM取款机。令人奇怪的是，Siruseri的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。

Banditji计划实施Siruseri有史以来最惊天动地的ATM抢劫。他将从市中心出发，沿着单向道路行驶，抢劫所有他途径的ATM机，最终他将在一个酒吧庆祝他的胜利。

使用高超的黑客技术，他获知了每个ATM机中可以掠取的现金数额。他希望你帮助他计算从市中心出发最后到达某个酒吧时最多能抢劫的现金总数。他可以经过同一路口或道路任意多次。但只要他抢劫过某个ATM机后，该ATM机里面就不会再有钱了。

例如，假设该城中有6个路口，道路的连接情况如下图所示：

市中心在路口1，由一个入口符号→来标识，那些有酒吧的路口用双圈来表示。每个ATM机中可取的钱数标在了路口的上方。在这个例子中，Banditji能抢劫的现金总数为47，实施的抢劫路线是：1-2-4-1-2-3-5。

输入格式

第一行包含两个整数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的单向的道路到达其中的至少一个酒吧。

题解

scc缩点+spfa
注意扩栈！

代码

#include<iostream>  #include<algorithm>  #include<cstdio>  #include<cmath>  #include<cstring>  #include<queue>  #include<vector>  #include<map>  #include<stack>  #define inf 1e9  #define maxn1 500005  using namespace std;#define f(x) for(i=1;i<=x;i++)  int n,m,s,p;  int Len[maxn1];  int Last[maxn1],Next[maxn1],End[maxn1],St[maxn1];  int be[maxn1],v[maxn1],nv[maxn1];  int dfn[maxn1],low[maxn1];    bool nb[maxn1],instack[maxn1];  bool bar[maxn1];  int vis,scc,cnt;  stack<int>q;  void tj(int x)  {          int i,j,k;          dfn[x]=low[x]=++vis;          q.push(x);instack[x]=true;          for(i=Last[x];i;i=Next[i])          {                  int en=End[i];                  if(dfn[en]==0){                          tj(en);                          low[x]=min(low[x],low[en]);                  }                  else if(instack[en])                  {                      low[x]=min(dfn[en],low[x]);                      }          }          if(low[x]==dfn[x])          {                  scc++;                  do{                          k=q.top();                          q.pop();                          be[k]=scc;                          nv[scc]+=v[k];                          instack[k]=false;                          if(bar[k]) nb[scc]=true;                  }while(k!=x);          }  }int dis[maxn1];bool mark[maxn1];queue<int>qq;void spfa(){    dis[s]=nv[s];    qq.push(s);    mark[s]=true;    while(qq.size())    {        int x=qq.front();        qq.pop();        mark[x]=false;        for(int i=Last[x];i;i=Next[i])        {            int en=End[i];            if(dis[en]<dis[x]+Len[i])            {                dis[en]=dis[x]+Len[i];                if(mark[en]==false) {                mark[en]=true;                qq.push(en);}            }        }    }}  void insert(int x,int y)  {          cnt++;        St[cnt]=x;          Next[cnt]=Last[x];          Last[x]=cnt;          End[cnt]=y;  }  void insert2(int x,int y,int z)  {          cnt++;          Next[cnt]=Last[x];          Last[x]=cnt;          End[cnt]=y;        Len[cnt]=z;  }  int main_main()  {          int i,j;       scanf("%d%d",&n,&m);       for(i=1;i<=m;i++)       {               int x,y;               scanf("%d%d",&x,&y);               insert(x,y);       }       for(i=1;i<=n;i++) scanf("%d",&v[i]);       scanf("%d%d",&s,&p);       for(i=1;i<=p;i++)       {               int x;               scanf("%d",&x);               bar[x]=true;       }       for(i=1;i<=n;i++)       {               if(dfn[i]==0) tj(i);       }       s=be[s];//   cout<<s<<" ";     memset(Last,0,sizeof(Last));     cnt=0;     for(i=1;i<=m;i++)       {                      int a=St[i],b=End[i];                      if(be[a]!=be[b])                      {                              insert2(be[a],be[b],nv[be[b]]);//                            cout<<be[a]<<"to"<<be[b];             }       }       spfa();     int maxx=0;     for(i=1;i<=scc;i++)     {        if(nb[i]) maxx=max(dis[i],maxx);     }//   f(scc) cout<<dis[scc]<<" ";     cout<<maxx;//    f(scc) cout<<nv[i]<<" ";}const int main_stack=16;  char my_stack[128<<20];int main() {  __asm__("movl %%esp, (%%eax);\n"::"a"(my_stack):"memory");  __asm__("movl %%eax, %%esp;\n"::"a"(my_stack+sizeof(my_stack)-main_stack):"%esp");  main_main();  __asm__("movl (%%eax), %%esp;\n"::"a"(my_stack):"%esp");  return 0;}