左程云算法通关041 最大公约数、同余原理

最新推荐文章于 2024-07-06 15:43:38 发布

温柔の风与诗

最新推荐文章于 2024-07-06 15:43:38 发布

阅读量627

点赞数 26

文章标签：算法

本文链接：https://blog.csdn.net/weixin_74776795/article/details/135887363

版权

本文介绍了欧几里得算法及其在求解最大公约数的应用，通过编程实现快速求解。同时探讨了最小公倍数的计算方法以及在LeetCode问题中的实际应用，涉及二分搜索和容斥原理。还提到了针对大整数运算的Stein算法和同余原理，以提高效率。

摘要由CSDN通过智能技术生成

最大公约数

1.欧几里得算法——辗转相除法

公式为gcd(a,b) = gcd(b,a mod b)

例如要求50,30最大公因数，过程为

50%30=20，余数不为0 ，继续递归gcd(30,20)

30%20=10，余数不为0 ，继续递归gcd(20,10)

20%10=0，余数为0，停止，gcd(10,0)此时10即是最大公因数

2.一行代码即可求出最大公约数

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

求最小公倍数

int lcm(int a,int b){
    return (long)a/gcd(a,b)*b;
}

3.证明过程

假设a,b是两个正整数，u是二者的公因子，并且a%b=r;

也就是我们要证明gcd(a,b)=gcd(b,r)

首先我们有a=k*b+r,k=0,1,2...

r=a-k*b,k=0,1,2...

并且我们知道a=s*u,b=t*u;代入可得

r=s*u-k*t*u=u*(s-kt)；这说明u如果是a和b的公因子，那么u也是r的因子。

我们假设b和r的公因子为m,b=m*v,r=n*v;

a=k*m*v+n*v=v*(k*m+n)；这说明m如果是b和r的因子，那么也是a的因子

因此我们得出结果，a和b的公因数集合=b和r的公因数集合。

4.时间复杂度

若a>b,时间复杂度大致为O(log a)的三次方.

5.经典习题

力扣（LeetCode）官网 - 全球极客挚爱的技术成长平台

LeetCode878 第N个神奇数字

一个正整数如果能被 a 或 b 整除，那么它是神奇的。

给定三个整数 n , a , b ，返回第 n 个神奇的数字。因为答案可能很大，所以返回答案对 10^9+ 7取模后的值。

本题思路：利用二分搜索和容斥原理。

粗略定出第n个数字的范围为（1，n*a）。因为假设没有b的加入，那么第n个神奇的数字必然是n*a,但是由于有b加入，所以搜索范围会缩短。

我们假设有一个函数f(x)可以计算出从1到x一共有多少个神奇的数字，那么我们就可以利用二分答案法来计算出第n个神奇的数字在哪里。那么核心点就在于如何找到这个f(x)

f(x)：在1-x范围上，能被a整除的有x/a这么多数字，能被b整除的有x/b这么多数字，有些数字既能被a又能被b整除，所以我们计算重复了，应该减去这一部分，也就是x/lcm(a,b)。这和离散数学里的容斥原理是一样的。

class Solution {
    public int nthMagicalNumber(int n, int a, int b) {
        long lcm=lcm(a,b);
        long answer=0;
        for(long l=0,r=(long)n*Math.min(a,b),m=0;l<=r;){
            m=(l+r)/2;
            if(m/a+m/b-m/lcm >=n){
                r=m-1;
                answer=m;
            }
            else{
                l=m+1;
            }
        }
        return (int)(answer%1000000007);
    }
    public long gcd(long a,long b){
        return b==0?a:gcd(b,a%b);
    }
    public long lcm(long a,long b){
        return (long)a*b/gcd(a,b);
    }
}