Codeforces Round #717 (Div. 2)-D. Cut-（双指针+倍增）

D. Cut

双指针 动态规划 倍增

---

分析题意可知，满足条件的分块中，所有元素均互质，即：gcd（a[L] ~ a[R] ）= 1。

先考虑怎么分块，从暴力法入手，可以枚举1~n为起点，往右遍历，每添加一个元素，这个元素的所有质因子都添加进容器表示已存在，当下一个要添加的元素的质因子已经出现了，那就不能添加，起点的最长分块也结束。

考虑用双指针，当要添加的元素，的任意一个质因子已经出现过，代表这个元素要加进来就不能保证gcd全为1了，
则需要把前方左指针对应元素的质因子全部删掉，且左指针++，直到这个要添加的元素的质因子都未出现过。

双指针尺取伪代码模板：

int L=1,R=1;
while(R<=n)
{
    add(a[R]);
    while(当前状态不符合)
    {
        delete(a[L]);
        L++;
    }
    记录(a[L]~a[R])对应答案;
    R++;
}

得出所有点跳一次能到的最远距离，图解：

后面要求找出询问下的【Li，Ri】的最小分块数，
基于我们已经预处理好了每个点一步能跳到的最远距离，显然暴力从区间内起点开始，并贪心的按双指针预处理好的最远距离往后跳，跳了多少次到达右边界答案就是多少，
但是这样是会T的，因为最坏情况下，1e5次查询，每次查询都是【1,1e5】，时间复杂度就炸了。

考虑用倍增思想进行优化，见状态转移方程代码：

最后每次查询的时候也是倍增的连跳但不碰边界，属于倍增的模板：

---

代码：

int n,m,k,q,a[maxn],b[maxn];
int dp[maxn][33],cnt[112345];
//dp[i][j]从i开始,跳 1<<(j-1) 次的最右可达距离+1

void del(int x)
{
    for(int i=2;i*i<=x;i++)
        while(x%i==0)
            x/=i,cnt[i]--;
    if(x>1) cnt[x]--;
}
void solve()
{
    cin>>n>>q;
    Qin(a);
    int L=1,R=1;
    while(R<=n)
    {
        //add--探测区间[l,R-1]添加a[R]
        int x=a[R];
        for(int i=2;i*i<=x;i++)
        {
            if(x%i==0)
            {
                while(cnt[i])//del左侧
                {
                    dp[L][1]=R;
                    del(a[L]);
                    L++;
                }
                while(x%i==0)
                    x/=i,cnt[i]++;
            }
        }
        if(x>1)
        {
            while(cnt[x])//del左侧
            {
                dp[L][1]=R;
                del(a[L]);
                L++;
            }
            cnt[x]++;
        }
        R++;
    }
    for(int i=L;i<R;i++)    dp[i][1]=R; //最后一段向右刷新(此R=n+1
    for(int i=0;i<=30;i++)  dp[n+1][i]=n+1;//防止越界
    
    //dp[i][j]从i开始,跳 1<<(j-1) 次的最右可达距离+1
    for(int i=n;i>=1;i--)
    {
        for(int j=2;j<=30;j++)
            dp[i][j] = dp[min(n+1,dp[i][j-1])][j-1];
        //起始点st跳2^j次能到的位置,
        //等于从st跳2^(j-1)次到达的st'开始,再跳2^(j-1)次能到的位置.
    }
    
    int x,y,ans;
    while(q--)
    {
        cin>>x>>y;
        ans=0;
        for(int i=30;i>=1;i--)
        {
            if(dp[x][i]<=y)//跳1<<(j-1)次的最右可达距离+1
            {
                ans += (1<<(i-1));
                x = dp[x][i];
            }
        }
        printf("%d\n",ans+1);
    }
}
int main()
{
	IOS;
//	int t;cin>>t;while(t--)
	solve();
	return 0;
}