最大公约数、欧几里得算法、拓展欧几里得算法、同余方程

最大公约数(Greatest Common Divisors)。

欧几里得算法

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

拓展欧几里得算法

int exgcd(int a, int b, int &x, int &y) {
    if (!b) { x = 1; y = 0; return a; }
    int result = exgcd(b, a % b, x, y);
    int temp = x; x = y; y = temp - a / b * y;
    return result;
}

同余方程

应用拓展欧几里得算法。如洛谷P1082:

/* Luogu 1082 同余方程
 * Au: GG
 * ax ≡ 1 (mod b)
 * ax + by = 1
 */
#include <bits/stdc++.h>
using namespace std;

int exgcd(int a, int b, int &x, int &y) {
    if (!b) { x = 1; y = 0; return a; }
    int result = exgcd(b, a % b, x, y);
    int temp = x; x = y; y = temp - a / b * y;
    return result;
}

int main() {
    int a, b, x, y;
    scanf("%d%d", &a, &b);
    exgcd(a, b, x, y);
    printf("%d\n", (x + b) % b);
    return 0;
}
posted @ 2017-08-05 16:57  greyqz  阅读(160)  评论(0编辑  收藏  举报