【JZOJ 4771】爬山

Description

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

Solution

很明显，做这题的第一件事就是用tarjan缩环，
缩完环以后，就直接从S点跑SPFA最大代价即可，

注意：本题用递归版tarjan会爆栈！！！

复杂度：$O(n+m)$

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define efo(i,q) for(int i=A[q];i;i=B[i][0])
using namespace std;
typedef long long LL;
const int N=500500,maxlongint=2147483640;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n,S;
LL ans,b[N],d1[N];
int A[N],B[2*N][2],B0,HE[N];
int g[N];
int dfn[N],low[N],za[N];
bool z1[N];
int d[N],z[N];
int ZA[N][2];
void link(int q,int w){if(!A[q])HE[q]=B0+1;B[++B0][0]=A[q],A[q]=B0,B[B0][1]=w;}
void tarjan(int q)
{
    int I=1;bool OK;
    ZA[I][0]=q;ZA[I][1]=0;
    while(I)
    {
        OK=1;q=ZA[I][0];
        if(!ZA[I][1])z[q]=z1[q]=1,za[dfn[q]=low[q]=++za[0]]=q,ZA[I][1]=A[q];
            else low[q]=min(low[q],low[B[ZA[I][1]][1]]),ZA[I][1]=B[ZA[I][1]][0];
        for(int i=ZA[I][1];i;i=B[i][0])
            if(z[B[i][1]])low[q]=min(low[q],low[B[i][1]]);
            else if(!z1[B[i][1]])
            {
                ZA[I][1]=i;
                ZA[++I][0]=B[i][1];
                ZA[I][1]=0;
                OK=0;break;
            }
        if(OK)
        {
            if(low[q]==dfn[q])
            {
                int t;
                while((t=za[za[0]--])!=q)z[t]=0,g[t]=q;
                g[q]=q;z[q]=0;
            }
            I--;
        }
    }
}
int gf(int q){return g[q]=(g[q]==q?q:gf(g[q]));}
int main()
{
    int q,w;
    read(n),read(m);
    fo(i,1,m)read(q),read(w),link(q,w);
    fo(i,1,n)b[i]=read(q);
    fo(i,1,n)if(!z1[i])tarjan(i);
    fo(i,1,n)if((q=gf(i))!=i)
    {
        b[q]+=b[i];b[i]=0;
        B[HE[i]][0]=A[q];A[q]=A[i];A[i]=0;
    }
    fo(i,1,n)efo(j,i)B[j][1]=g[B[j][1]];
    S=g[read(S)],read(q);
    fo(i,1,q)read(w),z1[w]=0;
    memset(z,0,sizeof(z));
    int i=1,j=1;d[i]=S;z[S]=1;
    while(i<=j)
    {
        q=d[i++];
        efo(k,q)if(d1[B[k][1]]<d1[q]+b[q]&&B[k][1]!=q)
        {
            d1[B[k][1]]=d1[q]+b[q];
            if(!z[B[k][1]])d[++j]=B[k][1],z[B[k][1]]=1;
        }
        z[q]=0;
    }
    fo(i,1,n)if(!z1[i])ans=max(ans,b[g[i]]+d1[g[i]]);
    printf("%lld\n",ans);
    return 0;
}