线段树区间更新+区间求和模板(数组实现)洛谷p3372,p3373
模板题目1:https://www.luogu.org/problemnew/show/P3372
懒惰标记讲解:https://www.cnblogs.com/wushengyang/p/11194456.html
这题目只需要用一个懒惰标记,不用考虑顺序。
AC代码:
#include <stdio.h>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <string>
#include <cstring>
#include <algorithm>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <ctime>
#include <vector>
#include <fstream>
#include <list>
#include <iomanip>
#include <numeric>
using namespace std;
typedef long long int ll;
typedef unsigned long long int ull;
const int inf = 0x3f3f3f3f;
int dir[][]={{,},{,},{,},{,-},{-,},{-,-},{,-},{-,}};
#define pi acos(-1)
#define me0(s) memset(s,0,sizeof(s))
#define me1(s) memset(s,1,sizeof(s))
#define mef(s) memset(s,-1,sizeof(s))
#define meinf(s) memset(s,inf,sizeof(s))
const int N=;
#define ls rt<<1
#define rs rt<<1|1
ll a[N],sum[N*],lazy[N*];
void push_up(ll rt){
sum[rt]=sum[ls]+sum[rs]; //更新和
}
void push_down(ll m,ll rt){ //m为总区间长度
if(lazy[rt]){ //如果懒惰标记存在
lazy[ls]+=lazy[rt]; //左子树懒惰标记等于其本身加父亲节点的值
lazy[rs]+=lazy[rt];
sum[ls]+=lazy[rt]*(m-(m>>)); //懒惰标记乘个数
sum[rs]+=lazy[rt]*(m>>);
lazy[rt]=;
}
}
void build(ll l,ll r,ll rt){ //建树
lazy[rt]=; //可以建树直接赋值0
if(l==r){
sum[rt]=a[l];
return ;
}
ll m=(l+r)>>;
build(l,m,ls);
build(m+,r,rs);
push_up(rt);
}
void update(ll L,ll R,ll c,ll l,ll r,ll rt){//区间更新[l,r],[L,R]代表总区间
if(l<=L&&r>=R){
lazy[rt]+=c; //更新懒惰标记
sum[rt]+=c*(R-L+); //更新节点存的sum值
return ;
}
push_down(R-L+,rt); //所有的节点都更新一遍lazy标记
ll m=(L+R)>>; //总区间的一半
if(l<=m) update(L,m,c,l,r,ls); //注意这里递归的是总区间为变量
if(r>m) update(m+,R,c,l,r,rs);
push_up(rt);
}
ll query(ll L,ll R,ll l,ll r,ll rt){ //[l,r]区间
if(l<=L&&r>=R) return sum[rt];
push_down(R-L+,rt);
ll m=(L+R)>>; //总区间的一半
ll ans=;
if(l<=m) ans+=query(L,m,l,r,ls); //这个递归的也是总区间变量
if(r>m) ans+=query(m+,R,l,r,rs);
return ans;
}
int main(int argc, char * argv[])
{
std::ios::sync_with_stdio(false);
int n,m;
cin>>n>>m;
for(int i=;i<=n;i++) cin>>a[i];
build(,n,);
while(m--){
ll op,x,y,k;
cin>>op;
if(op==){
cin>>x>>y>>k;
update(,n,k,x,y,);
}
if(op==){
cin>>x>>y;
cout<<query(,n,x,y,)<<endl;
}
}
return ;
}
模板题目2:https://www.luogu.org/problemnew/show/P3373
这道题目由于存在多种区间更改操作,所以懒惰标记之间有关系,比如如果先求加法的话给他加上一个懒惰标记,再做乘法的话以前的懒惰标记加的那一部分也同时需要做一次乘法,这样才能完整。
#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef unsigned long long int ull;
const int inf = 0x3f3f3f3f;
int dir[][]={{,},{,},{,},{,-},{-,},{-,-},{,-},{-,}};
#define pi acos(-1)
#define ls rt<<1
#define rs rt<<1|1
#define me0(s) memset(s,0,sizeof(s))
#define me1(s) memset(s,1,sizeof(s))
#define mef(s) memset(s,-1,sizeof(s))
#define meinf(s) memset(s,inf,sizeof(s))
const int N=;
ll mod;
ll a[N],sum[N*],lazya[N*],lazym[N*]; //lazya存的加的数,lazym存乘
void push_up(int rt){
sum[rt]=(sum[ls]+sum[rs])%mod;
}
void push_down(ll m,ll rt){
lazya[ls]=(lazya[rt]+lazya[ls]*lazym[rt])%mod;//乘法的分配率
lazya[rs]=(lazya[rt]+lazya[rs]*lazym[rt])%mod;
lazym[ls]=(lazym[rt]*lazym[ls])%mod; //懒惰标记直接乘
lazym[rs]=(lazym[rt]*lazym[rs])%mod;
sum[ls]=(lazya[rt]*(m-(m>>))+sum[ls]*lazym[rt])%mod;
sum[rs]=(lazya[rt]*(m>>)+sum[rs]*lazym[rt])%mod;
lazya[rt]=;
lazym[rt]=;
}
void build(ll l,ll r,ll rt){
lazym[rt]=; //初始化乘值
if(l==r){
sum[rt]=a[l];
return ;
}
ll m=(l+r)>>;
build(l,m,ls);
build(m+,r,rs);
push_up(rt);
}
void update_a(ll L,ll R,ll c,ll l,ll r,ll rt){//做加法
if(l<=L&&r>=R){
sum[rt]=(sum[rt]+(R-L+)*c)%mod;
lazya[rt]=(c+lazya[rt])%mod;
return ;
}
push_down(R-L+,rt);
ll m=(L+R)>>;
if(l<=m) update_a(L,m,c,l,r,ls);
if(r>m) update_a(m+,R,c,l,r,rs);
push_up(rt);
}
void update_m(ll L,ll R,ll c,ll l,ll r,ll rt){//做乘法
if(l<=L&&r>=R){
sum[rt]=(sum[rt]*c)%mod;
lazya[rt]=(c*lazya[rt])%mod;
lazym[rt]=(c*lazym[rt])%mod;
return ;
}
push_down(R-L+,rt);
ll m=(L+R)>>;
if(l<=m) update_m(L,m,c,l,r,ls);
if(r>m) update_m(m+,R,c,l,r,rs);
push_up(rt);
}
ll query(ll L,ll R,ll l,ll r,ll rt){
if(l<=L&&r>=R){
return sum[rt];
}
push_down(R-L+,rt);
ll m=(L+R)>>;
ll ans=;
if(l<=m) ans=(ans+query(L,m,l,r,ls))%mod;
if(r>m) ans=(ans+query(m+,R,l,r,rs))%mod;
return ans%mod;
}
int main(int argc, char * argv[])
{
std::ios::sync_with_stdio(false);
ll n,m,x,y,k,op;
cin>>n>>m>>mod;
for(int i=;i<=n;i++) cin>>a[i];
build(,n,);
while(m--){
cin>>op;
if(op==){
cin>>x>>y>>k;
update_m(,n,k,x,y,);
}
if(op==){
cin>>x>>y>>k;
update_a(,n,k,x,y,);
}
if(op==){
cin>>x>>y;
cout<<query(,n,x,y,)<<endl;
}
}
return ;
}
线段树区间更新+区间求和模板(数组实现)洛谷p3372,p3373的更多相关文章
- POJ.2299 Ultra-QuickSort (线段树 单点更新 区间求和 逆序对 离散化)
POJ.2299 Ultra-QuickSort (线段树 单点更新 区间求和 逆序对 离散化) 题意分析 前置技能 线段树求逆序对 离散化 线段树求逆序对已经说过了,具体方法请看这里 离散化 有些数 ...
- POJ.3321 Apple Tree ( DFS序 线段树 单点更新 区间求和)
POJ.3321 Apple Tree ( DFS序 线段树 单点更新 区间求和) 题意分析 卡卡屋前有一株苹果树,每年秋天,树上长了许多苹果.卡卡很喜欢苹果.树上有N个节点,卡卡给他们编号1到N,根 ...
- HDU.1394 Minimum Inversion Number (线段树 单点更新 区间求和 逆序对)
HDU.1394 Minimum Inversion Number (线段树 单点更新 区间求和 逆序对) 题意分析 给出n个数的序列,a1,a2,a3--an,ai∈[0,n-1],求环序列中逆序对 ...
- hdu 1166线段树 单点更新 区间求和
敌兵布阵 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submis ...
- hdu1166(线段树单点更新&区间求和模板)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1166 题意:中文题诶- 思路:线段树单点更新,区间求和模板 代码: #include <iost ...
- hdu2795(线段树单点更新&区间最值)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2795 题意:有一个 h * w 的板子,要在上面贴 n 条 1 * x 的广告,在贴第 i 条广告时要 ...
- hdu1394(枚举/树状数组/线段树单点更新&区间求和)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1394 题意:给出一个循环数组,求其逆序对最少为多少: 思路:对于逆序对: 交换两个相邻数,逆序数 +1 ...
- HDU 1166 敌兵布阵(线段树点更新区间求和裸题)
Problem Description C国的死对头A国这段时间正在进行军事演习,所以C国间谍头子Derek和他手下Tidy又开始忙乎了.A国在海岸线沿直线布置了N个工兵营地,Derek和Tidy的任 ...
- POJ 2892 Tunnel Warfare(线段树单点更新区间合并)
Tunnel Warfare Time Limit: 1000MS Memory Limit: 131072K Total Submissions: 7876 Accepted: 3259 D ...
- HDU 3308 LCIS(线段树单点更新区间合并)
LCIS Given n integers. You have two operations: U A B: replace the Ath number by B. (index counting ...
随机推荐
- mysql重点,表查询操作和多表查询
表单查询 1. 完整的查询语句语法 select distinct(* or 字段名 or 四则运算 )from 表名 where 条件 group by 条件 having 条件 order by ...
- Vue.js Ajax(axios)
Vue.js 2.0 版本推荐使用 axios 来完成 ajax 请求. Axios 是一个基于 Promise 的 HTTP 库,可以用在浏览器和 node.js 中. Github开源地址: ht ...
- 正则表达式r和re
# coding:utf-8 import re print 'a\ws' print r'a\nb' # r'': 一般用在正则表达式中,称为原始字符串,作用是将Python语法中的反斜杠转义给 取 ...
- HDU-1501-Zipper-字符串的dfs
Given three strings, you are to determine whether the third string can be formed by combining the ch ...
- vue-router路由跳转判断用户是否存在
router.beforeEach((to, from, next) => { //console.log("to:", (to)); //console.log(" ...
- Web开发之Tomcat&Servlet
<!doctype html>01 - JavaEE - Tomcat&Servlet figure:first-child { margin-top: -20px; } #wri ...
- springboot整合jpa和mybatis实现主从复制
百度多方参考终于配出我自己的了,以下仅供参考 参考https://www.cnblogs.com/cjsblog/p/9712457.html 代码 首先数据源配置 spring.datasource ...
- CVE-2016-0095提权漏洞分析
1 前言 瞻仰了k0shl和鹏哥 的漏洞分析,感慨万千,任重而道远. 2 系统环境和工具 windows 7 32旗舰版 windbg 3 poc 3.1poc复现 首先k0shl大佬给出的poc() ...
- OpenGL键盘交互响应事件
GLUT允许我们编写程序,在里面加入键盘输入控制,包括了普通键,和其他特殊键(如F1,UP).在这一章里我们将学习如何去检测哪个键被按下,可以从GLUT里得到些什么信息,和如何处理键盘输入. 处理 ...
- C++ 系列:Boost Thread 编程指南
转载自:http://www.cppblog.com/shaker/archive/2011/11/30/33583.html 作者: dozbC++ Boost Thread 编程指南0 前言1 创 ...