cf 1154G Minimum Possible LCM

...这题关键在他的a[i]<=1e7

那么我们知道lcm(a,b)=a*b/gcd(a,b)；

那么我们只要枚举每一个因数d，不管他是不是gcd

然后找出能被这个d整除的最小的两个数字a,b

那么对于这个因数d，tmp=a*b/d,ans=min(ans,tmp)

由于我们枚举了1-1e7所有的质因子，所以就算a*b/d不是lcm，但之后总会枚举到a*b/gcd(a,b)

而使用桶来装数字并查找d的倍数的的最小两个，复杂度为 1e7/d

那么总复杂度就是调和级数和  nlnn ，最后我程序跑出来是1.4s, 而时限是4s

 

#include<bits/stdc++.h>
#define maxl 10000010
using namespace std;

int n,l,r;
int a[maxl];
int num[maxl][2],cnt[maxl];
long long ans;

inline void prework()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		++cnt[a[i]];
		if(cnt[a[i]]==1)
			num[a[i]][0]=i;
		if(cnt[a[i]]==2)
			num[a[i]][1]=i;
	}
}

inline long long gcd(long long a,long long b)
{
	return b==0 ? a:gcd(b,a%b); 
}

inline void mainwork()
{
	ans=1ll*a[1]/gcd(a[1],a[2])*a[2];
	l=1;r=2;int sum,ll,rr;long long tmp;
	for(int d=1;d<maxl;d++)
	{
		sum=0;tmp=1;
		for(int j=d;j<maxl;j+=d)
		{
			if(cnt[j]>=2)
			{
				if(sum==1)
					sum++,rr=num[j][0];
				if(sum==0)
					sum+=2,ll=num[j][0],rr=num[j][1];
			}
			else if(cnt[j]==1)
			{
				sum++;
				if(sum==1) ll=num[j][0];
				if(sum==2) rr=num[j][0];
			}
			if(sum==2)
			{
				tmp=1ll*a[ll]*a[rr]/d;
				if(tmp<ans)
				{
					ans=tmp;
					l=ll;r=rr;
				}
				break;
			}
		}
	}
}

inline void print()
{
	if(l>r) swap(l,r);
	printf("%d %d",l,r);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}