码迷,mamicode.com
首页 > 编程语言 > 详细

公约数算法

时间:2014-10-31 23:26:05      阅读:232      评论:0      收藏:0      [点我收藏+]

标签:blog   io   os   sp   数据   div   log   bs   amp   

/*
对于已知的两个自然数m, n,假设m>n
计算m除以n,将得到的余数记做r
如果r=0,则此时的n为求得的最大公约数。否则,将n的值保存在m中,将r的值保存在n中,
重复执行下去。
*/

//欧几里得->辗转相除法


#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>
#include <iostream>
#include <iomanip>
#include <algorithm>

using namespace std;

int gcd(int a, int b)
{
    int m, n, r; //m大 n小
    if(a>b)
    {
        m=a; n=b;
    }
    else
    {
        m=b; n=a;
    }
    r=m%n;
    while(r!=0)
    {
        m=n;
        n=r;
        r=m%n;
    }
    return n;
}

int main()
{
    int n, m, dd;
    while(scanf("%d %d", &n, &m)!=EOF)
    {
        dd=gcd(n, m);
        printf("%d\n", dd);
    }
    return 0;
}

 

/*
欧几里得算法的执行效率很高,但是如果参与运算的数据非常大的话就暴露了缺点。
计算机中的整数最多是64位,如果低于64位,取模运算比较简单,直接用%就行。
但是如果对于计算两个超过64位的整数的模,gcd算法效率十分低,这是便有了Stein算法。
*/


#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>
#include <iostream>
#include <iomanip>
#include <algorithm>

using namespace std;

int Stein(int a, int b)
{
    int m, n, r;
    if(a>b)
    {
        m=a; n=b;
    }
    else
    {
        m=b; n=a;
    }
    if(n==0)
      return m;
    if(m%2==0 && n%2==0 )
    {
        return 2*Stein(m/2, n/2);
    }
    if(m%2==0)
      return Stein(m/2, n);
    if(n%2==0)
      return Stein(m, n/2);
    return Stein((m+n)/2, (m-n)/2);
}

int main()
{
    int n, m, dd;
    while(scanf("%d %d", &n, &m)!=EOF)
    {
        dd=Stein(n, m);
        printf("%d\n", dd);
    }
    return 0;
}

 

公约数算法

标签:blog   io   os   sp   数据   div   log   bs   amp   

原文地址:http://www.cnblogs.com/yspworld/p/4066037.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!