【bzoj 1798】[AHOI2009]维护序列
2016-08-21 02:07
218 查看
题目描述
老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成。 有长为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
题解:经典线段树题目,同时有两个标记,一个加法标记,一个乘法标记,每个标记维护的意义为:下面的子树中,要先把每一项都乘以乘法标记,再加上加法标记。
设序列A = {a1,a2,a3,…,an},如果每一项先乘以p1,则序列变为{p1*a1,p1*a2,p1*a3,...,p1*an},再加上p2,则序列变为{p1*a1+p2,p1*a2+p2,p1*a3+p2,...,p1*an+p2},再乘以p3,则序列变为{p1*p3*a1+p2*p3,p1*p3*a2+p2*p3,p1*p3*a3+p2*p3,...,p1*p3*an+p2*p3}。
由此可见,在添加标记或者下放标记合并时,若新加乘法标记,则原有的乘法标记和加法标记都乘以新加的乘法标记,若新加加法标记,则与前面的乘法标记无关,直接加在加法标记上,根据乘法结合律,此做法成立。
#include <bits/stdc++.h> using namespace std; #define lch (now<<1) #define rch ((now<<1)+1) #define mid ((l+r)>>1) #define fa (now>>1) #define len (r-l+1) const int maxn = 100010; typedef long long ll; struct node{ ll sum, oldsum, add, mul; }; node tree[maxn<<2]; int n, p, m; void update(int now){ if(now == 0) return; tree[fa].oldsum = tree[fa].sum; tree[fa].sum += (tree[now].sum - tree[now].oldsum); tree[fa].sum = ((tree[fa].sum%p)+p)%p; update(fa); } void build(int now, int l, int r){ if(l == r){ scanf("%lld", &tree[now].sum); tree[now].sum %= p; tree[now].mul = 1, tree[now].add = 0, tree[now].oldsum = 0; return; } build(lch, l, mid); build(rch, mid+1, r); tree[now].sum = tree[lch].sum + tree[rch].sum; tree[now].sum %= p; tree[now].add = 0, tree[now].mul = 1, tree[now].oldsum = 0; } void down(int now, int l, int r){ if(tree[now].add == 0 && tree[now].mul == 1) return; tree[lch].oldsum = tree[lch].sum; tree[rch].oldsum = tree[rch].sum; ll mul = tree[now].mul, add = tree[now].add; tree[lch].add *= mul, tree[rch].add *= mul; tree[lch].add %= p, tree[rch].add %= p; tree[lch].add += add, tree[rch].add += add; tree[lch].add %= p, tree[rch].add %= p; tree[lch].mul *= mul, tree[rch].mul *= mul; tree[lch].mul %= p, tree[rch].mul %= p; tree[lch].sum = ((tree[lch].sum * mul)%p + add*(mid-l+1)%p), tree[rch].sum = ((tree[rch].sum * mul)%p + add*(r-mid)%p); tree[lch].sum %= p, tree[rch].sum %= p; tree[now].add = 0, tree[now].mul = 1; } void modify(int now, int l, int r, int pos1, int pos2, int change, bool is){ if(l == pos1 && r == pos2){ if(is == 0){ tree[now].add += change; tree[now].add %= p; tree[now].oldsum = tree[now].sum; tree[now].sum += len*change; tree[now].sum %= p; update(now); } else{ tree[now].add *= change, tree[now].mul *= change; tree[now].add %= p, tree[now].mul %= p; tree[now].oldsum = tree[now].sum; tree[now].sum *= change; tree[now].sum %= p; update(now); } return; } down(now,l,r); if(pos2 <= mid) modify(lch,l,mid,pos1,pos2,change,is); else if(pos1 >= mid+1) modify(rch,mid+1,r,pos1,pos2,change,is); else{ modify(lch,l,mid,pos1,mid,change,is); modify(rch,mid+1,r,mid+1,pos2,change,is); } return; } int que(int now, int l, int r, int pos1, int pos2){ if(l == pos1 && r == pos2){ return tree[now].sum; } down(now,l,r); if(pos2 <= mid) return que(lch, l, mid, pos1, pos2) % p; else if(pos1 >= mid+1) return que(rch, mid+1, r, pos1, pos2) % p; else return (que(lch, l, mid, pos1, mid) + que(rch, mid+1, r, mid+1, pos2)) % p; } int main(){ scanf("%d%d", &n, &p); build(1,1,n); int q; scanf("%d", &q); for(int i = 1; i <= q; i ++){ int t; scanf("%d", &t); if(t == 1){ int pos1, pos2, change; scanf("%d%d%d", &pos1, &pos2, &change); change %= p; modify(1,1,n,pos1, pos2, change, 1); }else if(t == 2){ int pos1, pos2, change; scanf("%d%d%d", &pos1, &pos2, &change); change %= p; modify(1,1,n,pos1, pos2, change, 0); }else{ int pos1, pos2; scanf("%d%d", &pos1, &pos2); printf("%d\n", ((que(1,1,n,pos1,pos2)%p)+p)%p); } } return 0; }
20160821 by yhf_2015
相关文章推荐
- bzoj 1798: [Ahoi2009]Seq 维护序列seq 线段树
- BZOJ 1798: [Ahoi2009]Seq 维护序列seq (线段树乘法加法的混合操作)
- _bzoj1798 [Ahoi2009]Seq 维护序列seq【线段树 lazy tag】
- [BZOJ1798] [Ahoi2009]Seq 维护序列seq
- BZOJ 1798 [Ahoi2009]维护序列seq (线段树)
- [BZOJ 1798][Ahoi2009]Seq 维护序列seq
- 【线段树】【BZOJ1798】【AHOI2009】维护序列
- bzoj1798: [Ahoi2009]Seq 维护序列seq(线段树)
- 线段树 BZOJ1798 [Ahoi2009]Seq 维护序列seq
- 【AHOI2009】bzoj1798 维护序列
- BZOJ1798 [Ahoi2009]Seq 维护序列seq
- [BZOJ]1798: [Ahoi2009]Seq 维护序列seq
- bzoj 1798 [Ahoi2009]Seq 维护序列seq
- 【bzoj1798】[Ahoi2009]Seq 维护序列seq
- bzoj 1798: [Ahoi2009]Seq 维护序列seq 双标记线段树
- 【BZOJ1798】【Ahoi2009】维护序列seq(线段树)
- [bzoj1798][Ahoi2009]Seq 维护序列seq
- BZOJ1798:[Ahoi2009]Seq 维护序列 线段树
- bzoj1798: [Ahoi2009]Seq 维护序列seq
- BZOJ1798 【AHOI2009】 seq维护序列 <线段树>