Codeforces Round #717 (Div. 2) D. Cut

题目：D. Cut

题意:给定一个n个数的数组，q次询问，问每次[l,r]区间需要分成多少个区间，每个区间的乘积是区间所有的的最小公倍数。

首先分析区间所有数的乘积是所有数的最小公倍数，那么成立每一个数之间gcd都是1，每个数都是互质的，唯一分解可以证明。n是1e5的，先进行分解质因数，比如到第i个数，取他的左边第一个与他不互质的数，那么以i为点的这个数就分割区间完了，之后再往左递推就可以计算出需要分割成几个区间。考虑一下极端情况，递推过去每一个数都是不互质的，复杂度为O(n)，超时。考虑化简利用倍增。

#include<cmath>
#include<iostream>
#include<sstream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<istream>
#include<vector>
#include<stack>
#include<set>
//#include<unordered_map>
#include<map>
#include<algorithm>
#include<queue>
#define mmp make_pair
#define inf 0x3f3f3f3f
#define llinf 0x7fffffffffffffff
using namespace std;
typedef long long  ll;
typedef pair<int,int> PP;
typedef double lld;
const ll mod=1e9+7;
const int MAXN=2e5+2;
int a[MAXN];
int first[MAXN],dp[22][MAXN],f[MAXN];
//f:最左边不互质的数
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) {
        f[i]=f[i-1];
        for(int j=2;j*j<=a[i];++j) {
            if(a[i]%j==0) {
                f[i]=max(f[i],first[j]);
                first[j]=i;
                while(a[i]%j==0) a[i]/=j;
            }
        }
        if(a[i]^1) {
            f[i]=max(first[a[i]],f[i]);
            first[a[i]]=i;
        }
        dp[0][i]=f[i];//以i为右区间端点 向左分割2^0个 的最左区间
    }
    for(int i=1;i<=20;++i) {
        for(int j=1;j<=n;++j) {
            dp[i][j]=dp[i-1][dp[i-1][j]];
        }
    }
    while(q--) {
        int l,r;
        scanf("%d %d",&l,&r);
        int ans=0;
        for(int i=20;i!=-1;--i) {
            if(dp[i][r]>=l) {
                ans+=(1<<i);
                r=dp[i][r];
            }
        }
        ++ans;
        printf("%d\n",ans);
    }
    return 0;
}