【笔记】Binary Extended Euclidean Algorithm

最新推荐文章于 2021-12-18 20:17:59 发布

xueyan0096

最新推荐文章于 2021-12-18 20:17:59 发布

阅读量405

点赞数

本文链接：https://blog.csdn.net/xueyan0096/article/details/83716647

版权

本文详细介绍了扩展欧几里得算法，包括除法版本、二进制版本以及无符号和有符号版本。该算法用于求解非负整数a和b的线性组合，得到最大公约数gcd(a, b)。文章通过代码示例和逻辑分析阐述了算法的实现过程，并特别讨论了在二进制移位操作中如何处理偶数和奇数的情况。" 5772579,895694,Windows环境下Xapian搜索引擎的编译与配置教程,"['Windows开发', 'C++', '搜索引擎', 'Xapian', '编译工具']

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

扩展欧几里得算法

给定非负整数a, b，求解向量(u1, u2, u3)，使得au1 + bu2 = u3 = gcd(a, b)。

扩展欧几里得算法的除法版本

引入辅助向量(v1, v2, v3)，使得av1 + bv2 = v3 代码如下：

/// <summary>
/// 返回{x, y, gcd}, 使得<paramref name="a"/>x + <paramref name="b"/>y = gcd(<paramref name="a"/>, <paramref name="b"/>)
/// </summary>
public static int[] ExtendedEuclid(int a, int b)
{
    if (a < 0 || b < 0)
        throw new ArgumentOutOfRangeException("a and b should all be nonegative");

    if (a == 0) return new int[] { 0, 1, b };
    if (b == 0) return new int[] { 1, 0, a };

    var aZeros = a.NumberOfTrailingZeros();
    var bZeros = b.NumberOfTrailingZeros();
    var k = Math.Min(aZeros, bZeros);
    a >>= k; b >>= k;

    var u = new int[] { 1, 0, a };
    var v = new int[] { 0, 1, b };
    var t = new int[3];
    while(v[2] != 0)
    {
        var q = u[2] / v[2];
        t[0] = u[0] - v[0] * q;
        t[1] = u[1] - v[1] * q;
        t[2] = u[2] - v[2] * q;

        u[0] = v[0]; u[1] = v[1]; u[2] = v[2];

        v[0] = t[0]; v[1] = t[1]; v[2] = t[2];
    }

    u[2] <<= k;

    return u;
}