QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#737611 | #9622. 有限小数 | Zxyoul# | WA | 238ms | 3860kb | C++23 | 2.1kb | 2024-11-12 16:25:13 | 2024-11-12 16:25:16 |
Judging History
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();
}
}
詳細信息
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)