“求公约数”和“求公倍数”之类问题的“逆问题”

题目：

Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。现在，刚刚放学回家的Hankson正在思考一个有趣的问题。

今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在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

源代码：

#include<stdio.h>

int divisor (int a,int b)    /*自定义函数求两数的最大公约数*/

{

  int  temp;          /*定义整型变量*/

  if(a<b)             /*通过比较求出两个数中的最大值和最小值*/

    { temp=a;a=b;b=temp;} /*设置中间变量进行两数交换*/

   while(b!=0)           /*通过循环求两数的余数，直到余数为0*/

    {

      temp=a%b;

      a=b;              /*变量数值交换*/

      b=temp;

    }

  return (a);            /*返回最大公约数到调用函数处*/ 

}

int multiple (int a,int b)  /*自定义函数求两数的最小公倍数*/

{

  int divisor (int a,int b); /*自定义函数返回值类型*/

  int temp;

  temp=divisor(a,b);  /*再次调用自定义函数，求出最大公约数*/

  return  (a*b/temp); /*返回最小公倍数到主调函数处进行输出*/

}

 

void main()

{

int a[100][4];

int i,n,j,m,t=0,max;

printf("输入数据的组数：");

scanf("%d",&n);

for(j=0;j<n;j++)

{

      for(i=0;i<4;i++)

  {

          scanf("%d",&a[j][i]);

  }

}//输入n组数据

 

for(j=0;j<n;j++)

{

max=a[j][3]/a[j][2];

t=0;

for(m=max;m<=a[j][3];m++)//题目中x也就是程序中的m的范围是最小公倍数除以另一个因子到最小公倍数这个值

{

 

if(a[j][1]==divisor (a[j][0],m)&&a[j][3]==multiple (a[j][2],m) )

{

            t=t+1;

}

}//如果符合要求，就给t加一，算出的t就是满足了要求的数字个数

printf("%d\n",t) ;

}

}

        或者可以将外部的求最大公约数和最小公倍数的函数删除，而直接利用它们的数学关系可以做出来，这个程序较多的使用了for循环，也有一部分多位数组的知识，这部分可以很好地复习到之前学到的知识，再有一点就是在做这个程序时，x的取值范围必须自己好好斟酌一下，取值范围的选取直接影响了最终数值的正确与否。