HDU 1573（X问题）

关注

转载最新推荐文章于 2026-03-15 12:50:26 发布 278 阅读 AI 写同款· GEO 优化 ›

X问题 HDU - 1573(中国剩余定理水题) 求在小于等于N的正整数中有多少个X满足：X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 < a[i] <= 10)。 Input 输入数据的第一行为一个正整数T，表示有T组测试数据。每组测试数据的第一行为两个正整数N，M (0 < N <= 1000... 阅读详情

#include <iostream>
using namespace std;
const int MAXM = 15;

int A[MAXM], B[MAXM];

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

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

int inv(int a, int b)
{
    int x = 0, y = 0;
    exgcd(a, b, x, y);
    x = (x % b + b) % b;
    if (!x)
        x += b;
    return x;
}

int main()
{
    int T;
    cin >> T;
    while (T--)
    {
        int N, M;
        cin >> N >> M;
        bool flag = true;
        for (int i = 1; i <= M; i++)
            cin >> A[i];
        for (int i = 1; i <= M; i++)
        {
            cin >> B[i];
            B[i] %= A[i];
        }

        int m1, m2, c1, c2, t, ans;
        for (int i = 2; i <= M; i++)
        {
            m1 = A[i - 1]; m2 = A[i]; c1 = B[i - 1]; c2 = B[i];
            t = gcd(m1, m2);
            if ((c2 - c1) % t != 0)
            { 
                flag = false;
                break; 
            }
            A[i] = m1 * m2 / t;
            B[i] = ((inv(m1 / t, m2 / t) * ((c2 - c1) / t) % (m2 / t) * m1 + c1) % A[i] + A[i]) % A[i];
        }
        if (!flag)
            cout << 0 << endl;
        else
        {
            if (N < B[M])
            { 
                cout << 0 << endl;
                continue; 
            }
            ans = (N - B[M]) / A[M] + 1;
            if (!B[M])
                --ans;
            cout << ans << endl;
        }
    }
    return 0;
}