数据结构——ST表

ST表的定义

ST表,又名稀疏表,是一种基于倍增思想,用于解决可重复贡献问题的数据结构

倍增思想

这里列举一个去寻找一个区间内的最大值的例子

因为每次会将将区间增大一倍,所以才被称之为倍增思想 ,这种思想十分好用,建议友友们平常多多练习

ST表的适用范围

如果A区间和B区间可能有重叠的部分

但是并不影响A+B区间的答案,能通过 A区间答案 和 B区间答案 就加工出来

那么对应的区间询问,就是一个可重复贡献问题

例如:区间最大值,区间最小值、区间公约数等,但是区间求和就不符合这个要求

再例如:区间按位与、区间按位或,ST表都能高效地解决

ST表的优势和劣势

RMQ问题(Range Maximum/Minimum Query)可以用ST表维护,也可以用线段树等结构维护
ST表的优势:构建过程时间复杂度0(n*logn),单次查询时间复杂度0(1),代码量较小
ST表的劣势:需要空间较大,能维护的信息非常有限,不支持修改操作

只能维护静态区间内的可重复贡献类问题,不能维护动态的

模版

预处理模版:

cpp 复制代码
//st[i][j]表示的是以i为起点,长度为2^j的区间的最值
//外层循环遍历的是长度的指数,内层循环遍历的是起点 
for(int j=1;j<20;j++)
{
	for(int i=1;i+(1<<j)-1<=n;i++)
	{
		st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
	}
} 

查询模版:

cpp 复制代码
//maxn是统计给定区间的最值
//k是用来判断区间的长度
//l,r是给定的区间范围 
for(int i=1;i<=m;i++)
{
	cin>>l>>r;
	int k=log2(r-l+1);
	maxn=max(dp[l][k],dp[r-(1<<k)+1][k]);
} 

例题:

P3865 【模板】ST 表 && RMQ 问题

思路:十分板正的一个ST表板题,就是去求区间内的一个最大值,纯板题

cpp 复制代码
#include<bits/stdc++.h>  
using namespace std;  
#define int long long  

int n, m;  
int l, r;  
int st[100005][20];  

signed main() {  
    ios::sync_with_stdio(0);  
    cin.tie(0);  
    cout.tie(0);  
    cin >> n >> m;  
    
    for (int i = 1; i <= n; i++) {  
        cin >> st[i][0];  
    }  
    for (int j = 1; j < 20; j++) {  
        for (int i = 1; i + (1 << j) - 1 <= n; i++) {  
            st[i][j] = max(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]);  
        }  
    }  
    for (int i = 1; i <= m; i++) {  
        cin >> l >> r;  
        int k = log2(r - l + 1);  
        cout << max(st[l][k], st[r - (1 << k) + 1][k]) << "\n";  
    }  

    return 0;  
}

P2880 [USACO07JAN] Balanced Lineup G

思路:也是ST表的板题,只不过是去求一遍最大值的ST表,然后再去求一遍最小值的ST表,然后减去就可以了

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,q;
int f_max[50005][17];
int f_min[50005][17];
int l,r;
signed main()
{
	cin>>n>>q;
	for(int i=1;i<=n;i++)
	{
		cin>>f_max[i][0];
		f_min[i][0]=f_max[i][0];
	}
	for(int j=1;j<=16;j++)
	{
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			f_max[i][j]=max(f_max[i][j-1],f_max[i+(1<<(j-1))][j-1]);
			f_min[i][j]=min(f_min[i][j-1],f_min[i+(1<<(j-1))][j-1]);
		}
	}
	for(int i=1;i<=q;i++)
	{
		cin>>l>>r;
		int k=log2(r-l+1);
		cout<<max(f_max[l][k],f_max[r-(1<<k)+1][k])-min(f_min[l][k],f_min[r-(1<<k)+1][k])<<"\n";
	}
}

P2251 质量检测

思路1:滑动窗口问题,可以直接单调队列解决,去寻找m区间长度的窗口内的最小值,但是和今天的主题不相符,这种写法就靠看官自己去琢磨了

思路2:ST表,这题是求区间内的最小值,但是却是给你规定了区间大小,让你找出整个区间的这么大小的最小值,直接预处理一遍,然后遍历即可,但是如果m不是2^k,就找到最大的k, 然后用m-(2^k)+1的k区间长度,去找最小值即可

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m;
int st[1000005][25]; 
signed main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>st[i][0];
	}
	for(int j=1;j<=21;j++)
	{
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			st[i][j]=min(st[i][j-1],st[i+(1<<(j-1))][j-1]);
		}
	}
	int k=log2(m);
	for(int i=1;i+m-1<=n;i++)
	{
		cout<<min(st[i][k],st[i+m-(1<<k)][k])<<"\n";
	}
	return 0;
}

P1890 gcd区间

思路:区间gcd也是一种可重复贡献问题,因此也可以用ST表,我们直接预处理一遍,然后直接跑即可

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m;
int st[1005][15];
int l,r;
int gcd(int a,int b)
{
	if(b==0)
	return a;
	return gcd(b,a%b);
}
signed main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>st[i][0];
	}
	for(int j=1;j<=14;j++)
	{
		for(int i=1;i+(1<<j)-1<=n;i++)
		{
			st[i][j]=gcd(st[i][j-1],st[i+(1<<(j-1))][j-1]);
		}
	}
	for(int i=1;i<=m;i++)
	{
		cin>>l>>r;
		int k=log2(r-l+1);
		cout<<gcd(st[l][k],st[r-(1<<k)+1][k])<<"\n";
	}
	return 0;
}

P4155 [SCOI2015] 国旗计划

思路:这题其实一开始我也没看出来是ST表的题目,后来想了好久我才恍然大悟,我们首先来看题,每个士兵都有自己负责的一个区域,为L和R,然后呢,我们有m个边防站,然后我们要遍历第i个士兵必须参与的情况下,最少需要多少个哨兵

我们首先需要拆环为链,我们之间将这个数组拷贝一遍,再延伸出来一个长度,然后我们去处理士兵负责的范围,加入右区间的端点,小于左区间端点,那么我们就直接将右区间的端点+m即可,然后我们再去复制每个哨兵一次,其区间的值都+m

然后我们就跑ST表,表示从i开始往后2^j次方能够达到什么位置,如果是当前这个人的左区间+m的位置,那么就说明覆盖了整个哨所,那就是可以了,我们需要去处理一下

然后直接写即可

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,m;
int l,r;
struct node{
	int l;
	int r;
	int pos;
}a[400005];
int st[400005][25];
int sum[400004];
bool cmp(node a,node b)
{
	return a.l<b.l;
}
signed main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>l>>r;
		if(l>r)
		{
			r+=m;
		}
		a[i].l=l;
		a[i].r=r;
		a[i].pos=i;
	}
	sort(a+1,a+1+n,cmp);
	for(int i=1;i<=n;i++)
	{
		a[i+n].l=a[i].l+m;
		a[i+n].r=a[i].r+m;
	}
	for(int i=1,j=1;i<=2*n;i++)
	{
		while(j<=2*n&&a[j].l<=a[i].r)
		{
			j++;
		}
		st[i][0]=j-1;
	}
	for(int j=1;j<=20;j++)
	{
		for(int i=1;i+(1<<j)-1<=2*n;i++)
		{
			st[i][j]=st[st[i][j-1]][j-1];
		}
	}
	for(int i=1;i<=n;i++)
	{
		int flag=a[i].l+m;
		int ans=0;
		int p=i;
		for(int j=20;j>=0;j--)
		{
			if(st[p][j]&&a[st[p][j]].r<flag)
			{
				ans+=1<<j;
				p=st[p][j];
			}
		}
		sum[a[i].pos]=ans+2;
	}
	for(int i=1;i<=n;i++)
	{
		cout<<sum[i]<<" ";
	}
	return 0;
}

Frequent values

思路:我们应当用到题目说的有序,我们可以将每一种数视为一个桶,然后将每个桶内的个数去求出来,然后去记录每个桶的左端点和右端点,然后我们在询问的时候,然后我们去预处理桶的数目的ST表,表示从第i个桶开始长度为2^j次方的桶内的最大值,然后我们再去找两边的最大值,然后去处理三部分最大值的最大值即可

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,q;
int a;
int l[500005];
int r[500005];
int st[500005][24];
int num[500005];
int x,y;
void ini()
{
	memset(l,0,sizeof(l));
	memset(r,0,sizeof(r));
	memset(st,0,sizeof(st));
	memset(num,0,sizeof(num));
}
signed main()
{
	while(cin>>n)
	{
		if(n==0)
		{
			break;
		}
		ini();
		cin>>q;
		int len=0;
		int flag=-0x3f3f3f3f;
		for(int i=1;i<=n;i++)
		{
			cin>>a;
			if(a!=flag)
			{
				r[len]=i-1;
				len++;
				l[len]=i;
				flag=a;
			}
			st[len][0]++;
			num[i]=len;
		}
		r[len]=n;
		for(int j=1;j<=20;j++)
		{
			for(int i=1;i+(1<<j)-1<=n;i++)
			{
				st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
			}
		}
		for(int i=1;i<=q;i++)
		{
			cin>>x>>y;
			if(x>y)
			swap(x,y);
			if(num[x]==num[y])
			{
				cout<<y-x+1<<"\n";
			}
			else
			{
				int ans=0;
				int L=num[x]+1;
				int R=num[y]-1;
				if(L<=R)
				{
					int k=0;  
				    while((1<<(k+1))<=R-L+1) k++;
					ans=max(ans,max(st[L][k],st[R-(1<<k)+1][k]));
				}
				ans=max(ans,max(y-l[num[y]]+1,r[num[x]]-x+1));
				cout<<ans<<"\n";
			}
		}
	}
	return 0;
}
相关推荐
太空漫步1110 分钟前
逆序对个数
java·算法·排序算法
23级二本计科16 分钟前
分治_归并_归并排序(逆序对)
数据结构·算法·排序算法
CHENWENFEIc17 分钟前
基础排序算法详解:冒泡排序、选择排序与插入排序
c语言·数据结构·算法·排序算法·学习方法·冒泡排序
养佳工具人19 分钟前
Java各种排序
java·算法·排序算法
都适、隶仁ミ20 分钟前
【密码学】SM4算法
linux·运维·服务器·算法·网络安全·密码学·网络攻击模型
yangpipi-1 小时前
数据结构(C语言版)-4.树与二叉树
c语言·开发语言·数据结构
C++oj1 小时前
普及组集训--图论最短路径设分层图
数据结构·算法·图论·最短路径算法
In 20291 小时前
图论【Lecode_HOT100】
算法·图论
ruleslol1 小时前
java基础概念49-数据结构2
java·数据结构
荒古前1 小时前
小发现,如何高级的顺序输出,逆序输出整数的每一位(栈,队列)
数据结构·c++·算法