QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#744159#9622. 有限小数RosmontispesWA 64ms3792kbC++202.1kb2024-11-13 21:01:472024-11-13 21:01:53

Judging History

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

  • [2024-11-13 21:01:53]
  • 评测
  • 测评结果:WA
  • 用时:64ms
  • 内存:3792kb
  • [2024-11-13 21:01:47]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
const int inf = 1e9;
long long exgcd(long long a,long long b,long long &x,long long &y)
{
    if(!b)
    {
        x = 1,y = 0;
        return a;
    }
    else
    {
        long long g = exgcd(b,a % b,y,x);
        y -= a / b * x;
        return g;
    }
}



void solve()
{
    long long a,b;
    cin>>a>>b;
    long long w = 1;
    long long p = b;
    while(p % 2 == 0)
        p /= 2,w *= 2;
    while(p % 5 == 0)
        p /= 5,w *= 5;
    if(p == 1){
        cout<<0<<" "<<1<<"\n";
        return;
    }
    long long x,y;
    long long g = exgcd(p,a,x,y);
    long long c = g;
    long long d = c / g;
    //x大于等于0的最小解,在这基础上x + b/g,y - a/g为一个新解
    x = ((__int128)(x * d) % (a / g) + a / g) % (a / g);
    y = (c - p * x) / a;
    if(y >= 0){
        long long k = y / (p / g);
        y -= k * (p / g);
        x += k * (a / g);
        if(y >= 0){
            y -= p / g;
            x += a / g;
        }
    }
    y = -y;

    // cout<<p<<" "<<x<<" "<<a<<" "<<y<<" "<<p / g<<" "<<a / g<<"\n";
    // cout<<3 * g<<" "<<b * (-y)<<"\n";
    long long ans = 1e9,bot = 0;
    for(long long i = 1;i <= inf;i *= 5){
        for(long long j = 1;i * j <= inf;j *= 2){
            long long cc = i * j;

            long long aa = y,bb = p / g;
            long long xx,yy;
            long long gg = exgcd(aa,bb,xx,yy);
            //y + k * p / g = cc
            if(cc % gg != 0)
                continue;
            long long dd = cc / gg;
            xx = ((__int128)(xx * dd) % (bb / gg) + bb / gg) % (bb / gg);
            yy = (cc - aa * xx) / bb;
            // if(cc == 4){
                long long ansc = (xx * x + (a / g) * yy) * p - cc * a;
                if(ansc < ans){
                    ans = ansc;
                    bot = b * cc;
                }
            // }
        }
    }
    cout<<ans<<" "<<bot<<"\n";

}
int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    
    int T;
    std::cin >> 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: 3748kb

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: 64ms
memory: 3792kb

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 350
1 231680000
23 3936
1 36096000
3 2086400000
1 63360
0 1
1 31232
0 1
1 4093640704
1 44032
1 757760
1 11714560
1 14227079168
1 2944
1 31
1 6
1 289600000
1 455000
1 58880
1 51840
0 1
1 304
0 1
1 11140071424
1 19328000
1 20054016000
1 3892314112
1 608
1 72192
3 24800
1 192
1 7462...

result:

wrong answer Integer 2086400000 violates the range [1, 10^9]