【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;}