zoukankan      html  css  js  c++  java
  • 欧几里得算法/欧几里得扩展算法-python

    说在开头。

    出于对欧几里得的尊重,先简单介(cou)绍(ge)一(zi)下(shu).。

    欧几里得,古希腊人,数学家。他活跃于托勒密一世时期的亚历山大里亚,被称为“几何之父”。

    他最著名著作几何原本》是欧洲数学的基础,提出五大公设,欧几里得几何,被广泛的认为是历史上最成功的教科书。

    欧几里得也写了一些关于透视圆锥曲线球面几何学数论的作品。(https://baike.baidu.com/item/欧几里得/182343?fr=aladdin)

    ----------------------------------华丽的分割线-------------------------------------

    以下都是自己在看书的时候自己想的一些思路,总结。

    如有雷同纯属意外,如有异议可以py(仅限妹纸)。

     

    欧几里得算法,也就是所谓的辗转相除法。人话就是,用来求最大公约数的方法。

    证明:gcd(a,b)= gcd(b,a%b)    

    (gcd(a,b),ab的最大公约数;a%b,求模运算,也就是求余数运算。比如7%2=1。)

    证明过程:

    设正整数a,b(a>b)。

    设:gcd(a,b)= c  (我就要设c)

    则有:c|a ,c|b  (|,表示可以整除)

    设: a = bx + a%b  (x为整数,因为这是求余运算,emmm,这样说应该够平易近人了吧)

    则有:a%b = a - bx  (别说恒等变换不知道。233333)

    因为:c|b  则 c|bx  (因为c|b=k(k为整数),则kx为整数,这样说够清楚了吧。emmmm)

    则:(a%b)/c = (a - bx)/ c    ====》  (a%b)/c = a/c - bx/c 

    因为:c|a , c|bx

    故:a/c,bx/c都为整数,所以a/c - bx/c 也是整数。所以(a%b)也可以被c整数。

    即:c |(a%b)

    又因为:c|b

    故:gcd(b,a%b)= c
    则:gcd(a,b)= gcd(b,a%b) 得证。(这证明过程,我就不信全网还有比这写的更清楚的。)

    证明完这个我们就可以通过迭代,反复相除(辗转相除)来求ab的最大公约数了。

    emmm,为什么就可以了呢。因为gcd(a,b)= gcd(b,a%b)  就是一个反复的过程。

    比如我可以继续写:gcd(a,b)= gcd(b,a%b) = gcd(a%b,b%(a%b))=gcd(r1,r2)=·······=gcd(rn-1,rn)(rn=rn-2%rn-1

    这样是不是可以更清楚点了。。

    代码实现(递归实现

    1、递归操作:辗转相除

    2、递归结束操作:余数为0

    """

    这只是一个简单版本。

    比如,对没有最大公约数的情况并没有做判断。

    """

    def gcd(a,b):
    if a < b:
    a, b = b, a
    if a % b != 0:
    return gcd(b,a%b)
    return b

    ------------------------------分割线-------------------------------------
    回家睡觉。明天再说。困死了。
    先给出扩展欧几里得扩展算法实现。有空再说。
    -----------------------------分割线---------------------------------------
    好吧,之前太懒了。。
    万恶的欧几里得扩展:
    简单来说就是求ax+by=z的通解。
    下面就先来说说什么是通解。

    首先,设数字a,b。(a,b>0),并且ab存在最大公约数(没有公约数还玩什么欧几里得)
    设最大公约数为c,即gcd(a,b)= c
    然后有 ax+by = c,求出的这个x,y就是所谓的ax+by=z的通解。
    为什么这么说呢。这里就用简单的办法来说明一下。(可能不严谨)
    拿上面的ax+by=c,两边同时乘一个k,即a(xk)+b(yk)= ck
    同时令ck=z,则有a(xk)+b(yk)= z。刚才已经把x,y求出来了。
    那么这里都乘一个k就可以算出ax+by=z中的x和y了。
    (自认为虽然不严谨,但很清楚了。嗯嗯)

    然后,接下来就来说一说。这个ax+by=gcd(a,b)的解,x,y怎么算吧。
    (不严谨的证明开始了~~)

    ----------------没时间了,明天再说-----------------------

    
    
  • 相关阅读:
    ●BZOJ 3926 [Zjoi2015]诸神眷顾的幻想乡
    ●BZOJ 1396 识别子串
    ●UVA 1608 Non-boring sequences
    ●SPOJ 8222 NSUBSTR–Substrings
    ●SPOJ 7258 Lexicographical Substring Search
    ●CodeForces 429D Trick_Function
    ●BZOJ 2555 SubString
    ●洛谷P2664 树上游戏
    ●洛谷P3168 [CQOI2015]任务查询系统
    【UVA1057】Routing
  • 原文地址:https://www.cnblogs.com/sbxlqswl/p/9839714.html
Copyright © 2011-2022 走看看