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

D. Cut

感觉这场的题思维上都非常妙，一开始建模想着用网络流跑二分匹配，但是觉得码量过大不太可能是cf正解，看题解果然，贪心就行，另外用到了倍增，感觉很不chuo。
官方题解地址

题意

给出$n$个数和$q$个询问。

对于每个询问，给出一个区间，问你最少把这个区间分成多少份（连续），使得每个区间内的数都两两不互质。

想法

可以证明，用贪心的方法，对于每个开始数找一个最大区间作为划分，接着这个区间后的第一个数最为新区间的第一个数即可得到最优解（最少划分）。

预处理方法为倒序处理，定义数组$go[i]$为下标$i$跳转到的下一个区间第一个数的下标，当处理第$i$个数的时候，遍历这个数的所有因子（质），找到有这个因子的下一个位置，更新为$go[i]=min(go[i],nxt[j])$，同时更新$nxt[j]$。

有了以上跳转数组之后，考虑$log(n)$的方法去跳转遍历区间，用到倍增方法（有一定单调性）。倍增处理方式参考lca等经典思想，直接见代码，如下。

代码

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>

using namespace std;
typedef long long ll;
const int maxn=1e5+5;
const int inf=0x3f3f3f3f;

int a[maxn];
int go[maxn];
int nxt[maxn];
int dp[maxn][35];
int n,q;

void init(){
	for(int i=n;i>=1;--i){
		dp[i][0]=go[i];
		for(int j=1;j<32;++j){
			if(dp[i][j-1]>n) break;
			dp[i][j]=dp[dp[i][j-1]][j-1];
		}
	}
}

int solve(int l,int r){
	int ans=1;
	for(int i=30;i>=0;--i){
		if(dp[l][i]<=r) {
			l=dp[l][i];
			ans+=1<<i;
		}
	}
	return ans;
}

int main(){
	
	scanf("%d%d",&n,&q);
	memset(nxt,inf,sizeof(nxt));
	memset(dp,inf,sizeof(dp));
	for(int i=1;i<=n;++i) scanf("%d",a+i),go[i]=n+1;
	go[n+1]=n+1;
	
	for(int i=n;i>=1;--i){
		int tmp=a[i];
		for(int j=2;j*j<=tmp;++j){
			if(tmp%j) continue;
			go[i]=min(go[i],nxt[j]);
			while(tmp%j==0) tmp/=j;
			nxt[j]=i;
		}
		if(tmp>1){
			go[i]=min(go[i],nxt[tmp]);
			nxt[tmp]=i;
		}
		go[i]=min(go[i],go[i+1]);
	}

	init();
	
	int l,r;
	while(q--){
		scanf("%d%d",&l,&r);
		printf("%d\n",solve(l,r));
	}
	
	return 0;
}