高级ST表（区间最长不下降子序列 & 区间GCD

已于 2022-04-12 16:16:38 修改
阅读量196 收藏
点赞数
分类专栏： 练习 文章标签： c++ c语言 算法
于 2022-03-09 09:28:44 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45419976/article/details/123368205
版权

文章目录

一、区间最长不下降子序列（不一定连续

题面

在这里插入图片描述
在这里插入图片描述

Tips：

  • 子序列没说连续！！
  • Hint：最长上升子序列只需要考虑区间内是否有0,1{0, 1}0,1子序列，考虑维护两个数组LastOne和FirstZero，为了使得该hint含蓄一些，数组含义请自己脑补qwq。最长不下降子序列只有三种情况{000…000},{111…111},{000…111}，前两种情况利用前缀和维护即可，第三种情况可以枚举0,1的分界点，总复杂度为O(NM)，会超时。考虑用ST表优化这一过程，做到O(1)查询，最终复杂度为O(NlogN+M)。

T两个点代码：

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const int N = 1000005;
int n,m,opt,l,r;
int a[N],end0[N][20] = {0},end1[N][20] = {0};
int last1[N][20] = {0},first0[N][20] = {0};
int sum[N] = {0};

void get_inform(){
	int tmp = log(n)/log(2);
	for(int i=1;i<=tmp;i++)
		for(int j=1;j+(1<<i)-1<=n;j++){
			if(last1[j+(1<<(i-1))][i-1] < j+(1<<(i-1))) 
				last1[j][i] = last1[j][i-1];
			else last1[j][i] = max(last1[j][i-1] , last1[j+(1<<(i-1))][i-1]);
			if(first0[j][i-1] >= j+(1<<(i-1))) 
				first0[j][i] = first0[j+(1<<(i-1))][i-1];
			else first0[j][i] = min(first0[j][i-1],first0[j+(1<<(i-1))][i-1]);
			end0[j][i] = end0[j][i-1]+end0[j+(1<<(i-1))][i-1];
			end1[j][i] = max(end0[j][i-1]+end1[j+(1<<(i-1))][i-1]  ,
								end1[j][i-1]+sum[j+(1<<i)-1]-sum[j+(1<<(i-1))-1]);
			
		}
}

int findattack(int L,int R){
	int tmp = log(R-L+1)/log(2);
	int maxx;
	if((1<<tmp) == R-L+1) {
		return max(end0[L][tmp],end1[L][tmp]);
	}
	else maxx = max(end0[L][tmp]+findattack(L+(1<<tmp),R),end1[L][tmp]+sum[R]-sum[L+(1<<tmp)-1]);
	return maxx;
}

int findgreat(int L,int R){
	int tmp = log(R-L+1)/log(2);
	int last ,first;
	if(last1[R-(1<<tmp)+1][tmp] <= R-(1<<tmp)) 
		last = last1[L][tmp];
	else last = max(last1[L][tmp], last1[R-(1<<tmp)+1][tmp]);
	if(first0[L][tmp] >= L+(1<<tmp)) 
		first = first0[R-(1<<tmp)+1][tmp];
	else first = min(first0[L][tmp], first0[R-(1<<tmp)+1][tmp]);
	
	if(first > last) return 1;
	return 2;
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		sum[i] = sum[i-1]+a[i];
		if(a[i]==0) {
			last1[i][0] = i-1;
			first0[i][0] = i;
			end0[i][0] = 1;
			end1[i][0] = 0;
		}
		else {
			last1[i][0] = i;
			first0[i][0] = i+1;
			end0[i][0] = 0;
			end1[i][0] = 1;
		}
	}
	get_inform();//last0 & first1 
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&opt,&l,&r);
		if(opt == 1) printf("%d\n",findattack(l,r));
		else printf("%d\n",findgreat(l,r));
	}
}

是考虑用end0和end1（两个st表）来维护每个长为2k的区间以0和1结尾的最长子序列长度，这样做的好处是在预处理的时候可以很方便地处理好数据，缺点是在L,R查询的时候需要一定程度的递推（[L+(1<<tmp),R]的最长不下降子序列长度不能直接得出来，需要递推）。是从区间相互连接的端点进行考虑的。

T一个点代码：

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const int N = 1000005;
int n,m,opt,l,r;
int a[N],st[N][20] = {0};
int last1[N][20] = {0},first0[N][20] = {0};
int sum[N] = {0};

void get_inform(){
	int tmp = log(n)/log(2);
	for(int i=1;i<=tmp;i++)
		for(int j=1;j+(1<<i)-1<=n;j++){
			if(last1[j+(1<<(i-1))][i-1] < j+(1<<(i-1))) 
				last1[j][i] = last1[j][i-1];
			else last1[j][i] = max(last1[j][i-1] , last1[j+(1<<(i-1))][i-1]);
			if(first0[j][i-1] >= j+(1<<(i-1))) 
				first0[j][i] = first0[j+(1<<(i-1))][i-1];
			else first0[j][i] = min(first0[j][i-1],first0[j+(1<<(i-1))][i-1]);
			
			st[j][i] = max(st[j][i-1]+sum[j+(1<<i)-1]-sum[j+(1<<(i-1))-1]  ,
						st[j+(1<<(i-1))][i-1]+(1<<(i-1))-(sum[j+(1<<(i-1))-1]-sum[j-1]));
			
		}
}

int findattack(int L,int R){
	int tmp = log(R-L+1)/log(2);
	int maxx;
	maxx = max(st[L][tmp]+sum[R]-sum[L+(1<<tmp)-1] , 
				st[R-(1<<tmp)+1][tmp]+(R-(1<<tmp)-L+1)-(sum[R-(1<<tmp)]-sum[L-1]));
	return maxx;
}

int findgreat(int L,int R){
	int tmp = log(R-L+1)/log(2);
	int last ,first;
	if(last1[R-(1<<tmp)+1][tmp] <= R-(1<<tmp)) 
		last = last1[L][tmp];
	else last = max(last1[L][tmp], last1[R-(1<<tmp)+1][tmp]);
	if(first0[L][tmp] >= L+(1<<tmp)) 
		first = first0[R-(1<<tmp)+1][tmp];
	else first = min(first0[L][tmp], first0[R-(1<<tmp)+1][tmp]);
	
	if(first > last) return 1;
	return 2;
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		sum[i] = sum[i-1]+a[i];
		if(a[i]==0) {
			last1[i][0] = i-1;
			first0[i][0] = i;
		}
		else {
			last1[i][0] = i;
			first0[i][0] = i+1;
		}
		st[i][0] = 1;
	}
	get_inform();//last0 & first1 
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&opt,&l,&r);
		if(opt == 1) printf("%d\n",findattack(l,r));
		else printf("%d\n",findgreat(l,r));
	}
}

昨天半夜想到的……闭上眼睛开始睡觉，突然意识到应该怎样O(1)查询：st[i][j]存的是[i,i+(1<<j)-1]区间最长不下降子序列长度，在合并或者查询的时候，对于大的区间分成两部分（可能有重叠区域）则0和1的分界点一定属于其中某一个区间，所以用max(st[左区间]+右边1数,st[右区间+左边0数])可以完美考虑到所有的情况。是从0，1分界点的角度出发进行考虑的。
可惜还是T掉了一个点，咋整呢
经过一些测试证明是常数大了，要减小常数

AC代码：

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const int N = 1000005;
int n,m,opt,l,r;
int a[N],st[N][20] = {0};
int last1[N][20] = {0},first0[N][20] = {0};
int sum[N] = {0};

void get_inform(){
	int tmp = log(n)/log(2);
	for(int i=1;i<=tmp;i++)
		for(int j=1;j+(1<<i)-1<=n;j++){
			if(last1[j+(1<<(i-1))][i-1] < j+(1<<(i-1))) 
				last1[j][i] = last1[j][i-1];
			else last1[j][i] = max(last1[j][i-1] , last1[j+(1<<(i-1))][i-1]);
			if(first0[j][i-1] >= j+(1<<(i-1))) 
				first0[j][i] = first0[j+(1<<(i-1))][i-1];
			else first0[j][i] = min(first0[j][i-1],first0[j+(1<<(i-1))][i-1]);
			
			st[j][i] = max(st[j][i-1]+sum[j+(1<<i)-1]-sum[j+(1<<(i-1))-1]  ,
						st[j+(1<<(i-1))][i-1]+(1<<(i-1))-(sum[j+(1<<(i-1))-1]-sum[j-1]));
			
		}
}

void findattack(){
	int tmp = log(r-l+1)/log(2);
	int maxx = max(st[l][tmp]+sum[r]-sum[l+(1<<tmp)-1] , 
				st[r-(1<<tmp)+1][tmp]+(r-(1<<tmp)-l+1)-(sum[r-(1<<tmp)]-sum[l-1])); //声明时赋值应该会更快 
	printf("%d\n",maxx);
}

void findgreat(){
	int tmp = log(r-l+1)/log(2);
	int last ,first;
	if(last1[r-(1<<tmp)+1][tmp] <= r-(1<<tmp)) 
		last = last1[l][tmp];
	else last = max(last1[l][tmp], last1[r-(1<<tmp)+1][tmp]);
	if(first0[l][tmp] >= l+(1<<tmp)) 
		first = first0[r-(1<<tmp)+1][tmp];
	else first = min(first0[l][tmp], first0[r-(1<<tmp)+1][tmp]);

	if(first > last) printf("1\n");
	else printf("2\n");
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		sum[i] = sum[i-1]+a[i];
		if(a[i]==0) {
			last1[i][0] = i-1;
			first0[i][0] = i;

		}
		else {
			last1[i][0] = i;
			first0[i][0] = i+1;
		}
		st[i][0] = 1;
	}
	get_inform();//last0 & first1 
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&opt,&l,&r);
		if(opt == 1) findattack();//函数传参全部删了，还改成了void类型直接输出，减少参数传来传去的复杂度 
		else findgreat();
	}
	return 0;
}

修改点：

  • 可以声明时赋值的都声明即赋值了
  • 查询时用到的函数全部变成了void类型，直接在内部输出，进去的参数也删掉了，全部用全局变量，减少了传参的复杂度，竟然过了hhh

二、区间GCD

题面

在这里插入图片描述

GCD怎么用st表呢？因为GCD是可以合并的啊……数学不好一时间没想出来……
害 下面放代码

AC代码：

#include<iostream>
#include<cmath>
using namespace std;
const int N = 500005;
int n,m,L,R;
int a[N] = {0};
int st[N][20] = {0};

int gcd(int x,int y){
	int r = x%y;
	while(r){
		x = y;
		y = r;
		r = x%y;
	}
	return y;
}

void getst(){
	int tmp = log(n)/log(2);
	for(int i=1;i<=tmp;i++){
		for(int j=1;j+(1<<i)-1<=n;j++){
			st[j][i] = gcd(st[j][i-1],st[j+(1<<(i-1))][i-1]);
		}
	}
	return ;
}

int binary1(int l,int r,int target){//找最大值 
	int mid;
	while(l<=r){
		mid = (l+r)>>1;
		int tmp = log(mid-L+1)/log(2);
		if(gcd(st[L][tmp],st[mid-(1<<tmp)+1][tmp]) >= target)
			l = mid+1;
		else r = mid-1;
	}
	return l-1;
}

int binary2(int l,int r,int target){//找最小值 
	int mid;
	while(l<=r){
		mid = (l+r)>>1;
		int tmp = log(mid-L+1)/log(2);
		if(gcd(st[L][tmp],st[mid-(1<<tmp)+1][tmp]) <= target)
			r = mid-1;
		else l = mid+1;
	}
	return r+1;
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		st[i][0] = a[i];
	} 
	getst();
	for(int i=1;i<=m;i++){
		scanf("%d%d",&L,&R);
		int tmp = log(R-L+1)/log(2);
		int ans = gcd(st[L][tmp],st[R-(1<<tmp)+1][tmp]);
		printf("%d %d\n",ans,binary1(L,n,ans)-binary2(L,n,ans));
	}
	return 0; 
}

思路

在这里插入图片描述
在这里插入图片描述

矿矿不想吃饭
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
ACM.algorithm.rar_GCD矩阵_匈牙利_区间匹配_最短路 三维_模板元
09-14
各种算法模板（二分图最大匹配匈牙利算法、最小生成树prime和kruskal算法、Dijkstra算法、两点最短路径负权值边SPFA算法、图任意两点最短路径Floy算法、网络最大流SAP算法、网络最大流最小费用算法、乘法逆元gcd扩展...
A. A Twisty Movement(翻转区间求最长不下降子序列)
代码笔记本
09-12 381
https://codeforces.com/problemset/problem/933/A 题意翻译 给定一个序列 A，你可以翻转其中的一个区间内的数，求翻转后的序列的最长不下降子序列的长度。（|A|\le 2000,1\le a_i \le 2∣A∣≤2000,1≤ai​≤2） 感谢@touristWang 提供的翻译 输入输出样例 输入 #1复制 4 1 2 1 2 输出 #1复制 4 输入 #2复制 10 1 1 2 2 2 1 1 2 2 1 输出 #2复制 .
【线性动态规划】最长不下降子序列（LIS）
DC_Fan的博客
01-25 570
最长不下降子序列（线性动态规划）    我们知道，递推就是按照递推公式不断往前求解的一个过程，对于当前，只有一个状态描述当前的值。若某些问题并非由一个状态而是由多个状态不断推导，那么这种方法就是动态规划，简称DP。动态规划是运筹学的一个分支，是将问题分解成各个阶段，由相邻的两个阶段根据状态转移方程推到求解的一种方法。    所谓的线性动态规划，就是该问题模型是线性的，数据结构表现为线性表的形式。 题目 题干：    最长不下降子序列。给定一个长度为n的序列a，求出这个序列中最长不下降子序列，所谓最长不下降子
求序列最长不下降子序列_最长不下降子序列nlogn算法详解
weixin_39667626的博客
12-19 477
今天花了很长时间终于弄懂了这个算法……毕竟找一个好的讲解真的太难了，所以励志我要自己写一个好的讲解QAQ这篇文章是在懂了这个问题n^2解决方案的基础上学习。解决的问题：给定一个序列，求最长不下降子序列的长度(nlogn的算法没法求出具体的序列是什么)定义：a[1..n]为原始序列，d[k]表示长度为k的不下降子序列末尾元素的最小值，len表示当前已知的最长子序列的长度。初始化：d[1]=a[1];...
小黄的刷题之路(十八)——ST表求解区间gcd问题
m0_56543373的博客
11-10 1105
ST表，稀疏表，求解区间gcd问题，事实上ST表可以求解所有可重复贡献问题
子序列-反转区间求最长不下降子序列
PushyTao的博客
07-25 413
题目描述 小Z有一个01序列A=(A1,A2,A3,…,An)。他可以进行一次操作：选择一个区间L,R将其反转。 例如，对于序列A=(1,0,0,1,1)，当L=2,R=4时，原序列将变为(1,1,0,0,1)。 小Z希望：通过这一次反转操作，使得这个序列的最长不下降子序列的长度尽可能的大，并想让你输出这个最大值。 一个序列的不下降子序列定义为：对于一个序列(p1,p2,…,pk)满足1≤p1<p2<…<pk≤n(n≤819200)且Ap1≤Ap2≤…≤Apk。则序列(Ap1,Ap2,…,
Gcd&&Lcm(1).pptx
02-09
这是一个关于学习最大公约数与最小公倍数的PPT，里面包含了程序编写以及各类知识点。知识点简练易懂。同时还可以到本人博客（寒假培训——GCD，LCM）中查询有关的例题链接并配有AC代码。
GCD常用和不常用API说明GCD-Demo-master.zip
最新发布
04-12
GCD常用和不常用API说明GCD_Demo-master.zip
GCD的高级用法
05-12
GCD的一些高级用法，是不容错过的多线程开发学习材料
计算GCD的简单java程序
07-12
创建最常用的除数计算器；。。。GCD计算器允许您快速查找一组数字的最大公约数；。。。当两个数的最大公约数为1时，这两个数称为互质或相对素数。
Codeforces Round #736 (Div. 2) D ST表维护区间gcd，二分
qq_46264636的博客
08-02 130
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=2e5+10,mod=1e9+7; #define PII pair<int,int> // #define x first // #define y second #define PB push_back bool vis[N]; ll ST[N][100]; ll b[N],a[N]; int n; void ST_in
最长不下降子序列（上升同理）
qpfxgg的博客
05-11 406
最长不下降就是一条个数最多的，不下降（可以相等）的序列 比如 1 1 1 1 最长不下降子序列就是4 13，7，9，16，38，24，37，18，44，19，21，22，63，15 最长不下降子序列就是7，9，16，18，19，21，22，63 下面我们来分析一下思路 我们使用一个二维数组s[10][3],每一行的第一个代表个数，第几个，第二个代表他自身数据。第三行代表从他开始的不下降子序列的长度，第四个代表从他开始的不下降子序列的下一个元素的坐标 下面我们使用上面的例子来详细说明一下算法过程 有这样一组
【网络流24题】 最长不下降子序列问题
hipamp的博客
02-02 198
题意： 给定正整数序列x1,...,xnx_1,...,x_nx1​,...,xn​ （1）计算其最长不下降子序列的长度 sss。 （2）计算从给定的序列中最多可取出多少个长度为s的不下降子序列。 （3）如果允许在取出的序列中多次使用 x1x_1x1​ 和 xnx_nxn​，则从给定序列中最多可取出多少个长度为 sss 的不下降子序列。 分析 第一问，普及组问题。 第二，三问，要用网络流解决。 第...
静态区间 / gcd区间【ST表】
//(^ o ^)//
07-06 256
>Link ybtoj静态区间 luogu P1890 >解题思路 其实就是ST表的模板套上一个gcd就行了 算了一下时间复杂度，大概O(nlognlogai)O(nlognloga_i)O(nlognlogai​)（?）差不多可以过的 >代码 #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #define N 50010 usin
最长不下降子序列（可以改成最长上升子序列）
lyc1635566ty的博客
07-29 678
#include #define maxn 101111 using namespace std; int a[maxn],b[maxn]; int Search(int num,int low,int high) { int mid; while(low<=high) { mid=(low+high)/2; if(num>=b[mid]) low=
luoguP1890 gcd区间 [st表][gcd]
dirushi6322的博客
08-01 102
题目描述 给定一行n个正整数a[1]..a[n]。 m次询问，每次询问给定一个区间[L,R]，输出a[L]..a[R]的最大公因数。 输入输出格式 输入格式： 第一行两个整数n，m。 第二行n个整数表示a[1]..a[n]。 以下m行，每行2个整数表示询问区间的左右端点。 保证输入数据合法。 输出格式： 共m行，每行表示一个询问的答案。 输入输出样例 输入样...
【GCD】【ST表】静态区间
blog
07-06 262
给出n个数a1,a2,...an m个询问，每次给出l,r，输出al,al+1,...ar-1,ar的最大公因数
GCD （ST表，二分求区间查询）
yhn19951008的博客
07-22 1949
HUD 5726 GCD    给一个序列，多次查询区间的最大公约数，并求出同样是这个最大公约数的区间有多少个。 区间查询采用ST表，第二问查询利用区间向右延伸最大公约数递减的规律可通过二分快速找到右边界。把第一问的答案先求出来，表示要查询的公约数加入map。 然后枚举左端点二分右端点，把已经从存在map的值加上找到的区间数，不存在的就不用管了。 #include #inclu
gcd区间 (ST表)(爱思创算法四)
HackerYY的博客
12-04 451
爱思创GCD区间(ST表) 内附AC代码
C++在一个区间里生成两个数的最大公约数
05-14
可以使用欧几里得算法（辗转相除法）来求解最大公约数。代码如下： ```c++ #include <iostream> using namespace std; int gcd(int a, int b) { if (b == 0) { return a; } return gcd(b, a % b); } int main() { int l, r; cin >> l >> r; for (int i = l; i <= r; i++) { for (int j = i + 1; j <= r; j++) { int g = gcd(i, j); cout << "gcd(" << i << ", " << j << ") = " << g << endl; } } return 0; } ``` 这里我们使用了两个循环来遍历区间内的所有数对，并计算它们的最大公约数。其中，`gcd` 函数使用递归实现欧几里得算法。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值