P1029 [NOIP2001 普及组] 最大公约数和最小公倍数问题题解

题目

输入两个正整数x0​,y0​，求出满足下列条件的P,Q的个数：

1. P,Q是正整数。

2. 要求P,Q以x0​为最大公约数，以y0​为最小公倍数。

试求：满足条件的所有可能的P,Q的个数。

输入输出格式

输入格式

一行两个正整数x0​,y0​。

输出格式

一行一个数，表示求出满足条件的P,Q的个数。

输入输出样例

输入样例

3 60

输出样例

4

解析

这里先补充一些最大公约数和最小公倍数的知识。设a1,a2两个整数，如果d|a1,d|a2，那么d就称为a1和a2的公约数，其中最大的称为a1和a2的最大公约数，记作(a1,a2)。一般地，可以类似地定义k个整数a1,a2,……,ak的公约数和最大公约数，后者记作(a1,……,ak).

设a1,a2是两个整数，如果a1|l,a2|l，那么l就称为a1和a2的公倍数，其中最小的称为a1和a2的最小公倍数，记作[a1,a2]。一般地，可以类似地定义k个整数a1,a2……,ak的公倍数和最小公倍数，后者记作[a1,a2……,ak].

下面给出四个有关最大公约数和最小公倍数的常见性质和结论。

性质1 对任意整数m，m(a1,a2,……,ak)=(ma1,ma2,……,mak)，即整数同时成倍放大，最大公约数也放大相同倍数。该性质同样适用于最小公倍数情况。

性质2 对任意整数x，(a1,a2)=(a1,a2+a1x)，即一个整数加上另外一个整数的任意倍数，它们的最大公约数不变。该性质不适用于最小公倍数情况。

性质3 (a1,a2,……,ak)=((a1,a2),……,ak)以及一个显然的推论(a1,a2,……,ak+r)=((a1,a2,……,ak),(ak+1,……ak+r))。

性质4 [a1,a2](a1,a2)=a1a2，最大公约数*最小公倍数=原来两个数的乘积。

根据性质2可以给出我们一个高效的求两数最大公约数的算法：每次让较大的数对较小的数取模可以缩小问题规模而保持最大公约数不变，然后重复这个步骤。递归边界使某数变成0.而此时另一个数即为所求答案。于是得到以下代码：

int gcd(int x,int y){
	if(y==0){
		return x;
	}
	else{
		return gcd(y,x%y);
	}
}

这种利用两数相除求最大公约数的方法叫做辗转相除法或Euclid算法。

很简单利用性质4，用两数之积除以它们的最大公约数，代码如下：

int lcm(int x,int y){
	return x/gcd(x,y)*y;
}

针对这个题目利用性质4可以得到PQ=xy，而又因为P|y，所以不妨枚举y的所有约数，假设为P，同时计算出Q，再检验这对P和Q是否满足条件。这里运用的是枚举约数的方法，时间复杂度为

#include<cstdio>
using namespace std;
int x,y,p,q,cnt;
int gcd(int x,int y){
	return y?gcd(y,x%y):x;
}
int main(){
	scanf("%d%d",&x,&y);
	for(int k=1;k<=y/k;k++){
		if(y%k==0){
			if(gcd(k,y/k*x)==x){//一组是P=k，Q=x*y/k
				cnt++;
			}
			if(y/k!=k){//注意重复判断
				if(gcd(y/k,k*x)==x){//另一组是P=y/k,Q=k*x
					cnt++;
				}
			}
		}
	}
	printf("%d",cnt);
	return 0;
}