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

最新推荐文章于 2022-11-22 23:15:14 发布

学不会dp QAQ

最新推荐文章于 2022-11-22 23:15:14 发布

阅读量93

点赞数

分类专栏：数据结构思维文章标签： acm竞赛

本文链接：https://blog.csdn.net/qq_50554695/article/details/120916214

版权

数据结构同时被 2 个专栏收录

25 篇文章 1 订阅

订阅专栏

思维

9 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何运用主席树解决区间查询中选出连续子数组满足特定条件的最大值问题。通过排序、构建主席树以及二分查找，实现了高效求解动态区间最值的算法。文章详细阐述了代码实现过程，并给出了样例输入和输出。

摘要由CSDN通过智能技术生成

题目链接

题意

有 $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;
	} 
}