Codeforces 1154 G. Minimum Possible LCM (贪心+简单数论)

题目：给定数组A，找到下标I和J使得A[i]和A[j]的最小公倍数最小，数据范围是2<=n<=1e6,A[i]<=1e7

思路：如果是按照正常做法求解GCD和LCM，那么是O(n^2)复杂度，题目要求是NlgN内求解，所以转换思路每次枚举GCD按照当前的GCD去选出满足条件的答案。需要注意的是对于同一个GCD：d而言，两个数均为d的倍数，那么这两个数最小即可，没有必要继续枚举下面的数，否则结果也会更大：2,4与4,8都是2的倍数，但是明显第一个是满足条件的。

AC code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e7+5;

ll n,cnt[maxn],a[maxn];
vector<ll>res;
ll ans=1e17+5;

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

int main()
{
    ll limit=0,value1=0,value2=0;
    scanf("%lld",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]);
        cnt[a[i]]++;
        limit=max(limit,a[i]);
    }
    for(int i=1;i<=limit;i++)
    {
        int j;res.clear();
        for (j=i;j<=limit;j+=i)
        {
            if (!cnt[j]) continue;
            int tmp=cnt[j];
            while ((tmp) && (res.size()<2))
            {
                res.push_back(j);
                tmp--;
            }
            if (res.size()==2) break;
        }
        if (res.size()!=2) continue;
        ll now=res[0]*res[1]/gcd(res[0],res[1]);
        if (now<ans)
            ans=now,value1=res[0],value2=res[1];
    }
    //cout<<value1<<" "<<value2<<endl;
    ll index1=0,index2=0;
    for(int i=1;i<=n;i++)
    {
        if(index1==0&&a[i]==value1)
            index1=i;
        else if(index2==0&&a[i]==value2)
            index2=i;
    }
    if(index1>index2)
        swap(index1,index2);
    printf("%lld %lld\n",index1,index2);
    return 0;
}