zoukankan      html  css  js  c++  java
  • ACM数论之旅3---最大公约数gcd和最小公倍数lcm(苦海无边,回头是岸( ̄∀ ̄))

    gcd(a, b),就是求a和b的最大公约数

    lcm(a, b),就是求a和b的最小公倍数

    然后有个公式

    a*b = gcd * lcm     ( gcd就是gcd(a, b), ( •̀∀•́ ) 简写你懂吗)

    解释(不想看就跳过){

      首先,求一个gcd,然后。。。

      a / gcd 和 b / gcd 这两个数互质了,也就是 gcd(   a / gcd ,b / gcd  )  =  1,然后。。。

      lcm = gcd *  (a / gcd) * (b / gcd)

      lcm = (a * b) / gcd

      所以。。a*b = gcd * lcm

    }

    所以要求lcm,先求gcd

    辣么,问题来了,gcd怎么求

    辗转相除法

    while循环

    1 LL gcd(LL a, LL b){
    2     LL t;
    3     while(b){
    4         t = b;
    5         b = a % b;
    6         a = t;
    7     }
    8     return a;
    9 }

    还有一个递归写法

    1 LL gcd(LL a, LL b){
    2     if(b == 0) return a;
    3     else return gcd(b, a%b);
    4 }
    5 
    6 LL gcd(LL a, LL b){
    7     return b ? gcd(b, a%b) : a;
    8 }
    9 //两种都可以

    辣么,lcm = a * b / gcd

    (注意,这样写法有可能会错,因为a * b可能因为太大  超出int  或者 超出 longlong)

    所以推荐写成 : lcm = a / gcd * b

    然后几个公式自己证明一下

    gcd(ka, kb) = k * gcd(a, b)

    lcm(ka, kb) = k * lcm(a, b)

    上次做题碰到这个公式

    lcm(S/a, S/b) = S/gcd(a, b)

    S = 9,a = 4,b = 6,小数不会lcm,只好保留分数形式去通分约分。

    当我看到右边那个公式。。。。

    (╯°Д°)╯┻━┻

    这TM我怎么想的到,给我证明倒是会证。 T_T

  • 相关阅读:
    出错处理函数abort、exit、atexit、strerror. . .
    linux查询系统信息命令
    [转载]比google和百度强十倍的搜索类网站
    陶  朱  商  经
    ip的划分,超详细.【网管常识】
    linux的hostname修改详解
    勤于寻找谈话资料
    Windows常用命令集
    C语言中printf格式
    How to disable SELinux
  • 原文地址:https://www.cnblogs.com/linyujun/p/5167914.html
Copyright © 2011-2022 走看看