poj2891

难得见中级中的水题。。

一开始做，想着存下来，如下：

#include <iostream>

using namespace std;

#define LL long long 
#define MAXN 100000000

LL ext_gcd(LL a, LL b, LL& x, LL& y)
{
    LL t, ret;
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    ret = ext_gcd(b, a%b, x, y);
    t = x, x = y, y = t - a / b*y;
    return ret;
}

LL modular_linear_system(LL b [], LL w [], LL k)
{
    LL d, x, y, a = 0, m, n = 1, i;
    for (i = 0; i < k; i++)
        n *= w[i];
    for (i = 0; i < k; i++)
    {
        m = n / w[i];
        d = ext_gcd(w[i], m, x, y);
        a = (a + y*m*b[i])%n;
    }
    return (a + n)%n;
}

int main()
{
    int t;
    while (cin >> t)
    {
        LL * r = new LL[MAXN];
        LL * a = new LL[MAXN];
        for (int i = 0; i < t; i++)
        {
            cin >> r[i] >> a[i];
        }
        LL ans = modular_linear_system(a, r, t);
        cout << ans << endl;
    }
}

就和这个主题的名字一样。。too young too naive...各种RE，根本没有告诉k的范围，乱猜哪行哇~

然后改了改，不存，逐个算就是了

#include <iostream>
#include <algorithm>

using namespace std;

#define LL long long
#define MAXN 100000000

LL ext_gcd(LL a, LL b, LL& x, LL& y)
{
    LL t, ret;
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    ret = ext_gcd(b, a%b, x, y);
    t = x, x = y, y = t - a / b*y;
    return ret;
}
LL modular_linear(LL a, LL b, LL m)
{
    LL d, x, y;
    d = ext_gcd(a, m, x, y);
    if (b % d != 0) return -1;
    return (x * (b / d) % m + m) % m;
}

int main()
{
    LL a1, a2, m1, m2, k, y;
    while (cin>>k)
    {
        bool flag = false;
        cin >> m1 >> a1;
        if (k == 1) y = a1;

        for (int i = 1; i < k; i++)
        {
            cin >> m2 >> a2;
            if (flag) continue;
            LL y = modular_linear(m2, a1 - a2, m1);
            if (y == -1) { flag = true; continue; }
            a1 = a2 + m2 * y;
            m1 = m1*m2 / (__gcd(m1, m2));
            a1 = (a1 % m1 + m1) % m1;
        }
        if (flag) cout << -1 << endl;
        else cout << a1 << endl;
    }
    return 0;
}