QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#748956 | #9622. 有限小数 | snow_miku | WA | 167ms | 3648kb | C++20 | 2.5kb | 2024-11-14 22:08:48 | 2024-11-14 22:08:56 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll mod = 998244353;
const ll INF = 1e9 + 7;
const int N = 35;
ll gcd(ll a, ll b) {
if (!b) return a;
return gcd(b, a % b);
}
ll exgcd(ll a, ll b, ll &x, ll &y) {
if (!b) {
x = 1;
y = 0;
return a;
}
ll d = exgcd(b, a % b, x, y);
ll t = x;
x = y;
y = t - (a / b) * y;
return d;
}
struct info {
ll p2 = 1, p5 = 1, p = 1;
};
ll pow2[N], pow5[N];
vector<pair<ll, ll>> p25;
void init() {
pow2[0] = 1;
pow5[0] = 1;
for (int i = 1;i <= 30;i++) {
pow2[i] = pow2[i - 1] * 2LL;
}
for (int i = 1;i <= 15;i++) {
pow5[i] = pow5[i - 1] * 5LL;
}
for (int i = 0;i <= 30;i++) {
if (pow2[i] <= INF) {
for (int j = 0;j <= 15;j++) {
if (pow5[j] * pow2[i] <= INF) p25.push_back(make_pair(pow2[i], pow5[j]));
else break;
}
}
else break;
}
}
info check(ll x) {
info res;
while (x % 2 == 0) {
res.p2 *= 2LL;
x /= 2;
}
while (x % 5 == 0) {
res.p5 *= 5LL;
x /= 5;
}
res.p = x;
return res;
}
void solve()
{
ll a, b; cin >> a >> b;
a %= b;
info T = check(b);
ll P = T.p;
// ll x = T.p2 * T.p5;
// cout << P << " " << T.p2 << " " << T.p5 << '\n';
if (P == 1) {
cout << 0 << " " << 1 << '\n';
return;
}
// vector<ll> fct;
// for (int i = 3;i * i <= P;i++) {
// if (P % i == 0) {
// fct.push_back(i);
// if (P % i != i) fct.push_back(P / i);
// }
// }
// for (int q : fct) {
// }
ll c = INF;
ll d = 1;
for (auto [t2, t5] : p25) {
// if (t2 < T.p2 || t5 < T.p5) continue;
ll x = 0, k = 0;
exgcd(P, -T.p2 * T.p5, k, x);
// cout << x << " " << k << '\n';
x *= abs(t2 * t5 * a - P);
ll y = t2 * t5 * T.p2 * T.p5 * P;
ll g = gcd(x, y);
x /= g, y /= g;
if (c > x) {
c = x, d = y;
// cout << t2 << " " << t5 << " " << P << '\n';
}
}
cout << c << " " << d << '\n';
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0);
init();
// for (auto [x, y] : p25) {
// cout << x << " " << y << '\n';
// }
int t; cin >> t;
while(t--) solve();
return 0;
};
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 1ms
memory: 3648kb
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: 167ms
memory: 3572kb
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:
2 3 3 106 7 30 3 7 11 1810 37 246 5 282 5 326 28 99 0 1 5 122 0 1 3 244 21 344 36 37 15 143 13 53 5 92 1 31 1 3 9 362 18 91 3 46 10 81 0 1 45 76 0 1 17 166 22 151 23 306 28 29 3 19 11 1128 128 31 5 12 3 278 4 295 3 109 259 76 9 43 14 495 7 2850 1 1272 0 1 7 1168 71 17 7 179 49 13 18 17 2 215 36 13 1...
result:
wrong answer The result is not terminating.(Testcase 1)