扩展欧几里得算法

标签：

首先、扩展欧几里得定理：对于两个不全为0的整数a、b，必存在一组解x,y，使得ax+by==gcd(a,b);

int gcd(int a,int b)
{
    int t,d;
    if(b==0)
    {
        x=1;
        y=0;   //不明处1
return a;
    }
    d=gcd(b,a%b);
    t=x;
    x=y;
    y=t-(a/b)*y;  //不明处2
return d;
}

x y 用全局变量表示

不明处1：由扩展欧几里得定理：ax+by==gcd(a,b)---式1，而此时b==0，也就是说gcd(a,0)==a。原式变为ax+by==a --> x==1,y==0。应该够清楚了吧

不明处2：这里先说明一下我的一些规则，x,y表示第一次递归时的值，x1,y1表示第二次递归时的值。那么

gcd(a,b)==gcd(b,a%b)，同时都代入式1，有ax+by==b*x1+(a%b)*y1。将右边变形一下

b*x1+(a%b)*y1==b*x1+(a-(a/b)*b)*y1==a*y1+b*(x1-(a/b)*y1)，最终得到ax+by==a*y1+b*(x1-(a/b)*y1)

也就是说，上一深度的x等于下一深度的y1，上一深度的y等于下一深度的x1-(a/b)*y1。    需要注意，上面推导时用的除法都是整型除法

那么对于一般的不定式ax+by==c，它的解应该是什么呢。很简单，x1=x*(c/gcd(a,b)),y1=y*(c/gcd(a,b))。

    ax+by=gcd(a,b)=gcd(b,a%b)=bx+(a-(int)a/b*b)y=ay+b(x-(a-(int)a/b*y)

 1 int extended_gcd(int a, int b, int &x, int &y)
 2 {
 3   int ret, tmp；
 4         if (!b) {
 5             x = 1; y = 0;  return a;
 6         }
 7         ret = extended_gcd(b, a % b, x, y);
 8         tmp = x;
 9         x = y;
10         y = tmp - a / b * y;
11        return ret;
12 }

      

bool modularLinearEquation(int a,int b,int n)
     {
        int x,y,x0,i;
        int d=Extended_Euclid(a,n,x,y);     //ax=b (mod n) 等价于ax+ny=b
        if(b%d) 
            return false;
        x0=x*(b/d)%n;
        for(i=1;i<=d;i++)
           printf("%d\n",(x0+i*(n/d))%n);
        return true;
     }

扩展欧几里得算法

标签：

原文地址：http://www.cnblogs.com/linliu/p/4914556.html