luogu P1890ybtojRMQ问题课堂过关例题2静态区间 & gcd区间

Posted 2021-06-27 SSL_ZZL

【例题2】静态区间 & gcd区间

• Link
• 题目大意
• 解题思路
• Code

---

Link

luogu P1890 gcd区间
ybtoj【RMQ问题课堂过关】【例题2】静态区间
题面//因为不知道侵不侵权所以就是题面是私密的，有账号的直接看转送门就可了

---

题目大意

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

---

解题思路

类似这种两个区间的ans值再取ans，就是这个大区间得ans值，都可以用st表（大概吧…）

---

Code

#include <iostream>
#include <cstdio>

using namespace std;

int n, m, x, y, lg[100100], f[100100][30];

void log() {
	int root = 2, i = 1, now = 2;
	for(; now * 2 <= n; i++) {
		while(root < now *2)
			lg[root++] = i;
		now *= 2;
	}
	while(root <= n)
		lg[root++] = i;
}

int gcd(int x,int y) {  //把max改成了gcd就过了......
	int r = x % y;
	while(r) {
		x = y, y = r;
		r = x % y;
	}
	return y;
}

int main() {
	scanf("%d %d", &n, &m);
	log();
	for(int i = 1; i <= n; i++)
		scanf("%d", &f[i][0]);
	for(int j = 1; j <= lg[n]; j++)
		for(int i = 1; i <= n; ++i)
			if (i + (1 << (j - 1)) > n) f[i][j] = f[i][j - 1];
				else f[i][j] = gcd(f[i][j - 1], f[i + (1 << (j - 1))][j - 1]);
	for(int i = 1; i <= m; i++) {
		scanf("%d %d", &x, &y);
		int len = lg[y - x + 1];
		printf("%d\\n", gcd(f[x][len],f[y - (1 << len) + 1][len])); 
	}
}