QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#784906#9622. 有限小数Nyanya_AC ✓34ms3832kbC++231.3kb2024-11-26 16:17:372024-11-26 16:17:38

Judging History

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

  • [2024-11-26 16:17:38]
  • 评测
  • 测评结果:AC
  • 用时:34ms
  • 内存:3832kb
  • [2024-11-26 16:17:37]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int T;
	cin >> T;
	while (T--)
	{
		long long a, b;
		cin >> a >> b;
		long long z = b;
		while (z % 2 == 0)
			z /= 2;
		while (z % 5 == 0)
			z /= 5;
		vector<long long> p2(30);
		vector<long long> p5(13);
		p2[0] = p5[0] = 1;
		for (int i = 1; i <= 29; ++i)
			p2[i] = p2[i - 1] * 2;
		for (int i = 1; i <= 12; ++i)
			p5[i] = p5[i - 1] * 5;
		auto exgcd = [&](auto self, long long a, long long b, long long &x, long long &y) -> long long
		{
			if (b == 0)
			{
				x = 1;
				y = 0;
				return a;
			}
			long long d = self(self, b, a % b, y, x);
			y -= a / b * x;
			return d;
		};
		long long ansc = 1e9, ansd;
		for (int i = 0; i <= 29; ++i)
			for (int j = 0; j <= 12; ++j)
			{
				if (p2[i] * p5[j] * z > 1e9)
					break;
				long long k = p2[i] * p5[j];
				long long d = k * z;
				long long g = gcd(b, d);
				long long A = b / g;
				long long B = z;
				long long C = -a * (d / g);
				long long c, y;
				long long D = exgcd(exgcd, A, B, c, y);
				c = ((c * (C / D)) % (B / D) + (B / D)) % (B / D);
				if (c < ansc)
					ansc = c, ansd = d;
			}
		cout << ansc << " " << ansd << "\n";
	}
}

/*
4
1 2
2 3
3 7
19 79
*/

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 0ms
memory: 3696kb

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: 34ms
memory: 3488kb

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: 32ms
memory: 3832kb

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: 11ms
memory: 3608kb

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)