Description 给出一个长度在 100 000 以内的正整数序列,大小不超过 10^12。
求一个连续子序列,使得在所有的连续子序列中,它们的GCD值乘以它们的长度最大。
对于固定的右端点r,区间(l..r)的gcd单调不增,且每次变小一定变为原来的因数,所以“有用的位置”【gcd相同的取最左端】不超过logai个。每次r右移,可能会使一些原来gcd不同的位置变得相同,所以需要O(logai)维护gcd,删除元素并更新最优解。总时间复杂度O(nlogai)。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define LL long long
int f[2][100010],top[2];
LL a[100010],g[2][100010];
LL gcd(LL x,LL y)
{
return y?gcd(y,x%y):x;
}
int main()
{
int i,j,k,m,n,p,q,T;
LL ans,x,y,z;
scanf("%d",&T);
while (T--)
{
scanf("%d",&n);
for (i=1;i<=n;i++)
scanf("%lld",&a[i]);
top[0]=top[1]=ans=0;
for (i=1;i<=n;i++)
{
top[i&1]=0;
for (j=1;j<=top[i&1^1];j++)
{
x=gcd(g[i&1^1][j],a[i]);
if (top[i&1]&&g[i&1][top[i&1]]==x)
continue;
top[i&1]++;
f[i&1][top[i&1]]=f[i&1^1][j];
g[i&1][top[i&1]]=x;
ans=max(ans,x*(i-f[i&1][top[i&1]]+1));
}
if (a[i]==g[i&1][top[i&1]])
continue;
top[i&1]++;
f[i&1][top[i&1]]=i;
g[i&1][top[i&1]]=a[i];
ans=max(ans,a[i]);
}
printf("%lld\n",ans);
}
}