[Codeforces Round #717 (Div. 2)] D.Cut (倍增)

题目链接：D. Cut

题意

给一个长度为n的数组，有q次查询，每次查询给定一个区间[l,r]，问该区间能至少划分成多少个连续子区间使得子区间里的数两两互质。

题解

我们首先可以预处理得到以每个点x作为右端点向左划分一个子区间时的左端点的下标f[x]。

用$dp[i][j]$表示以j为区间右端点向左划分2^i个子区间后的左端点下标
初始化：$dp[0][x]=f[x]$
转移方程：$dp[i][j]=dp[i-1][dp[i-1][j]-1]$

（注意这里转移方程里dp[i-1][j]是以j为右端点向左划分2^(i-1)个区间后的左端点，所以转移时需要左端点-1作为新的右端点进行递推）

得到了dp的值，我们可以固定查询范围最右边r，贪心地从最长的划分区间$2^x$选取

1. 如果左端点≥左区间，说明划分的这些区间都在范围内的，答案+区间数（$2^x$）
2. 否则存在一些划分区间不在范围内，我们只需缩短区间长度。

这个思想有点像把答案转为二进制形式，从最大的往小加，一定能得到正确的结果。

小思考，为什么ans一开始赋值1？😃

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<cassert>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<deque>
#include<iomanip>
#include<list>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#include<unordered_set>
#include<unordered_map>
using namespace std;
//extern "C"{void *__dso_handle=0;}
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define pii pair<int,int>
#define lowbit(x) x&-x

const double PI=acos(-1.0);
const double eps=1e-6;
const ll mod=1e9+7;
const int inf=0x3f3f3f3f;
const int maxn=1e5+10;
const int maxm=100+10;
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);

int f[maxn],dp[25][maxn],a[maxn],last[maxn];

int main()
{
	int n,q; scanf("%d%d",&n,&q);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
	{
		int x=a[i],cnt=f[i-1];
		for(int j=2;j*j<=x;j++)
		{
			if(x%j==0)
			{
				cnt=max(cnt,last[j]+1);
				last[j]=max(last[j],i);
				while (x%j==0) x/=j;
			}
		}
		if(x^1)
		{
			cnt=max(cnt,last[x]+1);
			last[x]=max(last[x],i);
		}
		f[i]=cnt;
	}
	for(int i=1;i<=n;i++) dp[0][i]=f[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]-1];
	while (q--) {
		int l,r; scanf("%d%d",&l,&r);
		int ans=1;
		for(int i=20;i>=0;i--)
		{
			if(dp[i][r]-1>=l)
			{
				ans+=(1<<i);
				r=dp[i][r]-1;
			}
		}
		printf("%d\n",ans);
	}
}