Array Stabilization (GCD version)

题目链接
题意：给我们一个有n个元素的数组a成环形，然后b数组为 b[i] = gcd(a[i], a[i+1]),且b[n] = gcd(a[n],a[1])。然后再把b数组赋给a数组，问这样几次以后可以使a数组当中所有的元素相等。
题解：易证gcd(gcd(a,b),c) == gcd(a,b,c)。有了这个结论后我们就好写题了。我们可以知道要最后元素相等的话，那么肯定是a[i] = gcd(a[1],a[2]…a[n])（1<=i<=n)。且第一次为a[i] = gcd(a[i],a[i+1]),第二次为，a[i] =gcd(gcd(a[i],a[i+1)),gcd(a[i+1],a[i+2])) = gcd(a[i],a[i+1],a[i+2])。
但我们有了这个性质以后暴力这个题还是会超时的，我们要明白加入这个数在第k次后为 gcd(a[1],a[2]…a[n])，那么他再进行操作是不会改变的，由此我们可以得到单调性，我们可以用二分来进行优化，我们肯定也不能每次暴力进行gcd，考虑到这是个区间问题，所以我们选用比较简单的st表来存储区间gcd，这个题就可以完成了。

#include<bits/stdc++.h>
#define ll long long
#define pr pair<ll,ll>
#define ios ios::sync_with_stdio(false)
#define CRL(a) memset(a,0,sizeof a)
#define endl "\n"

using namespace std;
const int maxn = 4e5 + 5;
int st[maxn][32];
int a[maxn];
int query(int l,int r)
{
    int k = log2(r - l + 1);
    return __gcd(st[l][k], st[r - (1 << k) + 1][k]);
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n;
        cin >> n;
        for (int i = 1; i <= n;i++)
        {
            cin >> a[i];
            st[i][0] = a[i + n] = st[i + n][0] = a[i];
        }
        int GCD = a[1];
        bool flag = 1;
        for (int i = 2; i <= n;i++)
        {
            GCD = __gcd(GCD, a[i]);
            if(a[i]!=a[1])
                flag = 0;
        }
        if(flag)
        {
            cout << 0 << endl;
            continue;
        }
        for (int j = 1; (1 << j) <= n + n;j++)
        {
            for (int i = 1; i + (1 << j) - 1 <= n + n;i++)
            {
                st[i][j] = __gcd(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]);
            }
        }
        int ans = 0;
        for (int i = 1; i <= n;i++)
        {
            int l = i, r = n + i - 1;
            while(l<=r)
            {
                int mid = l + r >> 1;
                if(query(i,mid)!=GCD)
                {
                    l = mid + 1;
                }
                else 
                {
                    r = mid - 1;
                }
            }
            ans = max(ans, l - i);
        }
        cout << ans << endl;
    }
}