【hiho一下-95】 扩展欧几里得算法

      例题是hiho的一道题：

      题目

      这道题主要用到了欧几里得扩展算法，并且求一下最小正整数解就好了，没什么特别的。不过，再次学习了一下欧几里得扩展算法，有了更深的理解。欧几里得扩展算法可以有两种形式：

      1. ax + by = c = bx +(a%b)y 证明很容易,根据欧几里得算法，设t=(a,b)，a=k1*t, b=k2*t, (a%b)=k3*t，那么c也能表示成t的倍数，只要修改x和y，就一定能相等。

      2. a*t1(k-1) + b*t2(k-1) = t3(k-1)

          a*t1(k) + b*t2(k) = t3(k)

          a*t1(k+1) + b*t2(k+1) = t3(k+1)

          t3(k+1) = t3(k-1)%t3(k)

      只要一直迭代下去，直到t3(k+1)=0，那么t3(k)就是最小公约数。有人可能会迷惑，这个公约数是谁和谁的呢？那么如果我们设置t3(0)=a，t3(1)=b，那么就是(a,b)，由此可以得到最初的两个式子：

      a*1 + b*0 = a

      a*0 + b*1 = b

      只要推理一下，就能得到递推公式：

      t3(k+1) = t3(k-1) % t3(k)

      t1(k+1) = t1(k-1) - t1(k)*tmp

      t2(k+1) = t2(k-1) - t2(k)*tmp

      tmp = t3(k-1)/t3(k)

      证明方法请自行百度一下：证明

      下面是题目的代码：

//
//  hiho.cpp
//  hiho
//
//  Created by lhq on 16/4/25.
//  Copyright © 2016年 lhq. All rights reserved.
//

#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;

/*
 ax + by = gcd = bx + a%by = bx + (a - a/b * b)y = ay + (x - a/b*y)b
 
 v1 > v2
 s1 + v1*t = s2+ v2*t - km
 (v1-v2)*t + mk = s2 - s1
 
 v2 > v1
 s2 + v2*t = s1 + v1*t -km
 (v2-v1)*t + km = s1 - s2
 
 at + mk = gcd
 求t,k   t,k > 0 且为整数
 
 ax + by = gcd
 a*(x + b/gcd * u)
 b*(y - a/gcd * u)
 
 求最小的x
 
 */

// ax + by = c  a = b , b = a%b 直到b = 0 主动改变a,b
long long gcd_ex(long long& x, long long& y, long long a, long long b){
    if(b == 0){
        x = 1;
        y = 0;
        return a;
    }
    long long tmp = gcd_ex(x, y, b, a%b);
    long long tmpy = y;
    y = x - (a/b)*y;
    x = tmpy;
    return tmp;
}

// at1 + bt2 = t3   t3 = t1%t2 直到t3=1 主动改变x,y  注意：初始值设置为，t3(k-1)=a, t3(k)=b，所以肯定能得到t3=(a,b)
long long gcd_ex_2(long long& x, long long& y, long long a, long long b){
    long long to_1,to_2,to_3;
    long long tm_1,tm_2,tm_3;
    long long tn_1,tn_2,tn_3;
    //不存在最大公约数
    if(a == 0 || b == 0){
        x = y = 0;
        return 0;
    }
    to_1 = 1; to_2 = 0; to_3 = a;
    tm_1 = 0; tm_2 = 1; tm_3 = b;
    for(tn_3 = to_3%tm_3; tn_3!=0; tn_3=to_3%tm_3){
        long long tmp = to_3 / tm_3;
        tn_1 = to_1 - tmp * tm_1;
        tn_2 = to_2 - tmp * tm_2;
        
        to_1 = tm_1; to_2 = tm_2; to_3 = tm_3;
        tm_1 = tn_1; tm_2 = tn_2; tm_3 = tn_3;
    }
    x = (tm_1%(b/tm_3) + b/tm_3)%(b/tm_3);   //最小正整数解
    y = (tm_2%(b/tm_3) + a/tm_3)%(a/tm_3);   //最小正整数解
    return tm_3;
}

int main(){
    long long s1,s2,v1,v2,m;
    long long x,y,a,b,c,gcd;
    while(cin>>s1>>s2>>v1>>v2>>m){
        a = (v1 - v2);
        b = m;
        c = (s2 - s1);
        if(a<0){
            a = -a; c = -c;
        }
        if(c<0){
            c += m;
        }
        gcd = gcd_ex_2(x, y, a, b);
//        gcd = gcd_ex(x, y, a, b);
        if(c % gcd){
            cout<<-1<<endl;
        }else{
            b = b/gcd;
            x = (x*c/gcd)%b;
            if(x>0){
                cout<<x<<endl;
            }else{
                while(x<0){
                    x += b;
                }
                cout<<x<<endl;
            }
        }
    }
    return 0;
}