链接: 题目链接.
先分块,把n个数分成根号n个块,每个块根号n个数,(如果不是完全平方数区间个数就+1)
num数组记录每个数
lef数组维护当前这个数所在区间的左边界
rig数组维护当前这个数所在区间的右边界
blo数组维护当前这个数属于哪一个块
lazy数组维护每个块增加的数
思路:
修改
1.对于每次修改,找到L所在的块和R所在的块
2.如果L和R在相邻块或者同一块,直接枚举区间num[i]+C;
3.否则,枚举第L和R中间的块,lazy[i]+C,并且枚举L所在的块和R所在的块的num[i]+C;
查询
第i个数的值就是num[i]+lazy[blo[i]]
#include <bits/stdc++.h>
#define ll long long
#define T int T;scanf("%d", &T);while(T--)
using namespace std;
const int mod = 998244353;
const int maxn = 1e6+10;
int n;
int lef[maxn];
int rig[maxn];
int blo[maxn];
int num[maxn];
int lazy[maxn];
void update(int l,int r,int c) {
int L = blo[l];//l所在块的左边界
int R = blo[r];//r所在块的右边界
if(L+1>=R) {//相邻
for(int i = l; i <= r; i++) {
num[i] += c;
}
} else {
for(int i = L+1; i < R; i++) {//L+1到R-1这些块中所有数+C
lazy[i] += c;
}
update(l,rig[L],c);//l到l所在块的右边界所有数+C
update(lef[R],r,c);//r到r所在块的左边界所有数+C
}
}
void solve() {
int block = sqrt(n);
for(int i = 1; i <= n; i++) {
scanf("%d", num+i);
}
for(int i = 1; i <= block; i++) {
lef[i] = (i-1)*block+1;//维护第i个块的左边界
rig[i] = i*block;//维护第i个块的右边界
}
if(rig[block]!=n) {//如果不能刚好分成根号n个块,就加一个块
block++;
lef[block] = rig[block-1]+1;
rig[block] = n;
}
for(int i = 1; i <= block; i++) {
for(int j = lef[i]; j <= rig[i]; j++) {
lazy[i] = 0;
blo[j] = i;//维护每个数属于哪个块
}
}
for(int i = 0; i < n; i++){
int op,l,r,c;
scanf("%d %d %d %d", &op, &l, &r, &c);
if(!op) {
update(l,r,c);
} else {
printf("%d\n", num[r]+lazy[blo[r]]);
}
}
}
int main() {
scanf("%d", &n);
solve();
return 0;
}
题目链接.
思路:跟上一题相比,增加一个vector用于记录每一块的数,并排序,方便二分查找
修改:
1.修改一整块的时候,只要lazy+C
2.修改一块的某一段,枚举修改每个点,并且重新维护每一块的vector
查询:
1.查询一整块的时候用二分查出这一块小于c*c的个数
2.查询某一块的某一段枚举查询
#include <bits/stdc++.h>
#define ll long long
#define T int T;scanf("%d", &T);while(T--)
using namespace std;
const int mod = 998244353;
const int maxn = 1e6+10;
int n;
int lef[maxn];
int rig[maxn];
int blo[maxn];
ll num[maxn];
ll lazy[maxn];
vector<ll> a[maxn];
void update(int l,int r,ll c) {
int L = blo[l];
int R = blo[r];
if(L+1>=R) {
for(int i = l; i <= r; i++) {
num[i] += c;
}
a[L].clear();//更改某一块的某一段,重新维护当前块的vector
for(int i = lef[L]; i <= rig[L]; i++) {
a[L].push_back(num[i]);
}
sort(a[L].begin(),a[L].end());
a[R].clear();
for(int i = lef[R]; i <= rig[R]; i++) {
a[R].push_back(num[i]);
}
sort(a[R].begin(),a[R].end());
// cout << "as" << endl;
} else {
for(int i = L+1; i < R; i++) { // 更新某一块,直接lazy+C
lazy[i] += c;
}
update(l,rig[L],c);
update(lef[R],r,c);
}
}
ll query(int l,int r,ll c) {
int L = blo[l];
int R = blo[r];
ll ans = 0;
if(L+1>=R) {
for(int i = l; i <= r; i++) { // 当前这个数是num[i]+lazy[blo[i]],和上一题查找的数相同
if(num[i]+lazy[blo[i]] < c) ans++;
}
} else {
for(int i = L+1; i < R; i++) {//x+lazy<c 所以x<c-lazy
ans += lower_bound(a[i].begin(),a[i].end(),c-lazy[i])-a[i].begin();
}
ans += query(l,rig[L],c) + query(lef[R],r,c);
}
return ans;
}
void solve() {
int block = sqrt(n);
for(int i = 1; i <= n; i++) {
scanf("%lld", num+i);
}
for(int i = 1; i <= block; i++) {
lef[i] = (i-1)*block+1;
rig[i] = i*block;
}
if(rig[block]!=n) {
block++;
lef[block] = rig[block-1]+1;
rig[block] = n;
}
for(int i = 1; i <= block; i++) {
for(int j = lef[i]; j <= rig[i]; j++) {
lazy[i] = 0;
blo[j] = i;
a[i].push_back(num[j]);
}
sort(a[i].begin(),a[i].end());
}
for(int i = 0; i < n; i++){
int op,l,r;
ll c;
scanf("%d %d %d %lld", &op, &l, &r, &c);
if(!op) {
update(l,r,c);
} else {
printf("%lld\n", query(l,r,c*c));
}
}
}
int main() {
scanf("%d", &n);
solve();
return 0;
}
题目链接.
思路:二分查找比c小的最大值
修改:和上一题一样
1.修改一整块的时候,只要lazy+C
2.修改一块的某一段,枚举修改每个点,并且重新维护每一块的vector
查询:
1.查询一整块的时候用二分查出这一块小于c的位置pos,pos-1就是小于c的最大值
2.查询某一块的某一段枚举维护小于c的最大值
#include <bits/stdc++.h>
#define ll long long
#define T int T;scanf("%d", &T);while(T--)
using namespace std;
const int mod = 998244353;
const int maxn = 1e6+10;
int n;
int lef[maxn];
int rig[maxn];
int blo[maxn];
ll maxx[maxn];
ll num[maxn];
ll lazy[maxn];
vector<ll> a[maxn];
void update(int l,int r,ll c) {
int L = blo[l];
int R = blo[r];
if(L+1>=R) {
for(int i = l; i <= r; i++) {
num[i] += c;
}
a[L].clear();
for(int i = lef[L]; i <= rig[L]; i++) {
a[L].push_back(num[i]);
}
sort(a[L].begin(),a[L].end());
a[R].clear();
for(int i = lef[R]; i <= rig[R]; i++) {
a[R].push_back(num[i]);
}
sort(a[R].begin(),a[R].end());
} else {
for(int i = L+1; i < R; i++) {
lazy[i] += c;
}
update(l,rig[L],c);
update(lef[R],r,c);
}
}
ll ans;
int f;
void query(int l,int r,ll c) {
int L = blo[l];
int R = blo[r];
if(L+1>=R) {
for(int i = l; i <= r; i++) {
if(num[i]+lazy[blo[i]] < c) { //更新小于c的最大值
ans = max(ans,num[i]+lazy[blo[i]]);
f = 1;
}
}
} else {
for(int i = L+1; i < R; i++) {//二分小于c的最大值
int pos = lower_bound(a[i].begin(),a[i].end(),c-lazy[i])-a[i].begin();
if(pos!=0 && a[i][pos-1]+lazy[i] < c) {
ans = max(ans,a[i][pos-1]+lazy[i]);
}
}
query(l,rig[L],c);
query(lef[R],r,c);
}
}
void solve() {
int block = sqrt(n);
for(int i = 1; i <= n; i++) {
scanf("%lld", num+i);
}
for(int i = 1; i <= block; i++) {
lef[i] = (i-1)*block+1;
rig[i] = i*block;
}
if(rig[block]!=n) {
block++;
lef[block] = rig[block-1]+1;
rig[block] = n;
}
for(int i = 1; i <= block; i++) {
for(int j = lef[i]; j <= rig[i]; j++) {
lazy[i] = 0;
blo[j] = i;
a[i].push_back(num[j]);
}
sort(a[i].begin(),a[i].end());
}
for(int i = 0; i < n; i++){
int op,l,r;
ll c;
scanf("%d %d %d %lld", &op, &l, &r, &c);
ans = INT_MIN;
f = 0;
if(!op) {
update(l,r,c);
} else {
query(l,r,c);
if(ans==INT_MIN) {
printf("-1\n");
} else {
printf("%lld\n", ans);
}
}
}
}
int main() {
scanf("%d", &n);
solve();
return 0;
}