QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#760479#9622. 有限小数1205691775WA 35ms4080kbC++232.9kb2024-11-18 17:06:582024-11-18 17:06:59

Judging History

你现在查看的是最新测评结果

  • [2024-11-18 17:06:59]
  • 评测
  • 测评结果:WA
  • 用时:35ms
  • 内存:4080kb
  • [2024-11-18 17:06:58]
  • 提交

answer

#pragma optimize GCC("Ofast")
#include <bits/stdc++.h>

using i64 = long long;
constexpr int fk_max = 1e9;


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

std::array<i64, 3> exgcd(i64 a, i64 b) {
    if (b == 0) {
        return {a, 1, 0};
    }
    auto [g, x, y] = exgcd(b, a % b);
    return {g, y, x - a / b * y};
}

i64 inline floor(i64 a, i64 b) {
    return a / b - ((a ^ b) < 0 && a % b);
}
i64 inline ceil(i64 a, i64 b) {
    return -floor(-a, b);
}

std::array<i64, 2> inline valid(i64 a, i64 b, i64 c, i64 d, i64 x, i64 y) {
    // auto [d, x, y] = exgcd(a, b);
    if (c % d) return {0, -1};
    x *= c / d;
    y *= c / d;
    i64 dx = b / d;
    i64 dy = -a / d;
    // std::cerr << a << " " << b << " " << c << " " << x << " " << y << " " << "\n";
    i64 k = dx > 0 ? ceil(-x, dx) : floor(-x, dx);
    x += k * dx, y += k * dy;

    // std::cerr << a << " " << b << " " << c << " " << x << " " << y << " " << "\n";
    return {y >= 0, x};

}
using namespace std;
map<pair<int,int>,pair<int,int>> ans;
void solve(){
    i64 a, b;
    std::cin >> a >> b;
    if (ans.count({a, b})) {
        auto [c,d] = ans[{a,b}];
        cout<< c<< ' '<<d<< '\n';
        return;
    }
    i64 remb = b;
    while(remb % 2 == 0){
        remb /= 2;
    }
    while(remb % 5 == 0) {
        remb /= 5;
    }

    std::vector<i64> zy;
    for (int i = 1; 1ll * i * i <= remb; ++i) {
        if (remb % i == 0) {
            zy.push_back(i);
        }
        if (1ll * i * i == remb) continue;

        zy.push_back(remb / i);
    }
    std::sort(std::begin(zy), std::end(zy));

    i64 c = 1e18, d = 1e18;

    __int128 cx_a = a;
    if (cx_a % remb == 0) {
        std::cout << "0 1\n";
        return;
    }

    // ad + bc + xbd = 0
    // for (auto zz : zy) {
    auto [g, x, y] = exgcd(b, -remb * b);
    // i64 cnt = 0;
    for(i64 i = 1; i <= fk_max; i = i * 2) {
        for (i64 j = i; j <= fk_max; j *= 5) {
            // cnt++;
            // continue;
            i64 dd = j  * b;
            if (dd > fk_max) break;
            auto [ok, cc] = valid(b, - remb *  b , - a * dd, g, x, y);
            if (!ok) continue;
            // std::cerr << ": " << dd  << "\n";
            i64 gg = std::gcd(cc, dd);
            cc /= gg;
            dd /= gg;
            if (cc < c) {
                c = cc;
                d = dd;
                if (c == 1) break;
            } else if (cc == c) {
                d = std::min(dd, d);
            }
        }
    }
        // std::cout<< cnt<< '\n';
        // break;
    // }
    std::cout << c << " ";
    std::cout << d;
    std::cout <<  "\n";
    ans[{a,b}] = {c,d};
}

int main() {
    std::cin.tie(nullptr);
    std::ios::sync_with_stdio(false);

    int t;
    std::cin >> t;
    while(t--){
        solve();
    }

    return 0;
}

詳細信息

Test #1:

score: 100
Accepted
time: 1ms
memory: 3516kb

input:

4
1 2
2 3
3 7
19 79

output:

0 1
1 3
1 14
3 316

result:

ok 4 case(s)

Test #2:

score: -100
Wrong Answer
time: 35ms
memory: 4080kb

input:

10000
11 12
28 53
17 60
2 35
17 181
80 123
68 141
79 163
71 99
13 64
33 61
15 32
16 61
11 86
33 74
128 143
40 53
7 23
30 31
5 6
86 181
73 91
13 23
71 81
1 2
7 38
117 160
33 83
129 151
88 153
25 58
16 19
19 141
95 124
43 96
71 139
11 59
106 109
93 152
34 43
17 99
1 57
20 159
16 25
5 73
159 170
172 17...

output:

1 12
1 54272
1 60
1 7
1 231680000
23 3936
1 36096000
5 326
1 63360
0 1
1 31232
0 1
1 4880
1 10750
1 18500
1 11714560
1 331250
1 2944
1 31
1 6
1 289600000
1 455000
1 58880
1 51840
0 1
1 304
0 1
1 415
1 19328000
1 765000000
1 4640
1 608
1 72192
3 775
1 48
3 347500
1 944
1 43600
1 76
1 430000
1 6336
1 ...

result:

wrong answer Jury found better answer than participant's 1 < 2 (Testcase 8812)