1179: [Apio2009]Atm
Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 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
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跑一下最长路即可
1 #include <bits/stdc++.h> 2 #define ll long long 3 #define inf 1000000 4 #define eps 1e-7 5 using namespace std; 6 inline int read(){ 7 int x=0;int f=1;char ch=getchar(); 8 while(!isdigit(ch)) {if(ch==‘-‘) f=-1;ch=getchar();} 9 while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();} 10 return x*f; 11 } 12 const int MAXN=1e6+10; 13 struct node{ 14 int y,next; 15 }e[MAXN]; 16 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]; 17 inline void insert(int xx,int yy){ 18 e[++len].y=yy;e[len].next=linkk[xx];linkk[xx]=len; 19 } 20 void spfa(){ 21 memset(vis,0,sizeof(vis)); 22 dis[ine[s]]=sum[ine[s]];vis[ine[s]]=0; 23 q[++tail]=ine[s]; 24 while(head!=tail){ 25 if(head==MAXN-1) head=0; 26 int tn=q[++head]; 27 for(int i=linkk[tn];i;i=e[i].next){ 28 if(dis[e[i].y]<dis[tn]+sum[e[i].y]){ 29 dis[e[i].y]=dis[tn]+sum[e[i].y]; 30 if(!vis[e[i].y]){ 31 if(tail==MAXN-1) tail=0; 32 q[++tail]=e[i].y;vis[e[i].y]=1; 33 } 34 } 35 } 36 vis[tn]=0; 37 } 38 } 39 void rebuild(){ 40 memset(linkk,0,sizeof(linkk));len=0; 41 for(int i=1;i<=m;i++){ 42 if(ine[x[i]]!=ine[y[i]]){ 43 insert(ine[x[i]],ine[y[i]]); 44 } 45 } 46 } 47 inline void tarjin(int st){ 48 dfn[st]=low[st]=++dfs_clock; 49 stark[++top]=st;vis[st]=1; 50 for(int i=linkk[st];i;i=e[i].next){ 51 if(!dfn[e[i].y]){ 52 tarjin(e[i].y); 53 low[st]=min(low[st],low[e[i].y]); 54 } 55 else if(vis[e[i].y]) low[st]=min(low[st],dfn[e[i].y]); 56 } 57 if(low[st]==dfn[st]){ 58 int k;tot++; 59 do{ 60 k=stark[top--]; 61 ine[k]=tot; 62 vis[k]=0; 63 sum[tot]+=a[k]; 64 }while(k!=st); 65 } 66 } 67 int main(){ 68 n=read();m=read(); 69 for(int i=1;i<=m;i++){ 70 x[i]=read();y[i]=read(); 71 insert(x[i],y[i]); 72 } 73 for(int i=1;i<=n;i++){ 74 a[i]=read(); 75 } 76 for(int i=1;i<=n;i++){ 77 if(!dfn[i]) tarjin(i); 78 } 79 s=read();p=read(); 80 rebuild(); 81 spfa(); 82 int mn=-inf; 83 for(int i=1;i<=p;i++){ 84 int xx=read(); 85 mn=max(mn,dis[ine[xx]]); 86 } 87 cout<<mn<<endl; 88 return 0; 89 }