QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#815842 | #9622. 有限小数 | TosakaUCW | AC ✓ | 26ms | 3636kb | C++20 | 2.5kb | 2024-12-15 17:45:47 | 2024-12-15 17:45:47 |
Judging History
answer
#include <bits/stdc++.h>
using i64 = long long;
#define int i64
#define pb push_back
#define ep emplace
#define eb emplace_back
using std::cerr;
// using namespace std::views;
// using namespace std::ranges;
using std::max, std::min, std::swap, std::array;
using std::cin, std::cout, std::string, std::vector;
using std::ostream;
template <class T1, class T2> ostream &operator<<(ostream &os, const std::pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << '\n'; }
using pii = std::pair<int, int>;
#define fi first
#define se second
double eps = 1e-7;
const int inf = 1e9;
int p2[31], p5[14];
i64 exgcd(i64 a, i64 b, i64 &x, i64 &y) {
if (b == 0) {
x = 1;
y = 0;
return a;
}
i64 g = exgcd(b, a % b, y, x);
y -= a / b * x;
return g;
}
// ax + b = 0 (mod m)
std::pair<i64, i64> sol(i64 a, i64 b, i64 m) {
assert(m > 0);
b *= -1;
i64 x, y;
i64 g = exgcd(a, m, x, y);
if (g < 0) {
g *= -1;
x *= -1;
y *= -1;
}
if (b % g != 0) {
return {-1, -1};
}
x = x * (b / g) % (m / g);
if (x < 0) {
x += m / g;
}
return {x, m / g};
}
void solve() {
int a; cin >> a;
int b; cin >> b;
int p = b;
while (p % 2 == 0) p /= 2;
while (p % 5 == 0) p /= 5;
pii ans = {inf, inf};
for (int x : p2) {
for (int y : p5) {
int d = p * x * y;
if (d > 1'000'000'000) break;
auto [c, g] = sol(b, a * d, p * p);
// cout << c << ' ' << d << '\n';
// cout << "bc = " << b * c;
// cout << ", ad = " << a * d;
// cout << ", p = " << p << '-';
// cout << c << ' ' << d << '\n';
if (c < ans.fi) ans = {c, d};
}
}
cout << ans.fi << ' ' << ans.se << '\n';
}
signed main() {
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
p2[0] = p5[0] = 1;
for (int i = 1; i <= 30; i++) p2[i] = p2[i - 1] * 2LL;
for (int i = 1; i <= 13; i++) p5[i] = p5[i - 1] * 5LL;
int T; cin >> T;
for (; T--; solve());
// solve();
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 3604kb
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: 0
Accepted
time: 26ms
memory: 3632kb
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 3 1 828125000 1 15 1 7 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 15 1 289600000 1 455000 1 115000000 1 1265625 0 1 1 1484375 0 1 1 415 1 235937500 1 765000000 1 181250 1 2968750 1 4406250 3 775 1 3 3 347500 1 944...
result:
ok 10000 case(s)
Test #3:
score: 0
Accepted
time: 24ms
memory: 3552kb
input:
10000 649 915 163 504 53 235 130 131 109 475 214 757 719 788 479 975 515 811 367 972 69 221 21 44 53 157 77 398 227 332 38 87 145 976 396 403 203 381 675 838 157 309 827 962 17 56 455 654 377 809 605 908 907 929 631 978 286 621 289 462 425 428 139 144 259 368 623 642 767 773 32 281 235 624 88 871 39...
output:
1 11712 1 630 1 11750000 1 131 1 95 9 19379200 1 15760 7 9750 33 129760000 1 199065600 1 3536 1 1375 2 98125 9 31840 1 66400 1 54375 1 38125000 4 31484375 1 952500000 9 838000 1 98880000 3 4925440 1 4375 1 638671875 17 2071040 1 2270000 4 2903125 7 1222500 31 12420000 17 288750 1 107000 1 45 1 115 1...
result:
ok 10000 case(s)
Test #4:
score: 0
Accepted
time: 8ms
memory: 3636kb
input:
10000 92916 121039 44812 49081 648307 724965 758916 780077 99256 205723 326746 382303 288096 996145 218860 951883 12801 35014 128665 169196 201524 298553 448151 688549 68962 180601 834591 924280 492975 701087 54245 64901 335012 440951 264723 555529 178065 246602 48529 77438 549077 578801 406123 4154...
output:
2207 774649600 2046 766890625 8323 1449930 21161 780077 2139 82289200 2301 244673920 4299 398458 7309 95188300 37 10941875 491 845980000 36866 186595625 94399 688549000 8011 288961600 239 1478848 36723 876358750 1047 8307328 2769 11023775 2725 71107712 2144 385315625 3413 619504000 15679 11576020 92...
result:
ok 10000 case(s)