题面戳此处:https://www.lydsy.com/JudgeOnline/problem.php?id=1179
题意大概就是给一个图,其中每个点都有固定的权值,求从一个起点到固定的几个终点能够获得的最大收获。点和边都可以重复经过,但是一个点的权值只能在第一次经过它时被获取一次。
因为是有向图,首先想到把spfa最短路算法瞎搞成最长路,但是和求最短路时不能有负权环一样,用spfa求最长路也得先把正权环整掉。这就必须用到tarjan将这个有向图的每个强连通分量都缩成点,重新构图。(听起来十分复杂其实实现起来很简单,毕竟tarjan是个伟人,tarjan算法也是个优雅而实用的算法)。
首先构图,正常tarjan,dfn是时间戳,表示图中此点是第几个被遍历到的(其实用意只在于给一个乱七八糟的图一个序,使得先遍历到的点和的序比后遍历到的点的序小),low[u]在程序中表示u在栈中最早的祖先的dfn值(听起来难以理解是记这个是用来干什么的?),其实low[u]在图中的意义就是如果存在从u指向u的祖先(比u早进栈而且遍历到u时还没有被弹出栈的节点)的一条边,那么low[u]将会记录这些祖先的dfn值的最小的那个(实现就是在遍历到还在栈中的点时将自己的low值和这位祖先的dfn值取min,然后千万别忘了递归回溯时不断将low值和自己的儿子的low值取min,因为自己儿子能到达的点自己也能通过儿子到达)。
事实上,这个在上述过程中具有最小dfn值的祖先就是u所在的强连通分量的“根”,最后将这一强连通分量找出来的线索就是这个根(把根提出来,底下根连的节点也都被提出来了,实现时的操作是将栈从栈顶将根都弹出栈,因为此时从栈顶元素到根都是一个强连通分量的了)。
那么,怎么找到一个强连通分量的根呢?那就很简单了啊,当一个节点的dfn值等于他的low值时,它就是它所在的强连通分量的根喽。
找到了强连通分量后,我们要为之后的缩点做准备,那就是把一个强联通分量中的点都赋予同一个新节点标号,并累加这个分量中的所有点的权值在新节点的身上。
接着,我们就可以快乐地缩点构图了!这里很好理解,详见下面代码中的gt函数。
spfa将 < 改为 >,走起。
最后记得查询的是所有酒馆缩完点之后的节点标号的dis,还有要注意的就是原图的节点标号等一系列东西不要和新建的图的东西搞混了!
代码:
#include<cstdio> #include<iostream> #include<cstring> using namespace std; #define maxn 500005 int n,m,g[maxn],tmp,s,val[maxn]; int sta[maxn],col[maxn],top; int p,pub[maxn],gr[maxn],nn,ral[maxn]; int dfn[maxn],low[maxn],tim,now; bool v[maxn]; int dis[maxn],q[maxn*4]; int ans; struct la{ int c,to; }e[maxn],ede[maxn]; inline void tu(int x,int y){e[++tmp].c=y; e[tmp].to=g[x]; g[x]=tmp;} inline void tu2(int x,int y){ede[++tmp].c=y; ede[tmp].to=gr[x]; gr[x]=tmp;} inline int mmin(int x,int y){return x<y?x:y;} inline int mmax(int x,int y){return x>y?x:y;} void tarjan(int u){ dfn[u]=low[u]=++tim; v[u]=1; sta[++top]=u; for(int i=g[u];i;i=e[i].to){ int d=e[i].c; if(!dfn[d]){ tarjan(d); low[u]=mmin(low[u],low[d]); } else if(v[d]){//still in the sta[] //which mean it is not in a known block low[u]=mmin(low[u],dfn[d]); } } if(dfn[u]==low[u]){// it is the root of this block nn++; while(sta[top]!=u){ now=sta[top]; top--; ral[nn]+=val[now]; col[now]=nn; v[now]=0; } ral[nn]+=val[u]; col[u]=nn; v[u]=0; top--; } } inline void gt(){ tmp=0; for(int i=1;i<=n;i++){ for(int j=g[i];j;j=e[j].to){ int d=e[j].c; if(col[d]!=col[i]){ tu2(col[i],col[d]); } } } } inline void spfa(int s){ memset(dis,0,sizeof(dis)); int h=1,r=1; memset(v,0,sizeof(v)); q[h]=col[s];v[col[s]]=1; dis[col[s]]=ral[col[s]]; while(h<=r){ int d=q[h];v[d]=0;h++; for(int i=gr[d];i;i=ede[i].to){ int u=ede[i].c; if(ral[u]+dis[d]>dis[u]){ dis[u]=dis[d]+ral[u]; if(!v[u]){ q[++r]=u; v[u]=1; } } } } } int main(){ scanf("%d%d",&n,&m); for(int i=1,a,b;i<=m;i++){ scanf("%d%d",&a,&b); tu(a,b); } for(int i=1;i<=n;i++){ scanf("%d",&val[i]); } scanf("%d%d",&s,&p); for(int i=1;i<=p;i++){ scanf("%d",&pub[i]); } tarjan(s); gt(); spfa(s); for(int i=1;i<=p;i++){ ans=mmax(ans,dis[col[pub[i]]]); } printf("%d",ans); }如对代码或者讲解有疑问,欢迎私信或评论交流!