QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#736533#9622. 有限小数xiangguiWA 32ms3912kbC++142.0kb2024-11-12 11:40:112024-11-12 11:40:11

Judging History

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

  • [2024-11-12 11:40:11]
  • 评测
  • 测评结果:WA
  • 用时:32ms
  • 内存:3912kb
  • [2024-11-12 11:40:11]
  • 提交

answer

#include <bits/stdc++.h>
#define int long long

using namespace std;

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

// bool getans(int cur2) {
//     int x, y;
//     exgcd()
// }

void solve() {
    int a, b;
    scanf("%lld%lld", &a, &b);
    int tmp = 1;
    int tb = b;
    while (tb % 2 == 0) {
        tb /= 2;
        tmp *= 2;
    }
    while (tb % 5 == 0) {
        tb /= 5;
        tmp *= 5;
    }
    // cerr << b << ' ' << tmp << '\n';
    int cur = 1;
    int ans = 1000000000, ans2 = 0;
    // int 
    auto getans = [&tb, &b, &a](int cur2, int d) -> int {
        int x, y;
        // cerr << b << ' ' << 1 << '\n';
        // if (cur2 == 5) cerr << b << ' ' << d << '\n';
        int g = exgcd(tb, d, x, y);
        y = -y;
        // if (cur2 == 5)
        // cerr << x << ' 'a << y << '\n';
        x *= cur2 * a;
        y *= cur2 * a;
        y %= tb / g;
        y += tb / g;
        y %= tb / g;
        return y;
    };
    for (int i = 0; i <= 20; ++i) {
        // int cu
        int cur2 = cur;
        for (int j = 0; j <= 20; ++j) {
            // cerr << cur2 * tmp << '\n';
            int tmp2 = getans(cur2, __gcd(tmp, cur2));
            // if (tmp2 == 0) cerr << cur2 << '\n';
            // if (cur2 )
            // if (cur2 == 2) cerr << cur2 * tmp << ' ' << tmp2 << '\n';
            // if (tmp2 == 1) cerr << cur2 << '\n';
            if (tmp2 < ans) {
                ans = tmp2;
                ans2 = cur2 * b / __gcd(tmp, cur2);
            }
            if (cur2 > 1000000000 / (5 * b)) break;
            cur2 *= 5;
        }
        if (cur > 1000000000 / (2 * b)) break;
        cur *= 2;
    }
    printf("%lld %lld\n", ans, ans2);
}

signed main() {
    int T;
    scanf("%lld", &T);
    while (T--) {
        solve();
    }
    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 3912kb

input:

4
1 2
2 3
3 7
19 79

output:

0 2
1 3
1 4375
3 316

result:

ok 4 case(s)

Test #2:

score: -100
Wrong Answer
time: 32ms
memory: 3904kb

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 828125000
1 60
1 109375
1 231680000
23 960937500
1 36096000
5 326
1 63360
0 64
1 61000
0 32
1 4880
1 10750
1 18500
1 11714560
1 331250
1 898437500
1 31
1 6
1 289600000
1 455000
1 115000000
1 1265625
0 2
1 29687500
0 160
1 415
1 235937500
1 765000000
1 181250
1 2968750
1 4406250
3 24800
1 480
...

result:

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