p3627[APIO2009] 抢掠计划

Posted 2021-01-12 顾z

Description

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

Banditji 计划实施 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 个有酒吧的路口的编号。

Output

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

(Tarjan)缩点跑最长路,

在求强联通分量的时候维护每个联通分量的(val)(即当前强联通分量中所有(atm)的钱数.

不能再次建边的时候再求.(会出锅,但是原理我不太知道。

再对强联通分量建边跑最长路即可.

注意：(Dijkstra)的贪心策略不可跑最长路。

这里数组可以重复使用,不过这题没有卡数组,就没有重复使用了 qwq.

代码

#include<cstdio>
#include<cctype>
#include<queue>
#define N 500008
#define R register
using namespace std;
inline void in(int &x)
{
    int f=1;x=0;char s=getchar();
    while(!isdigit(s)){if(s=='-')f=-1;s=getchar();}
    while(isdigit(s)){x=x*10+s-'0';s=getchar();}
    x*=f;
}
int n,m,head[N],tot,low[N],dfn[N],col,belong[N],idx,stk[N],pos;
int val[N],v[N],h[N],ttt,p,t[N],dis[N],top,s,ans;
struct cod{int u,v;}edge[N<<2],e[N<<2];
bool inq[N],vis[N];
inline void add(int x,int y)
{
    edge[++tot].u=head[x];
    edge[tot].v=y;
    head[x]=tot;
}
inline void ado(int x,int y)
{
    e[++ttt].u=h[x];
    e[ttt].v=y;
    h[x]=ttt;
}
inline void spfa(int s)
{
    for(R int i=1;i<=col;i++)dis[i]=-214748364;
    queue<int>q;
    q.push(s);vis[s]=true;dis[s]=val[s];
    while(!q.empty())
    {
        int u=q.front();q.pop();vis[u]=false;
        for(R int i=h[u];i;i=e[i].u)
        {
            if(dis[e[i].v]<dis[u]+val[e[i].v])
            {
                dis[e[i].v]=dis[u]+val[e[i].v];
                if(!vis[e[i].v])
                {
                    vis[e[i].v]=true;
                    q.push(e[i].v);
                }
            }
        }
    }
    for(R int i=1;i<=p;i++)
        ans=max(ans,dis[belong[t[i]]]);
    printf("%d",ans);
}
void tarjan(int x)
{
    dfn[x]=low[x]=++idx;
    stk[++top]=x;inq[x]=true;
    for(R int i=head[x];i;i=edge[i].u)
    {
        if(!dfn[edge[i].v])
        {
            tarjan(edge[i].v);
            low[x]=min(low[x],low[edge[i].v]);
        }
        else if(inq[edge[i].v])
            low[x]=min(low[x],dfn[edge[i].v]);
    }
    if(low[x]==dfn[x])
    {
        col++;
        int now=-1;
        while(now!=x)
        {
            now=stk[top--];
            inq[now]=false;
            belong[now]=col;
            val[col]+=v[now];
        }
    }
}
int main()
{
    in(n),in(m);
    for(R int i=1,x,y;i<=m;i++)
    {
        in(x),in(y);
        add(x,y);
    }
    for(R int i=1;i<=n;i++)in(v[i]);
    for(R int i=1;i<=n;i++)
        if(!dfn[i])tarjan(i);
    for(R int i=1;i<=n;i++)
        for(R int j=head[i];j;j=edge[j].u)
            if(belong[edge[j].v]!=belong[i])
                ado(belong[i],belong[edge[j].v]);
    in(s);in(p);
    for(R int i=1;i<=p;i++)in(t[i]);
    spfa(belong[s]);
}