GCD & 拓展GCD 入门——HDU - 2504，UVA - 10673

---

gcd，欧几里得算法or辗转相除法，给定两个数a.b，gcd（a,b）用于计算a和b的最大公因数，而拓展欧几里得算法则可以用于求解a*x+b*y=gcd（a,b）等一类问题的解x和y。有很多的问题也可以抽象转化为gcd或者拓展gcd。

---

1.GCD：

int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}

很简单的一行代码就可以实现gcd，也可以从中看出其辗转相除的特点。

====例题：又见GCD HDU - 2504====
这道题给出了gcd(a,b)=c等式中的a和c，要你求出满足（！=c）的最小的b。这里我们可以利用c是a和b 的最大公因数，所以a和b都是c的倍数，通过枚举倍数，检查gcd（a，i*c）是否等于c就可以得到最小的b。

#include<iostream>
#include<algorithm>
#include<string.h>
#include<stdio.h>
using namespace std;
int a,b,c;

int gcd(int a,int b){
    return b==0?a:gcd(b,a%b);
}

int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d",&a,&c);
        for(int i=2;i*c<a;i++){
            b=i*c;
            if(gcd(a,b)==c)break;
        }
        printf("%d\n",b);

    }
    return 0;
}

2.拓展GCD

int extgcd（int a,int b,int &x,int &y）{
    int d=a;
    if(b!=0){
        d=extgcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else{
    x=1;y=0;
    }
    return d;
}

在我看来这是一个递归然后回溯的过程，先不断的辗转相除，递归到最后得到a’x+0*y=gcd(a,0)=a;解出x=1，y=0，然后回溯的过程求出原式的x和y。对于某一层递归有bx’+(a%b)y’=gcd(a,b)，因为a%b=a-(a/b)* b;带入得到ay’+b(x’-(a/b)y’)=gcd(a,b)也就解释了代码。
当c=gcd（a，b）的倍数时将x和y乘上相应倍数就能得出答案。
例题：Play with Floor and Ceil UVA - 10673

#include<iostream>
#include<algorithm>
#include<string.h>
#include<stdio.h>
using namespace std;
long long x,k;
long long a,b;

void ex_gcd(long long a,long long b,long long &g,long long &x,long long &y){
    if(b==0){g=a;x=1;y=0;}
    else{ex_gcd(b,a%b,g,y,x);y-=(a/b)*x;}
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld",&x,&k);
        long long p,q;
        if(x%k==0){p=1;q=k-1;}
        else{
        a=x/k;
        b=x/k+1;
        long long c=1;
        ex_gcd(a,b,c,p,q);
        p=p*x;q=q*x;
        }
        printf("%lld %lld\n",p,q);
    }
    return 0;
}