【CERC2013】bzoj4052 Magical GCD

Description 给出一个长度在 100 000 以内的正整数序列，大小不超过 10^12。
求一个连续子序列，使得在所有的连续子序列中，它们的GCD值乘以它们的长度最大。

对于固定的右端点r，区间(l..r)的gcd单调不增，且每次变小一定变为原来的因数，所以“有用的位置”【gcd相同的取最左端】不超过logai个。每次r右移，可能会使一些原来gcd不同的位置变得相同，所以需要O(logai)维护gcd，删除元素并更新最优解。总时间复杂度O(nlogai)。

#include
#include
#include
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);
    }
}