[洛谷P2023] [AHOI2009]维护序列
洛谷题目链接:[AHOI2009]维护序列
题目描述
老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成。 有长为N的数列,不妨设为a1,a2,…,aN 。有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2)把数列中的一段数全部加一个值; (3)询问数列中的一段数的和,由于答案可能很大,你只需输出这个数模P的值。
输入输出格式
输入格式:
第一行两个整数N和P(1≤P≤1000000000)。 第二行含有N个非负整数,从左到右依次为a1,a2,…,aN, (0≤ai≤1000000000,1≤i≤N)。 第三行有一个整数M,表示操作总数。 从第四行开始每行描述一个操作,输入的操作有以下三种形式: 操作1:“1 t g c”(不含双引号)。表示把所有满足t≤i≤g的ai改为ai×c (1≤t≤g≤N,0≤c≤1000000000)。 操作2:“2 t g c”(不含双引号)。表示把所有满足t≤i≤g的ai改为ai+c (1≤t≤g≤N,0≤c≤1000000000)。 操作3:“3 t g”(不含双引号)。询问所有满足t≤i≤g的ai的和模P的值 (1≤t≤g≤N)。 同一行相邻两数之间用一个空格隔开,每行开头和末尾没有多余空格。
输出格式:
对每个操作3,按照它在输入中出现的顺序,依次输出一行一个整数表示询问结果。
输入输出样例
输入样例#1:
7 43
1 2 3 4 5 6 7
5
1 2 5 5
3 2 4
2 3 7 9
3 1 3
3 4 7
输出样例#1:
2
35
8
说明
【样例说明】
初始时数列为(1,2,3,4,5,6,7)。
经过第1次操作后,数列为(1,10,15,20,25,6,7)。
对第2次操作,和为10+15+20=45,模43的结果是2。
经过第3次操作后,数列为(1,10,24,29,34,15,16}
对第4次操作,和为1+10+24=35,模43的结果是35。
对第5次操作,和为29+34+15+16=94,模43的结果是8。
测试数据规模如下表所示
数据编号 1 2 3 4 5 6 7 8 9 10
N= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000
M= 10 1000 1000 10000 60000 70000 80000 90000 100000 100000
Source: Ahoi 2009
一句话题意: 写一中数据结构支持区间加法/乘法,以及区间求和.
题解: 维护方式多种多样,这里使用的是无旋treap.当然这题主要是来练习标记下放的.可以维护一个乘法标记,一个加法标记,在下放的时候要先下放乘法标记. 为什么呢?假设先有一个加法标记打在这个区间上,那么直接将加法标记也乘上这个修改的值就可以了.假设先有乘法标记,那么也会先将乘法标记下放,所以是没有问题的.
还有就是修改的时候下放到左右儿子的标记的值都只需要用这个点的标记值来修改,比如说乘法:
应该是
(t[ls].sum *= t[x].mul) %= mod; (t[rs].sum *= t[x].mul) %= mod;
而不是
(t[ls].sum *= t[ls].mul) %= mod; (t[rs].sum *= t[ls].mul) %= mod;
因为之前如果有标记打在这个区间上,那么它就已经对这个区间的sum,val等信息都已经修改过了.
然后就是每个地方取mod都处理好,开个long long就可以了.
#include<bits/stdc++.h>
#define debug out(root), printf("\n")
using namespace std;
typedef long long lol;
const int N=100000+5;
lol n, m, mod, a[N], root, cnt = 0, r1, r2, r3;
struct treap{
lol ch[2], rd, size, val, sum, tag, mul;
}t[N];
inline lol gi(){
lol ans = 0, f = 1; char i = getchar();
while(i<'0'||i>'9'){if(i=='-')f=-1;i=getchar();}
while(i>='0'&&i<='9'){ans=ans*10+i-'0';i=getchar();}
return ans * f;
}
inline lol newnode(lol val){
t[++cnt].rd = rand(); t[cnt].size = t[cnt].mul = 1;
t[cnt].val = t[cnt].sum = val; t[cnt].tag = 0;
return cnt;
}
inline void up(lol x){
lol ls = t[x].ch[0], rs = t[x].ch[1];
t[x].sum = (t[ls].sum+t[x].val+t[rs].sum)%mod;
t[x].size = t[ls].size+t[rs].size+1;
}
inline void pushdown(lol x){
lol ls = t[x].ch[0], rs = t[x].ch[1];
if(t[x].mul != 1){
(t[ls].mul *= t[x].mul) %= mod; (t[rs].mul *= t[x].mul) %= mod;
(t[ls].tag *= t[x].mul) %= mod; (t[rs].tag *= t[x].mul) %= mod;
(t[ls].val *= t[x].mul) %= mod; (t[rs].val *= t[x].mul) %= mod;
(t[ls].sum *= t[x].mul) %= mod; (t[rs].sum *= t[x].mul) %= mod;
t[x].mul = 1;
}
if(t[x].tag){
(t[ls].tag += t[x].tag) %= mod; (t[rs].tag += t[x].tag) %= mod;
(t[ls].val += t[x].tag) %= mod; (t[rs].val += t[x].tag) %= mod;
(t[ls].sum += (t[x].tag*t[ls].size)%mod) %= mod;
(t[rs].sum += (t[x].tag*t[rs].size)%mod) %= mod;
t[x].tag = 0;
}
}
inline void split(lol x, lol k, lol &a, lol &b){
if(x == 0){a = b = 0; return;} pushdown(x);
if(k <= t[t[x].ch[0]].size) b = x, split(t[x].ch[0], k, a, t[x].ch[0]);
else a = x, split(t[x].ch[1], k-t[t[x].ch[0]].size-1, t[x].ch[1], b); up(x);
}
inline lol merge(lol x, lol y){
if(x == 0 || y == 0) return x+y;
pushdown(x), pushdown(y);
if(t[x].rd < t[y].rd){
t[x].ch[1] = merge(t[x].ch[1], y);
up(x); return x;
} else {
t[y].ch[0] = merge(x, t[y].ch[0]);
up(y); return y;
}
}
inline void insert(lol val, lol pos){
split(root, pos, r1, r2);
root = merge(r1, merge(newnode(val), r2));
}
inline lol query(lol x, lol y){
lol res;
split(root, y, r2, r3); split(r2, x-1, r1, r2);
res = t[r2].sum; root = merge(r1, merge(r2, r3));
return res;
}
inline void update(lol x, lol y, lol z){
z %= mod; split(root, y, r2, r3); split(r2, x-1, r1, r2);
(t[r2].tag += z) %= mod; (t[r2].val += z) %= mod;
(t[r2].sum += z*t[r2].size) %= mod;
pushdown(r2); root = merge(r1, merge(r2, r3));
}
inline void multiply(lol x, lol y, lol z){
z %= mod; split(root, y, r2, r3); split(r2, x-1, r1, r2);
(t[r2].mul *= z) %= mod; (t[r2].tag *= z) %= mod;
(t[r2].val *= z) %= mod; (t[r2].sum *= z) %= mod;
pushdown(r2); root = merge(r1, merge(r2, r3));
}
int main(){
lol opt, x, y, z; n = gi(), mod = gi();
for(lol i=1; i<=n; i++) x = gi(), insert(x%mod, i-1);
m = gi();
for(lol i=1; i<=m; i++){
opt = gi(), x = gi(), y = gi();
if(opt == 1) z = gi(), multiply(x, y, z);
if(opt == 2) z = gi(), update(x, y, z);
if(opt == 3) printf("%lld\n", query(x, y));
}
return 0;
}
[洛谷P2023] [AHOI2009]维护序列的更多相关文章
- 洛谷P2023 [AHOI2009]维护序列(线段树区间更新,区间查询)
洛谷P2023 [AHOI2009]维护序列 区间修改 当我们要修改一个区间时,要保证 \(ax+b\) 的形式,即先乘后加的形式.当将区间乘以一个数 \(k\) 时,原来的区间和为 \(ax+b\) ...
- 洛谷 P2023 [AHOI2009]维护序列
P2023 [AHOI2009]维护序列 题目描述 老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成. 有长为N的数列,不妨设为a1,a2,…,aN .有如下三种操作形式: (1)把数列中 ...
- 洛谷 P2023 [AHOI2009]维护序列 题解
P2023 [AHOI2009]维护序列 题目描述 老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成. 有长为N的数列,不妨设为a1,a2,-,aN .有如下三种操作形式: (1)把数列中 ...
- 【题解】洛谷P2023 [AHOI2009] 维护序列(线段树)
洛谷P2023:https://www.luogu.org/problemnew/show/P2023 思路 需要2个Lazy-Tag 一个表示加的 一个表示乘的 需要先计算乘法 再计算加法 来自你谷 ...
- 洛谷 P2023 [AHOI2009]维护序列 || 线段树加法和乘法运算
原理倒是非常简单.设原数为x,加法的lazytag为b,乘法的lazytag为a,操作数为c,那么原式为ax+b,乘上c后(ax+b)c=(ac)*x+b*c,加上c后(ax+b)+c=ax+(b+c ...
- 洛谷 2023 [AHOI2009]维护序列
洛谷 2023 [AHOI2009]维护序列 洛谷原题传送门 这个题也是一道经典的线段树模版(其实洛谷的模版二改一下输入顺序就能AC),其中包括区间乘法修改.区间加法修改.区间查询三个操作. 线段树的 ...
- 洛谷P3373 【模板】线段树 2 && P2023 [AHOI2009]维护序列——题解
题目传送: P3373 [模板]线段树 2 P2023 [AHOI2009]维护序列 该题较传统线段树模板相比多了一个区间乘的操作.一提到线段树的区间维护问题,就自然想到了“懒标记”:为了降低时间复 ...
- P2023 [AHOI2009]维护序列 题解(线段树)
题目链接 P2023 [AHOI2009]维护序列 解题思路 线段树板子.不难,但是...有坑.坑有多深?一页\(WA\). 由于乘法可能乘\(k=0\),我这种做法可能会使结果产生负数.于是就有了这 ...
- [P2023][AHOI2009]维护序列(线段树)
题目描述 老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成. 有长为N的数列,不妨设为a1,a2,…,aN .有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2)把数列中的一 ...
随机推荐
- BZOJ 3166 HEOI2013 ALO 可持久化trie+st表
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3166(洛谷上也有) 题意概述: 给出一个序列,对于一个区间,其权值为区间中的次大值亦或区 ...
- VS2013运行C++报错:This function or variable may be unsafe. Consider using fopen_s instead.
在vs2013中运行时发生的关于方法调用的不安全错误. 1.更换方法,但是有些方法更改后参数不变,所以可能比较麻烦. 2.添加一条预处理器定义: 点击项目——>属性——>c/c++——&g ...
- 总结python 元组和列表的区别
python的基本类型中有元组和列表这么俩个,但是这哥俩却比较难于区分,今天就来用简单的实例说明两者的不同. 列表:1.使用中括号([ ])包裹,元素值和个数可变 实例: aaa = ['sitena ...
- 对Objective-C中runtime的理解
Objective-C是面向runtime(运行时)的语言,在应用程序运行的时候来决定函数内部实现什么以及做出其它决定的语言.程序员可以在程序运行时创建,检 查,修改类,对象和它们的方法,Object ...
- PagedDataSource数据绑定控件和AspNetPager分页控件结合使用列表分页
1.引用AspNetPager.dll. 2.放置Repeater数据绑定控件. <asp:Repeater ID="Repeater1" runat="serve ...
- WITH REPLACE 含义
RESTORE DATABASE db_CSharp from disk='backup.bak' WITH REPLACE WITH REPLACE后面是限定条件,with replace意思是替换 ...
- SSM整合步骤
第一步:mybatis和spring整合 mybatis-spring-1.2.2:是mybatis官方出的包: mybatis的包: mybatis和spring的整合包: spring及sprin ...
- BZOJ 1179 Atm(强连通分量缩点+DP)
题目说可以通过一条边多次,且点权是非负的,所以如果走到图中的一个强连通分量,那么一定可以拿完这个强连通分量上的money. 所以缩点已经很明显了.缩完点之后图就是一个DAG,对于DAG可以用DP来求出 ...
- bzoj2676 Contra
题意: 给定N,R,Q,S 有N个关卡,初始有Q条命,且任意时刻最多只能有Q条命 每通过一个关卡,会得到u分和1条命,其中u=min(最近一次连续通过的关数,R) 若没有通过这个关卡,将失去一条命,并 ...
- P1349 广义斐波那契数列
题目描述 广义的斐波那契数列是指形如an=p*an-1+q*an-2的数列.今给定数列的两系数p和q,以及数列的最前两项a1和a2,另给出两个整数n和m,试求数列的第n项an除以m的余数. 输入输出格 ...