@codeforces

目录

• @description@
• @solution@
• @accepted code@
• @details@

---

@description@

给定若干个三维空间的点 (xi, yi, zi)，求一个坐标都为整数的点 P，使得 P 到这些点的最大曼哈顿距离最小。

原题传送门。

@solution@

显然三分套三分套三分。

看到最大值，把绝对值 |x| 拆成 max(x, -x)。接着二分最大距离 d，则 max(...) ≤ d。

因此得到如下不等式组：

[egin{cases} l_1 leq x + y + z leq r_1 \ l_2 leq x + y - z leq r_2 \ l_3 leq x - y + z leq r_3 \ l_4 leq - x + y + z leq r_4 \ end{cases} ]

仿照二维情况将曼哈顿距离转切比雪夫距离的方式，作代换 (a = x + y - z, b = x - y + z, c = - x + y + z)。

则有：(x = frac{a + b}{2}, y = frac{a + c}{2}, z = frac{b + c}{2}, x + y + z = a + b + c)。

当 (x, y, z) 都是整数时，(a, b, c) 同奇同偶。不妨先枚举奇偶性，则可把原不等式变形为如下形式：

[egin{cases} l_1' leq a' + b' + c' leq r_1' \ l_2' leq a' leq r_2' \ l_3' leq b' leq r_3' \ l_4' leq c' leq r_4' \ end{cases} ]

这样做的好处是，我们只留下了一个 (a', b', c') 互相制约的不等式。

剩下的只需要贪心地把 (a', b', c') 先设置为最小值，然后往上调整即可。时间复杂度 (O(nlog A))。

@accepted code@

#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;

typedef long long ll;

const ll INF = ll(3E18);
const int dx[4] = {1, 1, 1, -1};
const int dy[4] = {1, 1, -1, 1};
const int dz[4] = {1, -1, 1, 1};

ll le[4], ri[4];
ll a1, b1, c1;
bool get() {
	for(int i=0;i<4;i++)
		if( le[i] > ri[i] ) return false;
	
	a1 = le[1], b1 = le[2], c1 = le[3];
	if( a1 + b1 + c1 > ri[0] ) return false;
	else {
		if( a1 + b1 + c1 < le[0] ) {
			if( ri[1] + b1 + c1 >= le[0] ) {
				a1 = le[0] - b1 - c1;
				return true;
			} else {
				a1 = ri[1];
				if( a1 + ri[2] + c1 >= le[0] ) {
					b1 = le[0] - a1 - c1;
					return true;
				} else {
					b1 = ri[2];
					if( a1 + b1 + ri[3] >= le[0] ) {
						c1 = le[0] - a1 - b1;
						return true;
					} else return false;
				}
			}
		} else return true;
	}
}

ll lb[4], ub[4];
ll ansx, ansy, ansz;
bool check(ll d) {
	for(int o=0;o<=1;o++) {
		le[0] = (lb[0] - d) - 3*o, ri[0] = (ub[0] + d) - 3*o;
		for(int i=1;i<4;i++) le[i] = (lb[i] - d) - o, ri[i] = (ub[i] + d) - o;
		for(int i=0;i<4;i++) le[i] = ceil((long double)le[i] / 2), ri[i] = floor((long double)ri[i] / 2);
		if( get() ) {
			ll a = 2*a1 + o, b = 2*b1 + o, c = 2*c1 + o;
			ansx = (a + b) / 2, ansy = (a + c) / 2, ansz = (b + c) / 2;
			return true;
		}
	}
	return false;
}
void solve() {
	int n; scanf("%d", &n);
	for(int i=0;i<4;i++) lb[i] = -INF, ub[i] = INF;
	for(int i=1;i<=n;i++) {
		ll x, y, z; scanf("%lld%lld%lld", &x, &y, &z);
		for(int j=0;j<4;j++) {
			lb[j] = max(lb[j], dx[j]*x + dy[j]*y + dz[j]*z);
			ub[j] = min(ub[j], dx[j]*x + dy[j]*y + dz[j]*z);
		}
	}
	
	ll l = 0, r = INF;
	while( l < r ) {
		ll m = (l + r) >> 1;
		if( check(m) ) r = m;
		else l = m + 1;
	}
	check(r); printf("%lld %lld %lld
", ansx, ansy, ansz);
}

int main() {
	int T; scanf("%d", &T);
	while( T-- ) solve();
}

@details@

一开始本来想转类切比雪夫距离结果发现好像二维三维不一样。

然后尝试从立体几何入手想象，发现我完全没学过立几。

果然这是一道数学题啊。数学题好难。