洛谷P2934 [USACO09JAN]安全出行Safe Travel

首先读懂题意，这道题目的意思是
对于每一个点$x$，断去$1->x$最短路径上的最后一条边，求此时的最短路径
感觉一下子除了暴力没有思路
那么我们反向思考一下，一条不在最短路径上的边会有什么贡献呢
我们假设有一颗最短路径树，那么这一条不在最短路径上的边，一定连接了两个不为父子关系的节点

假设有一颗树是这样的
虚线为一条非树边，我们断去u上方的边，就变成了了这样

$$di{s}_u=min(di{s}_u,di{s}_v+w)$$
可以更新一次答案了，但是对于每一个点我们如果都遍历所有边的话，仍然会超时，我们考虑一次性求完所有贡献，能求出多少呢，$u$到$LCA(u,v)$,如何实现呢，存一个父亲即可
然而还是会T，因为每条边都这么求的话重复遍历也会爆炸，我们再考虑，如何跳过哪些已经求出最优值的点，我们可以在跳的过程中进行并查集缩点，当然，这时候我们的边就需要提前排序了
这里通过骚操作实现，我们存储w的权值变为$di{s}_u+di{s}_v+w$，这样我们更新的时候用新权值减去更新点的深度即可，而且因为减去的距离是一个定值，我们可以直接排序，保证先遍历的边尽量最优

并且对于$u$的祖先，上式仍然成立，因为$di{s}_u-di{s}_{find(u)}$实际上就是祖先到u的距离，再加上$di{s}_v+w$，就是新路径的长度，这样我们就完成了一个类似于最小生成树式的更新

讲的可能不是很清楚，但是画个图对着代码就能理解了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 2e5+7;
const int INF = 2147483647;
struct nod
{
	int from,to,w;
}way[maxn*2];
struct node
{
	int to,next,w;
}edge[maxn*2];
int cnt,head[maxn];
void add(int from,int to,int w)
{
	edge[++cnt].to=to;
	edge[cnt].w=w;
	edge[cnt].next=head[from];
	head[from]=cnt;
}
int n,m,f[maxn];
int find(int x){return f[x]==x?x:f[x]=find(f[x]);}
int ans[maxn],dis[maxn],fa[maxn];
void dij()
{
	priority_queue< pair<int,int> >q;
	for(int i=1;i<=n;i++)dis[i]=INF;
	dis[1]=0;
	q.push(make_pair(-dis[1],1));
	while(!q.empty())
	{
		int f1=-q.top().first;
		int f2=q.top().second;
		q.pop();
		if(f1!=dis[f2])continue;
		for(int i=head[f2];i;i=edge[i].next)
		{
			int to=edge[i].to;
			int w=edge[i].w;
			if(f1+w<dis[to])
			{
				dis[to]=f1+w;
				fa[to]=f2;
				q.push(make_pair(-dis[to],to));
			}
		}
	}
}

bool cmp1(nod a,nod b){return a.w<b.w;}

int num;
void work(nod a)
{
	int u=a.from,v=a.to,w=a.w;
	while(find(u)!=find(v))
	{
		num++;
		if(dis[find(u)]<dis[find(v)])swap(u,v);
		ans[find(u)]=min(ans[find(u)],w-dis[find(u)]);
		u=f[find(u)]=fa[find(u)];
	}
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1,x,y,z;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z);
		add(y,x,z);
	}
	dij();
	int tot=0;
	for(int i=2;i<=cnt;i+=2)
	{
		int u=edge[i].to,v=edge[i-1].to;
		if(fa[u]!=v&&fa[v]!=u)
		{
			way[++tot].from=u;
			way[tot].to=v;
			way[tot].w=dis[u]+dis[v]+edge[i].w;
		}
	}
	for(int i=1;i<=n;i++)f[i]=i,ans[i]=INF;
    sort(way+1,way+1+tot,cmp1);
	for(int i=1;i<=tot&&num<n-1;i++)
		work(way[i]);
	int no=-1;
	for(int i=2;i<=n;i++)
	{
		printf("%d\n",ans[i]==INF?no:ans[i]);
	}
}