QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#137947#6725. Pick UpYarema#WA 296ms6480kbC++171.8kb2023-08-10 19:32:282023-08-10 19:32:29

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-08-10 19:32:29]
  • 评测
  • 测评结果:WA
  • 用时:296ms
  • 内存:6480kb
  • [2023-08-10 19:32:28]
  • 提交

answer


#include <bits/stdc++.h>
using namespace std;

#define SZ(a) (int)a.size()
#define ALL(a) a.begin(), a.end()
#define FOR(i, a, b) for (int i = (a); i<(b); ++i)
#define RFOR(i, b, a) for (int i = (b)-1; i>=(a); --i)
#define MP make_pair
#define PB push_back
#define F first
#define S second
#define FILL(a, b) memset(a, b, sizeof(a))

typedef long long LL;
typedef pair<int, int> PII;
typedef vector<int> VI;

struct point
{
	int x, y;
	
	point() {}
	point(int a, int b): x(a), y(b) {}
	
	LL dis(point p)
	{
		return LL(abs(x - p.x)) + abs(y - p.y);
	}
};

point A, B, C;
int dx, dy;
int p = 0;
const int N = 1 << 18;
LL ans[N];
int a[N], b[N];

LL f(point X)
{
	LL d1 = A.dis(X) * b[p];
	LL d2 = B.dis(X) * a[p];
	LL d3 = C.dis(X) * a[p];
	return max(d1, d2) + d3;
}

LL g(int st)
{
	int cx = (C.x - A.x);
	if (dx) cx /= dx;
	int cy = (C.y - A.y);
	if (dy) cy /= dy;
	
	point pt(A.x, A.y);
	pt.x += cx * min(st, dx);
	pt.y += cy * min(max(0, st - dx), dy);
	LL a1 = f(pt);
	
	pt = point(A.x, A.y);
	pt.y += cy * min(st, dy);
	pt.x += cx * min(max(0, st - dy), dx);
	LL res = min(a1, f(pt));
	ans[p] = min(res, ans[p]);
	return res;
}


void solve()
{
	cin >> a[p] >> b[p];
	cin >> A.x >> A.y;
	cin >> B.x >> B.y;
	cin >> C.x >> C.y;
	ans[p] = A.dis(C) * b[p];
	
	dx = abs(A.x - C.x);
	dy = abs(A.y - C.y);
	
	LL lx = 0, rx = dx + dy;
	while (lx + 47 < rx)
	{
		LL d = (rx - lx) / 3;
		LL m1 = lx + d;
		LL m2 = rx - d;
		if (g(m1) < g(m2))
			rx = m2;
		else
			lx = m1;
	}
	FOR (i, lx, rx + 1)
		g(i);
	p++;
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	
	int t;
	cin >> t;
	//scanf("%d", &t);
	FOR (tt, 0, t)
	{
		solve();
	}
	cout << fixed << setprecision(6);
	FOR (i, 0, t) cout << double(ans[i]) / (a[i] * b[i]) << '\n';
	
	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

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

input:

3
1 2
0 2 1 0 2 2
1 3
1 1 0 1 3 1
1 2
0 0 100 100 1 1

output:

1.500000
1.000000
2.000000

result:

ok 3 numbers

Test #2:

score: -100
Wrong Answer
time: 296ms
memory: 6480kb

input:

100000
158260522 877914575
602436426 24979445 861648772 623690081 433933447 476190629
731963982 822804784
450968417 430302156 982631932 161735902 880895728 923078537
317063340 492686568
773361868 125660016 650287940 839296263 462224593 492601449
576823355 782177068
404011431 818008580 954291757 1604...

output:

-192892619.817197
357225688.693232
458102224.121226
179827368.403294
116356303.767098
332260084.022387
-1050979906.589531
-627918246.871372
-569142933.760094
668430178.066735
436914678.144082
-1240483100.711399
243210140.395646
-239515133.426356
195435029.774080
20467141.159575
-3980017391.055724
-9...

result:

wrong answer 1st numbers differ - expected: '1.3846072', found: '-192892619.8171970', error = '139312156.1372279'