CF786B Legacy 线段树优化建图

最新推荐文章于 2022-03-23 21:46:51 发布
最新推荐文章于 2022-03-23 21:46:51 发布
阅读量136 收藏 1
点赞数
分类专栏： 数据结构 图论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_51068403/article/details/117308186
版权

传送门

文章目录

题意：

实现如下连边后跑最短路。
在这里插入图片描述

思路：

优化建图板子题，优化思路就是将区间分割成若干个线段树上的线段，与线段树分治有点类似，由于有点向区间也有区间向点的边，那么我们需要建两颗线段树，点向区间的树是自顶向下连边，区间向点的是自底向上连边，最后两棵树的最后一层需要与原来的 n n n个点连双向(未经说明都是单向边)，比如下面这个图(盗用了日报的图)
在这里插入图片描述
让后我们就按照图来建边就好啦，我这里第一颗编号是 [ 1 , 4 ∗ n ] [1,4*n] [1,4∗n]，第二棵编号是 [ 4 ∗ n + 1 , 8 ∗ n ] [4*n+1,8*n] [4∗n+1,8∗n]，绿色点是 [ 8 ∗ n + 1 , 9 ∗ n ] [8*n+1,9*n] [8∗n+1,9∗n]。
当然可以舍去绿色的点，直接以线段树叶子节点为绿色的点，下面也给出代码了，只需要记一下 l e a f [ i ] leaf[i] leaf[i]即可。

// Problem: B. Legacy
// Contest: Codeforces - Codeforces Round #406 (Div. 1)
// URL: https://codeforces.com/problemset/problem/786/B
// Memory Limit: 256 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<LL,int> PII;

const int N=900200,M=N*40,mod=1e9+7,INF=0x3f3f3f3f;
const LL inf=0x3f3f3f3f3f3f3f3f;
const double eps=1e-6;

int n,q,s,b1,b2;
int e[M],ne[M],w[M],h[N],idx;
LL dis[N];
bool st[N];
//第一颗线段树1-4*n 第二棵4*n+1-8*n 点8*n+1-9*n

void add(int a,int b,int c) {
	e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

void build(int u,int l,int r) {
	if(l==r) {
		add(u,l+8*n,0); add(l+8*n,u,0);
		add(u+4*n,l+8*n,0); add(l+8*n,u+4*n,0);
		return;
	}
	add(u,u*2,0); add(u,u*2+1,0);
	add(u*2+4*n,u+4*n,0); add(u*2+1+4*n,u+4*n,0);
	int mid=(l+r)>>1;
	build(u<<1,l,mid); build(u<<1|1,mid+1,r);
}

void change1(int u,int l,int r,int ql,int qr,int st,int cs) {
	if(ql<=l&&qr>=r) {
		add(st+8*n,u,cs);
		return;
	}
	int mid=(l+r)>>1;
	if(ql<=mid) change1(u<<1,l,mid,ql,qr,st,cs);
	if(qr>mid) change1(u<<1|1,mid+1,r,ql,qr,st,cs);
}

void change2(int u,int l,int r,int ql,int qr,int st,int cs) {
	if(ql<=l&&qr>=r) {
		add(u+4*n,st+8*n,cs);
		return;
	}
	int mid=(l+r)>>1;
	if(ql<=mid) change2(u<<1,l,mid,ql,qr,st,cs);
	if(qr>mid) change2(u<<1|1,mid+1,r,ql,qr,st,cs);
}

void dijkstra() {
	priority_queue<PII,vector<PII>,greater<PII> >q; q.push({0ll,s+8*n});
	memset(dis,0x3f3f3f3f,sizeof(dis));
	dis[s+8*n]=0;
	while(q.size()) {
		PII u=q.top(); q.pop();
		if(st[u.Y]) continue;
		st[u.Y]=1;
		for(int i=h[u.Y];~i;i=ne[i]) {
			int j=e[i];
			if(dis[j]>dis[u.Y]+w[i]) {
				dis[j]=dis[u.Y]+w[i];
				q.push({dis[j],j});
			}
		}
	}
	for(int i=1;i<=n;i++) printf("%lld ",dis[i+8*n]==inf? -1:dis[i+8*n]);
}

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);

	memset(h,-1,sizeof(h));
	cin>>n>>q>>s;
	build(1,1,n);
	while(q--) {
		int op,v,l,r,c;
		scanf("%d",&op);
		if(op==1) {
			scanf("%d%d%d",&l,&r,&c);
			add(l+8*n,r+8*n,c);
		}
		else if(op==2) {
			scanf("%d%d%d%d",&v,&l,&r,&c);
			change1(1,1,n,l,r,v,c);
		}
		else {
			scanf("%d%d%d%d",&v,&l,&r,&c);
			change2(1,1,n,l,r,v,c);
		}
	}
	dijkstra();

	




	return 0;
}
/*

*/










// Problem: B. Legacy
// Contest: Codeforces - Codeforces Round #406 (Div. 1)
// URL: https://codeforces.com/problemset/problem/786/B
// Memory Limit: 256 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<LL,int> PII;

const int N=800200,M=N*20,mod=1e9+7,INF=0x3f3f3f3f;
const LL inf=0x3f3f3f3f3f3f3f3f;
const double eps=1e-6;

int n,q,s,b1,b2;
int e[M],ne[M],w[M],h[N],idx;
int leaf[N];
LL dis[N];
bool st[N];
//第一颗线段树1-4*n 第二棵4*n+1-8*n 

void add(int a,int b,int c) {
	e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}

void build(int u,int l,int r) {
	if(l==r) {
		leaf[l]=u;
		return;
	}
	int mid=(l+r)>>1;
	add(u,u*2,0); add(u,u*2+1,0);
	add(u*2+4*n,u+4*n,0); add(u*2+1+4*n,u+4*n,0);
	build(u<<1,l,mid); build(u<<1|1,mid+1,r);
}

void change1(int u,int l,int r,int ql,int qr,int st,int w) {
	if(ql<=l&&qr>=r) {
		add(st,u,w);
		return;
	}
	int mid=(l+r)>>1;
	if(ql<=mid) change1(u<<1,l,mid,ql,qr,st,w);
	if(qr>mid) change1(u<<1|1,mid+1,r,ql,qr,st,w);
}

void change2(int u,int l,int r,int ql,int qr,int st,int w) {
	if(ql<=l&&qr>=r) {
		add(u+4*n,st,w);
		return;
	}
	int mid=(l+r)>>1;
	if(ql<=mid) change2(u<<1,l,mid,ql,qr,st,w);
	if(qr>mid) change2(u<<1|1,mid+1,r,ql,qr,st,w);
}

void dijkstra() {
	memset(dis,0x3f3f3f3f,sizeof(dis));
	priority_queue<PII,vector<PII>,greater<PII>>q;
	q.push({s,leaf[s]}); dis[leaf[s]]=0;
	while(q.size()) {
		PII u=q.top(); q.pop();
		if(st[u.Y]) continue;
		st[u.Y]=1;
		for(int i=h[u.Y];~i;i=ne[i]) {
			int j=e[i];
			if(dis[j]>dis[u.Y]+w[i]) {
				dis[j]=dis[u.Y]+w[i];
				q.push({dis[j],j});
			}
		}
	}
	for(int i=1;i<=n;i++) printf("%lld ",dis[leaf[i]]==inf? -1:dis[leaf[i]]);
}

int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);

	memset(h,-1,sizeof(h));
	cin>>n>>q>>s;
	build(1,1,n);
	for(int i=1;i<=n;i++) add(leaf[i],leaf[i]+4*n,0),add(leaf[i]+4*n,leaf[i],0);
	while(q--) {
		int op,v,l,r,c;
		scanf("%d",&op);
		if(op==1) {
			scanf("%d%d%d",&l,&r,&c);
			add(leaf[l],leaf[r],c);
		}
		else if(op==2) {
			scanf("%d%d%d%d",&v,&l,&r,&c);
			change1(1,1,n,l,r,leaf[v],c);
		}
		else {
			scanf("%d%d%d%d",&v,&l,&r,&c);
			change2(1,1,n,l,r,leaf[v],c);
		}
	}
	dijkstra();
	




	return 0;
}
/*

*/
YB Lin
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
CF786B Legacy
weixin_33938733的博客
05-09 70
CF786B Legacy 给定一张有向图，有三种边： \(u\to v\) ，边权为 \(w\) \(u\to [l,\ r]\) ，边权为 \(w\) \([l,\ r]\to u\) ，边权为 \(w\) \(n,\ m\leq10^5,\ w\leq10^9\) 线段树优化建图，最短路 线段树优化建图板子题，注意空间大小 由于有 \(n\log n\) 条边， 且普通 dijks...
线段树优化建图 || CF786B Legacy
weixin_30678821的博客
08-05 102
题面：786B - Legacy 代码： 1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 #include<queue> 5 #define ll long long 6 using namespace std; ...
[CF786B]Legacy
StaroForgin的博客
08-23 139
Legacy 题解 看到题目应该是容易想到最短路的，但是由于区间到单点与单点到区间的边我们需要想些办法来维护。 于是，我们就想到了通过虚点来进行维护。但总不能每个区间都建一个点，只能利用线段树来对各个区间的虚点进行维护。 针对一个区间需要建两个点，一个入点，一个出点。将大区间与小区间连上边，出边从下往上，入边从上往下，其边权都是0。 由于大小区间内已经连了边，再加上最短路的性质，懒标记什么的就没必要了。最后建完边后再去跑一边Dijkstra即可，注意需要加堆优化。 由于，边数应该不超过，总时间复
Luogu CF786B Legacy
YYHS_WSF的博客
11-08 816
https://www.luogu.org/problemnew/show/CF786B  CF786B Legacy 104通过 404提交 题目来源CodeForces 786B 评测方式RemoteJudge 标签 难度省选/NOI- 时空限制2000ms / 256MB  提交  题解    提示：收藏到任务计划后，可在首页查看。 最新讨论显示 推荐的相关题目...
CF786B Legacy（线段树优化建图）
weixin_30607659的博客
09-18 153
传送门 ovo这题该怎么做呢？我们首先考虑暴力建图，但是因为建图的操作太多直接就会MLE，所以这个就别想了…… 我们考虑如何优化建图。因为发现一个点可以向一个区间连边，一个区间也可以向一个点连边，想到区间很容易想到线段树……所以我们可以使用一个点来代替一段区间进行建图。 具体的操作可以这样，我们建立两棵线段树，一棵用于存所有的出边（就是区间向点连边），每个点向其父亲建一条边权为0的边，另一棵...
Codeforces-786B-Legacy （线段树+最短路）
On the way
03-24 1700
题目链接：Codeforces-786B-Legacy看到区间很容易就能想到线段树。 对于操作2，我们可以建立一棵线段树，线段树的某一个结点表示区间 [l,r][l,r] ，那么在图中建立一个结点 uu 表示区间 [l,r][l,r] ，新建点 uu 指向 l,l+1,l+2...rl,l+1,l+2...r 的边，权值为 00 。 当得到一个操作 (v,l,r,w)(v,l,r,w) 时，查
CF786B Legacy 线段树优化Dijkstra
weixin_43908980的博客
08-24 250
文章目录一.题目二.题解三.CodeThanks ！ 一.题目 传送门 翻译： 二.题解 如果这道题目没有后面两个操作，那么就是一个Dijkstra跑过，但是当加上后两个操作，如果vvv向每个区间内的点连一条边将会连n2n^2n2条边，存不下也跑不过。 这时就得使用线段树来优化Dijkstra。 我们不如直接将点vvv连到区间l−rl-rl−r上或将区间l−rl-rl−r连到点vvv上，具体操作下面我将线段树的打法解释一下： 首先要建两颗线段树treein[Index],treeout[Index]tre
codeforces786B - Legacy
juruo? juruo!
01-04 323
题面在这里题意：三种操作： 1.给a->b加一条边 2.给a->[b,c]分别加一条权值相同的边 3.给[a,b]->c分别加一条权值相同的边 然后问一个点到其他点的最短路做法：考虑在线段树上加边 一个点对应一个区间，就把那个点和那个区间所对应的线段树节点编号加边 然后线段树上父亲和儿子再分别连一条权值为0的边 然后就可以直接跑最短路了 好妙啊。代码：/***************
CodeForces - 786B Legacy【最短路+线段树优化建边】
niiick
11-08 411
Time limit 2000 ms Memory limit 262144 kB Rick and his co-workers have made a new radioactive formula and a lot of bad guys are after them. So Rick wants to give his legacy to Morty before bad guys ca...
【题解】codeforces786B Legacy
ZJZ的博客
03-31 1153
题目链接 题意：给定一个n个点m条带权值的“边”的有向图，其中“边”分为3类。第1类“边”是点u到点v的边，第2类“边”是点u到区间[l,r]中的点的边（即r-l+1条边），第3类“边”是区间[l,r]中的点到点u的边。求以点s为源点的最短路。 分析：和atcoder2336类似。考虑给图加一些辅助点和边。对第2类“边”按如下方式处理：对点集[1,n]建立一棵线段树，线段[l,r]维护dis[
P3085 [USACO13OPEN]Yin and Yang G 点分治
m0_51068403的博客
02-19 2756
文章目录题意：思路： 传送门 题意： 给你一颗nnn个点的树，每条边为黑色或者白色，问满足以下条件的路径条数：路径上存在一个不是端点的点，使得两端点到该点的两条路径上两种颜色的边数相等。 1≤n≤1000001\le n\le 1000001≤n≤100000 思路： 统计树上路径问题显然需要用到点分治了，这个题维护的信息比较麻烦，想明白了思路还需要考虑如何做才能使代码变得简便好写。 考虑点分治每个步骤需要算的贡献，在选出重心之后，之后的子树中就不会包含重心这个点，所以贡献需要算每颗包含重心的子树内的路径以
The 2021 ICPC Asia Taipei Regional F. What a Colorful Wall 扫描线 + 并查集
m0_51068403的博客
03-23 1159
文章目录题意:思路 传送门 题意: 给你平面nnn个矩形，每个矩形有一种颜色，依次给出矩形以及其的颜色，后面的矩形会覆盖前面的矩形，问最终有多少种颜色。 1≤n≤4000,0≤x1<x2<228,0≤y1<y2<228,1≤c≤n1\le n\le 4000,0\le x_1<x_2< 2^{28},0\le y_1<y_2< 2^{28},1\le c\le n1≤n≤4000,0≤x1​<x2​<228,0≤y1​<y2​<228,
可持久化Splay 学习笔记
m0_51068403的博客
05-12 693
可持久化Splay是怎么回事呢？Splay相信大家都很熟悉，但是可持久化Splay是怎么回事呢，下面就让小编带大家一起了解吧。 　　可持久化Splay，其实就是将Splay持久化一下，大家可能会很惊讶Splay怎么可以持久化呢？但事实就是这样，小编也感到非常惊讶。 　　这就是关于可持久化Splay的事情了，大家有什么想法呢，欢迎在评论区告诉小编一起讨论哦！ ...
XXI Open Cup. Grand Prix of Korea I. Query On A Tree 17 树剖 + 二分 + 树带权重心
m0_51068403的博客
08-23 597
传送门 文章目录题意：思路： 题意： 给你一棵树，每棵树初始权值都为000，现在给你两个操作： (1)(1)(1)将uuu的子树权值全部加111。 (2)(2)(2)将(u,v)(u,v)(u,v)路径上的点权值都加111。 每次输出一个点xxx,满足∑y=1Na[y]∗dis(x,y)\sum_{y=1}^Na[y]*dis(x,y)∑y=1N​a[y]∗dis(x,y)最小，如果多个点相同，输出深度最小的点。 思路： 观察这个式子，发现就是一个带权重心的问题，如果不在重心显然向重心走更优，需要支持修改。
2020EC-final
m0_51068403的博客
11-22 565
传送门 文章目录B - Rectangle Flip 2题意：思路：A - Namomo Subsequence题意：思路：D - City Brain题意:思路: B - Rectangle Flip 2 题意： 给你一个n∗mn*mn∗m的矩阵，接下来n∗mn*mn∗m秒每秒都会消失一个格子，问每个时刻矩阵中构成的矩形有多少个。 n,m≤500n,m\le 500n,m≤500 思路： 先说一下复杂度n4n^4n4但是跑不满的算法，对于每个删去的点，枚举左边以及右边每个位置，维护其能到的上下界，比如上界
from torchtext.legacy import data, 用SGD优化器的代码
最新发布
02-06
上述代码是使用 torchtext.legacy 中的 TabularDataset 和 Iterator 类来加载和处理数据，然后使用 PyTorch 中的 nn.Module 和 optim.SGD 类来构建模型和优化器，完成训练。 注意, torchtext.legacy 为过时的版本，...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值