p3627&bzoj1179 抢掠计划（ATM)

传送门（洛谷）

传送门（bzoj）

题目

Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定， 在每个路口都设立了一个 Siruser

i 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。Bandit

ji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心 出发，沿着单向道路行驶，抢劫所有他

途径的 ATM 机，最终他将在一个酒吧庆 祝他的胜利。使用高超的黑客技术，他获知了每个 ATM 机中可以掠取的

现金数额。他希 望你帮助他计算从市中心出发最后到达某个酒吧时最多能抢劫的现金总数。他可 以经过同一路口

或道路任意多次。但只要他抢劫过某个 ATM 机后，该 ATM 机 里面就不会再有钱了。 例如，假设该城中有 6 个

路口，道路的连接情况如下图所示：

市中心在路口 1，由一个入口符号→来标识，那些有酒吧的路口用双圈来表示。每个 ATM 机中可取的钱数标在了

路口的上方。在这个例子中，Banditji 能抢 劫的现金总数为 47，实施的抢劫路线是：1-2-4-1-2-3-5。

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。

接下来M行，每行两个整数，这两个整数都在1到N之间，

第i+1行的两个整数表示第i条道路的起点和终点的路口编号。

接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。

接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。

接下来的一行中有P个整数，表示P个有酒吧的路口的编号

N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。

输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

Output

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

Sample Input

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

Sample Output

47

分析

Tarjan缩点的好题，先用Tarjan缩点，求出每个环的总钱数和这个环中是否有酒吧，然后重新建图跑spfa求最大点权和，但由于要在酒吧结束，所以最终答案只能在有酒吧的环中取。

插曲

作为一个最短路只会写folyd和dijstar的蒟蒻，我自己yy了一个spfa，然而它居然对了（QAQ）

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
int bar[510000],mey[510000],belong[510000],cnt,tmy[510000],tbr[510000];
int dfn[510000],low[510000],vis[510000],sum;
int d[510000],ist[510000],maxn,iqu[510000];
vector<int>v[500010];
vector<int>nv[500010];
stack<int>a;
queue<pair<int,int> >q;
int n,m,s,p;
inline void read(int &x){
      int f=1;x=0;
      char s=getchar();
      while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
      while(s>='0'&&s<='9'){x=x*10+(s-'0');s=getchar();}
      x*=f;
}
inline void tarjan(int x){
      low[x]=dfn[x]=++cnt;
      a.push(x);
      ist[x]=1;
      int i,j,k;
      for(i=0;i<v[x].size();i++)
         if(!dfn[v[x][i]]){
             tarjan(v[x][i]);
             low[x]=min(low[x],low[v[x][i]]);
         }else if(ist[v[x][i]]){
             low[x]=min(low[x],dfn[v[x][i]]);
         }
      if(dfn[x]==low[x]){
          sum++;
          int ok=0;
          while(1){
              int u=a.top();
              ist[u]=0;
              if(bar[u])ok=1;
              belong[u]=sum;
              tmy[sum]+=mey[u];
              a.pop();
              if(u==x)break;
          }
          if(ok)tbr[sum]=1;
      }
      return;
}
int main()
{     int i,j,k,x,y;
      read(n),read(m);
      for(i=1;i<=m;i++){
          read(x),read(y);
          v[x].push_back(y);
      }
      for(i=1;i<=n;i++){
          read(mey[i]);
      }
      scanf("%d%d",&s,&p);
      for(i=1;i<=p;i++){
          read(x);
          bar[x]=1;
      }
      for(i=1;i<=n;i++)
         if(!dfn[i])
           tarjan(i);
      for(i=1;i<=n;i++)
         for(j=0;j<v[i].size();j++)
            if(belong[i]!=belong[v[i][j]]){
                nv[belong[i]].push_back(belong[v[i][j]]);
            }
      d[belong[s]]=tmy[belong[s]];
      q.push(make_pair(belong[s],tmy[belong[s]]));
      iqu[belong[s]]=1;
      while(!q.empty()){
          int u=q.front().first;
          q.pop();
          iqu[u]=0;
          if(tbr[u])maxn=max(maxn,d[u]);
          for(i=0;i<nv[u].size();i++){
            d[nv[u][i]]=max(d[nv[u][i]],d[u]+tmy[nv[u][i]]);
            if(!iqu[nv[u][i]]){
              q.push(make_pair(nv[u][i],d[nv[u][i]]));
              iqu[nv[u][i]]=1;
            }
          }
      }
      cout<<maxn<<endl;
      return 0;
}