题目大意
定义一个区间是友好的,当m>=2的时候区间l,r之间%mod值相等。给出一个数组,求最大友好区间长度。
思路
//gcd的性质
gcd(a,b) = gcd(b,a-b)
gcd(a,b,c) = gcd(b,b-a,c-b)
由性质可以知道,区间gcd就是差分数组的gcd,再使用线段树去维护就可以了。区间gcd的维护方法有两种,一种是st表,一种是线段树,我使用的线段树。
知道区间gcd之后,题目所求就变成了区间gcd(l,r)>=2的条件下,l,r差的最大值。
由于如果右边gcd为1了,左端点到右端点的gcd肯定也是1,所以可以使用双指针维护来求最大区间。给定右端点,寻找最远的左端点,每次更新答案。
AC代码
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define ll long long
#define endl "\n"
const int max_n = 2e5 + 100;
int buf[max_n];
int arr[max_n];
struct Node
{
int l, r;
int v;
} tr[max_n * 4];
int gcd(int a, int b)
{
if (!b)
return a;
return gcd(b, a % b);
}
void pushup(int u)
{
tr[u].v = gcd(tr[u << 1].v, tr[u << 1 | 1].v);
}
void build(int u, int l, int r)
{
tr[u] = {l, r};
if (l == r)
{
tr[u].v = arr[r];
return;
}
int mid = l + r >> 1;
build(u << 1, l, mid);
build(u << 1 | 1, mid + 1, r);
pushup(u);
}
int query(int u, int l, int r)
{
if (tr[u].l >= l && tr[u].r <= r)
{
return tr[u].v;
}
int mid = tr[u].l + tr[u].r >> 1;
if (l > mid)
return query(u << 1 | 1, l, r); //全在右边
if (r <= mid)
return query(u << 1, l, r); //全在左边
//两边都有
return gcd(query(u << 1, l, r), query(u << 1 | 1, l, r));
}
signed main()
{
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
int T;
cin >> T;
while (T--)
{
int n;
cin >> n;
for (int i = 1; i <= n; i++)
cin >> buf[i];
for (int i = 1; i <= n - 1; i++)
{
arr[i] = abs(buf[i + 1] - buf[i]);
}
if (n == 1)
{
cout << 1 << endl;
continue;
}
build(1, 1, n - 1);
int ans = 0;
int l = 1, r = 1;
while (r < n)
{
while (l <= r && query(1, l, r) == 1)
l++;
ans = max(ans, r - l + 2);
r++;
}
cout << ans << endl;
}
return 0;
}