bzoj 1179 [Apio2009]Atm 缩点+最短路

[Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 4290  Solved: 1893
[Submit][Status][Discuss]

Description

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

HINT

题解：Tarjan缩点+最短路很显然的吧

 

#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<queue>

#define N 500007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,S,p;
int dis[N],w[N],val[N];
int cnt,hed[N],rea[N<<1],nxt[N<<1];
int cnt1,hed1[N],rea1[N<<1],nxt1[N<<1];
int tim,tot,scc,bel[N],stk[N],dfn[N],low[N];
bool flag[N],instk[N],flag1[N];

void add(int u,int v)
{
    nxt[++cnt]=hed[u];
    hed[u]=cnt;
    rea[cnt]=v;
}
void add1(int u,int v)
{
    nxt1[++cnt1]=hed1[u];
    hed1[u]=cnt1;
    rea1[cnt1]=v;
}
void Tarjan(int u)
{
    dfn[u]=low[u]=++tim;
    stk[++tot]=u,instk[u]=true;
    for (int i=hed[u];i!=-1;i=nxt[i])
    {
        int v=rea[i];
        if (!dfn[v])
        {
            Tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if (instk[v]) low[u]=min(low[u],dfn[v]);
    }
    if (low[u]==dfn[u])
    {
        int x=-1;scc++;
        while (x!=u)
        {
            x=stk[tot--];
            bel[x]=scc;
            instk[x]=false;
        }
    }
}
void Spfa(int st)
{
    queue<int>q;q.push(bel[st]);bool boo[N];
    memset(boo,0,sizeof(boo));boo[bel[st]]=true;
    for (int i=1;i<=scc;i++) dis[i]=0;
    dis[bel[st]]=val[bel[st]];
    while (!q.empty())
    {
        int u=q.front();q.pop();
        for (int i=hed1[u];i!=-1;i=nxt1[i])
        {
            int v=rea1[i];
            if (dis[u]+val[v]>dis[v])
            {
                dis[v]=dis[u]+val[v];
                if (!boo[v]) q.push(v),boo[v]=true;
            }
        }
        boo[u]=false;
    }
}
int main()
{
    freopen("fzy.in","r",stdin);
    freopen("fzy.out","w",stdout);
    
    memset(hed,-1,sizeof(hed));
    memset(hed1,-1,sizeof(hed1));
    n=read(),m=read();
    for (int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        add(x,y);
    }
    for (int i=1;i<=n;i++) w[i]=read();
    S=read(),p=read();
    for (int i=1;i<=p;i++)
        flag[read()]=true;
    for (int i=1;i<=n;i++)
        if (!dfn[i]) Tarjan(i);
    for (int u=1;u<=n;u++)
    {
        for (int i=hed[u];i!=-1;i=nxt[i])
        {
            int v=rea[i];
            if (bel[u]!=bel[v]) add1(bel[u],bel[v]);
        }
        flag1[bel[u]]|=flag[u];
        val[bel[u]]+=w[u];
    }
    Spfa(S);
    int ans=0;
    for (int i=1;i<=scc;i++)
        ans=max(ans,dis[i]*flag1[i]);
    printf("%d\n",ans);
}