A Simple Problem with Integers
| Time Limit: 5000MS | Memory Limit: 131072K | |
| Total Submissions: 73163 | Accepted: 22585 | |
| Case Time Limit: 2000MS | ||
Description
You have N integers, A1, A2, ... , AN. You need to deal with two kinds of operations. One type of operation is to add some given number to each number in a given interval. The other is to ask for the sum of numbers in a given interval.
Input
The first line contains two numbers N and Q. 1 ≤ N,Q ≤ 100000.
The second line contains N numbers, the initial values of A1, A2, ... , AN. -1000000000 ≤ Ai ≤ 1000000000.
Each of the next Q lines represents an operation.
"C a b c" means adding c to each of Aa, Aa+1, ... , Ab. -10000 ≤ c ≤ 10000.
"Q a b" means querying the sum of Aa, Aa+1, ... , Ab.
Output
You need to answer all Q commands in order. One answer in a line.
Sample Input
10 5 1 2 3 4 5 6 7 8 9 10 Q 4 4 Q 1 10 Q 2 4 C 3 6 3 Q 2 4
Sample Output
4 55 9 15
Hint
The sums may exceed the range of 32-bit integers.
Source
POJ Monthly--2007.11.25, Yang Yi
最裸的线段树lazy标记
1 #include <iostream>
2 #include <cstdio>
3 using namespace std;
4 const int N = 100005;
5 typedef long long LL;
6 LL sum[N<<2]; //sum用来存储每个节点的子节点数值的总和
7 LL add[N<<2];//add用来记录该节点的每个数值应该加多少
8 struct Node{
9 int l,r;//表示改点的左右区间
10 int mid(){//结构体函数
11 return (l+r)>>1;
12 }
13 } tree[N<<2];
14
15 void PushUp(int rt){//算某一节点的左右孩子值的和
16 sum[rt] = sum[rt<<1] + sum[rt<<1|1];
17 }
18
19 void PushDown(int rt,int m){//rt当前节点 m此节点的区间长度
20 if(add[rt]!=0){//如果当前节点lazy标记不为 0
21 add[rt<<1] += add[rt];//左右孩子lazy累加
22 add[rt<<1|1] += add[rt];
23 sum[rt<<1] += add[rt]*(m-(m>>1));//更新计算左右孩子
24 sum[rt<<1|1] += add[rt] * (m>>1);
25 add[rt] = 0;//小细节，但很重要，不能忘记lazy用过后清零
26 }
27 }
28
29 void build(int l,int r,int rt){
30 tree[rt].l = l;
31 tree[rt].r = r;
32 add[rt] = 0;//lazy标记记为0
33 if(l == r){
34 scanf("%lld",&sum[rt]);//把子节点信息记录在sum里
35 return ;
36 }
37 int m = tree[rt].mid();
38 build(l,m,rt<<1);//建立左右孩子，这时编号是按照类似“宽搜”来编的
39 build(m+1,r,rt<<1|1);
40
41 PushUp(rt);//算出此节点的权值
42 }
43
44 void update(int c,int l,int r,int rt){//当前节点rt，在l和r区间上加上c
45 if(l<=tree[rt].l&&tree[rt].r<=r){//当前节点所表示的区间完全被所要更新的区间包含
46 add[rt]+=c;//lazy累加,add[i]数组是针对i左右孩子的
47 sum[rt]+=(LL)c*(tree[rt].r-tree[rt].l+1);//这句话很重要和下面的PushUP()类似
48 return;
49 }
50 PushDown(rt,tree[rt].r - tree[rt].l + 1);//用rt的lazy更新其子节点
51 int m = tree[rt].mid();
52 if(l<=m) update(c,l,r,rt<<1);
53 if(m+1<=r) update(c,l,r,rt<<1|1);
54 PushUp(rt);
55 }
56
57 LL query(int l,int r,int rt){//当前节点为rt，求l，到r的和
58 if(l<=tree[rt].l&&tree[rt].r<=r){//区间完全包含，可以直接返回
59 return sum[rt];
60 }
61 //因为此时用到rt节点，所以才更新lazy
62 PushDown(rt,tree[rt].r-tree[rt].l + 1);
63
64 int m = tree[rt].mid();
65 LL res = 0;
66
67 if(l<= m)//要求的区间有一部分在mid的左边
68 res += query(l,r,rt<<1);
69 if(m+1<=r)
70 res += query(l,r,rt<<1|1);
71 return res;
72 }
73
74 int main(){
75
76 int n,m;
77
78 scanf("%d %d",&n,&m);
79 build(1,n,1);
80
81 while(m--){
82 char ch[2];
83 scanf("%s",ch);
84 int a,b,c;
85 if(ch[0] == ‘Q‘){
86 scanf("%d %d", &a,&b);
87 printf("%lld\n",query(a,b,1));
88 }
89 else{
90 scanf("%d %d %d",&a,&b,&c);
91 cin>>a>>b>>c;
92 update(c,a,b,1);
93 }
94 }
95
96 return 0;
97 }
还有一种较好理解的方法：
1 #include<iostream>
2 #include<cstdio>
3 #include<cstdlib>
4 #include<cmath>
5 #include<cstring>
6 #include<algorithm>
7 using namespace std;
8 typedef long long LL;
9 const LL NN=100005;
10 LL a[NN];
11 LL N,tot,T,M;
12 struct node{
13 LL num;
14 LL lc,rc;
15 LL l,r;
16 LL sum;
17 LL lazy;
18 };
19
20 node segtree[2000000];
21
22 void Build(LL root,LL l,LL r){
23 segtree[root].l=l;
24 segtree[root].r=r;
25 segtree[root].lazy=0;
26 segtree[root].num=tot;
27 if(l==r){
28 segtree[root].sum=a[l];
29 return ;
30 }
31
32 LL mid=(l+r)>>1;
33
34 segtree[root].lc=++tot;
35 segtree[segtree[root].lc].num=tot;
36 Build(tot,l,mid);
37
38 segtree[root].rc=++tot;
39 segtree[segtree[root].rc].num=tot;
40 Build(tot,mid+1,r);
41 segtree[root].sum=segtree[segtree[root].lc].sum+
42 segtree[segtree[root].rc].sum;
43
44 }
45
46 void updatason(LL root){
47 LL d=segtree[root].lazy;
48 if(d!=0){
49 segtree[segtree[root].lc].lazy+=(LL)d;
50 segtree[segtree[root].rc].lazy+=(LL)d;
51
52 segtree[segtree[root].lc].sum+=(LL)d*(segtree[segtree[root].lc].r
53 -segtree[segtree[root].lc].l+1);
54 segtree[segtree[root].rc].sum+=(LL)d*(segtree[segtree[root].rc].r
55 -segtree[segtree[root].rc].l+1);
56 segtree[root].lazy=0;
57 }
58 }
59
60 LL query(LL root,LL l,LL r){
61 if(l<=segtree[root].l&&segtree[root].r<=r){
62 return segtree[root].sum;
63 }
64 updatason(root);
65 LL ans=0;
66 LL mid=(segtree[root].r+segtree[root].l)>>1;
67 if(l<=mid) ans+=(LL)query(segtree[root].lc,l,r);
68 if(mid+1<=r) ans+=(LL)query(segtree[root].rc,l,r);
69 return ans;
70 }
71
72 void changesect(LL root,LL l,LL r,LL d){
73 if(l<=segtree[root].l&&segtree[root].r<=r){
74 segtree[root].sum+=(LL)(segtree[root].r-segtree[root].l+1)*d;
75 segtree[root].lazy+=(LL)d;
76 return ;
77 }
78 updatason(root);
79 LL mid=(segtree[root].r+segtree[root].l)>>1;
80 if(l<=mid) changesect(segtree[root].lc,l,r,d);
81 if(mid+1<=r) changesect(segtree[root].rc,l,r,d);
82 segtree[root].sum=(LL)segtree[segtree[root].lc].sum+
83 (LL)segtree[segtree[root].rc].sum;
84 }
85 int main(){
86
87 scanf("%d%d",&N,&M);
88 for(LL i=1;i<=N;i++)
89 scanf("%lld",&a[i]);
90 tot++;
91 Build(1,1,N);
92 LL ll,rr,de;
93 char x[2];
94 for(LL i=1;i<=M;i++){
95 scanf("%s",x);
96 if(x[0]==‘C‘){
97 cin>>ll>>rr>>de;
98 changesect(1,ll,rr,de);
99 continue;
100 }
101 if(x[0]==‘Q‘){
102 cin>>ll>>rr;
103 printf("%lld\n",query(1,ll,rr));
104 }
105 }
106 return 0;
107 }
由于此题数据较大，所以要将int 都改为long long 。但这也给了我一个惨痛的教训，如果要在定义变量的时候该类型，千万不要忘记改变函数返回值的类型和输入改成“lld”，这个题耗了我一天加一上午，但也让我对线段树更加深入了解，试遍了各种错误，模板更加完善，，，也不错。。。。
时间： 2024-10-11 07:20:50