问题描述
已知正整数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;
}