2021-07-10 ST表专项

ST表专项

写在前面：

ST表，怎么说呢，不难，还较为简单，如果DP学的好，就完全不需要担心。

ST表的原理是依靠二进制做预处理，然后使用动态转移的方法快速接近数据的预处理，然后在位运算得出答案。

总之，记住下面这个动态转移方程公式，就可以AKST表。
$$f_{ij}=\max (f_{i{j}_{-1}},f_{i_{+1<<(j-1)}{j}_{-1}})$$
这个便是ST表的动态转移方程，不难看出ST表是基于预处理和DP的一种区间查询算法，速度极快，但是，只能用于查询，并不支持修改的操作，修改均需要提前做好。

ST表主要用于RMQ问题，可以来进行区间查询。（不会ST表的话，其实可以单调队列查询（小声哔哔））

只讲定理，毫无意义，直接给大家看一下原题。

T1:

题目描述：

2008年9月25日21点10分，酒泉卫星发射中心指控大厅里，随着指挥员一声令下，长征二号F型火箭在夜空下点火起飞，神舟七号飞船载着翟志刚、刘伯明、景海鹏3位航天员，在戈壁茫茫的深邃夜空中飞向太空，开始人类漫步太空之旅。第583秒，火箭以7.5公里／秒的速度，将飞船送到近地点200公里、远地点350公里的椭圆轨道入口。而此时，火箭的燃料也消耗殆尽，即将以悲壮的方式与飞船告别。这个过程，在短短不到10分钟时间内，翟志刚和他的两名战友体会到了从超重到失重的过程。

除了超重和失重的感觉之外，就是浩瀚的长空中璀璨的星星，和地面上看到的星星不同，在太空中看到的星星是成一条直线的，一共有N(1<=N<=100,000)颗星星，编号为1到N,每个星星有自己的体积，由于在飞船中很无聊，除了不停地玩弄手中失重的书和笔之外没有别的事可干，此时翟志刚说我们来玩游戏吧，一共玩了M轮(1<=M<=100,000)，每一轮都是给出两个整数L和R(1<=L<=R<=N)，询问第L到第R颗星星之间最大星星的体积，每次答对的人就可以多休息一段时间。

由于翟志刚还要进行太空漫步，所以他现在请你帮忙，你得到的回报就是太空饼干。

正解：

ok，妥妥的裸题，我们就用这道题来练练练手吧！

首先，如果你认真看题的话，显然，这是一道标准的区间查询RMQ问题，用线段树也可以做出来，不过，我们先讲一讲ST表吧。

首先，三部曲：初始化，动态转移方程，实现

我们设
$$f_{i,j}=表示在区间i到2的j次方中的最大值$$
还记得上文的方程吗？这里，我们便可直接套用
$$f_{ij}=\max (f_{i{j}_{-1}},f_{i_{+1<<(j-1)}{j}_{-1}})$$
还差做后一步，实现！其实没什么好说的，数组方程都好了（小声哔哔…）

来看看代码吧！

code：

#include<bits/stdc++.h>
using namespace std;
int x,y,z,m,n,a[500001],f[500001][20],k,sum[500001]; 
int main()
{
	scanf("%d%d",&m,&n);
	for(int i=1;i<=m;i++) scanf("%d",&a[i]);
	for(int i=1;i<=m;i++) f[i][0]=a[i];
	for(int j=1;j<=17;j++)
		for(int i=1;i<=m;i++)
			f[i][j]=max(f[i][j-1],f[i+ (1<< (j-1) ) ][j-1]);
	while(n--)
	{
		scanf("%d%d",&x,&y);
		z=0;
		while(1<<(z+1)<=y-x+1)z++;
		sum[++k]=max(f[x][z],f[y-(1<<z)+1][z]);
	}
	for(int i=1; i<=k; i++)
		printf("%d\n",sum[i]);
	return 0;
}

T2:

题目描述：

经过几个月辛勤的工作，FJ决定让奶牛放假。假期可以在1…N天内任意选择一段（需要连续），每一天都有一个享受指数W。但是奶牛的要求非常苛刻，假期不能短于P天，否则奶牛不能得到足够的休息；假期也不能超过Q天，否则奶牛会玩的腻烦。FJ想知道奶牛们能获得的最大享受指数。

正解：

（出现了，三大信息学公敌！——Farmer John 和 他的奶牛！）

首先，如果你认真看题的话，显然，这还是一道标准的区间查询RMQ问题，不过看题中…“每一天都有一个享受指数W”

emm…前缀和！！！！

首先，三部曲：初始化，动态转移方程，实现

我们设
$$f_{i,j}=表示在天数i到2的j次方中的最大值$$
还记得上文的方程吗？这里，我们便可稍微更改一番
$$f_{ij}=\max (f_{i{j}_{-1}},f_{i_{+d_{j-1}}{j}_{-1}})$$
还差做后一步，实现！其实没什么好说的，数组方程都好了（小声哔哔…）

来看看代码吧！

code：

#include<bits/stdc++.h>
using namespace std;
const int N=100005;
long long a[N],s[N],f[N][30],d[30];
long long stb(int l,int r)
{
    int x=log2(r-l+1);
    return max(f[l][x],f[r-d[x]+1][x]);
}
int main()
{
	d[0]=1;
	for(int i=1;i<=25;i++) d[i]=d[i-1]*2;
	int n,p,q;
	long long ans=-0x7ffffffffffff;
	scanf("%d%d%d",&n,&p,&q);
	for(int i=1;i<=n;i++) 
	{
		scanf("%lld",&a[i]);
		s[i]=s[i-1]+a[i];
		f[i][0]=s[i];
	}
	for(int j=1;d[j]<=n;j++)
	{
		for(int i=1;i+d[j]-1<=n;i++) 
		{
			f[i][j]=max(f[i][j-1],f[i+d[j-1]][j-1]);
		}
	}
	for(int i=1;i<=n-p+1;i++)
	{
		long long cnt=stb(i+p-1,min(n,i+q-1));
		ans=max(ans,cnt-s[i-1]);
	}
	printf("%lld",ans);
	return 0;
}

写在最后：

这篇是我今天更新的第二篇啦！

不过，还差一个专题：树状数组！

来！等一下！记得支持哦！

如有错误，请评论区提出！

谢谢你们！！