Hankson的“逆问题”

2. 问题

已知正整数a0,a1,b0,b1，设某未知正整数x满足：
1、 x和a0的最大公约数是a1；
2、 x和b0的最小公倍数是b1。
Hankson的“逆问题”就是求出满足条件的正整数x。这样的x并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的x的个数。

3.程序代码

#include<stdio.h>
int gcd(int a, int b)	//最大公约数函数
{
    return b == 0 ? a : gcd(b, a % b);
}

int main () {
    int n, sum=0;//n为输入数的组数，sum计数
	int a0, a1, b0, b1;
	printf("请输入准备数输入的组数n： ");
    scanf("%d",&n);
	while(n<=0)
	{
	printf("error,请输入准备数输入的组数n： ");
	scanf("%d",&n);
	}
    while(n--) 
	{
	
	printf("请输入4个数：  ");
	scanf("%d%d%d%d",&a0,&a1, &b0, &b1);    
  while((a0<=0 )||(a1<=0)||(b0<=0)||(b1<=0))
	{	printf("请重新输入4个数：  ");
		scanf("%d%d%d%d",&a0,&a1, &b0, &b1); 
	}
		for (int i = 1; i * i <= b1; i++) //i从1开始，i*i不超过b1,枚举
		{	
		
		//printf("%d\t",i);
		if (b1 % i == 0)
		{         
        if (gcd(a0, i) == a1 && (b0*i / gcd(b0, i)) == b1 ) //公约数为a1，公倍数为b1
		{sum++; 
		//printf(":%d",i);
		}											//满足条件计数一次
		else if (gcd(a0, (b1 / i)) == a1  &&  (b0*(b1 / i) / gcd(b0, (b1/i)) == b1))//(b1/i)
		{sum++;
		//printf(":%d",i);
		}											//满足条件计数一次	 
		}

		}
        printf("满足条件的数个数为：    %d\n\n",sum);//打印结果
		sum=0;
	 }
    return 0;
}

4.运行结果：