（C++）Hankson的“逆问题”-逆解最大公约数

问题描述

已知正整数a0,a1,b0,b1，设某未知正整数x满足：

1、  x和a0的最大公约数是a1；

2、  x和b0的最小公倍数是b1。

Hankson的“逆问题”就是求出满足条件的正整数x。但稍这样的x并不唯一，甚至可能不存在，编程求解满足条件的x的个数。

输入格式

输入第一行为一个正整数n，表示有n组输入数据。接下来的n行每行一组输入数据，为四个正整数a0，a1，b0，b1，每两个整数之间用一个空格隔开。输入数据保证a0能被a1整除，b1能被b0整除。

输出格式

输出共n行。每组输入数据的输出结果占一行，为一个整数。

对于每组数据：若不存在这样的x，请输出0；

若存在这样的x，请输出满足条件的x的个数；

样例输入

2

41 1 96 288

95 1 37 1776

样例输出

6

2

问题分析

根据已知条件很容易列出两个等式：

     gcd(x,a0)=a1......①

     lcm(x,b0)=b1......②

(注：gcd:greatest common divisor 最大公约数   lcm:lowest common multiple最小公倍数  )

我们又知道，最大公约数和最小公倍数之间有一个隐含的等式,即对于两个数a,b，有  :

     lcm(a,b)=(a*b)/gcd(a,b)......③

因此②可以转换为:

     (x*b0)/gcd(a,b0)=b1......④

网上有很多巧妙的的解法，但是本人小白一枚看不大懂......所以我还是采用了暴力枚举的方法。经同学指点，我知道了因为x和b0的最小公倍数一定比大，所以可以定义一个计数器，将x从1到b0枚举，满足等式①和②的计数器就加加一次。

伪代码描述

while(n>=0)

{

     for(x=1;x<=b1;x++)

      {

            if(gcd(x,a0)==a1 && (x*b0)/gcd(x,b0) ==b1)

                     cnt++;

        }

   输出i;

    n--;

}

运行结果

源代码

#include<iostream>
#include<math.h>
using namespace std;

//求最大公约数-辗转相除递归调用 
int gcd(int a,int b)
{
	return !(a%b)? b : gcd(b,a%b);
}

int SolveHankson(int &n)//求满足条件的解
{
	int a0,a1,b0,b1;
	int i=0;//i是计数器，统计满足条件的x个数 
	while(n--)
	{
		cin>>a0>>a1>>b0>>b1;
		for(int x=1;x<=b1;x++)
		{
			if(gcd(x,a0)==a1 && (x*b0)/gcd(x,b0) ==b1)
			{
				i++;
			}	
		}
		cout<<i<<endl;
		i=0;
	}
}
int main(void)
{
	int n;//n组输入数据
	cin>>n;
	SolveHankson(n);
	return 0;
}