[CF786B]Legacy

Legacy

题解

看到题目应该是容易想到最短路的,但是由于区间到单点与单点到区间的边我们需要想些办法来维护。

于是,我们就想到了通过虚点来进行维护。但总不能每个区间都建一个点,只能利用线段树来对各个区间的虚点进行维护。

针对一个区间需要建两个点,一个入点,一个出点。将大区间与小区间连上边,出边从下往上,入边从上往下,其边权都是0。

由于大小区间内已经连了边,再加上最短路的性质,懒标记什么的就没必要了。最后建完边后再去跑一边Dijkstra即可,注意需要加堆优化。

由于n,q\leq 10^5,边数m应该不超过(n+q)log_{n},总时间复杂度O\left( (n+q)\log^2 n\right ),能卡过去。

源码

#include<bits/stdc++.h>
using namespace std;
#define MAXN 100005
typedef long long LL;
typedef unsigned long long uLL;
const LL INF=0x7f7f7f7f7f7f;
const int mo=1e9+7;
typedef pair<int,LL> pii;
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
int n,m,s,tot,head[MAXN],cnt;
int trin[MAXN<<2],trout[MAXN<<2];
LL dis[MAXN*10];
bool vis[MAXN*10];
vector<pii>e[MAXN*10];
struct ming{
	int id;LL t;
	bool friend operator < (const ming &x,const ming &y){
		return x.t>y.t;
	}
};
priority_queue<ming> q;
void build(int rt,int l,int r){
	if(l==r)return (void)(trin[rt]=trout[rt]=l);
	int mid=l+r>>1;
	build(rt<<1,l,mid);
	build(rt<<1|1,mid+1,r);
	trout[rt]=++cnt;trin[rt]=++cnt;
	e[trout[rt<<1]].push_back(make_pair(trout[rt],0));
	e[trout[rt<<1|1]].push_back(make_pair(trout[rt],0));
	e[trin[rt]].push_back(make_pair(trin[rt<<1],0));
	e[trin[rt]].push_back(make_pair(trin[rt<<1|1],0));
}
void updatain(int rt,int l,int r,int al,int ar,int au,int aw){
	if(al<=l&&r<=ar){e[au].push_back(make_pair(trin[rt],aw));return ;}
	int mid=l+r>>1;
	if(al<=mid)updatain(rt<<1,l,mid,al,ar,au,aw);
	if(ar>mid)updatain(rt<<1|1,mid+1,r,al,ar,au,aw);
}
void updataout(int rt,int l,int r,int al,int ar,int au,int aw){
	if(al<=l&&r<=ar){e[trout[rt]].push_back(make_pair(au,aw));return ;}
	int mid=l+r>>1;
	if(al<=mid)updataout(rt<<1,l,mid,al,ar,au,aw);
	if(ar>mid)updataout(rt<<1|1,mid+1,r,al,ar,au,aw);
}
void dijkstra(){
	memset(dis,0x7f,sizeof(dis));
	dis[s]=0;q.push((ming){s,0});
	while(!q.empty()){
		ming t=q.top();q.pop();if(vis[t.id])continue;
		int siz=e[t.id].size();
		for(int i=0;i<siz;i++){
			pii a=e[t.id][i];if(!a.first)continue;
			if(1ll*a.second+dis[t.id]<dis[a.first])
				dis[a.first]=1ll*a.second+dis[t.id],
				q.push((ming){a.first,dis[a.first]});
		}
	}
}
signed main(){
	read(n);read(m);read(s);cnt=n;build(1,1,n);
	for(int i=1;i<=m;i++){
		int opt;read(opt);
		if(opt==1){
			int u,v,w;read(u);read(v);read(w);
			e[u].push_back(make_pair(v,w));
		}
		if(opt==2){
			int u,vl,vr,w;
			read(u);read(vl);read(vr);read(w);
			updatain(1,1,n,vl,vr,u,w);
		}
		if(opt==3){
			int v,ul,ur,w;
			read(v);read(ul);read(ur);read(w);
			updataout(1,1,n,ul,ur,v,w);
		}
	}
	dijkstra();
	for(int i=1;i<=n;i++)
		printf("%lld ",(dis[i]>INF-1)?-1:dis[i]);
	return 0;
}

谢谢!!!

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值