参考:https://www.cnblogs.com/H-Vking/p/4297973.html

题意:

虽然是比较裸的线段树,但是比较麻烦,并且有很多细节需要考虑,对着别人的ac代码debug了一个晚上。纪念一下

 这道题坑在有三种询问:set , add , mul。所以lazy标记要有三个,如果三个标记同时出现的处理方法——当更新set操作时,就把add标记和mul标记全部取消;当更新mul操作时,如果当前节点add标记存在,就把add标记改为:add * mul。这样的话就可以在PushDown()操作中先执行set,然后mul,最后add。

  麻烦在有三种询问:和 , 平方和 , 立方和。对于set和mul操作来说,这三种询问都比较好弄。

  对于add操作,和的话就比较好弄,按照正常方法就可以;

  平方和这样来推:(a + c)2 = a2 + c2 + 2ac  , 即sum2[rt] = sum2[rt] + (r - l + 1) * c * c + 2 * sum1[rt] * c;

  立方和这样推:(a + c)3 = a3 + c3 + 3a(a2 + ac) , 即sum3[rt] = sum3[rt] + (r - l + 1) * c * c * c + 3 * c * (sum2[rt] + sum1[rt] * c);

  几个注意点:add标记取消的时候是置0,mul标记取消的时候是置1;在PushDown()中也也要注意取消标记,如set操作中取消add和mul,mul操作中更新add; 在add操作中要注意sum3 , sum2 , sum1的先后顺序,一定是先sum3 , 然后sum2 , 最后sum1; int容易爆,还是用LL要保险一点; 最后就是运算较多,不要漏掉东西。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <list>
#include <iterator> using namespace std;
#define lson l , mid , rt << 1
#define rson mid + 1 , r , rt << 1 | 1
#define debug(x) cerr << #x << " = " << x << "\n"; typedef long long ll;
const int maxn = 1e5+;
const int mdd = ; struct node {
int l,r;
ll sum1,sum2,sum3;
ll lazy1,lazy2,lazy3;
}st[maxn<<]; void allmdd(int rt)
{
st[rt].sum1 = st[rt].sum1%mdd;
st[rt].sum2 = st[rt].sum2%mdd;
st[rt].sum3 = st[rt].sum3%mdd;
} void pushup(int rt)
{
st[rt].sum1 = (st[rt<<].sum1 + st[rt<<|].sum1)%mdd;
st[rt].sum2 = (st[rt<<].sum2 + st[rt<<|].sum2)%mdd;
st[rt].sum3 = (st[rt<<].sum3 + st[rt<<|].sum3)%mdd;
} void pushdown(int rt)
{
int len = st[rt].r - st[rt].l + ;
int llen = len - (len>>);
int rlen = len>>;
if(st[rt].lazy3)
{
ll c = st[rt].lazy3%mdd;
st[rt<<].lazy3 = st[rt<<|].lazy3 = st[rt].lazy3;
st[rt<<].lazy1 = st[rt<<|].lazy1 = ;
st[rt<<].lazy2 = st[rt<<|].lazy2 = ; st[rt<<].sum1 = (llen * (c%mdd))%mdd;
st[rt<<|].sum1 = (rlen * c)%mdd; st[rt<<].sum2 = (llen * (c * c)%mdd ) %mdd;
st[rt<<|].sum2 = ((rlen * ((c * c)%mdd))%mdd)%mdd; st[rt<<].sum3 = (llen * ((c * c)%mdd * (c % mdd))%mdd)%mdd;
st[rt<<|].sum3 = (rlen * ((c * c)%mdd * (c % mdd))%mdd)%mdd;
st[rt].lazy3 = ;
} if(st[rt].lazy2 != )
{
ll c = st[rt].lazy2 % mdd;
st[rt<<].lazy2 = (st[rt<<].lazy2 * c) % mdd;
st[rt<<|].lazy2 = (st[rt<<|].lazy2 * c) % mdd;
if(st[rt<<].lazy1){
st[rt<<].lazy1 = (st[rt<<].lazy1 * c)%mdd;
}
if(st[rt<<|].lazy1){
st[rt<<|].lazy1 = (st[rt<<|].lazy1 * c)%mdd;
}
st[rt<<].sum1 = (st[rt<<].sum1 * c)%mdd;
st[rt<<|].sum1 = (st[rt<<|].sum1 * c)%mdd; st[rt<<].sum2 = ((st[rt<<].sum2 * c%mdd)*c)%mdd;
st[rt<<|].sum2 = ((st[rt<<|].sum2 * c%mdd)*c)%mdd; st[rt<<].sum3 = ((((st[rt<<].sum3 * c)%mdd)*c%mdd)*c)%mdd;
st[rt<<|].sum3 = ((((st[rt<<|].sum3 * c)%mdd)*c%mdd)*c)%mdd;
st[rt].lazy2 = ;
} if(st[rt].lazy1)
{
ll c = st[rt].lazy1;
st[rt<<].lazy1 = (st[rt<<].lazy1 + c)%mdd;
st[rt<<|].lazy1 = (st[rt<<|].lazy1 +c)%mdd;
//注意sum3,sum2,sum1 的顺序;
st[rt<<].sum3 = st[rt<<].sum3 + ((c%mdd)*c%mdd)*c%mdd*llen + *c*((st[rt<<].sum2 + st[rt<<].sum1 * c)%mdd);
st[rt<<|].sum3 = st[rt<<|].sum3 + ((c%mdd)*c%mdd)*c%mdd*rlen +*c*((st[rt<<|].sum2 + st[rt<<|].sum1 * c)%mdd); st[rt<<].sum2 = st[rt<<].sum2 + (((llen *c)%mdd) * c)%mdd + *(st[rt<<].sum1)%mdd*c;
st[rt<<|].sum2 = st[rt<<|].sum2 + (((rlen *c)%mdd) * c)%mdd + *(st[rt<<|].sum1*c)%mdd; st[rt<<].sum1 = st[rt<<].sum1 + (llen*c)%mdd;
st[rt<<|].sum1 = st[rt<<|].sum1 + (rlen*c)%mdd; allmdd(rt<<);allmdd(rt<<|);
st[rt].lazy1 = ;
}
} void build(int l,int r,int rt)
{
st[rt].sum1 = st[rt].sum2 = st[rt].sum3 = ;
st[rt].lazy1 = st[rt].lazy3 = ;
st[rt].lazy2 = ;
st[rt].l = l,st[rt].r = r;
int mid = (l + r) >> ;
if(l == r)return;
build(lson);
build(rson);
pushup(rt);
}
void update(int l,int r,int rt,int L,int R,int c,int op)
{
if(l >= L&&r <= R)
{
if(op==)
{
st[rt].lazy1 += c;
st[rt].sum3 = (st[rt].sum3 + ((c*c%mdd)*c)%mdd*(r-l+)%mdd + (*c*(st[rt].sum2 + st[rt].sum1*c)%mdd))%mdd;
st[rt].sum2 = (st[rt].sum2 + (c*c)%mdd*(r-l+)%mdd + ( * st[rt].sum1 * c)%mdd)%mdd;
st[rt].sum1 = (st[rt].sum1 + (r-l+)*c%mdd)%mdd;
}
else if(op==)
{
st[rt].lazy2 = (st[rt].lazy2 * c)%mdd;
if(st[rt].lazy1)
st[rt].lazy1 = (st[rt].lazy1 * c)%mdd;
st[rt].sum1 = (st[rt].sum1 * c %mdd)%mdd;
st[rt].sum2 = (st[rt].sum2 * ((c * c) %mdd))%mdd;
st[rt].sum3 = (st[rt].sum3 * (((c * c)%mdd * c) % mdd))%mdd;
}
else if(op==)
{
st[rt].lazy3 = c;
st[rt].lazy1 = ;
st[rt].lazy2 = ;
st[rt].sum1 = (((r-l+)*c )% mdd);
st[rt].sum2 = (((((r-l+)* c)%mdd )* c %mdd))%mdd;
st[rt].sum3 = ((r-l+)* (((c * c)%mdd * c) % mdd));
}
return;
} int mid = (l+r)>>;
pushdown(rt);
if(L<=mid)update(lson,L,R,c,op);
if(R>mid)update(rson,L,R,c,op);
pushup(rt);
}
ll query(int l,int r,int rt,int L,int R,int c)
{
if(l>=L&&r<=R)
{
if(c == )return st[rt].sum1%mdd;
else if(c == )return st[rt].sum2%mdd;
else return st[rt].sum3%mdd;
}
int mid = (l+r)>>;
ll ans = ;
pushdown(rt);
if(L<=mid)ans = (ans + query(lson,L,R,c))%mdd;
if(R>mid)ans = (ans + query(rson,L,R,c))%mdd;
return ans;
}
int main(){
// freopen("input","r",stdin);
int n,m;
while(scanf("%d%d", &n, &m)==&&n+m)
{
// debug(n);
build(,n,);
while(m--)
{
int op;
int a,b,c;
scanf("%d%d%d%d", & op, & a,& b, & c);
if(op==){
printf("%lld\n", query(,n,,a,b,c)%mdd);
}
else {
update(,n,,a,b,c,op);
}
}
}
return ;
}

HDU4578 线段树(区间更新 + 多种操作)和平方,立方的更多相关文章

  1. HDU4578 线段树(区间更新 + 多种操作)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4578  , 线段树的区间更新 + 多种操作,好题. 虽然是比较裸的线段树,但是比较麻烦,并且有很多细节 ...

  2. hdu4578线段树区间更新

    /* 只有在区间中的数字不相同时才pushdown:往子区间传递数字再到子区间更新,同时该区间的flag置0 更新完左右子区间后进行pushup,如果左右子区间数字相同,那么把子区间合并,子区间数字置 ...

  3. hihoCoder 1080 : 更为复杂的买卖房屋姿势 线段树区间更新

    #1080 : 更为复杂的买卖房屋姿势 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho都是游戏迷,“模拟都市”是他们非常喜欢的一个游戏,在这个游戏里面他们 ...

  4. HDU 5023 A Corrupt Mayor's Performance Art(线段树区间更新)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5023 解题报告:一面墙长度为n,有N个单元,每个单元编号从1到n,墙的初始的颜色是2,一共有30种颜色 ...

  5. HDU 4902 Nice boat 2014杭电多校训练赛第四场F题(线段树区间更新)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4902 解题报告:输入一个序列,然后有q次操作,操作有两种,第一种是把区间 (l,r) 变成x,第二种是 ...

  6. HDU5039--Hilarity DFS序+线段树区间更新 14年北京网络赛

    题意:n个点的树,每个条边权值为0或者1, q次操作 Q 路径边权抑或和为1的点对数, (u, v)(v, u)算2个. M i修改第i条边的权值 如果是0则变成1, 否则变成0 作法: 我们可以求出 ...

  7. hihoCoder #1078 : 线段树的区间修改(线段树区间更新板子题)

    #1078 : 线段树的区间修改 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 对于小Ho表现出的对线段树的理解,小Hi表示挺满意的,但是满意就够了么?于是小Hi将问题 ...

  8. hdu 3966(树链剖分+线段树区间更新)

    传送门:Problem 3966 https://www.cnblogs.com/violet-acmer/p/9711441.html 学习资料: [1]线段树区间更新:https://blog.c ...

  9. POJ 3468:A Simple Problem with Integers(线段树区间更新模板)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 141093 ...

随机推荐

  1. 对Rust所有权、借用及生命周期的理解

    Rust的内存管理中涉及所有权.借用与生命周期这三个概念,下面是个人的一点粗浅理解. 一.从内存安全的角度理解Rust中的所有权.借用.生命周期 要理解这三个概念,你首要想的是这么做的出发点是什么-- ...

  2. python基础之条件控制

    行与缩进 python使用缩进来表示代码块,不需要使用大括号 {} ,同一个代码块的语句必须包含相同的缩进空格数. 条件控制 #条件控制语法: if 条件: 代码块 elif 条件: 代码块 else ...

  3. Scala集合(四)

    1. 集合 集合主要有三种: Sequence Map Set sequence是一种线性元素的集合,可能会是索引或者线性的(链表).map是包含键值对的集合,就像Java的Map,set是包含无重复 ...

  4. HTML/CSS:图片居中(水平居中和垂直居中)

    css图片居中(水平居中和垂直居中) css图片居中分css图片水平居中和垂直居中两种情况,有时候还需要图片同时水平垂直居中, 下面分几种居中情况分别介绍: css图片水平居中 1.利用margin: ...

  5. java实现发短信功能---腾讯云短信

    目录 java实现发短信功能 前言 开发环境 腾讯云 ---短信 代码 效果 结束语 java实现发短信功能 前言 如今发短信功能已经成为互联网公司的标配,本篇文章将一步步实现java发送短信 考察了 ...

  6. 世界十大OTA公司盘点

    世界十大OTA公司盘点 文/刘照慧(执惠旅游联合创始人,首发百度百家) 全球在线旅游公司(OTA)经过多年发展,已经形成较为成熟的商业模式,各大巨头跑马圈地,格局初现, 这两篇文章就梳理出全球按市值( ...

  7. 1.2模板templates

    一.模板使用 1. 配置模板目录 如果命令行创建的项目,需要手动配置模板文件目录,如果是Pycharm创建的项目,则无需配置 在项目根目录下创建模板目录,比如叫 templates,后续开发模板文件会 ...

  8. 电脑查询pico的mac

    配置好adb或者sdk后,  adb shell cat /sys/class/net/wlan0/address

  9. Powered by .NET Core 进展:第5次发布尝试(Windows部署)

    (图注:Windows 自带的性能监控,红色表示 CPU 占用,绿色表示 QPS) 今天中午 12:30 左右,我们进行了 .NET Core 博客站点的第5次发布(页脚会显示"Powere ...

  10. loging日志的使用

    2.日志: 记住怎么使用就好了 自己定义日志开始 import logging logger = logging.getLogger() # 创建一个logger fh = logging.FileH ...