Silver Cow Party

Silver Cow Party

给出一个n个点m条边的有向图，\(w[i][j]\)为i到j的边权，给出一个点X,询问从每个点出发到达X再回到起点的路径长度和最大值，\(n\leq 1000,m\leq 100000\)。

解

显然为单源最短路径的题目，考虑点为起点终点的取反，建反边，工具是bfs(边权为1),双端队列bfs(边权0,1),优先队列bfs(边权为正)。

显然从点X出发得到到达每一个点i的路径长度即该点i返回时路径长度，记做\(dis1[i]\),那么建反边以后，从X到达的点i的路径长度记做\(dis2[i]\)。

用优先队列bfs求出这两个数组后，取最大的\(dis1[i]+dis2[i]\)即可。

参考代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#define il inline
#define ri register
#define Size 1050
#define intmax 1061109567
using namespace std;
template<class free>
struct heap{
	free a[Size*Size];int n;
	il void push(free x){
		a[++n]=x;ri int p(n);
		while(p>1)
			if(a[p]<a[p>>1])
				swap(a[p],a[p>>1]),
					p>>=1;
			else break;
	}
	il void pop(){
		a[1]=a[n--];ri int p(1),s(2);
		while(s<=n){
			if(s<n&&a[s+1]<a[s])++s;
			if(a[s]<a[p])
				swap(a[s],a[p]),
					p=s,s<<=1;
			else break;
		}
	}
};
struct pi{
	int x,y;
	il bool operator<(const pi&a)const{
		return y<a.y;
	}
};
heap<pi>H;
struct point{
	point*next;int to,w;
}*head[Size],*Head[Size];
int dis1[Size],dis2[Size];
il void read(int&),link(int,int,int,point**),
					 dijistra(int,int*,point**);
int main(){
	int n,m,x;
	read(n),read(m),read(x);
	for(int i(1),u,v,w;i<=m;++i)
		read(u),read(v),read(w),
			link(u,v,w,head),
			link(v,u,w,Head);
	memset(dis1,-1,sizeof(dis1)),dijistra(x,dis1,head);
	memset(dis2,-1,sizeof(dis2)),dijistra(x,dis2,Head);
	int ans(-intmax);
	for(int i(1);i<=n;++i)
		ans=max(ans,dis1[i]+dis2[i]);
	printf("%d",ans);
	return 0;
}
il void dijistra(int s,int *dis,point**head){
	H.n=0,H.push({s,0});pi sour;
	while(H.n){
		sour=H.a[1],H.pop();
		if(dis[sour.x]>=0)continue;
		dis[sour.x]=sour.y;
		for(point*i(head[sour.x]);i!=NULL;i=i->next)
			if(dis[i->to]<0)
				H.push({i->to,i->w+sour.y});
	}
}
il void link(int u,int v,int w,point**head){
	head[u]=new point{head[u],v,w};
}
il void read(int &x){
	x^=x;ri char c;while(c=getchar(),c<'0'||c>'9');
	while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
}