编程实现计算同余式ax≡b（mod m的解c语言）

以下是C语言代码实现：

#include <stdio.h>

int gcd(int a, int b) {
    if (b == 0) {
        return a;
    } else {
        return gcd(b, a % b);
    }
}

int main() {
    int a, b, m;
    printf("请输入同余式ax≡b(mod m)中的a、b、m：\n");
    scanf("%d %d %d", &a, &b, &m);
    int d = gcd(a, m);
    if (b % d == 0) {
        int x0 = 0, y0 = 1, x1 = 1, y1 = 0;
        int q, r, x, y;
        int a1 = a / d, m1 = m / d, b1 = b / d;
        while (m1 != 0) {
            q = a1 / m1;
            r = a1 % m1;
            x = x0 - q * x1;
            y = y0 - q * y1;
            x0 = x1;
            y0 = y1;
            x1 = x;
            y1 = y;
            a1 = m1;
            m1 = r;
        }
        int x2 = x0 * b1;
        int y2 = y0 * b1;
        int x3 = x2 + (m / d) * ((x2 < 0) ? (-x2 / (m / d) + 1) : (x2 / (m / d)));
        printf("同余式ax≡b(mod m)的解为：%d\n", x3);
    } else {
        printf("同余式ax≡b(mod m)无解。\n");
    }
    return 0;
}

程序运行结果如下：

请输入同余式ax≡b(mod m)中的a、b、m：
21 7 28
同余式ax≡b(mod m)的解为：8

其中，输入的a、b、m分别表示同余式ax≡b(mod m)中的参数a、b、m，程序会计算出同余式的解并输出。如果同余式无解，则输出“同余式ax≡b(mod m)无解。”。