Submission #20140

# Submission time Handle Problem Language Result Execution time Memory
20140 2016-02-27T16:13:00 Z hongjun7 Min-cost GCD (GA9_mcg) C++14
57 / 100
1000 ms 1224 KB
#include <cstdio>
#include <vector>
#include <map>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
static char _buffer[1024];
static int _currentChar = 0;
static int _charsNumber = 0;
static inline int _read() {
	if (_charsNumber < 0) {
		exit(1);
	}
	if (!_charsNumber || _currentChar == _charsNumber) {
		_charsNumber = (int)fread(_buffer, sizeof(_buffer[0]), sizeof(_buffer), stdin);
		_currentChar = 0;
	}
	if (_charsNumber <= 0) {
		return -1;
	}
	return _buffer[_currentChar++];
}
static inline ll _readll() {
	int c, s; ll x;
	c = _read();
	while (c <= 32) c = _read();
	x = 0;
	s = 1;
	if (c == '-') {
		s = -1;
		c = _read();
	}
	while (c > 32) {
		x *= 10;
		x += c - '0';
		c = _read();
	}
	if (s < 0) x = -x;
	return x;
}
ll a, b, p, q;
map <pair<ll, ll>, ll> d;
map <pair<int, int>, ll> h;
int lim = 2e9; ll lim2 = 1e18;
inline ll g(int x, int y) {
	if (x == 0 || y == 0) return 0;
	if (h.count({ x, y })) return h[{x, y}];
	ll d1 = g(y, x%y) + p, d2 = lim2; int m;
	if (x >= y) {
		m = x / y;
		if (m <= d1 / q) d2 = g(x - m*y, y) + m*q;
	}
	else {
		m = y / x;
		if (m <= d1 / q) d2 = g(x, y - m*x) + m*q;
	}
	if (d1 > d2) d1 = d2;
	h[{x, y}] = d1;
	return d1;
}
inline ll f(ll x, ll y) {
	if (x == 0 || y == 0) return 0;
	if (x > lim || y > lim) {
		if (d.count({ x, y })) return d[{x, y}];
		ll d1 = f(y, x%y) + p, d2 = lim2, m;
		if (x >= y) {
			m = x / y;
			if (m <= d1 / q) d2 = f(x - m*y, y) + m*q;
		}
		else {
			m = y / x;
			if (m <= d1 / q) d2 = f(x, y - m*x) + m*q;
		}
		if (d1 > d2) d1 = d2;
		d[{x, y}] = d1;
		return d1;
	}
	return g((int)(x), (int)(y));
}
int main() {
	int T;
	for (scanf("%d", &T); T--; ) {
		a = _readll(); b = _readll(); p = _readll(); q = _readll();
		d.clear(); h.clear();
		printf("%lld\n", f(a, b));
	}
}
# Verdict Execution time Memory Grader output
1 Correct 0 ms 1224 KB Output is correct
2 Correct 0 ms 1224 KB Output is correct
3 Correct 0 ms 1224 KB Output is correct
4 Correct 0 ms 1224 KB Output is correct
5 Correct 0 ms 1224 KB Output is correct
6 Correct 0 ms 1224 KB Output is correct
7 Correct 0 ms 1224 KB Output is correct
8 Correct 0 ms 1224 KB Output is correct
9 Correct 0 ms 1224 KB Output is correct
10 Correct 0 ms 1224 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 81 ms 1224 KB Output is correct
2 Correct 149 ms 1224 KB Output is correct
3 Correct 157 ms 1224 KB Output is correct
4 Correct 172 ms 1224 KB Output is correct
5 Correct 80 ms 1224 KB Output is correct
6 Correct 155 ms 1224 KB Output is correct
7 Correct 164 ms 1224 KB Output is correct
8 Correct 95 ms 1224 KB Output is correct
9 Correct 151 ms 1224 KB Output is correct
10 Correct 115 ms 1224 KB Output is correct
11 Correct 145 ms 1224 KB Output is correct
12 Correct 84 ms 1224 KB Output is correct
13 Correct 138 ms 1224 KB Output is correct
14 Correct 2 ms 1224 KB Output is correct
15 Correct 0 ms 1224 KB Output is correct
16 Correct 0 ms 1224 KB Output is correct
17 Correct 1 ms 1224 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 4 ms 1224 KB Output is correct
2 Correct 3 ms 1224 KB Output is correct
3 Correct 4 ms 1224 KB Output is correct
4 Correct 4 ms 1224 KB Output is correct
5 Correct 4 ms 1224 KB Output is correct
6 Correct 3 ms 1224 KB Output is correct
7 Correct 5 ms 1224 KB Output is correct
8 Correct 4 ms 1224 KB Output is correct
9 Correct 3 ms 1224 KB Output is correct
10 Correct 4 ms 1224 KB Output is correct
11 Correct 0 ms 1224 KB Output is correct
12 Correct 3 ms 1224 KB Output is correct
13 Correct 4 ms 1224 KB Output is correct
14 Correct 0 ms 1224 KB Output is correct
15 Correct 0 ms 1224 KB Output is correct
16 Correct 0 ms 1224 KB Output is correct
17 Correct 0 ms 1224 KB Output is correct
18 Correct 0 ms 1224 KB Output is correct
19 Correct 5 ms 1224 KB Output is correct
20 Correct 6 ms 1224 KB Output is correct
21 Correct 6 ms 1224 KB Output is correct
22 Correct 3 ms 1224 KB Output is correct
23 Correct 5 ms 1224 KB Output is correct
24 Correct 6 ms 1224 KB Output is correct
25 Correct 7 ms 1224 KB Output is correct
26 Correct 6 ms 1224 KB Output is correct
27 Correct 6 ms 1224 KB Output is correct
28 Correct 4 ms 1224 KB Output is correct
29 Correct 4 ms 1224 KB Output is correct
30 Correct 5 ms 1224 KB Output is correct
31 Correct 5 ms 1224 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 808 ms 1224 KB Output is correct
2 Correct 791 ms 1224 KB Output is correct
3 Correct 804 ms 1224 KB Output is correct
4 Correct 783 ms 1224 KB Output is correct
5 Correct 781 ms 1224 KB Output is correct
6 Correct 788 ms 1224 KB Output is correct
7 Correct 792 ms 1224 KB Output is correct
8 Correct 798 ms 1224 KB Output is correct
9 Correct 811 ms 1224 KB Output is correct
10 Correct 791 ms 1224 KB Output is correct
11 Correct 793 ms 1224 KB Output is correct
12 Correct 823 ms 1224 KB Output is correct
13 Correct 806 ms 1224 KB Output is correct
14 Execution timed out 1000 ms 1224 KB Program timed out
15 Halted 0 ms 0 KB -