181021 gcd区间 题解

题目描述

给定一行n个正整数a[1]..a[n]。

m次询问，每次询问给定一个区间[L,R]，输出a[L]..a[R]的最大公因数。

输入格式

第一行两个整数n，m。

第二行n个整数表示a[1]..a[n]。

以下m行，每行2个整数表示询问区间的左右端点。

保证输入数据合法。

输出格式

共m行，每行表示一个询问的答案。

样例输入

5 3
4 12 3 6 7
1 3
2 3
5 5

样例输出

1
3
7

问题提示

对于30%的数据，n <= 100， m <= 10

对于60%的数据，m <= 1000

对于100%的数据，1 <= n <= 1000，1 <= m <= 1,000,000

思路

这是一道典型的RMQ问题，一般使用线段树来解决，但是此题输入的数无需改变，且线段树代码量太大，所以我们可以使用ST表加GCD来解决

先来讲解一下ST表：

该算法是基于倍增思想设计的在线算法，之所以称之为倍增，是因为算法记录从每个元素开始的连续的长度为2^k的区间内元素中的最值

类似于动态规划，st[i,j]表示从第i个数起连续2^j个数中的最值，即记录区间[i,i+2^j-1]中的最值

把st[i,j]平均分成两段（st[i,j]一定是偶数个数字），从i到i+2^(j-1)-1为一段，i+2^(j-1)到i+2^j-1为一段，长度均为2^(j-1)

于是得到状态转移方程为：

st[i,j]=max/min{st[i,j-1],st[i+2j-1,j-1]},其中st[i,0]=A[i]

由公式可以看出在写代码的时候应该先完成列，在计算行

接下来ST表模板

#include <cstdio>
#include <iostream>
#include <cmath>
using namespace std;
#define F(i,a,b) for(int i=(a);i<=(b);i++)
#define f(i,a,b) for(int i=(a);i>=(b);i--)
const int N=1e5+5;
int Log[N],st[N][20],n,m;
signed main()
{
    scanf("%d%d",&n,&m);
    F(i,1,n) scanf("%d",&st[i][0]);//动态规划边界
    F(i,2,n) Log[i]=Log[i/2]+1;
    F(j,1,Log[n])//描述长度(1>>j)<=n
        for(int i=1;i+(1<<j)-1<=n;i++)//描述下标位置 i~i+(1<<j-1)-1
            st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]);//状态转移方程
    return 0;
}

然后是GCD部分

GCD既是最大公约数，可以用辗转相除法实现，时间复杂度O(logn)

int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);//辗转相除法求最大公因数
}

最后将两部分结合即可

code:

#include <cstdio>
#include <iostream>
#include <cmath>
using namespace std;
#define F(i,a,b) for(int i=(a);i<=(b);i++)
#define f(i,a,b) for(int i=(a);i>=(b);i--)
const int N=1e5+5;
int Log[N],st[N][20],n,m;
int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);
}
signed main()
{
	scanf("%d%d",&n,&m);
	F(i,1,n) scanf("%d",&st[i][0]);
	F(i,2,n) Log[i]=Log[i/2]+1;
	F(j,1,Log[n])
		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]);//将求最大值变成求GCD 
    //处理输出
	int L,R;
	while(m--)
	{
		scanf("%d%d",&L,&R);
		int Lg=Log[R-L+1];
		printf("%d\n",gcd(st[L][Lg],st[R-(1<<Lg)+1][Lg]));//将求最大值变成求GCD 
	}
	return 0;
}