CF1549D Integers Have Friends(st表+二分/双指针)

题目传送门

题意

这道题让你找一个最长的区间，这个区间里面的每一个数，都要对于$m>=2$同余

思路

这道题用到一个同余的性质
如果 $x\equiv y(modp)$
那么$(x-y)$%$p=0$
那么对于一整个区间来说，我们只要满足这个区间的差分数组不互质就可以。
因此问题就转换成为了，维护一个差分数组，找到最长的不互质的子序列。
因此我们就要用到二分+st表，如果想优化你也可以双指针（100ms的优化）
同时这道题用scanf比cin快一倍。

二分

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
ll a[212100];
ll st[212100][30];
ll lg2[212100];ll n;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll query(ll l,ll r){return gcd(st[l][lg2[r-l+1]],st[r-(1<<lg2[r-l+1])+1][lg2[r-l+1]]);}
ll bsearch(ll l,ll r,ll p){

	ll mid;
	ll k=1;
	while(l<r){
		mid=l+r+1>>1;
		if(abs(query(p,mid))<2){
				r=mid-1;
		}else{
			l=mid;
		} 
	}
	if(abs(query(p,l))<2)return 1;
	return l;
}
int main(){
	int T;
	cin>>T;
	lg2[0]=-1;
    for(int i=2;i<=210000;i++)lg2[i]=lg2[i>>1]+1;
	while(T--){
		cin>>n;
		for(int i=1;i<=n;i++){
			cin>>a[i];
		}
		for(int i=n+1;i<=2*n&&i<=200000;i++)st[i][0]=0;
		for(int i=1;i<=n;i++){
			st[i][0]=a[i]-a[i-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]);
			}
		}
		ll p=1;
		for(ll i=1;i<=n;i++){
			ll u=bsearch(i+1,n,i+1);
			if(u==1)continue;
			p=max(p,u-i+1);
		}
		cout<<p<<endl; 
	}
}

双指针

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
ll a[212100];
ll st[212100][30];
ll lg2[212100];ll n;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll query(ll l,ll r){return gcd(st[l][lg2[r-l+1]],st[r-(1<<lg2[r-l+1])+1][lg2[r-l+1]]);}
int main(){
	int T;
	cin>>T;
	lg2[0]=-1;
    for(int i=2;i<=210000;i++)lg2[i]=lg2[i>>1]+1;
	while(T--){
		scanf("%lld",&n);
		//cin>>n;
		for(int i=1;i<=n;i++){
			scanf("%lld",&a[i]);
		//	cin>>a[i];
		}
		for(int i=n+1;i<=2*n&&i<=200000;i++)st[i][0]=0;
		for(int i=1;i<=n;i++){
			st[i][0]=a[i]-a[i-1];
			if(st[i][0]<0)st[i][0]=-st[i][0];
		}
		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]);
			}
		}
		ll p=1;
		ll ans=1;
		for(ll i=1,j=2;i<=n;i++){
			if(i==j)j=i+1;
			while(i<j&&j<=n){
				if(query(i+1,j)>=2){
					ans=max(ans,j-i+1);
					j++;
				}else{	
					break;
				}
			}
			
		}
		cout<<ans<<endl; 
	}
}