势能线段树

已于 2023-08-03 18:16:41 修改
阅读量120 收藏
点赞数
分类专栏： 数据结构 文章标签： 数据结构
于 2023-08-03 18:15:25 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_63761643/article/details/132088382
版权

目录

简单介绍

 题目

1. 上帝造题的七分钟 2 

2.SUM and REPLACE

3. And RMQ

总结 

简单介绍

 题目

1. 上帝造题的七分钟 2 

链接：https://www.luogu.com.cn/problem/P4145

 维护两种操作

1.区间开根号(下取整)

2.区间和询问

显然无法通过懒标记来计算区间开根号后的值，其由叶子结点本身的值决定。容易发现当一个数连续进行开根号操作会在很少的次数变为1，且值不再改变，1即为零势能点。因此我们可以维护区间max。一旦区间修改时发现此区间的max<=1时，我们不需要再次修改，直接return即可，否则向下递归修改。

Code:

#include<bits/stdc++.h>
using namespace std;
#define PII pair<int,int>
#define endl "\n"
#define int long long
const int N=1e5+10;
struct segment_tree {
    int a[N];
    struct node {
        int l,r;
        int mx,sum;
    }tr[N<<2];
    
    void build(int u,int l,int r) {
        tr[u].l=l,tr[u].r=r;
        if(l==r) {
            tr[u].mx=tr[u].sum=a[l];
            return ;
        }
        int mid=(l+r)>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    } 

    void modify(int u,int l,int r) {
        if(tr[u].mx<=1) return ;
        if(tr[u].l==tr[u].r) {
            tr[u].mx=sqrt(tr[u].mx);
            tr[u].sum=tr[u].mx;
            return ;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        if(l<=mid) modify(u<<1,l,r);
        if(r>mid) modify(u<<1|1,l,r);
        pushup(u);
    } 

    int query(int u,int l,int r) {
        if(tr[u].l>=l&&tr[u].r<=r) {
            return tr[u].sum;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        int res=0;
        if(l<=mid) res+=query(u<<1,l,r);
        if(r>mid) res+=query(u<<1|1,l,r);
        return res;
    }
    
    void pushup(int u) {
        tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;
        tr[u].mx=max(tr[u<<1].mx,tr[u<<1|1].mx);
    }
}ST;
signed main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int n,m;
    cin>>n;
    for(int i=1;i<=n;i++) cin>>ST.a[i];
    ST.build(1,1,n);
    cin>>m;
    while(m--) {
        int op,l,r;
        cin>>op>>l>>r;
        if(l>r) swap(l,r);
        if(op==0) ST.modify(1,l,r);
        else cout<<ST.query(1,l,r)<<endl;     
    }
}
2.SUM and REPLACE

链接：https://codeforces.com/contest/920/problem/F

定义f(x)为x因子的数量

维护三种操作

1.区间修改x=f(x)  

2.区间和查询

手动模拟f(x)可以发现，进行f(x)操作，数值单调不增，且x<=2时，其值不在改变，因此同上题一样维护区间最大值即可。f(x)可以O(nlogn)时间预处理出来。

#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define PII pair<int,int>
#define endl "\n"
#define int long long
const int N=3e5+10,M=1e6+10;
int d[M+1];
struct segment_tree {
    int a[N];
    struct node {
        int l,r;
        int mx,sum;
    }tr[N<<2];
    
    void build(int u,int l,int r) {
        tr[u].l=l,tr[u].r=r;
        if(l==r) {
            tr[u].mx=tr[u].sum=a[l];
            return ;
        }
        int mid=(l+r)>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    } 

    void modify(int u,int l,int r) {
        if(tr[u].mx<=2) return ;
        if(tr[u].l==tr[u].r) {
            tr[u].mx=d[tr[u].mx];
            tr[u].sum=tr[u].mx;
            return ;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        if(l<=mid) modify(u<<1,l,r);
        if(r>mid) modify(u<<1|1,l,r);
        pushup(u);
    } 

    int query(int u,int l,int r) {
        if(tr[u].l>=l&&tr[u].r<=r) {
            return tr[u].sum;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        int res=0;
        if(l<=mid) res+=query(u<<1,l,r);
        if(r>mid) res+=query(u<<1|1,l,r);
        return res;
    }

    void pushup(int u) {
        tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;
        tr[u].mx=max(tr[u<<1].mx,tr[u<<1|1].mx);
    }
};
signed main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    for(int i=1;i<=M;i++) {  //预处理
        for(int j=i;j<=M;j+=i) {
            d[j]++;
        }
    }
    int n,m;
    cin>>n>>m;
    segment_tree ST;
    for(int i=1;i<=n;i++) cin>>ST.a[i];
    ST.build(1,1,n);
    while(m--) {
        int op,l,r;
        cin>>op>>l>>r;
        if(op==1) ST.modify(1,l,r);
        else cout<<ST.query(1,l,r)<<endl;
    }
}
3. And RMQ

链接：

维护三个操作

1.区间按位与x

2.区间最大值

3.单点修改

这题零势能点藏得较深，我们考虑将x二进制展开，发现在x的二进制位为零的位置，区间所有数的二进制位也为零，则操作可以忽略。维护区间或   orsum=a_{l}|a_{l+2}|a_{l+2}|a_{r-1}|...|a_{r}  

当 orsum & x =orsum,则直接return 

Code:

#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define PII pair<int,int>
#define endl "\n"
const int N=4e5+10;
struct segment_tree {
    int a[N];
    struct node {
        int l,r;
        int mx,sum;
    }tr[N<<2];
    
    void build(int u,int l,int r) {
        tr[u].l=l,tr[u].r=r;
        if(l==r) {
            tr[u].mx=tr[u].sum=a[l];
            return ;
        }
        int mid=(l+r)>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    } 

    void modify(int u,int l,int r,int x) {
        if((tr[u].sum&x)==tr[u].sum) return ;
        if(tr[u].l==tr[u].r) {
            tr[u].mx=tr[u].mx&x;
            tr[u].sum=tr[u].mx;
            return ;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        if(l<=mid) modify(u<<1,l,r,x);
        if(r>mid) modify(u<<1|1,l,r,x);
        pushup(u);
    } 

    int query(int u,int l,int r) {
        if(tr[u].l>=l&&tr[u].r<=r) {
            return tr[u].mx;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        int res=0;
        if(l<=mid) res=max(res,query(u<<1,l,r));
        if(r>mid) res=max(res,query(u<<1|1,l,r));
        return res;
    }
    
    void update(int u,int k,int x) {
        if(tr[u].l==tr[u].r) {
            tr[u].mx=tr[u].sum=x;
            return ;
        }
        int mid=(tr[u].l+tr[u].r)>>1;
        if(k<=mid) update(u<<1,k,x);
        else update(u<<1|1,k,x);
        pushup(u);
    }

    void pushup(int u) {
        tr[u].sum=tr[u<<1].sum|tr[u<<1|1].sum;
        tr[u].mx=max(tr[u<<1].mx,tr[u<<1|1].mx);
    }
}ST;
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>ST.a[i];
    ST.build(1,1,n);
    while(m--) {
        int l,r,x;
        string op;
        cin>>op;
        if(op=="AND")  {
            cin>>l>>r>>x;
            ST.modify(1,l,r,x);
        }
        else if(op=="UPD") {
            cin>>l>>x;
            ST.update(1,l,x);
        }
        else {
            cin>>l>>r;
            cout<<ST.query(1,l,r)<<endl;
        }
    }
}

总结 

1.对于区间修改操作，修改操作会使得值在趋向零势能点前严格单调减少，在变为零势能点后不在变化。需要维护一个值来界定是否到达零势能

2.且题目不能出现其他非单调的区间修改操作，如区间加，区间乘等。如果有其他修改操作，可以通过构造形如 update1 ,update 2,update 1,update 2 的数据破坏单调性，从而使操作1复杂度变为暴力修改的O(nlogn)

Cambrain_
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
势能线段树专题
「 虚幻私塾」
04-11 578
Python微信订餐小程序课程视频 https://blog.csdn.net/m0_56069948/article/details/122285951 Python实战量化交易理财系统 https://blog.csdn.net/m0_56069948/article/details/122285941 势能线段树(吉司机线段树) 简单介绍和理解 我们知道传统的支持区间修改的线段树,我们都是靠lazylazylazy标记来节省开销的。可以使用lazylazylazy标记必须要满足下面两个条件： 区间节
势能线段树（均摊分析）
Edward The Bunny 的博客
03-23 987
对于每个区间维护区间内的 最大值MaxMaxMax 和 最小值MinMinMin，对于整除操作，如果有Max−⌊Maxd⌋=Min−⌊Mind⌋Max-\lfloor\frac{Max}{d}\rfloor = Min − \lfloor\frac{Min}{d}\rfloorMax−⌊dMax​⌋=Min−⌊dMin​⌋，就转化为区间减。否则直接向下递归。 考虑何时满足Max−⌊Maxd⌋=Min−⌊Mind⌋Max-\lfloor\frac{Max}{d}\rfloor = Min − \lfloor\
势能线段树（吉司机线段树）专题
jibinquan的博客
10-03 2671
势能线段树（吉司机线段树）专题 势能线段树在近期训练时遇到了好几次，但是由于本人太懒一直没补完，结果ICPC网络赛还真就出了一道势能线段树Orz……结果当然是没做出来……痛定思痛，这回把之前欠的一块儿补了。 简要介绍 ​ 我们知道，传统的支持区间修改的线段树是通过lazy标记来规避掉频繁大规模的单点修改来实现对时间开销的节省的。那么对于一种区间操作，其能使用lazy标记需要满足两个条件： 区间节点的值可以通过对当前结点lazy标记的计算来更新 多次不同的lazy标记可以实现就地的快速合并
势能线段树学习笔记
PLUS_CAT的博客
04-04 803
也叫吉司机线段树。
CF765F Souvenirs（势能线段树）
10-21 460
(⊙o⊙)？好强的势能线段树w(ﾟДﾟ)w
势能线段树模板题
weixin_45772744的博客
04-14 1001
题目链接: C-势能线段树模板题二_牛客竞赛数据结构专题班势能线段树、李超线段树 (nowcoder.com) 题意: 样例: 输入 10 3 10 10 10 10 10 10 10 10 10 10 1 1 5 2 1 10 5 3 1 10 输出 115 代码: #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N=2e5+10; typede.
势能线段树+弩蚊怒夏(详细思路)
m0_51780913的博客
09-03 512
势能线段树 势能线段树 势能： 信息学中，势能被用于计算某一个过程，或者某一类过程时间复杂度的总和 例如计算两个数之间的最大公倍数gcd的时间复杂度是O(log N)，而计算n个数之间的时间复杂度是O(n+log N) 总复杂度=数组循环遍历复杂度+gcd函数被调用的总复杂度 势能均摊复杂度 势能均摊复杂度是指把总的时间复杂度摊到操作次数或者循环次数上面 在上述n个数之间的gcd的时间复杂度为： O(N+logC)N=O(1+(log⁡CN))=O(1) \frac{O(N+logC)}{N}=O(1+(\
【题型总结】势能线段树
xhyu61的博客
06-22 341
那就好。我…我不擅长这种舞蹈，没有丢脸真是太好了。也就是心血来潮，想要任性一次…哈哈哈，这么说来，今天可真是任性的一天呀。 ——《原神》神里绫华
复杂度证明1
08-08
这里我们讨论的"复杂度证明1"涉及到的是一个特定的算法，它利用了数据结构——线段树，并着重于处理区间赋值和查询操作的复杂度。 线段树是一种自顶向下的平衡树结构，用于高效地支持区间查询和更新操作。在这个...
江苏省徐州市八年级物理下册8.2重力力的示意图学案无答案新版苏科版
09-09
通过画出带箭头的线段，我们可以清楚地看出物体受到的各个力的情况。例如，对于题目中的苹果，我们需要画出苹果受到的重力示意图，这通常包括一个指向下方的箭头，表示苹果受到的重力。 在质量与重力的关系上，我们...
贵州省遵义航天高级中学2017_2018学年高二物理上学期第一次9月月考试题文含解析
09-09
位移是从初位置到末位置的有向线段，是矢量。题中提到，坐高铁从杭州到南京，经过上海的路程大于直达南京的路程，但两次的位移是相同的。 3. 瞬时速度与平均速度：瞬时速度是物体在某一瞬间或某一位置的速度，而...
高一物理下册期末学业水平考试试题(卷）.doc
10-11
3. **位移概念**：位移是初位置到末位置的直线距离，可以是零、可以是线段的长度，也可以是两点间的弧长。例如，沿着半圆轨道从A点到B点的位移是半圆的直径长度2R。 4. **作用力与反作用力**：牛顿第三定律指出，...
数据结构(5.3_3)——由遍历序列构造二叉树
m0_65240792的博客
07-25 162
若只给出一棵二叉树的前/中/后/层 序遍历序列中的一种，不能唯一确定一棵二叉树 例：前序遍历序列：DAEFBCHGI 中序遍历序列：DAEFBCHGI
学习记录day19——数据结构 查找算法
最新发布
weixin_64918301的博客
07-25 213
在给定数据元素的某个值，在查找表中确定一个其关键字等于给定值的数据元素的操作，叫做查找。
数据结构---二叉树
weixin_47808575的博客
07-22 292
参考：二叉树遍历方法——前、中、后序遍历（图解）_二叉树遍历前序中序后序-CSDN博客
数据结构——二叉树性质
鲸落博客
07-23 412
这个性质很好记忆，观察一下图6-5-5。第一层是根结点，只有一个，所以2^(1-1)=2^0=1。第二层有两个，2^(2-1)=2=2。第三层有四个，2^(3-1)=2^2=4。第四层有八个，2^(4-1)=2^3=8。通过数据归纳法的论证，可以很容易得出在二叉树的第i层上至多有2^(i-1)个结点(i≥1)的结论。
人工势能场matlab
06-26
人工势能场通常被用于机器人路径规划、目标跟踪和机器人导航等领域，它将环境映射为一个势能场，机器人通过在势能场中移动来实现目标的达成。 人工势能场通常包括两部分：势能函数和机器人动力学模型。势能函数通常...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值