JZOJ 4771. 爬山

Problem

Description

国家一级爬山运动员h10今天获得了一张有着密密麻麻标记的地图，在好奇心的驱使下，他又踏上了去爬山的路。

对于爬山，h10有一个原则，那就是不走回头路，于是他把地图上的所有边都标记成了有向边。他决定从点S出发，每到达一个新的节点他就可以获得一定的成就值。同时h10又是一个很珍惜时间的运动员，他不希望这次爬山的成就值白白浪费，所以最后他一定要在一个存档点停下，保存自己的成就值。

请你计算出在此次爬山运动中h10能够得到的最大成就值。保证h10能走到存档点。

Input

第一行两个整数 N，M，表示点数和边数。

接下来 M 行，每行两个整数 u，v，表示u到v有一条有向边(没有自环)。

第 M+2 行 N 个正整数，表示每个点的成就值。

接下来一行两个整数 S，p，表示出发点和存档点个数。

下面一行 p 个整数，表示存档点。

Output

一个正整数，表示最大成就值。

Sample Input

5 7

5 1

3 1

2 5

3 5

4 3

4 2

4 5

7 6 3 2 2

4 3

1 5 2

Sample Output

17

Data Constraint

对于 30% 的数据， N,M≤1000，并且地图为有向无环图。

对于 100% 的数据， N,M≤500000。（数据有梯度，注意答案的大小）

Solution

我们用Tarjan找到所有强连通分量，并将每个强连通分量缩成一个点，因为强连通分量中所有的分数都可以获得。

然后再跑一遍Spfa最长路，然后取最高的存档点的值即可。

然而，链型数据可能会卡爆你的Tarjan!

所以我们要打人工栈。人工栈的原理与递归的原理相同，不过就是要再开一个数组。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define N 500010
#define LL long long
using namespace std;
int n,m,i,j,h,t,w,top,cnt,x,y,p,S,now,tot;
int g
,low
,dfn
,st
,head
,next
,go
,A
;
LL a
,d
,dis
,f
,ans;
bool b
,c
,e
,bz
;
inline int read()
{
int data=0;
char ch=0;
while (ch<'0' || ch>'9') ch=getchar();
while (ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();
return data;
}
void lb(int x,int y)
{
go[++tot]=y;
next[tot]=head[x];
head[x]=tot;
}
void tarjan(int x)
{
int T=1,i,now;bool o;
A[T]=x;
while (T)
{
if (dfn[A[T]]==-1)
{
h++;
dfn[A[T]]=low[A[T]]=h;
bz[A[T]]=1;
st[++top]=A[T];
}
o=0;
for (i=head[A[T]];i;i=next[i])
{
now=go[i];
if (dfn[now]==-1)
{
o=1;
break;
} else
if (bz[now]) low[A[T]]=min(low[A[T]],dfn[now]);
}
if (o) A[++T]=now;else
{
if (dfn[A[T]]==low[A[T]])
{
cnt++;
while (st[top+1]!=A[T])
{
bz[st[top]]=0;
g[st[top]]=cnt;
if (c[st[top]]) e[cnt]=1;
f[cnt]+=a[st[top]];
top--;
}
}
if (T>1) low[A[T-1]]=min(low[A[T-1]],low[A[T]]);
T--;
}
}
}
int main()
{
memset(dfn,255,sizeof(dfn));
n=read();m=read();
fo(i,1,m) x=read(),y=read(),lb(x,y);
fo(i,1,n) a[i]=read();
scanf("%d%d",&S,&p);
fo(i,1,p) x=read(),c[x]=1;
tot=0;
tarjan(S);
t=0,w=1;
d[1]=S;
dis[S]=f[g[S]];
while (t<w)
{
t++;
now=d[t];
for(i=head[now];i;i=next[i])
{
if (g[go[i]]!=g[now])
{
if (dis[now]+f[g[go[i]]]>dis[go[i]])
{
dis[go[i]]=dis[now]+f[g[go[i]]];
if (!b[go[i]])
{
b[go[i]]=1;
w++;
d[w]=go[i];
}
}
} else
{
if (dis[now]>dis[go[i]])
{
dis[go[i]]=dis[now];
if (!b[go[i]])
{
b[go[i]]=1;
w++;
d[w]=go[i];
}
}
}
}
b[now]=0;
}
fo(i,1,n) if (c[i]) ans=max(ans,dis[i]);
printf("%lld",ans);
}