[cf1473e]Minimum Path(分层图最短路)

是个分层图的问题。感觉知道的话似乎很容易就想到了,不知道的话……大概就不知道了qwq。

分层图就是说,如果在一个带权图上,我们对于每条路径都可以将其中k条边做一些改变,比如不计算其权值,那么我们就先建k+1层完全相同的初始图(权值为原值),然后层之间的连边情况不变,但权值为改变之后的。

比如如果图是1、2、3两两之间连一条权值为1的边,一条路径上可以有一条边不计算权值,求1到各点最短路。那么我们建图分两层,一层的1、2、3两两间连边权值为1,二层的1'、2'、3'两两间连边权值为1,然后1->2'=0,1->3'=0,2->1'=0,2->3'=0,3->1'=0,3->2'=0。然后以1为起点跑最短路,得到2'和3'dis即为答案。这样保证如果最后走到了第k+1层,一定进行过k次边权变化,且只要不走环,不会有相同的边变化了多次。

然后我们看这道题,他其实是相当于进行两次变化,一条边权值计为0,一条边权值*2。如何保证是最大边不计、最小边二倍呢?稍微思考一下,如果我们想要得到最短路,且一定经过这两次边权变化,那么对于同一条路径,不计最大值而计两边最小值必然是最优的结果。然后我们又发现一个问题,这两个变化类型不同,放在第1、2层之间的那个必然在2、3层之间的那个之前进行。所以我们要建两次图,一次假设最小边在最大边之前,另一次假设在后,得到的两个结果的最小值就是答案。还有就是,不要忘记考虑某些节点和1直接连的边就是1到它的最短路。

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define ll long long
using namespace std;
const int N=200010;
struct edge{
	int y,next,w;
}data[N*10];
struct node{
	int id;
	ll dis;
	bool vis;
	bool operator < (const node& x) const{
		return x.dis<dis;
	}
}a1[N*3],x;
int n,m,num,nn,h[N*3];
ll ans[N];
priority_queue<node,vector<node>,less<node> > q;
inline int read(){
	int x=0;char ch=getchar();
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x;
}
inline void addedge(int u,int v,int w){
	data[++num].y=v;data[num].w=w;data[num].next=h[u],h[u]=num;
}
inline ll min1(ll i,ll j){return i>j?j:i;}
int main(){
	n=read(),m=read();
	for(int i=2;i<=n;++i)ans[i]=0x3f3f3f3f3f3f3f3f;
	memset(h,-1,sizeof h);num=0;
	for(int i=1,u,v,w;i<=m;++i){
		u=read(),v=read(),w=read();
		if(u==1)ans[v]=min1(ans[v],w);
		if(v==1)ans[u]=min1(ans[u],w);
		addedge(u,v,w),addedge(v,u,w),
		addedge(u+n,v+n,w),addedge(v+n,u+n,w),
		addedge(u+n+n,v+n+n,w),addedge(v+n+n,u+n+n,w),
		addedge(u,v+n,0),addedge(u+n,v+n+n,w<<1),
		addedge(v,u+n,0),addedge(v+n,u+n+n,w<<1);
	}
	nn=n*3;
	a1[1].id=1,a1[1].dis=0,a1[1].vis=false;
	for(int i=2;i<=nn;++i)a1[i].id=i,a1[i].vis=false,a1[i].dis=0x3f3f3f3f3f3f3f3f;
	while(!q.empty())q.pop();
	q.push(a1[1]);
	while(!q.empty()){
		x=q.top();q.pop();
		a1[x.id].vis=false;
		for(int i=h[x.id],v;i!=-1;i=data[i].next){
			v=data[i].y;
			if(a1[x.id].dis+1ll*data[i].w<a1[v].dis){
				a1[v].dis=a1[x.id].dis+1ll*data[i].w;
				if(!a1[v].vis){q.push(a1[v]);a1[v].vis=true;}
			}
		}
	} 
	for(int i=2;i<=n;++i)ans[i]=min1(ans[i],a1[n*2+i].dis);
	for(int i=7;i<=num;i+=10){swap(data[i].w,data[i+1].w);swap(data[i+2].w,data[i+3].w);}
	a1[1].id=1,a1[1].dis=0,a1[1].vis=false;
	for(int i=2;i<=nn;++i)a1[i].id=i,a1[i].vis=false,a1[i].dis=0x3f3f3f3f3f3f3f3f;
	while(!q.empty())q.pop();
	q.push(a1[1]);
	while(!q.empty()){
		x=q.top();q.pop();
		a1[x.id].vis=false;
		for(int i=h[x.id],v;i!=-1;i=data[i].next){
			v=data[i].y;
			if(a1[x.id].dis+1ll*data[i].w<a1[v].dis){
				a1[v].dis=a1[x.id].dis+1ll*data[i].w;
				if(!a1[v].vis){q.push(a1[v]);a1[v].vis=true;}
			}
		}
	} 
	for(int i=2;i<=n;++i)ans[i]=min1(ans[i],a1[n*2+i].dis);
	for(int i=2;i<=n;++i)printf("%lld ",ans[i]);
	return 0;
} 

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值