bzoj 4488: [Jsoi2015]最大公约数

做这个题的时候蛋疼的想一个知道[l,r]的gcd那[l+1,r]的能不能快速算出来。

然后想了一下发现丝毫没有道理。

于是又扒了题解。

神奇的题解说，一个序列的gcd最多有log个，突然就GG了

诶，遇到这样的题还是多搞几组数据观察一下规律的好

#include <bits/stdc++.h>
#define LL long long
using namespace std;
inline LL ra()
{
    LL x=0; char ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x;
}

const int maxn=100005;

struct node{
    LL num; int x;
    bool operator < (const node &y) const {return num<y.num||num==y.num && x<y.x;}
}f[2][maxn];
LL a[maxn];
int n;

LL gcd(LL x, LL y) {return y==0?x:gcd(y,x%y);}

int main()
{
    n=ra(); 
    for (int i=1; i<=n; i++) a[i]=ra();
    int s1=0,s2=0;
    int p=0,q=1;
    LL ans=0; node d;
    for (int i=1; i<=n; i++)
    {
        for (int j=1; j<=s1; j++) f[p][j].num=gcd(a[i],f[p][j].num);
        d.num=a[i]; d.x=i; s1++; f[p][s1]=d;
        sort(f[p]+1,f[p]+s1+1);
        s2=0;
        for (int j=1; j<=s1; j++)
            if (f[p][j].num!=f[p][j-1].num)
                f[q][++s2]=f[p][j];
        for (int j=1; j<=s2; j++) ans=max(ans,f[q][j].num*(LL)(i-f[q][j].x+1));
        p^=1; q^=1; s1=s2;
    }
    printf("%lld\n",ans);
    return 0;
}