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