QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#737611#9622. 有限小数Zxyoul#WA 238ms3860kbC++232.1kb2024-11-12 16:25:132024-11-12 16:25:16

Judging History

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

  • [2024-11-12 16:25:16]
  • 评测
  • 测评结果:WA
  • 用时:238ms
  • 内存:3860kb
  • [2024-11-12 16:25:13]
  • 提交

answer

#include<bits/stdc++.h>
#define ll __int128
using namespace std;
const int maxn = 1e6 + 5;
const int B = 1000000000;
const ll INF = 1e18;
ll extend_gcd(ll a, ll b, ll &x, ll &y) {
    if(!b) {
        x = 1; y = 0;
        return a;
    }
    ll xt = 0, yt = 0;
    ll d = extend_gcd(b, a % b, xt, yt);
    x = yt;
    y = xt - yt * (a / b);
    return d;
}
void solve() {
    int aa, bb; cin >> aa >> bb;
    ll a = aa, b = bb;
    ll e = b;
    while(e % 2 == 0) { e /= 2; }
    while(e % 5 == 0) { e /= 5; }
    if(e == 1) {
        cout << 0 << " " << 1 << '\n';
        return;
    }
    vector< int > vec;
    for(int i = 1; i * i <= e; i++) {
        if(e % i == 0) {
            int c = i; vec.push_back(c);
            if(c * c == e) continue;
            c = e / i; vec.push_back(c);
        }
    }
    int ansc = b - a, ansd = b;
    // 枚举 d 的可能性
    for(auto item : vec) {
        for(ll i = 0, cur2 = 1; 1LL * cur2 * item <= B; i++, cur2 *= 2) {
            for(ll j = 0, cur5 = 1; 1LL * cur2 * item * cur5 <= B; j++, cur5 *= 5) {
                ll d = item * cur2 * cur5;
                ll down = b * d / __gcd(b, d);
                ll upa = a * down / b;
                ll upcbei = down / d;
                // c 最少是 upcbei 的倍数
                ll ext = e - upa % e;
                // upcbei * x % e == ext
                // 找个 x
                ll x, y;
                ll g = extend_gcd(upcbei, e, x, y);
                if(ext % g == 0) {
                    x *= ext / g;
                    int dx = e / g;
                    x = (x % dx + dx) % dx;
                    if(x == 0) x += dx;
                    ll c = x * upcbei;
                    if(c < ansc) {
                        ansc = c;
                        ansd = d;
                    }
                }
            }
        }
    }

    cout << ansc << " " << ansd << '\n';
}
int main() {
    ios::sync_with_stdio(0), cin.tie(0) ;
    int t = 1;
    cin >> t;
    while(t--) {
        solve();
    }
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

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

input:

4
1 2
2 3
3 7
19 79

output:

0 1
1 3
1 4375
3 316

result:

ok 4 case(s)

Test #2:

score: -100
Wrong Answer
time: 238ms
memory: 3860kb

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 1
1 61000
0 1
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 1
1 29687500
0 1
1 415
1 235937500
1 765000000
1 181250
1 2968750
1 4406250
3 24800
1 480
3 34...

result:

wrong answer The result is not terminating.(Testcase 321)