上机作业：4种算法求最大公约数

最新推荐文章于 2021-11-27 18:54:07 发布

浪味仙！

最新推荐文章于 2021-11-27 18:54:07 发布

阅读量343

点赞数

本文链接：https://blog.csdn.net/qq_44710356/article/details/88581517

版权

一.实验目的
1.明确算法的概念和特点。
2.通过对问题的分析，设计合理的算法解决问题；
二.实验内容
运行最大公约数的常用算法，并进行程序的调式与测试，要求程序设计风格良好，并添加异常处理模块（如输入非法等）。
1.辗转相除法
辗转相除法（又名欧几里德法）C语言中用于计算两个正整数a,b的最大公约数。
①函数嵌套调用
其算法过程为：前提：设两数为a,b设其中a 做被除数,b做除数，temp为余数
1、大数放a中、小数放b中；
2、求a/b的余数；
3、若temp=0则b为最大公约数；
4、如果temp!=0则把b的值给a、temp的值给a；
5、返回第二步；
②函数递归调用
2.穷举法（利用数学定义）
穷举法（也叫枚举法）穷举法求两个正整数的最大公约数的解题步骤：从两个数中较小数开始由大到小列举，直到找到公约数立即中断列举，得到的公约数便是最大公约数。
对两个正整数a,b如果能在区间[a,0]或[b,0]内能找到一个整数temp能同时被a和b所整除，则temp即为最大公约数。
3.更相减损法（等值算法）
更相减损术，是出自《九章算术》的一种求最大公约数的算法，它原本是为约分而设计的，但它适用于任何需要求最大公约数的场合。
翻译成现代语言如下：
（1）：任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简；若不是则执行第（2）步。
（2）：以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到所得的减数和差相等为止。
则第一步中约掉的若干个2与第（2）步中等数的乘积就是所求的最大公约数。
其中所说的“等数”，就是最大公约数。求“等数”的办法是“更相减损”法。
4.Stein算法
Stein算法由J. Stein 1961年提出，这个方法也是计算两个数的最大公约数。来研究一下最大公约数的性质，发现有 gcd( kx,ky ) = k*gcd( x,y ) 这么一个非常好的性质。
对两个正整数 x>y ：
1.均为偶数 gcd( x,y ) =2gcd( x/2,y/2 )；
2.均为奇数 gcd( x,y ) = gcd( (x+y)/2,(x-y)/2 )；
2.x奇y偶 gcd( x,y ) = gcd( x,y/2 )；
3.x偶y奇 gcd( x,y ) = gcd( x/2,y ) 或 gcd( x,y )=gcd( y,x/2 )；
现在已经有了递归式，还需要再找出一个退化情况。注意到 gcd( x,x ) = x ，就用这个。
①函数非递归调用
②函数递归调用
三.测试运行结果
在这里插入图片描述
四.经验归纳
五.源代码
1.辗转相除法
（1）嵌套调用

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.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 main()  
{
 int m,n,t,i; 
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=divisor(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组随机数的最大公约数所需时间为:%0.6f",i,duration);
}

（2）递归调用

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
int gcd (int a,int b)
{  if(a%b==0)
       return b;   
else  
       return gcd(b,a%b);
  }
int main()
{
 int m,n,t,i,p; /*定义整型变量*/
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=gcd(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组数据的最大公约数花费时间为:%0.6f",i,duration);
}

2.穷举法

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
int divisor (int a,int b) /*自定义函数求两数的最大公约数*/
{
    int  temp;          /*定义义整型变量*/
    temp=(a>b)?b:a;    /*采种条件运算表达式求出两个数中的最小值*/
    while(temp>0)     
    {
       if (a%temp==0&&b%temp==0) /*只要找到一个数能同时被a,b所整除，则中止循环*/
          break;    
       temp--;      /*如不满足if条件则变量自减，直到能被a,b所整除*/
    }
  return (temp); /*返回满足条件的数到主调函数处*/
}
main()
{
 int m,n,t,i; 
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=divisor(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组随机数的最大公约数所需时间为:%0.6f",i,duration);
}

3.更相减损法

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
int gcd(int m,int n)
{
	int i=0,temp,x;
	while(m%2==0 && n%2==0)  //判断m和n能被多少个2整除
	{
		m/=2;
		n/=2;
		i+=1;
	}
	if(m<n)     //m保存大的值
	{
		temp=m;
		m=n;
		n=temp;
	}
	while(x)
	{
		x=m-n;
		m=(n>x)?n:x;
		n=(n<x)?n:x;
		if(n==(m-n))
			break;
	}
	if(i==0)
		return n;
	else 
		return (int )pow(2,i)*n;
}
int main()
{
 int m,n,t,i; 
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=gcd(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组随机数的最大公约数所需时间为:%0.6f",i,duration);
}

4.Stein算法
（1）递归

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
int gcd(int u,int v)
{
    if (u == 0) return v;
    if (v == 0) return u;  // look for factors of 2
    if (~u & 1) // u is even
    {
        if (v & 1)  // v is odd
            return gcd(u >> 1, v);
        else  // both u and v are even
            return gcd(u >> 1, v >> 1) << 1;
    }
     if (~v & 1)  // u is odd, v is even
        return gcd(u, v >> 1); // reduce larger argument
    if (u > v)
        return gcd((u - v) >> 1, v);
     return gcd((v - u) >> 1, u);
}
int main()
{
 int m,n,t,i; 
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=gcd(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组随机数的最大公约数所需时间为:%0.6f",i,duration);
}

（2）非递归

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#include<math.h>
int Stein( unsigned int x, unsigned int y )  /*求最大公约数*/
{
        int factor = 0;
        int temp;
        if ( x < y )
        {
                temp = x;
                x = y;
                y = temp;
        }
        if ( 0 == y )
        {
                return 0;
        }
        while ( x != y )
        {
                if ( x & 0x1 )
                {/* when x is odd */
                        if ( y & 0x1 )
                        {/* when x and y are both odd */
                                y = ( x - y ) >> 1;
                                x -= y;
                        }
                        else
                        {/* when x is odd and y is even */
                                y >>= 1;
                        }
                }
                else
                {/* when x is even */
                        if ( y & 0x1 )
                        {/* when x is even and y is odd */
                                x >>= 1;
                                if ( x < y )
                                {
                                        temp = x;
                                        x = y;
                                        y = temp;
                                }
                        }
                        else
                        {/* when x and y are both even */
                                x >>= 1;
                                y >>= 1;
                                ++factor;
                        }
                }
        }
        return ( x << factor );
}
int main()
{
 int m,n,t,i; 
 double duration;  /*测量一个事件的持续时间*/
 int j[100],k[100];
 clock_t start,finish;
 srand((int)time(0));
 for(i=0;i<100;i++)  /*控制随机数的范围*/
	{
		j[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		k[i]=1+(int)(100.0*rand()/(RAND_MAX+1.0));
		printf("第%d组随机数为:%d\t%d\n",i+1,j[i],k[i]);
	}
 start=clock();
 for(i=0;i<100;i++)
	{
	 m=j[i];
	 n=k[i];
	 t=Stein(m,n);
	 printf("The higest common divisor in %d and %d is %d\n",j[i],k[i],t); /*输出最大公约数*/
	}
 finish=clock();
 duration=(double)(finish-start)/CLOCKS_PER_SEC;
 printf("计算%d组随机数的最大公约数所需时间为:%0.6f",i,duration);
}