BZOJ 1179[Apio2009]Atm

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

tarjan缩点，处理处每个点的点权，然后跑一个最长路即可。最后在所有能跑到的朋友家更新最大值。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<queue>
using namespace std;
const int maxn=500005;
struct dqs
{
int f,t;
}hh[maxn],ha[maxn];
int first[maxn],nxt[maxn],tot=0;
void build(int f,int t)
{
hh[++tot]=(dqs){f,t};
nxt[tot]=first[f];
first[f]=tot;
}
int tott=0,snum=0,cnt=0;
int dfn[maxn],low[maxn],stack[maxn],jlqlt[maxn];
bool in_stack[maxn];
void group(int x)
{
dfn[x]=low[x]=++tott;
stack[++snum]=x;
in_stack[x]=1;
for(int i=first[x];i;i=nxt[i])
{
int u=hh[i].t;
if(!dfn[u])
{
group(u);
low[x]=min(low[x],low[u]);
}
else if(in_stack[u])    low[x]=min(low[x],dfn[u]);
}
if(dfn[x]==low[x])
{
cnt++;
while(true)
{
jlqlt[stack[snum]]=cnt;
in_stack[stack[snum]]=0;
snum--;
if(stack[snum+1]==x) break;
}
}
}
int dian[maxn],cat[maxn],fir[maxn],nt[maxn],tot1,d[maxn];
bool used[maxn],fend[maxn],vis[maxn];
void add(int f,int t)
{
ha[++tot1]=(dqs){f,t};
nt[tot1]=fir[f];
fir[f]=tot1;
}
queue<int>q;
void spfa(int s)
{
q.push(s);
vis[s]=1;
while(!q.empty())
{
int x=q.front();
q.pop();
vis[x]=0;
for(int i=fir[x];i;i=nt[i])
{
int u=ha[i].t;
if(d[u]<d[x]+cat[u])
{
d[u]=d[x]+cat[u];
if(!vis[u])
{
q.push(u);
vis[u]=1;
}
}
}
}
}
int main()
{
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
build(x,y);
}
for(int i=1;i<=n;i++)   scanf("%d",&dian[i]);
for(int i=1;i<=n;i++)   if(!dfn[i]) group(i);
for(int i=1;i<=n;i++)
cat[jlqlt[i]]+=dian[i];
for(int i=1;i<=n;i++)
{
for(int j=first[i];j;j=nxt[j])
{
int u=hh[j].t;
if(jlqlt[i]!=jlqlt[u])  add(jlqlt[i],jlqlt[u]);
}
}
int s,tt;
scanf("%d%d",&s,&tt);
for(int i=1;i<=tt;i++)
{
int xx;
scanf("%d",&xx);
fend[jlqlt[xx]]=1;
}
for(int i=1;i<=cnt;i++) d[i]=cat[i];
spfa(jlqlt[s]);
int ans=0;
for(int i=1;i<=cnt;i++)
if(fend[i]==1&&d[i]>ans)
ans=max(ans,d[i]);
printf("%d\n",ans);
return 0;
}