NOI Oline 2022 入门组第二题 数学游戏

作者：hsez_yyh
链接：https://blog.csdn.net/yyh_getAC/article/details/123905285
来源：湖北省黄石二中竞赛组
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

今年NOI Oline真的是折磨人，差点给我劝退了……

入门组的第二题还是没写出来，只能打了个暴力，混了大概50分。比赛完了好好想了想，发现竟如此简单——真的好想去跳楼刷题。

下面我们来好好探讨下这道看似很难却极其简单的题目。

首先，我们先看看数据规模：t=500000次询问，如此，我们可以肯定，每次查询操作的时间复杂度应该不会超过log(t)。然而，这道恶心好题让人一看就没思路，各种数据结构和算法都感觉用不上——那就不用了，这题一看就要推公式啊。

我们先假设x,y的最大公约数为c，那么x，y可以写成x=ac，y=bc的形式，由于z=x*y*gcd(x，y)；那么z=a*b*(c^3);c=gcd(ac，bc)； 然后，我们就可以得到 c^2=gcd( a*(c^2),b*(c^2) )。此时有一个神奇的概念，因为c为x，y的最大公约数，那么a、b这两个数肯定互质，也就是说a、b分解质因数后，比如a=a1*a2*...*an ; b=b1*b2*...*bn ; a1到an中的任意一个数都不等于b1到bn中任意一个数。那么，我们可以肯定 a^2 绝对与 b 互质(有兴趣自己证明)。那么c^2=gcd( (a^2)*(c^2),b*(c^2) )。此时，我们好好观察一下前面几个等式，可以看到c^2=gcd( x^2，(z/x) )；而我们已知x和z，因此

我们就能很快的得到c=sqrt(gcd( x^2，(z/x) ) )；所以每次询问我们仅只用求出gcd而已。然后每次只用判断一下是否符合题目条件即可，复杂度完全可以接受。

下面是ac代码：

#include<bits/stdc++.h>
using namespace std;
int T;
long long xx,yy,zz;
long long get_gcd(long long a,long long b)
{
	if(a%b==0) return b;
	else return get_gcd(b,a%b);
}
int main()
{
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lld %lld",&xx,&zz);
		if(zz%xx!=0)
		{
			printf("-1\n");
			continue ;
		}
		long long vv=zz/xx;
		long long t=get_gcd(xx*xx,vv);
		long long d=sqrt(t);
		if(d*d!=t)
		{
			printf("-1\n");
			continue ;
		}
		if(zz%(xx*d)!=0)
		{
			printf("-1\n");
			continue ;
		}
		yy=zz/(xx*d);
		printf("%lld\n",yy);
	}
	return 0;
}