BZOJ 1179[Apio2009]Atm

Description

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号

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

50%的输入保证N, M<=3000。所有的输入保证N,
M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

tarjan缩点，处理处每个点的点权，然后跑一个最长路即可。最后在所有能跑到的朋友家更新最大值。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<queue>
using namespace std;
const int maxn=500005;
struct dqs
{
    int f,t;
}hh[maxn],ha[maxn];
int first[maxn],nxt[maxn],tot=0;
void build(int f,int t)
{
    hh[++tot]=(dqs){f,t};
    nxt[tot]=first[f];
    first[f]=tot;
}
int tott=0,snum=0,cnt=0;
int dfn[maxn],low[maxn],stack[maxn],jlqlt[maxn];
bool in_stack[maxn];
void group(int x)
{
    dfn[x]=low[x]=++tott;
    stack[++snum]=x;
    in_stack[x]=1;
    for(int i=first[x];i;i=nxt[i])
    {
        int u=hh[i].t;
        if(!dfn[u])
        {
            group(u);
            low[x]=min(low[x],low[u]);
        }
        else if(in_stack[u])    low[x]=min(low[x],dfn[u]);
    }
    if(dfn[x]==low[x])
    {
        cnt++;
        while(true)
        {
            jlqlt[stack[snum]]=cnt;
            in_stack[stack[snum]]=0;
            snum--;
            if(stack[snum+1]==x) break;
        }
    }
}
int dian[maxn],cat[maxn],fir[maxn],nt[maxn],tot1,d[maxn];
bool used[maxn],fend[maxn],vis[maxn];
void add(int f,int t)
{
    ha[++tot1]=(dqs){f,t};
    nt[tot1]=fir[f];
    fir[f]=tot1;
}
queue<int>q;
void spfa(int s)
{
    q.push(s);
    vis[s]=1;
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        vis[x]=0;
        for(int i=fir[x];i;i=nt[i])
        {
            int u=ha[i].t;
            if(d[u]<d[x]+cat[u])
            {
                d[u]=d[x]+cat[u];
                if(!vis[u])
                {
                    q.push(u);
                    vis[u]=1;
                }
            }
        }
    }
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        build(x,y);
    }
    for(int i=1;i<=n;i++)   scanf("%d",&dian[i]);
    for(int i=1;i<=n;i++)   if(!dfn[i]) group(i);
    for(int i=1;i<=n;i++)
        cat[jlqlt[i]]+=dian[i];
    for(int i=1;i<=n;i++)
    {
        for(int j=first[i];j;j=nxt[j])
        {
            int u=hh[j].t;
            if(jlqlt[i]!=jlqlt[u])  add(jlqlt[i],jlqlt[u]); 
        }
    }
    int s,tt;
    scanf("%d%d",&s,&tt);
    for(int i=1;i<=tt;i++)
    {
        int xx;
        scanf("%d",&xx);
        fend[jlqlt[xx]]=1;
    }
    for(int i=1;i<=cnt;i++) d[i]=cat[i];
    spfa(jlqlt[s]);
    int ans=0;
    for(int i=1;i<=cnt;i++)
        if(fend[i]==1&&d[i]>ans)
            ans=max(ans,d[i]);
    printf("%d\n",ans);
    return 0;
}