CF786B Legacy 线段树优化建图

传送门

文章目录

  • 题意:
  • 思路:

题意:

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

思路:

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315104.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【ZJOI2015】幻想乡 Wi-Fi 搭建计划【几何】【贪心】【dp】

传送门 题意&#xff1a;一个x∈(−∞,∞),y∈[0,R]x\in(-\infin,\infin),y\in[0,R]x∈(−∞,∞),y∈[0,R]的矩形中有nnn个点&#xff0c;矩形外有mmm个半径均为RRR的圆&#xff0c;有独立的代价cic_ici​。求覆盖最多的点所需的最小代价。 n,m≤100n,m\leq100n,m≤100 显然先…

.NET架构开发应知应会

.NET程序是基于.NET framework、.NET Core、Mono、UWP【.NET实现】开发和运行的 &#xff0c;定义以上【.NET实现】的标准规范称为.NET StandardL1&#xff1a;.NET Standard.NET标准是一组API集合&#xff0c;由上层三种【.NET实现】的Basic Class Library实现&#xff0c;更正…

几个冷门字符串算法的学习笔记(最小表示法,exKMP,Lyndon Word)

所有下标均从1开始 最小表示法 给定一个串&#xff0c;求字典序最小的循环同构。 我们把串复制一遍接在后面&#xff0c;然后求出[1,N][1,N][1,N]开始的长为NNN的子串中最小的 先设i1,j2i1,j2i1,j2 然后暴力找出iii和jjj往后匹配的第一个不同的位置&#xff0c;记为ikikik…

P6348 [PA2011]Journeys 线段树优化建图 区间连区间

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 每次连接[a,b][a,b][a,b]与[c,d][c,d][c,d]之间所有点&#xff0c;让后跑最短路。 思路&#xff1a; 比普通的优化建图能简单点&#xff0c;我们只需要加两个虚点之间边权为111&#xff0c;让后让某个点连…

.NET Core IdentityServer4实战 第Ⅴ章-单点登录

OiDc可以说是OAuth的改造版&#xff0c;在最初的OAuth中&#xff0c;我们需要先请求一下认证服务器获取下Access_token&#xff0c;然后根据Access_token去Get资源服务器, 况且OAuth1 和 2 完全不兼容&#xff0c;易用性差&#xff0c;而OIDC可以在登陆的时候就把信息返回给你&…

【CF594E】Cutting the Line 【贪心】【Lyndon Word】【扩展kmp】

传送门 题意&#xff1a;给一个字符串SSS和正整数kkk&#xff0c;将SSS分成最多kkk段&#xff0c;每段不变或翻转&#xff0c;使得最后的字典序最小。 ∣S∣≤5106|S|\leq5\times10^6∣S∣≤5106 发现不翻转可以看成拆成若干单字符分别翻转&#xff0c;所以先分析一下必须翻转…

一份好的工作总结才能帮你升职加薪

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「79」篇原创敬上最近有点忙&#xff0c;搬出之前攒的一篇文章来应急一下。一篇能助你挣更多钱的文章。好了&#xff0c;下面开始。我的读者们大部分是互联网行业的&#xff0c…

CF1385E Directing Edges 拓扑序

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个图和若干个边&#xff0c;有些是有向边&#xff0c;有些是无向边&#xff0c;让你给无向边定向&#xff0c;使得最终的图是DAGDAGDAG。 思路&#xff1a; 题目让构造DAGDAGDAG&#xff0c;比较容易…

【BZOJ3684】大朋友和多叉树【生成函数】【拉格朗日反演】【多项式幂函数】

传送门 题意&#xff1a;给定nnn和集合SSS&#xff0c;求含nnn个叶子结点、非叶子节点的儿子数在SSS内的树的个数 模 950009857(4532211)950009857(453\times2^{21}1)950009857(4532211)。结点无标号但儿子间有顺序。 n≤105n \leq 10^5n≤105 算是拉格朗日反演模版题了吧………

Codeforces Round #723 (Div. 2) D. Kill Anton 线段树 + 暴力

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个只有ANTOANTOANTO四个字母的字符串&#xff0c;你每次可以交换相邻两个&#xff0c;花费为111&#xff0c;让后让你打乱字符串&#xff0c;使得将打乱的字符串还原为原来的字符串的花费最小。 n≤1e…

腾讯开源软件镜像站上线

腾讯开源软件镜像站(Tencent Open Source Mirror Site)已于近日上线&#xff0c;其官方名称为「腾讯云软件源」&#xff0c;由腾讯云提供支持。官方表示搭建此开源镜像站的目的在于宣传自由软件的价值&#xff0c;提高自由软件社区文化氛围&#xff0c;推广自由软件在国内的应用…

【集训队作业2018】复读机【指数型生成函数】【单位根反演】【二项式定理】

传送门 单位根反演听着高级&#xff0c;其实没啥技术含量…… 本文是篇几乎没有证明的佛系讲解 单位根反演的式子长这样&#xff1a; 1n∑i0n−1ωnik[k∣n]\frac{1}{n}\sum_{i0}^{n-1}\omega_n^{ik}[k|n]n1​i0∑n−1​ωnik​[k∣n] 其实本质是IFFT 感觉懵的&#xff1f;…

ASP.NET Core on K8S学习初探(2)

“ [LOG] ASP.NET Core on K8S Starting...”在上一篇《单节点环境搭建》中&#xff0c;通过Docker for Windows在Windows开发机中搭建了一个单节点的K8S环境&#xff0c;接下来就是动人心弦的部署ASP.NET Core API到K8S了。但是&#xff0c;在部署之前&#xff0c;我还是把基本…

Educational Codeforces Round 96 E. String Reversa 线段树模拟序列交换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 与上一篇题解大同小异&#xff0c;无非就是不需要枚举排列了。 // Problem: E. String Reversal // Contest: Codeforces - Educational Codeforces Round 96 (Rated for Div. 2) // URL:…

【LOJ6363】「地底蔷薇」【点双】【指数型生成函数】【扩展拉格朗日反演】【多项式幂函数】

传送门 题意&#xff1a;给定nnn和集合SSS,求出nnn个点的「所有极大点双连通分量的大小都在SSS 内」的不同简单无向连通图的个数 模 998244353998244353998244353。 n,∑i∈Si≤105n,\sum_{i\in S}i \leq10^5n,∑i∈S​i≤105 道理我都懂&#xff0c;可为啥我百度搜地灵殿ex终…

ASP.NET Core on K8S学习初探(1)

“ [LOG] ASP.NET Core on K8S Starting...”01—写在之前当近期的一个App上线后&#xff0c;发现目前的docker实例&#xff08;应用服务BFF中台服务工具服务&#xff09;已经越来越多了&#xff0c;而我司目前没有专业的运维人员&#xff0c;发现运维的成本逐渐开始上来&#…

AtCoder Regular Contest 120 C - Swaps 2 线段树模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个序列a,ba,ba,b&#xff0c;每次可以执行一个操作&#xff1a;将a[i]a[i]a[i]与a[i1]a[i1]a[i1]交换&#xff0c;且让交换后的a[i]1,a[i1]−1a[i]1,a[i1]-1a[i]1,a[i1]−1&#xff0c;问将aaa变成bbb…

【BZOJ 4671】异或图 【斯特林反演】【线性基】【贝尔数复杂度】

传送门 题意&#xff1a;定义两个图的异或的边集为在两张图中恰出现一次的边。给sss张nnn个点的图的集合&#xff0c;求异或和为连通图的子集数。 s≤60,n≤10s \leq 60,n \leq 10s≤60,n≤10 设GiG_iGi​表示异或出iii个连通块的子集数&#xff0c;答案就是G1G_1G1​ GGG并不…

.NET Core IdentityServer4实战 第Ⅳ章-集成密码登陆模式

回顾下ClientCredentials模式&#xff0c;在ReSourceApi中定义了我们公开服务&#xff0c;第三方网站想要去访问ReSourceApi则需要在身份验证服务中获取toekn&#xff0c;根据token的内容&#xff0c;硬编码去访问公开服务&#xff08;ResApi&#xff09;,这个还是非常简单的&a…

Codeforces Round #656 (Div. 3) F. Removing Leaves 贪心 + 模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个贪心策略就是每次都找一个叶子节点最多的点&#xff0c;让后删掉他的kkk个叶子节点&#xff0c;现在我们就来考虑如何模拟这个过程。 我们整一个vector<set<int>>ve…