QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#738204 | #9622. 有限小数 | huo_hua_ya# | WA | 41ms | 3596kb | C++23 | 1.9kb | 2024-11-12 18:15:24 | 2024-11-12 18:15:24 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 1e6 + 7;
const int mod = 998244353;
#define ll __int128
//求解不定方程 ax+by=c 的最小整数解 x
ll exgcd(ll a, ll b, ll &x, ll &y)
{
if (!b) {x = 1; y = 0; return a;}
ll ans = exgcd(b, a % b, x, y);
ll tmp = x;
x = y;
y = tmp - a / b * y;
return ans;
}
ll x, y;
ll cal(ll a, ll b, ll c)
{
x = 0, y = 0;
ll ans = exgcd(a, b, x, y);
if (c % ans)return -1;
x *= c / ans; b /= ans;
// y *= c / ans;
if (b < 0)b = -b;
x = (x % b + b) % b;
return x;
}
void solve()
{
int a, b;
cin >> a >> b;
ll p = b, k = 1;
while (p % 2 == 0) {
p /= 2;
k *= 2;
}
while (p % 5 == 0) {
p /= 5;
k *= 5;
}
long long c = 1e18, t = 0;
auto check = [&](ll id) {
c = cal(k, a, p * id);
// cerr << c << "\n";
t = (id*p - 1ll * k * c) / a;
return t >= 1;
};
ll l = 1, r = 1e18;
while (l < r) {
ll mid = (l + r) >> 1;
if (check(mid))r = mid;
else l = mid + 1;
}
check(l);
auto cb = [&](ll x){
while(x % 2 == 0){
x /= 2;
}
while(x%5==0){
x/=5;
}
return x == 1;
};
// cerr << (long long)l <<'\n';
long long cx = c,tx = t;
for(int i=1;i<=3000;i++){
check(l + i);
// cerr <<c <<'\n';
if(c == 0)break;
if(cx > c&&t*p<=1000000000ll&&cb(t)){
cx = c;
tx = t;
}
}
long long d = tx * p;
// assert(d<=1000000000&&d>=1);
cout << cx << " " << d << '\n';
}
int main()
{
ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int _ = 1;
cin >> _;
while (_--)
solve();
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 3596kb
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: 41ms
memory: 3568kb
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 6 5 1696 0 27 1 7 11 1810 23 3936 5 282 5 326 5 3168 1 32 11 1220 7 6 13 183 5 129 2 37 7 1430 7 1325 2 69 1 31 1 6 9 362 15 1456 4 115 7 2025 0 1 1 304 2 20 1 415 22 151 16 19125 2 29 3 19 8 987 3 775 1 75 19 11120 1 944 3 109 1 76 2 215 1 6336 0 228 19 1113 3 1 3 1022 2 136 1 22912 3 65 2 119 7 ...
result:
wrong answer Jury found better answer than participant's 1 < 2 (Testcase 1)