这道题蛮简单的,但是也需要一点预备知识,两个数的最小公倍数和最大公约数的乘积等于这两个数的乘积.
证明·:x0 y0的最大公约数是A,最小公倍数是B,
x1=x0/A,y1=y0/A
那么x1 y1一定是互质的
所以x1 * y1 * A=B
展开一下就可以得到
A*B=x0 * y0
至于为什么x1y1A=B,可以去康一下最小公倍数的短除法的求法,一看就明白了.
知道这个定理以后题目就很简单了,满足条件的两个数的乘积就是x0y0, 现在就开始枚举,一个数x的因子我们只需要枚举1----sqrt(x),知道前面一半,后面一半也就知道了 就比如12 的一个因子是3 那么他的另一半因子就是4,所以我们枚举1—sqrt(x),假设1–sqrt(x)中一个数是t,那么另外一个数就是x0y0/t,所以我们要判断如果 x0y0%t == 0 满足gcd(x0y0/t,t)==x0, 那么可以判定x0*y0/t t这两个数是满足条件的数。
PS:有人可能会有疑问: x0y0%t == 0 gcd(x0y0/t,t)==x0 就可以判定满足条件,为啥判断条件不是满足最大公约数和最小公倍数?是这样的,我们之前已经说过了,如果两个数的最大公约数确定了,那么他的最小公倍数也就确定了…所以只需要判定满足最大公约数的条件就可以了.
用上面计算得出的结果*2就的到了百分之90正确答案,为啥是百分之90呢…这里还要注意,如果两个数的最大公约数和最小公倍数相同,那么这两个数一定是相同的,比如100 100,他的最大公约数是100,最小公倍数也是100,但是按照上面那个算法会得到结果是2,所以还要特别判断一下,如果两个数的最大公约数和最小公倍数相同,那么结果还要减去1
下面是AC代码~:
#include <iostream>
#include <cmath>
using namespace std;
int gcd(int a,int b); //求最大公约数
int main()
{
int n,m;
cin>>n>>m; //最大公约数 最小公倍数
int ans=0;
for(int i=1;i<=sqrt(n*m);i++)
{
if(n*m%i==0&&gcd(n*m/i,i)==n)
{
ans++;
}
}
ans*=2;
if(n==m)
{
ans--;
}
cout<<ans<<endl;
}
int gcd(int a,int b) //求最大公约数
{
if(a%b==0)
{
return b;
}
else
{
return gcd(b,a%b);
}
}