BZOJ1798: [Ahoi2009]Seq 维护序列seq
题解:线段树 维护一次函数即可
/**************************************************************
Problem: 1798
User: wang9897
Language: C++
Result: Accepted
Time:5780 ms
Memory:14572 kb
****************************************************************/
#include <bits/stdc++.h>
#define N 100005
#define ll long long
using namespace std;
int read(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return f*x;
}
typedef struct node{
int l,r;
ll sum,flag1,flag2;
}node;
node d[N<<2];ll a[N];
int n;ll lth;
void push(int x){
ll len=(d[x<<1].r-d[x<<1].l+1);
d[x<<1].sum=((d[x<<1].sum*d[x].flag1)%lth+(d[x].flag2*len)%lth)%lth;
len=(d[x<<1|1].r-d[x<<1|1].l+1);
d[x<<1|1].sum=((d[x<<1|1].sum*d[x].flag1)%lth+(d[x].flag2*len)%lth)%lth;
d[x<<1].flag1=(d[x<<1].flag1*d[x].flag1)%lth;d[x<<1].flag2=(d[x<<1].flag2*d[x].flag1)%lth;
d[x<<1|1].flag1=(d[x<<1|1].flag1*d[x].flag1)%lth;d[x<<1|1].flag2=(d[x<<1|1].flag2*d[x].flag1)%lth;
d[x<<1].flag2=(d[x<<1].flag2+d[x].flag2)%lth;d[x<<1|1].flag2=(d[x<<1|1].flag2+d[x].flag2)%lth;
d[x].flag1=1;d[x].flag2=0;
}
void up(int x){
d[x].sum=(d[x<<1].sum+d[x<<1|1].sum)%lth;
}
void built(int root,int l,int r){
if(l==r){
d[root].l=d[root].r=l;d[root].sum=a[l]%lth;d[root].flag1=1;d[root].flag2=0;
return ;
}
int mid=(l+r)>>1;
built(root<<1,l,mid);
built(root<<1|1,mid+1,r);
d[root].l=d[root<<1].l;d[root].r=d[root<<1|1].r;d[root].flag1=1;d[root].flag2=0;
up(root);
}
void update1(int root,int l,int r,ll c){
if(l<=d[root].l&&d[root].r<=r){
d[root].sum=(d[root].sum*c)%lth;d[root].flag1=(d[root].flag1*c)%lth;d[root].flag2=(d[root].flag2*c)%lth;
return ;
}
int mid=(d[root].l+d[root].r)>>1;
push(root);
if(l<=mid) update1(root<<1,l,r,c);
if(r>mid) update1(root<<1|1,l,r,c);
up(root);
//cout<<d[root].sum<<" "<<d[root].l<<" "<<d[root].r<<endl;
}
void update2(int root,int l,int r,ll c){
if(l<=d[root].l&&d[root].r<=r){
d[root].sum=(d[root].sum+(c*(d[root].r-d[root].l+1))%lth)%lth;d[root].flag2=(d[root].flag2+c)%lth;
// cout<<d[root].sum<<" "<<d[root].flag1<<" "<<d[root].flag2<<" "<<d[root].l<<" "<<d[root].r<<endl;
return ;
}
int mid=(d[root].l+d[root].r)>>1;
push(root);
if(l<=mid) update2(root<<1,l,r,c);
if(r>mid) update2(root<<1|1,l,r,c);
up(root);
}
ll ans;
void querty(int root,int l,int r){
if(l<=d[root].l&&d[root].r<=r){
ans=(ans+d[root].sum)%lth;
return ;
}
push(root);
int mid=(d[root].l+d[root].r)>>1;
if(l<=mid) querty(root<<1,l,r);
if(r>mid) querty(root<<1|1,l,r);
up(root);
}
int main(){
ios::sync_with_stdio(false);
scanf("%d%lld",&n,<h);
for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
built(1,1,n);
int op,l,r;ll c;
int q;scanf("%d",&q);
while(q--){
scanf("%d%d%d",&op,&l,&r);
if(op==1){
scanf("%lld",&c);c=c%lth;
update1(1,l,r,c);
}
else if(op==2){
scanf("%lld",&c);c=c%lth;
update2(1,l,r,c);
}
else{
ans=0;querty(1,l,r);
printf("%lld\n",ans);
}
}
return 0;
}
1798: [Ahoi2009]Seq 维护序列seq
Time Limit: 30 Sec Memory Limit: 64 MBSubmit: 8292 Solved: 3033
[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
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
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

浙公网安备 33010602011771号