关于最大公约数的三种解法之三（中学课本里的埃拉托色尼筛选法）

以前在中学里学习最大公约数的时候，课本里通常是这样解决问题，比如gcd（m,n），首先，找出m的所有质因数，第二步求出n的质因数，然后从第一步和第二步的质因数分解中找出所有的公因数，假设p是一个公因数，而且在m,n的质因数分解式分别出现过多次，最后将找出的质因数相乘，即为最大公约数，这个过程还是比其他的算法复杂的。

      比如60和24这两个数，我们得到：

                                                          60=2*2*3*5

                                                         24=2*2*3*2

                                                      gcd(60,24)=2*2*3=12

java  代码如下

package Sieve;

import java.util.Scanner;

public class sieve {
               public static int[] getPrime(int n)
               {
            	  int [] result1=new int[n];
            	  for(int i=2;i<n+1;i++)
            	  {
            		  result1[i-1]=i;
            		  
            	  }
            	 for(int i=2;i<n;i++)
            	 {
            		 for(int j=i+1;j<n+1;j++)
            		 {
            			 if(j%i==0)
            				 result1[j-1]=0;
            		 }
            		 
            	 }
            	 int []result2=getNoneZero(result1);
            	 return result2;
            	   
               }
               
               public static int[] getNoneZero(int []A)
               {
            	   int len=0;
            	   for(int i=0;i<A.length;i++)
            	   {
            		 if(A[i]!=0)
            			 len=len+1;
            		   
            		   
            	   }
            	   int[] result =new int[len];
            	   int k=0;
            	   for(int i=0;i<A.length;i++)
            	   {
            		   if(A[i]!=0)
            		   {
            			   result[k]=A[i];
            			   
            			   k++;
            		   }
            		   
            		   
            	   }
            	   return result;
            	   
               }
               
               public static int[] getNprime(int n)
               {
            	  int[] primes=getPrime(n);
            	  int[] result;
            	  int len=0;
            	  for(int i=0;i<primes.length;i++)
            	  {
            		  int temp=n;
            		  while(temp%primes[i]==0)
            		  {
            			  
            			  temp=temp/primes[i];
            			  len++;
            		  }
            	  }
            		  result=new int[len];
            		  int k=0;
            		  for(int i=0;i<primes.length;i++)
            		  {
            			  int temp=n;
            			  while(temp%primes[i]==0)
            			  {
            				  temp=temp/primes[i];
            				  result[k]=primes[i];
            				  k++;
            				
            			  }
            			
            	  }
            	  
            	   
            	   return result;
               }
               public static int[] getCommonPrime(int[] A , int[] B){
                   int[] result;
                   int lenA = A.length;
                   int lenB = B.length;
                   if(lenA < lenB){
                       result = new int[lenA];
                       for(int i = 0;i < lenA;i++){
                           int temp = A[i];
                           for(int j = 0;j < lenB;j++){
                               if(temp == B[j]){
                                   result[i] = A[i];
                                   B[j] = 0;
                                   break;
                               }
                           }
                       }
                   }
                   else{
                       result = new int[lenB];
                       for(int i = 0;i < lenB;i++){
                           int temp = B[i];
                           for(int j = 0;j < lenA;j++){
                               if(temp == A[j]){
                                   result[i] = B[i];
                                   A[j] = 0;
                                   break;
                               }
                           }
                       }
                   }
                   int[] result1 = getNoneZero(result);
                   return result1;
               }
               
               //求取数A和B的最大公约数
               public static void getMaxCommonDivisor(int A,int B){
                   int[] primesA =  getNprime(A);  //数A所有质因子
                   int[] primesB = getNprime(B);   //数B所有质因子
                   int[] resultPrime = getCommonPrime(primesA,primesB);  //数A和数B的公共质因数
                   int maxCommonDivisor = 1;
                   System.out.println(A+"和"+B+"的公共质因数为：");
                   for(int i = 0;i < resultPrime.length;i++){
                       maxCommonDivisor *= resultPrime[i];
                       System.out.print(resultPrime[i]+"\t");
                   }
                   System.out.println();
                   System.out.print(A+"和"+B+"的最大公约数为："+maxCommonDivisor);
               }
               
               public static void main(String[] args){
                   System.out.println("请输入数字A和数字B的值：");
                   Scanner in = new Scanner(System.in);
                   int a = in.nextInt();
                   int b = in.nextInt();
                   getMaxCommonDivisor(a,b);
               }
           }