基础算法——ST表(倍增思想)

文章介绍了如何在算法竞赛中使用ST表(SegmentTree)进行区间查询,包括RMQ(RangeMinimumQuery)的应用以及线段树的对比。通过实例展示了如何在P3865、P1440、P1816和P2880等问题中使用这些数据结构解决区间属性问题。
摘要由CSDN通过智能技术生成

原理

在这里插入图片描述
在这里插入图片描述
提示:只涉及查询时,可以用RMQ而不用线段树,线段树代码太长了,而涉及修改时,就得使用线段树啦。

应用

P3865(ST模板题)

P3865
在这里插入图片描述

快读模板

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
	return x*f;
}
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int f[N][22];
int n,m,l,r;

int main(){
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//还是要用scanf和printf最快 
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);//要从1开始,记数下标1 
	for(int j=1;j<=19;j++){
		for(int i=1;i+(1<<j)-1<=n;i++){
			f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
		}
	}
	while(m--){
		scanf("%d%d",&l,&r);
		int k=log2(r-l+1),ans=-1;
		printf("%d\n",max(f[l][k],f[r-(1<<k)+1][k]));
	}
}

P1440

P1440
在这里插入图片描述
法一:RMQ(ST表)
会MLE两个点
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
const int N=2e6+10;
int f[N][22];
int n,m,l,r;

int main(){
//	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//还是要用scanf和printf最快 
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);//要从1开始,记数下标1 
	for(int j=1;j<=21;j++){
		for(int i=1;i+(1<<j)-1<=n;i++){
			f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
		}
	}
	printf("0\n");
	if(n==1) return 0;
	for(int i=2;i<=n;i++){
		int l=max(1,i-m),r=i-1,len=r-l+1;
		int k=log2(len);
		printf("%d\n",min(f[l][k],f[r-(1<<k)+1][k]));
	}
	return 0;
}

总结归纳:每次查询一个区间属性时,先将区间两边端点和长度明确,再求k,最后求属性。l,r,len=r-l+1,k=log2(len)。

法二:单调队列

#include<bits/stdc++.h>
using namespace std;
const int N=2e6+10;
int q[N],tt=-1,hh,a[N];
int n,m;

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++){
		if(hh<=tt&&q[hh]<i-m) hh++;
		if(hh<=tt) printf("%d\n",a[q[hh]]);
		else printf("0\n");
		while(hh<=tt&&a[q[tt]]>=a[i]) tt--;
		q[++tt]=i;
	}
	return 0;
}

这个题目还可以用线段树维护。

P1816

P1816
查询区间最小值,当属范围在1e5范围内,RMQ表不会MLE。
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
const int N=2e6+10;
int f[N][22];
int n,m,l,r;

int main(){
//	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//还是要用scanf和printf最快 
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);//要从1开始,记数下标1 
	for(int j=1;j<=21;j++){
		for(int i=1;i+(1<<j)-1<=n;i++){
			f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
		}
	}
	while(m--){
		scanf("%d%d",&l,&r);
		int len=r-l+1,k=log2(len);
		printf("%d ",min(f[l][k],f[r-(1<<k)+1][k]));
	}
	return 0;
}

P2880

P2880
题意:求区间最大最小值的差值
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+10;
int f[N][22],g[N][22],a[N];
int n,m,l,r;

int main(){
//	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);//还是要用scanf和printf最快 
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&f[i][0]);//要从1开始,记数下标1 
	memcpy(g,f,sizeof g);
	for(int j=1;j<=21;j++){
		for(int i=1;i+(1<<j)-1<=n;i++){
			f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
			g[i][j]=max(g[i][j-1],g[i+(1<<(j-1))][j-1]);//f改变了 
		}
	}
	while(m--){
		scanf("%d%d",&l,&r);
		int len=r-l+1,k=log2(len);
		int mx=max(g[l][k],g[r-(1<<k)+1][k]);
		int mn=min(f[l][k],f[r-(1<<k)+1][k]);
//		cout<<mx<<" "<<mn<<" ";
		int ans=mx-mn;
		printf("%d\n",ans); 
	}
	return 0;
}

参考

董晓博客

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值