Loj 6280. 数列分块入门 4
题目描述
给出一个长为 nnn 的数列,以及 nnn 个操作,操作涉及区间加法,区间求和。
输入格式
第一行输入一个数字 nnn。
第二行输入 nnn 个数字,第 iii 个数字为 aia_iai,以空格隔开。
接下来输入 nnn 行询问,每行输入四个数字 opt\mathrm{opt}opt、lll、rrr、ccc,以空格隔开。
若 opt=0\mathrm{opt} = 0opt=0,表示将位于 [l,r][l, r][l,r] 的之间的数字都加 ccc。
若 opt=1\mathrm{opt} = 1opt=1,表示询问位于 [l,r][l, r][l,r] 的所有数字的和 mod (c+1)\bmod (c+1)mod(c+1)。
输出格式
对于每次询问,输出一行一个数字表示答案。
样例
样例输入
4
1 2 2 3
0 1 3 1
1 1 4 4
0 1 2 2
1 1 2 4
样例输出
1
4
数据范围与提示
对于 100% 100\%100% 的数据,1≤n≤50000,−231≤others 1 \leq n \leq 50000, -2^{31} \leq \mathrm{others}1≤n≤50000,−231≤others、ans≤231−1 \mathrm{ans} \leq 2^{31}-1ans≤231−1。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<bits/stdc++.h> 2 using namespace std; 3 4 #define F(i,a,b) for(int i=a;i<=b;i++) 5 #define D(i,a,b) for(int i=a;i>=b;i--) 6 #define ms(i,a) memset(a,i,sizeof(a)) 7 #define LL long long 8 #define st(x) ((x-1)*B+1) 9 #define ed(x) min(x*B,n) 10 #define bl(x) ((x-1)/B+1) 11 12 int inline read(){ 13 int x=0,w=0; char c=getchar(); 14 while (c<'0' || c>'9') w+=c=='-',c=getchar(); 15 while (c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar(); 16 return w? -x: x; 17 } 18 19 int const maxn=50003; 20 21 int n,B; 22 LL a[maxn],add[300],sum[300]; 23 24 void update(int l,int r,int LL z){ 25 int x=bl(l); 26 int y=bl(r); 27 if(x==y){ 28 F(i,l,r) a[i]+=z; 29 sum[x]+=(r-l+1)*z; 30 }else { 31 F(i,l,ed(x)) a[i]+=z; 32 sum[x]+=(ed(x)-l+1)*z; 33 F(i,st(y),r) a[i]+=z; 34 sum[y]+=(r-st(y)+1) * z; 35 F(i,x+1,y-1) add[i]+=z,sum[i]+=(ed(i)-st(i)+1)*z; 36 } 37 } 38 39 LL query(int l,int r,LL z){ 40 int x=bl(l); 41 int y=bl(r); 42 LL ans=0; 43 if(x==y){ 44 F(i,l,r) ans+=a[i]+add[x]; 45 }else { 46 F(i,l,ed(x)) ans+=a[i]+add[x]; 47 F(i,st(y),r) ans+=a[i]+add[y]; 48 F(i,x+1,y-1) ans+=sum[i]; 49 } 50 return ans % z; 51 } 52 int main(){ 53 n=read(); B=(int)sqrt(n); 54 F(i,1,n) a[i]=read(); 55 F(i,1,bl(n)) F(j,st(i),ed(i)) sum[i]+=a[j]; 56 57 F(i,1,n){ 58 int x=read(); 59 int l=read(); 60 int r=read(); 61 int c=read(); 62 if(x==0) update(l,r,c); 63 else printf("%lld\n",query(l,r,c+1)); 64 } 65 return 0; 66 }