BZOJ 1179 Apio2009 Atm

1179: [Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 4059  Solved: 1771
[Submit][Status][Discuss]

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的单向的道路到达其中的至少一个酒吧。

 

 

tarjin缩完点后spfa跑一下最长路即可

#include <bits/stdc++.h>
#define ll long long
#define inf 1000000
#define eps 1e-7
using namespace std;
inline int read(){
    int x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int MAXN=1e6+10;
struct node{
    int y,next;
}e[MAXN];
int linkk[MAXN],s,p,len=0,n,head,tail,m,x[MAXN],q[MAXN],y[MAXN],a[MAXN],stark[MAXN],top,ine[MAXN],tot,dfn[MAXN],low[MAXN],dfs_clock,sum[MAXN],dis[MAXN],vis[MAXN];
inline void insert(int xx,int yy){
    e[++len].y=yy;e[len].next=linkk[xx];linkk[xx]=len;
}
void spfa(){
    memset(vis,0,sizeof(vis));
    dis[ine[s]]=sum[ine[s]];vis[ine[s]]=0;
    q[++tail]=ine[s];
    while(head!=tail){
        if(head==MAXN-1) head=0;
        int tn=q[++head];
        for(int i=linkk[tn];i;i=e[i].next){
            if(dis[e[i].y]<dis[tn]+sum[e[i].y]){
                dis[e[i].y]=dis[tn]+sum[e[i].y];
                if(!vis[e[i].y]){
                    if(tail==MAXN-1) tail=0;
                    q[++tail]=e[i].y;vis[e[i].y]=1;
                }
            }
        }
        vis[tn]=0;
    }
}
void rebuild(){
    memset(linkk,0,sizeof(linkk));len=0;
    for(int i=1;i<=m;i++){
        if(ine[x[i]]!=ine[y[i]]){
            insert(ine[x[i]],ine[y[i]]);
        }
    }
}
inline void tarjin(int st){
    dfn[st]=low[st]=++dfs_clock;
    stark[++top]=st;vis[st]=1;
    for(int i=linkk[st];i;i=e[i].next){
        if(!dfn[e[i].y]){
            tarjin(e[i].y);
            low[st]=min(low[st],low[e[i].y]);
        }
        else if(vis[e[i].y]) low[st]=min(low[st],dfn[e[i].y]);
    }
    if(low[st]==dfn[st]){
        int k;tot++;
        do{
            k=stark[top--];
            ine[k]=tot;
            vis[k]=0;
            sum[tot]+=a[k];
        }while(k!=st);
    }
}
int main(){
    n=read();m=read();
    for(int i=1;i<=m;i++){
        x[i]=read();y[i]=read();
        insert(x[i],y[i]);
    }
    for(int i=1;i<=n;i++){
        a[i]=read();
    }
    for(int i=1;i<=n;i++){
        if(!dfn[i]) tarjin(i);
    }
    s=read();p=read();
    rebuild();
    spfa();
    int mn=-inf;
    for(int i=1;i<=p;i++){
        int xx=read();
        mn=max(mn,dis[ine[xx]]);
    }
    cout<<mn<<endl;
    return 0;
}