1798: [Ahoi2009]Seq 维护序列seq
Time Limit: 30 Sec Memory Limit: 64 MB
Submit: 6564 Solved: 2341
[Submit][Status][Discuss]
Description
老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成。 有长为N的数列,不妨设为a1,a2,…,aN 。有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2)把数列中的一段数全部加一个值; (3)询问数列中的一段数的和,由于答案可能很大,你只需输出这个数模P的值。
Input
第一行两个整数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)。 同一行相邻两数之间用一个空格隔开,每行开头和末尾没有多余空格。
Output
对每个操作3,按照它在输入中出现的顺序,依次输出一行一个整数表示询问结果。
Sample Input
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
Sample Output
2
35
8
HINT
【样例说明】
初始时数列为(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
分析:就是线段树的裸题。注意维护两个标记,一个是乘法的,一个是加法的。注意优先级的关系。注意下传乘法标记的时候加分标记也要乘一个乘法标记的那个值。和要分为乘法的部分和加法的部分。反正就是注意一下下传标记的时候标记之间的关系以及和之间的关系。最恶心的就是取模!!!太恶心了!!!
#include<cstdio>
#include<cstring>
#include<ctime>
#include<algorithm>
#define ll long long
using namespace std;
const int N = 1000010;
int n,p,m;
int a[N];
int root;
struct node{
ll sum;
int l,r;
ll flag;//加法标记
ll lazy;//乘法标记
}t[N];
void build(int root,int l,int r){
t[root].l=l,t[root].r=r,t[root].flag=0;t[root].lazy=1;
if(l==r){
t[root].sum=a[l],
t[root].flag=0;
t[root].lazy=1;
return ;
}
int mid=l+r>>1;
build(root<<1,l,mid);
build(root<<1|1,mid+1,r);
t[root].sum=(t[root<<1].sum+t[root<<1|1].sum)%p;
t[root].flag=0,t[root].lazy=1;
}
void pushdown(int root){
if(t[root].flag==0&&t[root].lazy==1) return ;
int l=t[root].l,r=t[root].r;
int mid=(l+r)>>1;
ll lazy=t[root].lazy,flag=t[root].flag;
lazy%=p,flag%=p;
t[root<<1].lazy=(t[root<<1].lazy*lazy)%p;
t[root<<1].flag=(t[root<<1].flag*lazy%p+flag)%p;
t[root<<1|1].lazy=(t[root<<1|1].lazy*lazy)%p;
t[root<<1|1].flag=(t[root<<1|1].flag*lazy%p+flag)%p;
t[root<<1].sum=(((t[root<<1].sum*lazy)%p+(flag*(mid-l+1)%p)%p)%p)%p;
t[root<<1|1].sum=(((t[root<<1|1].sum*lazy)%p+(flag*(r-mid)%p)%p)%p)%p;
t[root].flag=0;t[root].lazy=1;
}
void modify1(int root,int pos,int val,int delta){
int l=t[root].l,r=t[root].r;
if(pos<=l&&val>=r){
t[root].sum=(t[root].sum*delta)%p;
t[root].lazy=(t[root].lazy*delta)%p;
t[root].flag=(t[root].flag*delta)%p;
return ;
}
int mid=(l+r)>>1;
pushdown(root);
if(pos<=mid) modify1(root<<1,pos,val,delta);
if(val>mid) modify1(root<<1|1,pos,val,delta);
t[root].sum=(t[root<<1].sum%p+t[root<<1|1].sum%p)%p;
}
void modify2(int root,int pos,int val,int delta){
int l=t[root].l,r=t[root].r;
if(pos<=l&&val>=r){
t[root].sum=(t[root].sum+(delta*(r-l+1)%p)%p)%p;
t[root].flag=(t[root].flag+delta)%p;
return ;
}
int mid=(l+r)>>1;
pushdown(root);
if(pos<=mid) modify2(root<<1,pos,val,delta);
if(val>mid) modify2(root<<1|1,pos,val,delta);
t[root].sum=(t[root<<1].sum%p + t[root<<1|1].sum%p )%p;
}
ll query(int root,int pos,int val){
int l=t[root].l,r=t[root].r;
if(pos<=l&&val>=r) return t[root].sum%p;
ll ans=0;
int mid=(l+r)>>1;
pushdown(root);
if(pos<=mid) ans=(ans+query(root<<1,pos,val))%p;
if(val>mid) ans=(ans+query(root<<1|1,pos,val))%p;
t[root].sum=t[root<<1].sum%p+t[root<<1|1].sum%p;
return ans%p;
}
int main(){
n=1;
scanf("%d%d",&n,&p);
for(register int i=1;i<=n;i++) scanf("%d",&a[i]);
root=1;
build(root,1,n);
scanf("%d",&m);
while(m--){
int x;
scanf("%d",&x);
if(x==1){
int l,r,delta;
scanf("%d%d%d",&l,&r,&delta);
modify1(root,l,r,delta);
}
else if(x==2){
int l,r,delta;
scanf("%d%d%d",&l,&r,&delta);
modify2(root,l,r,delta);
}
else{
int l,r;
scanf("%d%d",&l,&r);
printf("%lld\n",query(root,l,r));
}
}
return 0;
}