【NOIP提高组】爬山

Description

国家一级爬山运动员h10今天获得了一张有着密密麻麻标记的地图，在好奇心的驱使下，他又踏上了去爬山的路。
对于爬山，h10有一个原则，那就是不走回头路，于是他把地图上的所有边都标记成了有向边。他决定从点S出发，每到达一个新的节点他就可以获得一定的成就值。同时h10又是一个很珍惜时间的运动员，他不希望这次爬山的成就值白白浪费，所以最后他一定要在一个存档点停下，保存自己的成就值。
请你计算出在此次爬山运动中h10能够得到的最大成就值。保证h10能走到存档点。

Solution

其实想要讲这道题并不是因为他有多难，而是这道题反映了我一些从未想到过的问题，主要在于自我反省。
首先的一个思想就是走最长路，但是对于一个环来说，无法用最短路算法实现，所以要用tarjan缩环（比赛的时候以为一条边只能走一次，所以压根没有想到缩环，主要是题目描述说不走回头路，让我以为只能走一遍）。改题的时候，发现自己的tarjan好像怎么都不对，回看我的其他tarjan，发现竟然没有一个是一样的……一脸懵逼！！！我竟然没有打过标准tarjan，赶紧学习。

Code

var
    last,next,data,t,u,v,zh,st,dfn,low,ch:array[0..500005] of longint;
    d:array[0..1500000] of longint;
    dis,a:array[0..500005] of int64;
    p:array[1..500000] of boolean;
    n,m,i,j,sum,tot,s,o:longint;
    ans:int64;
function max(x,y:int64):int64;
begin
    if x>y then exit(x) else exit(y);
end;
function min(x,y:longint):longint;
begin
    if xthen exit(x) else exit(y);
end;
procedure ins(x,y:longint);
begin
    inc(sum);
    data[sum]:=y;
    next[sum]:=last[x];
    last[x]:=sum;
end;
procedure tarjan(x:longint);
var i,k:longint;
begin
    zh[0]:=1;zh[1]:=x;st[1]:=last[x];
    while zh[0]<>0 do
    begin
        x:=zh[zh[0]];
        if dfn[x]=0 then
        begin
            inc(sum); dfn[x]:=sum;low[x]:=sum;
            p[x]:=true;inc(d[0]); d[d[0]]:=x;
        end;
        k:=st[zh[0]];
        while k<>0 do
        begin
            if dfn[data[k]]<>0 then
            begin
                if p[data[k]] then low[x]:=min(low[x],low[data[k]]);
            end
            else
            begin
                //tarjan(data[k]);
                //low[x]:=min(low[x],low[data[k]]);
                st[zh[0]]:=k;
                inc(zh[0]);
                zh[zh[0]]:=data[k];st[zh[0]]:=last[data[k]];
                break;
            end;
            k:=next[k];
        end;
        if k<>0 then continue;
        x:=zh[zh[0]];
        if dfn[x]=low[x] then
        begin
            inc(tot);
            while (d[0]<>0)and(dfn[d[d[0]]]<>low[d[d[0]]]) do
            begin
                p[d[d[0]]]:=false;
                ch[d[d[0]]]:=tot;
                dec(d[0]);
            end;
            if d[0]>0 then
            begin
                p[d[d[0]]]:=false;
                ch[d[d[0]]]:=tot;
                dec(d[0]);
            end;
        end;
        dec(zh[0]);
    end;
end;
procedure doit;
var x,k:longint;
begin
    i:=0;j:=1;d[1]:=s;dis[s]:=a[s];p[s]:=true;
    while ido
    begin
        inc(i);
        x:=d[i];k:=last[x];
        while k<>0 do
        begin
            if dis[x]+a[data[k]]>dis[data[k]] then
            begin
                dis[data[k]]:=dis[x]+a[data[k]];
                if p[data[k]]=false then
                begin
                    inc(j);
                    d[j]:=data[k];
                    p[data[k]]:=true;
                end;
            end;
            k:=next[k];
        end;
        p[x]:=false;
    end;
end;
procedure rebuild;
begin
    fillchar(last,sizeof(last),0);
    fillchar(next,sizeof(next),0);
    fillchar(data,sizeof(data),0);
    sum:=0;
    for i:=1 to n do a[ch[i]]:=a[ch[i]]+t[i];
    for i:=1 to m do
        if ch[u[i]]<>ch[v[i]] then ins(ch[u[i]],ch[v[i]]);
end;
begin
    readln(n,m);
    for i:=1 to m do
    begin
        readln(u[i],v[i]);
        ins(u[i],v[i]);
    end;
    for i:=1 to n do read(t[i]);
    readln(s,o);
    sum:=0;tot:=0;
    tarjan(s);
    rebuild;
    fillchar(p,sizeof(p),0);
    fillchar(dis,sizeof(dis),0);
    s:=ch[s];
    doit;
    for i:=1 to o do
    begin
        read(s);
        s:=ch[s];
        ans:=max(ans,dis[s]);
    end;
    writeln(ans);
end.