题目链接:
https://codeforces.com/contest/1549/problem/D
tle了无数次,甚至一开始上st表还没原始暴力过的样例多🔨,真的调st表及其查询调到心态崩了。(第一次用st表。啊不,主要是我太菜了😭)
题目大意:
求一个正整数序列中最长模某数(任意)同余序列(连续)长度。
解题思路:
1.a mod m = b mod m = c mod m,则有abs(a-b)%m==0&&(c-b)%m==0成立。因此模m同余可转化为差可被m整除。
2.st表的基本知识: st表适用于求解区间最值问题。st[i][j]表示从i位置到的区间最值,可由st[i][j-1]和st[
][j-1]生成。(可以简单的用数学方法进行证明)。也就是说st[i][0]就存放原始的第i个数。st表的创建时间复杂度为
,单次查询的时间复杂度为O(1)。
3.本题看作是在差数组中求解最大公约数>1的最长区间,可以使用st表存放最大公约数进行求解。
程序代码:
#include <bits/stdc++.h>
using namespace std;
long long int a[200005], c[200005]; //c存放差值,下标从0开始
long long int f[200005][30]; //st表
int t, n, len[200005];
long long int ans;
//最大公约数求解
long long int gcd(long long int a, long long int b)
{
long long int temp;
while(b!=0)
{
temp=a%b;
a=b;
b=temp;
}
return a;
}
//创建st表
void init()
{
len[1]=0;
for(int i=2; i<=n; i++) len[i]=len[i/2]+1;
int m=len[n]+1;
for(int i=1; i<n; i++) f[i][0]=c[i];
for(int j=1; j<m; j++)
for(int i=1; i+(1<<j)-1<n; i++)
f[i][j]=gcd(f[i][j-1], f[i+(1<<(j-1))][j-1]);
}
//单次区间查询
long long int query(int l, int r)
{
int k=len[r-l+1];
return gcd(f[l][k], f[r-(1<<k)+1][k]);
}
int main()
{
cin>>t;
while(t--)
{
ans=0;
cin>>n;
cin>>a[0];
for(int i=1; i<n; i++)
{
cin>>a[i];
c[i]=abs(a[i]-a[i-1]);
}
if(n==1)
{
cout<<1<<endl;
continue;
}
init();
for(int i=1; i+ans<n; i++)
while(i+ans<n&&query(i, i+ans)>1) ans++;
cout<<ans+1<<endl;
}
return 0;
}