Codeforces Round #717 (Div. 2) D. Cut（倍增）

题目链接：https://codeforces.com/problemset/problem/1516/D

从右往左，记录当前位置数最近的不互质数的位置，设当前位置为l，标记最近位置为go[l]，那么从go[l]到l这一块可以被分出来。（注意：go[l] >= go[l - 1]）

这样就可以每次询问时，从最右边开始跳，x = r, x = go[x]直到x<l，但是这么跳很可能会超时，比如整个序列所有的数都是一样的。这里用倍增的思想，dp[i][x]表示从x位置连续跳(1<<i)次到达的位置。
要实现这一点，转移方程为，dp[i][j]=dp[i - 1][dp[i - 1][j]]，那么O(logn)就可以求出一次询问的答案。

#include <bits/stdc++.h>
using namespace std;

using ll = long long;
const int N = 1e5 + 5;

int n, q, a[N], go[N], last[N], dp[30][N]; 
int cnt, vis[N], p[N], s[N];
// go[l]记录从l开始往左走第一个不与a[l]互质的位置, 但不可以小于go[l-1] 

void init() {
	s[1] = 1;
	for (int i = 2; i < N; i++) {
		if (!vis[i]) p[++cnt] = i, s[i] = i;
		for (int j = 1; j <= cnt; j++) {
			if (i * p[j] >= N) break;
			vis[i * p[j]] = 1;
			s[i * p[j]] = p[j];
			if (i % p[j] == 0)break;
		}
	}
}

int main(void) {
//	freopen("in.txt", "r", stdin);
	init();
	scanf("%d%d", &n, &q);
	for (int i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
		int x = a[i];
		go[i] = go[i - 1]; // go[i] >= go[i - 1]
		int bfprime = -1;
		while (s[x] != 1) {
			if (bfprime != s[x]) {
				go[i] = max(go[i], last[s[x]]);
			}
			last[s[x]] = i;
			bfprime = s[x];
			x /= s[x];
		}
		dp[0][i] = go[i]; 
	}
	for (int i = 1; i <= 20; i++)
		for (int j = 1; j <= n; j++)
			dp[i][j] = dp[i - 1][dp[i - 1][j]];
			// dp[i][j] 记录从位置j开始，连续跳 1 << i 次的位置 
	while (q--) {
		int l, r;
		scanf("%d%d", &l, &r);
		int ans = 0;
		for (int i = 20; i >= 0; i--)
			if (dp[i][r] >= l)
				ans += (1 << i), r = dp[i][r];
		printf("%d\n", ans + 1);
	}
	return 0;
}