【线段树+单调栈】Prime的把妹计划

题目

思路

考虑只做 $a_l\le a_r$ （然后再反过来做一遍）
首先预处理一个 $pre$ 表示 $a_i$ 左边第一个比 $a_i$ 大的数的位置。

对于每个 $a_i$ 我们分类讨论：

$a_i$ 大于等于栈顶元素

此时以 $pre$ 右边且在栈中的元素为左端点的区间最优右端点为 $i$ ，首先要满足在 $pre$ 右边，因为 $pre$ 就会成为最大值，占据右端点的位置，不合题意；其次要满足还在栈中，若不在栈中，则一定是因为它右边出现了比它小的元素才把它弹出的，那么它就不会是最小值了，也不合题意。
用一个下标与该单调栈一一对应的线段树维护这个过程，区间修改一段的右端点为 $i$ ，然后把 $a_i$ 压入栈顶。

$a_i$ 大于等于栈顶元素

不断弹栈，每弹出一个元素就把它加入用另一个与原序列下标一一对应的线段树，然后把第一棵线段树对应位置清空（不然会WA），
最后还是把 $a_i$ 压入栈顶。
考虑答案怎么求？
把询问按右端点排序，当 $a_i$ 被加入栈中之后，在两棵线段树上分别查询[ L , R ] [L,R][L,R]对应的区间最大值，更新答案。
注意 $L$ 在第一棵线段树中不一定存在对应的下标，取它右边最小的那个即可。

代码

#include<bits/stdc++.h>
#define rson k<<1|1
#define lson k<<1
using namespace std;
const int N=5e4+77;
struct query{int x,y,ans;}Q[N];
int a[N],le[N],ri[N],sta[N],mx[N<<2],mn[N<<2],f[235][235],Ls[235],n,q,siz,tot;
void build(int k,int l,int r)
{
	if(l==r) mx[k]=ri[l],mn[k]=le[l];
	else
	{
		int mid=l+r>>1;
		build(lson,l,mid);
		build(rson,mid+1,r);
		mx[k]=max(mx[lson],mx[rson]);
		mn[k]=min(mn[lson],mn[rson]);
	}
}
int querymax(int k,int l,int r,int x,int y,int z)
{
	if(mx[k]<=z) return 0;
	if(l==r) return l;
	int mid=l+r>>1,tmp;
	if(y<=mid) return querymax(lson,l,mid,x,y,z);
	if(x>mid) return querymax(rson,mid+1,r,x,y,z);
	tmp=querymax(lson,l,mid,x,mid,z);
	return tmp?tmp:querymax(rson,mid+1,r,mid+1,y,z);
}
int querymin(int k,int l,int r,int x,int y,int z)
{
	if(mn[k]>z) return 0;
	if(l==r) return l;
	int mid=l+r>>1,tmp;
	if(y<=mid) return querymin(lson,l,mid,x,y,z);
	if(x>mid) return querymin(rson,mid+1,r,x,y,z);
	tmp=querymin(rson,mid+1,r,mid+1,y,z);
	return tmp?tmp:querymin(lson,l,mid,x,mid,z);
}
void calc()
{
	build(1,1,n);
	memset(f,0,sizeof(f));
	for(int i=1; i<=tot; i++)
		for(int j=Ls[i]; j<Ls[i+1]; j++)
			f[i][i]=max(f[i][i],j-querymax(1,1,n,max(le[j]+1,Ls[i]),j,j)+1);
	for(int len=1; len<tot; len++)
		for(int i=1,j=len+1; j<=tot; i++,++j)
		{
			f[i][j]=f[i][j-1];
			for(int k=Ls[j]; k<Ls[j+1]; ++k)
				f[i][j]=max(f[i][j],k-querymax(1,1,n,max(le[k]+1,Ls[i]),k,k)+1);
		}
	for(int i=1,ans,l,r; i<=q; i++)
	{
		ans=0;
		l=(Q[i].x-1)/siz+1;
		r=(Q[i].y-1)/siz+1;
		if(Q[i].y<Ls[r+1]-1)
			for(int j=Ls[r--]; j<=Q[i].y; j++) ans=max(ans,j-querymax(1,1,n,max(le[j]+1,Q[i].x),j,j)+1);
		if(Q[i].x>Ls[l])
		{
			++l;
			if(f[l][r]>ans) ans=f[l][r];
			for(int j=Q[i].x;j<Ls[l];++j) ans=max(ans,querymin(1,1,n,j,min(ri[j]-1,Q[i].y),j)-j+1);
		}
		else if(f[l][r]>ans) ans=f[l][r];
		if(ans>Q[i].ans) Q[i].ans=ans;
	}
}
int main()
{
	int x,y,top;
	scanf("%d",&n),siz=sqrt(n);
	tot=n/siz; if(tot*siz<n) ++tot;
	for(int i=1; i<=tot; i++) Ls[i]=(i-1)*siz+1;
	Ls[tot+1]=n+1;
	for(int i=1; i<=n; i++) scanf("%d",&a[i]);
	scanf("%d",&q);
	for(int i=1; i<=q; i++) scanf("%d%d",&Q[i].x,&Q[i].y);
	for(int i=1,top=0; i<=n; i++)
	{
		while(top&&a[i]<=a[sta[top]]) --top;
		le[i]=sta[top],sta[++top]=i;
	}
	for(int i=n,top=0; i; i--)
	{
		while(top&&a[i]>=a[sta[top]]) --top;
		ri[i]=top?sta[top]:n+1,sta[++top]=i;
	}
	calc();
	for(int i=1,top=0; i<=n; i++)
	{
		while(top&&a[i]>=a[sta[top]]) --top;
		le[i]=sta[top],sta[++top]=i;
	}
	for(int i=n,top=0; i; i--)
	{
		while(top&&a[i]<=a[sta[top]]) --top;
		ri[i]=top?sta[top]:n+1,sta[++top]=i;
	}
	calc();
	for(int i=1; i<=q; i++) printf("%d\n",Q[i].ans);
}