BZOJ_1798_[AHOI2009]维护序列_线段树
BZOJ_1798_[AHOI2009]维护序列_线段树
题意:老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成。 有长为N的数列,不妨设为a1,a2,…,aN 。有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2)把数列中的一段数全部加一个值; (3)询问数列中的一段数的和,由于答案可能很大,你只需输出这个数模P的值。
分析:线段树上要打两个标记。要注意下传的顺序。显然先乘后加和先加后乘是不一样的。我们发现如果是先加后乘的话更改子树值的式子里会出现除法。不妨规定任何时候都先乘后加。推出的式子即为
t[lson]=(t[lson]*mul[pos]+add[pos]*(mid-l+1))%p;
mul[lson]=(mul[lson]*mul[pos])%p;
add[lson]=(add[lson]*mul[pos]+add[pos])%p;
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define ls p<<1
#define rs p<<1|1
#define LL long long
#define N 100050
LL add[N<<2],t[N<<2],mul[N<<2],mod;
int a[N],n,m;
void bt(int l,int r,int p){
mul[p]=1;
if(l==r){
scanf("%lld",&t[p]);return ;
}
int mid=l+r>>1;
bt(l,mid,ls);bt(mid+1,r,rs);
t[p]=(t[ls]+t[rs])%mod;
}
void pud(int l,int r,int p){
int mid=l+r>>1;
if(add[p]==0&&mul[p]==1)return ;
t[ls]=(t[ls]*mul[p]+add[p]*(mid-l+1))%mod;
t[rs]=(t[rs]*mul[p]+add[p]*(r-mid))%mod;
mul[ls]=mul[ls]*mul[p]%mod;
mul[rs]=mul[rs]*mul[p]%mod;
add[ls]=(add[ls]*mul[p]+add[p])%mod;
add[rs]=(add[rs]*mul[p]+add[p])%mod;
mul[p]=1;add[p]=0;
}
void upad(int l,int r,int x,int y,int c,int p){
if(x<=l&&y>=r){
add[p]=(add[p]+c)%mod;
t[p]+=1ll*(r-l+1)*c;t[p]%=mod;
return ;
}
int mid=l+r>>1;
pud(l,r,p);
if(x<=mid)upad(l,mid,x,y,c,ls);
if(y>mid)upad(mid+1,r,x,y,c,rs);
t[p]=(t[ls]+t[rs])%mod;
}
void upmu(int l,int r,int x,int y,int c,int p){
if(x<=l&&y>=r){
mul[p]=mul[p]*c%mod;
add[p]=add[p]*c%mod;
t[p]=t[p]*c%mod;
return ;
}
pud(l,r,p);
int mid=l+r>>1;
if(x<=mid)upmu(l,mid,x,y,c,ls);
if(y>mid)upmu(mid+1,r,x,y,c,rs);
t[p]=(t[ls]+t[rs])%mod;
}
LL query(int l,int r,int x,int y,int p){
if(x<=l&&y>=r)return t[p];
int mid=l+r>>1;
LL re=0;
pud(l,r,p);
if(x<=mid)re=(re+query(l,mid,x,y,ls))%mod;
if(y>mid)re=(re+query(mid+1,r,x,y,rs))%mod;
return re;
}
int main(){
scanf("%d%lld",&n,&mod);
bt(1,n,1);
scanf("%d",&m);
int op,x,y,z;
for(int i=1;i<=m;i++){
scanf("%d",&op);
if(op==1){
scanf("%d%d%d",&x,&y,&z);
upmu(1,n,x,y,z,1);
}else if(op==2){
scanf("%d%d%d",&x,&y,&z);
upad(1,n,x,y,z,1);
}else{
scanf("%d%d",&x,&y);
printf("%lld\n",query(1,n,x,y,1));
}
}
}
BZOJ_1798_[AHOI2009]维护序列_线段树的更多相关文章
- 洛谷P2023 [AHOI2009]维护序列(线段树区间更新,区间查询)
洛谷P2023 [AHOI2009]维护序列 区间修改 当我们要修改一个区间时,要保证 \(ax+b\) 的形式,即先乘后加的形式.当将区间乘以一个数 \(k\) 时,原来的区间和为 \(ax+b\) ...
- 1798. [AHOI2009]维护序列【线段树】
Description 老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成. 有长为N的数列,不妨设为a1,a2,…,aN .有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2 ...
- 【题解】洛谷P2023 [AHOI2009] 维护序列(线段树)
洛谷P2023:https://www.luogu.org/problemnew/show/P2023 思路 需要2个Lazy-Tag 一个表示加的 一个表示乘的 需要先计算乘法 再计算加法 来自你谷 ...
- BZOJ 1798: [Ahoi2009]Seq 维护序列seq( 线段树 )
线段树.. 打个 mul , add 的标记就好了.. 这个速度好像还挺快的...( 相比我其他代码 = = ) 好像是#35.. ---------------------------------- ...
- bzoj 1798: [Ahoi2009]Seq 维护序列seq 线段树 区间乘法区间加法 区间求和
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeO ...
- Bzoj 1798: [Ahoi2009]Seq 维护序列seq(线段树区间操作)
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MB Description 老师交给小可可一个维护数列的任务,现在小可 ...
- BZOJ1798: [Ahoi2009]Seq 维护序列seq[线段树]
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MBSubmit: 5504 Solved: 1937[Submit ...
- bzoj 1798: [Ahoi2009]Seq 维护序列seq (线段树 ,多重标记下放)
1798: [Ahoi2009]Seq 维护序列seq Time Limit: 30 Sec Memory Limit: 64 MBSubmit: 7773 Solved: 2792[Submit ...
- 【BZOJ】1798: [Ahoi2009]Seq 维护序列seq 线段树多标记(区间加+区间乘)
[题意]给定序列,支持区间加和区间乘,查询区间和取模.n<=10^5. [算法]线段树 [题解]线段树多重标记要考虑标记与标记之间的相互影响. 对于sum*b+a,+c直接加上即可. *c后就是 ...
随机推荐
- 高并发教程-基础篇-之nginx负载均衡的搭建
温馨提示:请不要盲目的进行横向扩展,优先考虑对单台服务器的性能优化,只有单台服务器的性能达到最优化之后,集群才会被最大的发挥作用. 一.架构图: 服务器准备:3台,ubuntu16.04系统maste ...
- 深入理解springMVC思想
转载:http://elf8848.iteye.com/blog/875830 深入理解Spring MVC 思想 目录 一.前言二.spring mvc 核心类与接口三.spring mvc ...
- Android Selector原理
android的selector对于android开发者而言再熟悉不过了,只要定义一个drawable目录下定义一个selector的xml文件,在布局文件中引用这个xml文件或者在代码中setBac ...
- Day 网络协议介绍 简单通信的实现
比如说实现两个手机之间的通信,需要做的几部: 服务端: 1,买手机 2,插卡 3,开机 4,等电话链接 5,基于建立的链接,收发协议 6,挂电话 7,关机 import socket #买手机 pho ...
- currval of sequence "follow_id_seq" is not yet defined in this session
postgresql上使用 select currval('follow_id_seq'); 报错: currval of sequence "follow_id_seq" is ...
- HttpDNS的坑以及一个针对安卓不太完善的测试方案
背景:单位因为域名劫持(具体表象是某个地区的用户ping不通域名或者因为DNS解析的ip跨网段导致访问速度很慢)需要运维经常去定位,于是提出了httpDNS方案. 想法是美好的,现实是残酷的.没引入这 ...
- C语言实现输出一组数字中的所有奇数
/*第二题*/ #include<stdio.h> //输入186732468 //输出173 //输入12345677 //输出13577 main(){ ;//输入的数字,数字的长度 ...
- BigDecimal常用的加减乘除算法、比较大小、保存两位小数点
项目中涉及到了BigDecimal的加.减.乘.比较大小.精确度的问题.所以在此总结一下,方便以后复习. //加法 BigDecimal coins = new BigDecimal("0& ...
- Flex 将默认日期格式转化成通用格式
flex 默认日期格式如:Wed Dec 16 00:00:00 GMT+0800 2015 想要得到的通用格式如:2015-12-16 转换方法如下: var sdate:String = &quo ...
- Visual Studio 和 c# 正则表达式
今天集中说说VS生产环境下的正则. Visual Sturdio 2012以上版本查找替换 对于VS的正则,准确说,是VS2012之后的IDE下VS的正则. VS的查找和替换功能支持基础的正则表达式, ...