【bzoj 1179】[Apio2009]Atm（Tarjan+spfa）

1179: [Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 3031  Solved: 1270
[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】

【由于本题中存在环，又是跑最长路，那么，只要遍历到环上的一个点，与它同属一个环的点也都会被遍历到。所以，把环缩成一个点处理，这个点的权值为环上所有点的权值之和。重新建图，跑spfa最长路即可】

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int a[500010],nxt[500010],p[500010],tot;
int d[500010],Nxt[500010],point[500010],tt;
int dis[500010],dft[500010],f[500010],size[500010],cnt,root;
int q[1000010],top,len[500010];
int n,m,s,t,num[500010],w[500010],val[500010],ans;
bool vis[500010];
inline void Add(int x,int y)
{
	tt++; d[tt]=y; Nxt[tt]=point[x]; point[x]=tt;
}
inline void add(int x,int y)
{
	tot++; a[tot]=y; nxt[tot]=p[x]; p[x]=tot;
}
void tarjan(int x)
{
	dft[x]=dis[x]=++cnt;
	q[++top]=x; vis[x]=1;
	for(int i=point[x];i!=-1;i=Nxt[i])
	 if(!dft[d[i]])
	  {
	  	tarjan(d[i]);
	  	dis[x]=min(dis[x],dis[d[i]]);
	  }
	  else
	   if(vis[d[i]]&&dis[x]>dft[d[i]]) dis[x]=dft[d[i]];
	if(dis[x]==dft[x])
	 {
	 	++root; int b;
	 	do{
	 		b=q[top--];
	 		size[root]++; 
	 		f[b]=root;
	 		vis[b]=0;
	 		val[root]+=w[b];
		 }while(b!=x);
	 }
}
inline void spfa(int s)
{
	memset(len,0,sizeof(len));
	memset(vis,0,sizeof(vis));
	queue<int>que;
	len[s]=val[s]; vis[s]=1; que.push(s);
	while(!que.empty())
	 {
	 	int u=que.front(); que.pop();
	 	vis[u]=0;
	 	for(int i=p[u];i!=-1;i=nxt[i])
	 	 if(len[a[i]]<len[u]+val[a[i]])
	 	  {
	 	  	len[a[i]]=len[u]+val[a[i]];
	 	  	if(!vis[a[i]])
	 	  	 {
	 	  	 	vis[a[i]]=1; 
	 	  	 	que.push(a[i]);
				}
		   }
	 }
	return;
}
int main()
{
	int i,j;
	memset(p,-1,sizeof(p));
	memset(nxt,-1,sizeof(nxt));
	memset(Nxt,-1,sizeof(Nxt));
	memset(point,-1,sizeof(point));
	scanf("%d%d",&n,&m);
	for(i=1;i<=m;++i)
	 {
	 	int x,y;
	 	scanf("%d%d",&x,&y);
	 	Add(x,y);
	 }
	for(i=1;i<=n;++i)scanf("%d",&w[i]);
	scanf("%d%d",&s,&t);
	for(i=1;i<=t;++i)scanf("%d",&num[i]);
	for(i=1;i<=n;++i)
	 if(!dft[i]) tarjan(i);
	for(i=1;i<=n;++i)
	 for(j=point[i];j!=-1;j=Nxt[j])
	  if(f[i]!=f[d[j]]) add(f[i],f[d[j]]);
	spfa(f[s]);
	for(i=1;i<=t;++i)
	 if(ans<len[f[num[i]]]) ans=len[f[num[i]]];
	printf("%d\n",ans);
	return 0;
}