UVA1642:魔法GCD

UVA1642:魔法GCD

質問の意味:

長さが与えられた\(N- \ 1当量5 ^ 10 \) 各番号\(a_iを\ leq10 12である^ {} \) 連続したシーケンスを探しているように、最大のサブ規則の長さの産物。

\(T \)データのセット。

アイデア:

間隔最大公約数テンプレートのタイトル。

列挙\((i、j)は\ ) 暴力、その後、時間の複雑さがある(O)2logn ^(N- \)\、確かにタイムアウトになります。

所定のシーケンス\(A \) 連続サブセグメント\(GCD \)有していて\(ログ(最大\ {a_iを \})\) 可能。

\(GCD(1、...、I)= GCD(GCD(1、...、I-1)、(I))\)

毎回だから、固定小数点の右、異なる見つけるために左に\(GCD \)の値を。

#include<bits/stdc++.h>
#define PLI pair<long long, int>
using namespace std;
typedef long long ll;
const int maxn = 1e5+10;
ll a[maxn];
int n;

ll gcd(ll a, ll b)
{
    if(b == 0) return a;
    return gcd(b, a%b);
}

//fir->gcd sec->右端点索引
vector<PLI> g[maxn];

void solve()
{
    scanf("%d", &n);
    for(int i = 1; i <= n; i++)
    {
        scanf("%lld", &a[i]);
        g[i].clear();
    }
    for(int i = 1; i <= n; i++)
    {
        ll x = a[i], y = i;
        g[i].push_back({x,y});
        for(int j = 0; j < g[i-1].size(); j++)
        {
            PLI p = g[i-1][j];
            ll t = gcd(x, p.first);
            if(t != x)
            {
                x = t; y = p.second;
                g[i].push_back({x, y});
            }
        }
    }

    ll ans = 0;
    for(int i = 1; i <= n; i++)
    {
        PLI p1, p2;
        for(int j = 0; j < g[i].size()-1; j++)
        {
            p1 = g[i][j];
            p2 = g[i][j+1];
            ans = max(ans, p1.first*(i-p2.second));
        }
        p1 = g[i][g[i].size()-1];
        ans = max(ans, (ll)(i)*p1.first);
    }

    cout << ans << endl;
}

int main()
{
    int T; scanf("%d", &T);
    while(T--) solve();
    return 0;
}

おすすめ

転載: www.cnblogs.com/zxytxdy/p/12306167.html