Bzoj2653: middle

题面

传送门

Sol

丽洁姐的题目还是棒棒的
考虑二分答案
(Check?)
把小于它的设为(-1),大于等于它的设为(1)
([a, b])求一个最大后缀子段和
([c, d])求一个最大前缀子段和
([b+1, c-1])求一个和
加起来如果大于等于(0),那么满足要求,这个数还可以变大,否则就只能缩小
每个数开一个线段树来做,空间开不下,用主席树即可

那么问题来了
会不会二分的答案不在这个区间内呢?
显然是不会的,如果区间外有个数满足要求,那么区间内一定会有个数大于等于它,显然区间内的那个数最优,而且也是满足要求的

那么做完了

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(1e5 + 5);

IL int Input(){
    RG int x = 0, z = 1; RG char c = getchar();
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

int rt[_], a[_], n, tot, q[4], id[_];
struct Mx{
	int lmax, rmax, sum, mx;

	IL void Init(){
		lmax = rmax = -_, sum = 0;
	}
} Ans;
struct HJT{
	int ls, rs;
	Mx v;
} T[_ * 20];

IL Mx Merge(RG Mx A, RG Mx B){
	RG Mx C;
	C.lmax = max(A.lmax, A.sum + B.lmax);
	C.rmax = max(B.rmax, A.rmax + B.sum);
	C.sum = A.sum + B.sum;
	return C;
}

IL void Build(RG int &x, RG int l, RG int r){
	x = ++tot;
	T[x].v.lmax = T[x].v.rmax = T[x].v.sum = r - l + 1;
	if(l == r) return;
	RG int mid = (l + r) >> 1;
	Build(T[x].ls, l, mid), Build(T[x].rs, mid + 1, r);
}


IL void Modify(RG int &x, RG int l, RG int r, RG int p){
	T[++tot] = T[x], x = tot;
	if(l == r){
		T[x].v.lmax = T[x].v.rmax = T[x].v.sum = -1;
		return;
	}
	RG int mid = (l + r) >> 1;
	if(p <= mid) Modify(T[x].ls, l, mid, p);
	else Modify(T[x].rs, mid + 1, r, p);
	T[x].v = Merge(T[T[x].ls].v, T[T[x].rs].v);
}

IL void Query(RG int x, RG int l, RG int r, RG int L, RG int R){
	if(L <= l && R >= r){
		Ans = Merge(Ans, T[x].v);
		return;
	}
	RG int mid = (l + r) >> 1;
	if(L <= mid) Query(T[x].ls, l, mid, L, R);
	if(R > mid) Query(T[x].rs, mid + 1, r, L, R);
}

IL int Check(RG int mid){
	RG int val = 0;
	if(q[1] + 1 <= q[2] - 1) Ans.Init(), Query(rt[mid], 1, n, q[1] + 1, q[2] - 1), val += Ans.sum;
	Ans.Init(), Query(rt[mid], 1, n, q[0], q[1]), val += Ans.rmax;
	Ans.Init(), Query(rt[mid], 1, n, q[2], q[3]), val += Ans.lmax;
	return val >= 0;
}

IL int Cmp(RG int x, RG int y){
	return a[x] < a[y];
}

int main(RG int argc, RG char* argv[]){
	n = Input(), Build(rt[1], 1, n), T[0].v.Init();
	for(RG int i = 1; i <= n; ++i) a[i] = Input(), id[i] = i;
	sort(id + 1, id + n + 1, Cmp);
	for(RG int i = 2; i <= n; ++i) rt[i] = rt[i - 1], Modify(rt[i], 1, n, id[i - 1]);
	for(RG int Q = Input(), ans = 0; Q; --Q){
		for(RG int i = 0; i < 4; ++i) q[i] = (Input() + ans) % n + 1;
		sort(q, q + 4);
		RG int l = 1, r = n;
		while(l <= r){
			RG int mid = (l + r) >> 1;
			if(Check(mid)) ans = a[id[mid]], l = mid + 1;
			else r = mid - 1;
		}
		printf("%d
", ans);
	}
    return 0;
}

原文地址:https://www.cnblogs.com/cjoieryl/p/8551784.html