实现N个数的最小公倍数和最大公约数

题目要求

• 基本要求： 求N个数的最大公约数和最小公倍数。用C或C++或java或python语言实现程序解决问题。
• 提高要求：Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。现在，刚刚放学回家的Hankson正在思考一个有趣的问题。
  今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在Hankson认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1，设某未知正整数x满足：
• x和a0的最大公约数是a1；
• 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

设计思路

• 主函数上采用以往实践较多的菜单功能来实现基本要求和提高要求。
• 本次程序设计并没有用到多高级的数据结构，但在保存每行输出结果的时候，采用一维数组num[N]保存（其中N表示行数）。另外在求N个数的最小公倍数和最大公约数，依然采用数组arr[N]保存这N个数（其中这里的N表示输入的N的个数）。
• 求最小公约数和最大公约数上采用辗转相除法的递归套用，尽管执行效率低于其他算法，但语言简洁。其中求最大公约数也就是先让最后两个数先求最大公约数，其结果在与之前倒数第三个数及需求最大公约数，以此类推，当进行到第一个数时即为这N个数的最大公约数。求最小公倍数，也是一样的道理，不过得利用到公式a*b/ gcd(a, b)先求出最大公约数再求出最小公倍数。
• 实现提高要求上，其实是根据数学知识对问题进行简化而得来的，题目强制规定了
  0整除a1，b0整除b1这两个条件。言外之意就是x需满足a1整除x且x整除b1。以最小公倍数x满足条件是[1,b1]，满足这个条件后，再判断x和a0最大公约是否为a1，是的话计数器flag变量加1，保存到数组num[N]中，最后输出结果。

流程图

源码展示

package com.divisor_and_multiciple;
/**
 * 该类封装和求最大公约数有关的函数
 * 最大公约数:先让最后两个数先求最最大公约数，其结果在与之前倒数第三个数及需求最大公约数，以此类推，
 * 求带第一个数便返回结果;
 * @author 万物甦醒
 *
 */
public class Divisor {
	//求两个数的最大公约数
	int gcd(int a,int b)
	{
		if (a<b) //保证a大于b
			{
				int temp=a;
				a=b;
				b=temp;
			}
		if (b==0)
			return a;
		else
			return gcd(b,a%b);
	}
	//求n个数的最小公倍数
	int ngcd(int[]arr,int n)
	{
		if(n==1)
			return arr[0];
		else
		return gcd(arr[n-1],ngcd(arr, n-1));//坐标n-1数与前面数的公约数 
	}
}

package com.divisor_and_multiciple;
/**
 * 该类封装和最小公倍数有关的几个函数
 * 思路原理和最大公约数大体相似
 * 其中求两个数的最小公倍数定理公式：a*b/gcd(a, b);
 * 因此，须先在该类中声明Divisor对象
 * @author 万物甦醒
 *
 */
public class Multiciple 
{
	//定义公约数类对象，为了后面调用求最大公约数
	Divisor divisor=new Divisor();
	//求两个数的最小公倍数
	int lcm(int a,int b)
	{
		return a*b/divisor.gcd(a, b);
	}
	
	//n个数的最小公倍数
	int nlcm(int []arr, int n)
	{
		if (n==1)
			return arr[0];
		else
			return lcm(arr[n-1],nlcm(arr,n-1));
	}
}

```java
package com.divisor_and_multiciple;

import java.util.Scanner;
/**
 * 求取N个数的最大公约数，和最小公倍数
 * 思路： 1.让用户输入N的大小，依次输入这N个数；
 *	    2.求最大公约数和最小公倍数采用辗转相除法gcd(a,b)。
 *		   最大公约数:先让最后两个数先求最最大公约数，其结果在与之前倒数第三个数及需求最大公约数，以此类推，
 *				    求带第一个数便返回结果;
 *		   最小公倍数和上面同理。
 * @author 万物甦醒
 *
 */
public class BasicHomework {
	void homeWork1()
	{
		int N; 
		System.out.print("请依次输入N的大小");
		Scanner scan=new Scanner(System.in);
		N=scan.nextInt(); 		//输入N的大小
		int []arr=new int[N];	//定义一个数组，存放这N个数

		//依次输入N个数
		System.out.print("请依次输入这"+N+"个数：");
		for(int i=0;i<N;i++)
		{
			arr[i]=scan.nextInt();
		}
		
		//定义公倍数和公约数两个对象
		Divisor div=new Divisor();
		Multiciple mul=new Multiciple();
		
		//输出结果
		System.out.println("这"+N+"个数的最大公约数为"+div.ngcd(arr, N));
		System.out.println("这"+N+"个数的最小公倍数为"+mul.nlcm(arr, N));
	}
}

```java
package com.divisor_and_multiciple;

import java.util.Scanner;
/**
 * 提高要求部分：
 * 大体思路：a0整除a1，b0整除b1这两个条件都是肯定的。x需满足 a1整除x,x整除b1
 * 由于最小公倍数x满足条件是[1,b1]，满足这个条件后，再判断x和a0最大公约是是否为a1，是的话
 * 计数器变量加1，最后输出结果；
 * @author 万物甦醒
 *
 */
public class HighHomework {
	void homeWork2()
	{
		int N=0; 
		int a0,a1,b0,b1;
		int x;
		int flag=0;//记录满足x的个数
		
		Multiciple multiciple=new Multiciple();
		Divisor divisor=new Divisor();
		

		System.out.println("请依次输入组数的大小");
		Scanner scan=new Scanner(System.in);
		N=scan.nextInt(); 		//输入N的大小

		int num[]=new int[N];//定义一维数组，保存每行的结果
		
		//依次执行每组组数
		for(int m=0;m<N;m++)
		{
			flag=0;//计数器，用于记录满足x的个数
			//System.out.print("依次正确输入a0 a1 b1 b2");
			//Scanner scan=new Scanner(System.in);
			a0=scan.nextInt();
			a1=scan.nextInt();
			b0=scan.nextInt();
			b1=scan.nextInt();
			//判断输入的异常性
			if(a1%a0!=0|b1%b0!=0)
			{
				System.out.println("a0没有整除a1或者b0没有整除b1");
				return;
			}
			//先判断符合最小公倍数的x所以情况，在判断这些情况那些是符合最大公约数条件的
			for(x=b1;x>=1;x--)
				//最小公倍数x满足条件是[1,b1]
			{   //是否满足x,b0的最小公倍数是b1	
				if(multiciple.lcm(x,b0)==b1)			
				{
					//满足上述条件的x是否满足和a0的最大公约数是a1
					if(divisor.gcd(x,a0)==a1) 
					{	
						//满足公约数和公倍数条件，加一
						flag++;     		
					}
				}
			}
			num[m]=flag;
		}
		
		for(int n=0;n<N;n++)
		{
			System.out.println(num[n]);	//依次输出每行满足x的个数
		}
		
		
	}
}

package com.divisor_and_multiciple;

import java.util.Scanner;

import javax.tools.Diagnostic;


public class Main {

	public static void main(String[] args)
	{
		int m=1;
		int type;//菜单类型
		BasicHomework h1=new BasicHomework();
		HighHomework h2=new HighHomework();
		while(m!=0)
		{
			System.out.println("请输入作业要求（基本OR提高）");
			System.out.println("1.基本要求");
			System.out.println("2.提高要求");
			System.out.println("3.退出");
			Scanner scan=new Scanner(System.in);
			type=scan.nextInt();        //获取类型
			switch(type)
			{
				case 1:h1.homeWork1();
						break;
				case 2:h2.homeWork2();
						break;
				case 3:m=0;
						break;
				default:
						{System.out.println("输入有误,重新输入");
						continue;
						}
			}
		}
	}

}

运行结果

心得

通过本次实践，有了以下的体会：
（1）编程细节上：定义二维数组时（最后源码改进时将其删去了），由于其特殊性，与一维数组的定义初始化不大相同。一开始错误得声明初始化为：
int arr[N][4]=new int[N][4];
其实正确的写法应是：
int arr[][]=new int[N][4];
（2）再处理异常输入的情况，目前我只通过if-else的嵌套调用语句判断，希望接下来抓紧提前学习JAVA的异常处理的知识。
（2）本次程序设计并没有用到多高级的数据结构，但在保存每行输出结果的时候，采用一维数组num[N]保存（其中N表示行数）。另外在求N个数的最小公倍数和最大公约数，依然采用数组arr[N]保存这N个数（其中这里的N表示输入的N的个数）。
（3）这个星期采用的时面向对象的程序设计，更加深刻理解了面向对象设计的思想，通过对象的调用函数，将一个大的问题简化了不少。