Codeforces Round #276 (Div. 1) E. Sign on Fence【主席树区间合并 + 二分】

题目链接

题意

有 $N$ 个栅栏，第 $i$ 个栅栏高为 $h_i$，宽为 $1$ ，$q$ 次询问区间 $[L,R]$ 中选 $w$ 个连续栅栏中最小的 $h$ 的最大值

题解

按 h 从大到小排序插入主席树，主席树 $X$ 维护插入第1~第X个栅栏，区间最大连续栅栏长度。

二分找$[L,R]$连续栅栏长度大于等于w的位置

#include<iostream>
#include<sstream>
#include<string>
#include<queue>
#include<map>
#include<unordered_map>
#include<set>
#include<vector>
#include<stack>
#include <utility>
#include<list>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<iomanip>
#include<time.h>
#include<random>
using namespace std;
#include<ext/pb_ds/priority_queue.hpp>
#include<ext/pb_ds/tree_policy.hpp>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
using namespace __gnu_pbds;
#include<ext/rope>
using namespace __gnu_cxx;

#define int long long
#define PI acos(-1.0)
#define eps 1e-9
#define lowbit(a) ((a)&-(a))
#define mid ((l+r)>>1)

const int mod = 1e9+7;
int qpow(int a,int b){
	int ans=1;
	while(b){
		if(b&1)ans=(ans*a)%mod;
		a=(a*a)%mod;
		b>>=1;
	}
	return ans;
}
const int INF = 0x3f3f3f3f;
const int N = 1e6+10;
struct node{
	int l,r;
	int llen,rlen,mlen;
}hjt[N*4];
struct Node{
	int w,id;
	bool friend operator<(const Node a,const Node b){
		return a.w>b.w;
	}
}a[N];
int root[N],cnt,n,m,res,ex;
void update(int &now,int pre,int l,int r,int pos){
	hjt[now=++cnt]=hjt[pre];
	if(l==r){
		hjt[now].llen=hjt[now].rlen=hjt[now].mlen=1;
		return;
	}
	if(pos<=mid)update(hjt[now].l,hjt[pre].l,l,mid,pos);
	else update(hjt[now].r,hjt[pre].r,mid+1,r,pos);
	int ls=hjt[now].l,rs=hjt[now].r;
	hjt[now].llen=hjt[ls].llen+(hjt[ls].llen==mid-l+1?hjt[rs].llen:0);
	hjt[now].rlen=hjt[rs].rlen+(hjt[rs].rlen==r-mid?hjt[ls].rlen:0);
	hjt[now].mlen=max(hjt[now].llen,hjt[now].rlen);
	hjt[now].mlen=max(hjt[now].mlen,hjt[ls].rlen+hjt[rs].llen);
	hjt[now].mlen=max(hjt[now].mlen,max(hjt[ls].mlen,hjt[rs].mlen));
}
void query(int now,int l,int r,int ql,int qr){
	if(ql<=l&&r<=qr){
		res=max(res,max(hjt[now].mlen,hjt[now].llen+ex));
		if(hjt[now].rlen==r-l+1)ex+=hjt[now].rlen;
		else ex=hjt[now].rlen;
		/*
		ex是右子树的右连续1长度，因为query是先递归右子树再左子树
		而右子树的右连续1长度可以连到左子树的左连续1上
		所以记录一下当前结点的右连续1(先递归到右子树，所以当前结点一定是相对右子树)
		*/
		return;
	}
	if(ql<=mid)query(hjt[now].l,l,mid,ql,qr);
	if(qr>mid)query(hjt[now].r,mid+1,r,ql,qr);
}
#define endl '\n'
signed main(){
	std::ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)cin>>a[i].w,a[i].id=i;
	sort(a+1,a+1+n);
	for(int i=1;i<=n;i++)update(root[i],root[i-1],1,n,a[i].id);
	cin>>m;
	while(m--){
		int x,y,w; cin>>x>>y>>w;
		int l=1,r=n,ans=0;
		while(l<=r){
			res=ex=0;
			query(root[mid],1,n,x,y);
			if(res>=w)ans=mid,r=mid-1;
			else l=mid+1;
		}
		cout<<a[ans].w<<endl;
	} 
}