思路:对于每一个询问的区间,如果要暴力解决的话,可以贪心的解决:从左端点到右端点取多个区间,每次直到不能融入同一个区间的时候才开一个新区间。
但是这样显然复杂度过不了,所以可以预处理每个点往右最多遍历到哪个位置可以放在同一个区间中。注意能遍历到哪个位置受到两个因素影响:1.该数自身不能遇到有相同因子的。2.该数右边的在区间内的数也不能遇到有相同因子的。
举个例子:2,3,9,4这个数列,考虑预处理2往右最多延伸到哪个数可以放进同一区间。2右边第一个有相同因子的数是4,所以往右的idx<4。同时,3和9有相同因子,所以往右的idx<3,两者取min,所以这个2往右最多延伸到idx=2,也就是只能取【2,3】。
我们可以用nxt数组记录往右延伸的最后一个位置的下一位置,比如上面的2,3,9,4就有nxt[1] = 3。
这样复杂度看起来好一些了,但是还是过不了题(实际上,这个做法能被轻易地卡成O(nq),显然不行)。
所以我们就要用到倍增大法。简单来说就是:走一次nxt是一次小跳,小跳太慢了,所以弄一些大跳:一次跳2步,4步,8步...这样就能变成O(qlogn)。
详见代码。(计算nxt的时候要处理一些细节内容,主要是质因数分解,仔细看看应该不难理解)
#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for (int i = (a); i <= (b); i++)
#define ROF(i, a, b) for (int i = (a); i >= (b); i--)
const int N = 1e5+5;
int n,q,a[N],nxt[N];
int f[30][N];
map<int,int> mp;
inline void init(){
nxt[n+1]=n+1; //末尾标记(用于nxt数组的单调性递推,用到性质nxt[i-1] <= nxt[i])
ROF(i,n,1){
int res = nxt[i+1];
int t=a[i];
FOR(j,2,sqrt(a[i]))
while(t%j==0){
if(mp[j] && mp[j]!=i) res=min(res,mp[j]);
mp[j] = i; //标记
t/=j;
}
if(t!=1 && mp[t] && mp[t]!=i) res=min(res,mp[t]);
mp[t] = i; //标记
nxt[i] = res;
}
FOR(i,0,20) f[i][n+1]= n+1; //末尾标记(用于ask函数中,跳到n+1位置后必定>r,此时需要返回)
FOR(i,1,n) f[0][i] = nxt[i];
FOR(i,1,20) FOR(j,1,n)
f[i][j] = f[i-1][f[i-1][j]];
}
int ask(int l,int r){
int cnt=1;
for(int i=20; i>=0; i--){
if(f[i][l] <= r){
l=f[i][l];
cnt+=(1<<i);
}
}
return cnt;
}
signed main(){
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin>>n>>q; FOR(i,1,n) cin>>a[i];
init();
int l,r;
while(q--){
cin>>l>>r;
cout<<ask(l,r)<<'\n';
}
}