谨以此文,纪念我逝去的这6个小时。原本昨天开开心心学学分块,但是入门2就卡住了,卡到生活都不能自理了。

题目分析

1

题目意思还是很清晰的。要求区间加法,区间查询符合条件的值。至于加法完全可以照搬前一道题目的方式,但是查询的操作着实有点耐人寻味了,因为区间查询意味着我们得采用分块的思想,不能暴力求解。首先我想到让每一个块变得有序,然后lower_bound查到第一个大于等于某个值的第一个元素的位置。首先对所有块排序的复杂度为√n * √n log_2(√n)也就是nlog_2(n)的复杂度,然后每一次查询最多是√n的复杂度,每一次区间加法是√n的复杂度,区间加法还要对固定两个残缺块进行重新排序(因为两边的残缺块可能会破坏有序结构,而中间的则会保持原来的顺序),这又需要2√nlog_2(√n)的复杂度。查询和加法一起是n√n的复杂度,所以整个算法就是n√n的复杂度,理论可行但是实践可惨死了。

首先,序列排序直接破坏了它的序列结构,这导致我在前一天20点-24点的提交一直不通过。后来我意识到了不能破坏它原有的序列,于是想到每个块用一个vector去保存。更新的时候clear再一个个push然后一个sort即可。然后又不行,发现是边界问题,因为我过于追求极限,因此会判断两边的块是否为完整块,完整则一起处理,不完整则先处理,这样导致如果没有处理,那么下面的循环条件没注意就不会被处理了。所以最好的方法就是,不管它,l那个点所在的块我就先处理,管你玩不完整的,右边的r也一样。

处理完了这个之后,我又被90分卡住了,最后发现是最后一个不完整的块它没有预先排序,但是其实是要的,虽然我也不知道为什么,因为我当时是感觉,它不可能会被当成一个完整的块,就基本没去管它了。但是事实就是这么魔幻:不管我就不让你过,就搞你心态。

然后来看看孩子的提交记录吧。

2

真的是太不容易了,下面给出我的AC代码。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
#include<bits/stdc++.h>
#define maxn 50005
using namespace std;
int a[maxn],bl[maxn],tag[maxn],block,n;
vector<int>v[maxn];


void update(int x){
v[x].clear();
for(int i=x*block;i<min((x+1)*block,n);i++){
v[x].push_back(a[i]);
}
sort(v[x].begin(),v[x].end());
}

void add(int l,int r,int c){
for(int i=l;i<=min(bl[l]*block+block-1,r);i++){
a[i]+=c;
}
update(bl[l]);
if(bl[l]==bl[r])return;
for(int i=bl[r]*block;i<=r;i++){
a[i]+=c;
}
update(bl[r]);

for(int i=bl[l]+1;i<bl[r];i++){
tag[i]+=c;
}

}
int query(int l,int r,int num){
int ans=0;
for(int i=l;i<=min(bl[l]*block+block-1,r);i++){
if(a[i]+tag[bl[i]]<num){
ans++;
}
}
if(bl[l]==bl[r])return ans;

for(int i=bl[r]*block;i<=r;i++){
if(a[i]+tag[bl[i]]<num){
ans++;
}
}

for(int i=bl[l]+1;i<bl[r];i++){
int k=lower_bound(v[i].begin(),v[i].end(),num-tag[i])-v[i].begin();
if(k>0)ans+=k;
}
return ans;
}

int main(){

cin>>n;
int q=n;
block=sqrt(n);
for(int i=0;i<n;i++){
cin>>a[i];
bl[i]=i/block;
}
for(int i=0;i<bl[n-1];i++){
update(i);
}
while(q--){
int opt,l,r,c;
cin>>opt>>l>>r>>c;
l--,r--;
if(opt==0)add(l,r,c);
else cout<<query(l,r,c*c)<<endl;
}


}