bzoj 1179: [Apio2009]Atm （spfa+tarjan）

1179: [Apio2009]Atm

Time Limit: 15 Sec   Memory Limit: 162 MB
Submit: 3026   Solved: 1268
[ 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的单向的道路到达其中的至少一个酒吧。

Source

[ Submit][ Status][ Discuss]

题解：tarjan+spfa

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define N 1000003
using namespace std;
int n,m,s,t;
int point[N],v[N],next[N],x[N],y[N],dis[N],can[N],val[N];
int belong[N],size[N],st[N],top,tot,cnt,sz,dfsn[N],ins[N],low[N];
void add(int x,int y)
{
	tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y;
	//cout<<x<<" "<<y<<endl;
}
void tarjan(int x)
{
	dfsn[x]=low[x]=++sz; st[++top]=x;
	ins[x]=1;
	for (int i=point[x];i;i=next[i])
	 if (!dfsn[v[i]]) {
	 	tarjan(v[i]);
	 	low[x]=min(low[x],low[v[i]]);
	 }else if (ins[v[i]])  low[x]=min(low[x],dfsn[v[i]]);
	if (dfsn[x]==low[x]) {
		++cnt; int j;
		do{
			j=st[top--];
			belong[j]=cnt;
			size[cnt]+=val[j]; ins[j]=0;
		}while (j!=x);
	}
}
void spfa()
{
	memset(dis,0,sizeof(dis));
	memset(can,0,sizeof(can));
	int t=belong[s];
    dis[t]=size[t]; can[t]=1;
    queue<int> p;
    p.push(t);
    while (!p.empty()) {
    	int now=p.front(); p.pop();
    	for (int i=point[now];i;i=next[i])
    	 if (dis[v[i]]<dis[now]+size[v[i]])
    	 {
    	 	dis[v[i]]=dis[now]+size[v[i]];
    	 	if (!can[v[i]]) {
    	 		can[v[i]]=1;
    	 		p.push(v[i]);
			 }
		 }
		can[now]=0;
	}
}
int main()
{
	freopen("a.in","r",stdin);
	//freopen("my.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=m;i++) {
		scanf("%d%d",&x[i],&y[i]);
		add(x[i],y[i]);
	}
	for (int i=1;i<=n;i++) scanf("%d",&val[i]);
	for (int i=1;i<=n;i++)
	 if (!dfsn[i]) tarjan(i);
	//cout<<endl;
	tot=0;
	memset(point,0,sizeof(point));
	memset(next,0,sizeof(next));
	for (int i=1;i<=m;i++)
	 if (belong[x[i]]!=belong[y[i]])   add(belong[x[i]],belong[y[i]]);
	scanf("%d%d",&s,&t);
	spfa();
	int ans=0;
	for (int i=1;i<=t;i++) {
		int x1; scanf("%d",&x1);
		ans=max(ans,dis[belong[x1]]);
	}
	printf("%d\n",ans);
}