数列分块整理..

链接: 题目链接.

先分块,把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;
}