CF1547F Array Stabilization (GCD version)（ST算法+数论）

题目传送门

思路

$ST$算法

$ST$算法能在$O(nlogn)$时间处理以后（动态规划思想），然后$O(1)$时间来查询的一种及其高效的做法。
最大公约数的一个性质：
$gcd(i,j,k)=gcd(gcd(i,j),k)$
因此设$st[i][j]$表示数列中在下标子区间$[i,i+2^j-1]$内的最大公约数，也就是从$i$开始的$2^j$个数的最大公约数，因此我们需要预处理出来所有的这样的区间最大公约数。
那么我们怎么样去实现这个东西呢？就是靠递推（动态规划的思想）
在递推时，我们将这个母区间分成两个子区间，进行递推：
$st[i][j]=gcd(st[i][j-1],st[i+2^{j-1}][j-1])$
这样我们就可以把长度为$2^j$的子区间的最大公约数转换成为两个半长度为$2^{j-1}$的左右子区间的最大公约数。

• 进行预处理的算法模板：

for(int i=1;(1<<i)<=n;i++)
    	for(int j=1;j+(1<<i)-1<=n;j++)
     	   st[j][i]=gcd(st[j][i-1],st[j+(1<<(i-1))][i-1]);

• 进行查询的算法模板

int check(int i,int mid){
	int k=log2(mid-i+1);
	return gcd(st[i][k],st[mid-(1<<k)+1][k]);
}

本题分析

这道题一看这是一个环，环我们不是很好处理，因此可以“断环成链”，开一个二倍大小的数组，来存这个链。
假设我们现在的数组是
$a_1,a_2,a_3,a_4......a_n$
我们处理一次了以后就变成了
$gcd(a_1,a_2),gcd(a_2,a_3),gcd(a_3,a_4)......gcd(a_{n-1},a_n),gcd(a_n,a_1)$
处理了第二次
$gcd(gcd(a_1,a_2),gcd(a_2,a_3)),gcd(gcd(a_2,a_3),gcd(a_3,a_4))......gcd(gcd(a_{n-1},a_n),gcd(a_n,a_1))$

现在我们发现了第二次的一些性质：
$gcd(a_1,a_2,a_3),gcd(a_2,a_3,a_4)......gcd(a_{n-1},a_n,a_1)$
因此只要我们找到每一个$i$,都能使$gcd(a_i,a_{i+1}......a_{i+k-1})$的最小的$k$,也就是这一个$k$能满足每一个$i$.

#include <iostream>
#include<cmath>
using namespace std;

int gcd(int a,int b){
	return b?gcd(b,a%b):a;
}
int a[2102100];
int st[2102100][20];
int _gcd;int n;
int ans=-2;
int check(int i,int mid){
	int k=log2(mid-i+1);
	return gcd(st[i][k],st[mid-(1<<k)+1][k]);
}
void bsearch(int l,int r,int i){
	int mid;
	while(l<r){
		mid=l+r>>1;
		if(check(i,mid)!=_gcd)l=mid+1;
		else r=mid;
		//cout<<mid<<endl;
	}
	
	ans=max(ans,l-i+1);
	//cout<<ans<<endl;
}
int main()
{
	int T;
	cin>>T;
	while(T--){
		ans=-1;
		cin>>n;
	int ju=0;
	for(int i=1;i<=n;i++){
		cin>>a[i];
		a[i+n]=a[i];
		if(a[i]!=a[1]){
			ju=1;
		}
		st[i+n][0]=st[i][0]=a[i];
		if(i==1)_gcd=a[i];
		else _gcd=gcd(_gcd,a[i]);
	}
	if(!ju){
		cout<<"0"<<endl;
		continue;
	}
	for(int i=1;1<<i<=2*n;i++)
    	for(int j=1;j+(1<<i)-1<=2*n;j++)
     	   st[j][i]=gcd(st[j][i-1],st[j+(1<<(i-1))][i-1]);
    for(int i=1;i<=n;i++){
    	bsearch(i,i+n-1,i);
    }
    cout<<ans-1<<endl;
}
}