HDU 1573 Romantic

题目链接 HDOJ 1573-Romantic

这个是我学数论的第一道题，扩展欧几里得算法的入门题了。下面我尝试着给出我对扩展欧几里得的解释：

首先是欧几里得算法，为了同扩展欧几里得区别开，我们用前者的另一个名字称呼之——辗转相除法。

辗转相除法在c++代码中的实现是一个递归的思路，它基于一个基本的数论结论：（记 gcd 为两个数a，b的最大公约数）

gcd ( a , b )  =  gcd ( b , a%b )

证明如下：

一方面，平凡的，存在k,r∈Z，使得a=kb+r，即a mod b = r

设d为a,b的公约数， 即d|a , d|b,由于 r=a-kb , 显然d|r

因此 d 是b和a mod b的公约数

另一方面，若d是b和a mod b的公约数，即d|b , d|r,

由于 a=kb+r, 同理有 d|a.

因此 d 是 a 和b 的公约数.

综上，（a,b）和（b， a mod b）有完全相同的公约数，自然也有相同的最大公约数，即gcd

这就是辗转相除法的理论依据，由此有个粗暴的递归求gcd的方法：

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

最后跳出递归的条件即 a%b=0，之前看到过一个文章，貌似这个方法的复杂度是 O(logN)

接下来是扩展欧几里得，它基于一个数论上的定理：对于给定的a,b ，方程 ax+by = gcd ( a , b )必有一组 (x,y ) 解，这组解由一个通式导出。

扩展欧几里得算法能够得到gcd，同时得到这组解。

推导：

不妨设 a>b ,两者相等的情况是平凡的

1. 当b=0 时，显然gcd (a,b)=a, 这时x=1,y=0

2. 基于递归：设 ax1+by1=gcd(a,b)

再设 bx2+(a mod b)y2=gcd(a,a mod b)

由于 gcd(a,b)=gcd(a mod b)

有 ax1+by1=bx2+(a mod b)y2

即 ax1+by1= bx2+ [a-(a/b)*b]y2         //这里是整数意义上的除法，求余的另一种写法而已

进而 ax1+by1 = bx2+ay2-(a/b)*b*y2

因此 ax1+by1 = ax2+b[x2-(a/b*y2) ]

考虑到a和b是给出的，上式为恒等式，对应项相等，因此有：

x1=x2   y1=x2-(a/b)*y2

这就得到了解的递推式，这样的递归是有限次的，可以终止的，因为总会有一次递归出现 a/b=0 的情况，这样自底向上得到解

进而还能有一个推论 ：关于a,b整数方程 ax+by=c,若 gcd(a,b)|c , 该方程存在整数解；否则无整数解。

且这组解的通式为：

x = x1 + b/gcd(a , b)*k

y = y1  - a/gcd(a , b)*k      k∈Z 

其中x1,y1是扩展欧几里得实现中得到的一组特解。

有了以上的知识，这个问题就可以解决了，题意即寻找方程 ax+by=1 中x为最小正整数的解

这里就相当于令c=1，于是用扩展欧几里得求gcd，如果不为一，必然无解，输出sorry即可。否则循环求出符合要求的特解，具体见代码

#include<cstdio>

using namespace std;
int ex_eucild(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    {
        int gcd=ex_eucild(b,a%b,x,y);
        int temp=x;
        x=y;
        y=temp-a/b*y;
        return gcd;
    }
}
int m,n;
int main()
{
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        int x,y;
        int ans=ex_eucild(m,n,x,y);
        if(ans==1)
        {
            while(x<0)
            {
                x+=n;
                y-=m;
            }
            printf("%d %d\n",x,y);
        }
        else
        {
            printf("sorry\n");

        }
    }
}

数论难就难在推导证明了。。。好歹自己是数学系的，也别太懒了。。