HDU-3579-Hello Kiki （利用拓展欧几里得求同余方程组）

#include "bits/stdc++.h"
using namespace std;
int a[105], b[105];
// 拓展欧几里得C++模板 
int ex_gcd(int a, int b, int &x, int &y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    int ans = ex_gcd(b, a % b, y, x);
    y -= a / b * x;
    return ans;
}
int solve(int n) {
    // 任何数对1取余都是0,所以初始化ans = 0, lcm = 1; 
    int x, y, ans = 0, lcm = 1;
    for (int i = 0; i < n; i++) {
        int gcd = ex_gcd(lcm, a[i], x, y);
        if ((b[i] - ans) % gcd) {
            return -1;
        }
        // 拓展欧几里得求得的x和y是对于gcd而言的。乘完之后才是对于 (b[i] - ans) 的 x 
        x *= (b[i] - ans) / gcd;
        a[i] /= gcd;
        // 使 x 成为最小非负数解 
        x = (x % a[i] + a[i]) % a[i];
        // 更新ans 
        ans += x * lcm;
        // 更新lcm 
        lcm *= a[i];
    } 
    // 本题要求最小正整数解，如果ans是0，要加一个最小公倍数也就是lcm 
    return ans ? ans : lcm;
}
int main() {
    int t, n, ca = 1;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &n);
        for (int i = 0; i < n; i++) {
            scanf("%d", &a[i]);
        }
        for (int i = 0; i < n; i++) {
            scanf("%d", &b[i]);
        }
        printf("Case %d: %d\n", ca++, solve(n));
    }
    return 0;
}