洛谷 P3964 松鼠聚会

题目链接

参考链接

题解：

这道题目的距离其实就是切比雪夫距离，又称棋盘距离，也就是$D=max(|x_2-x_1|,|y_2-y_1|)$。与之有一定联系的距离是曼哈顿距离，假设$A(x_1,y_1),B(x_2,y_2)$，那么$|AB|=|x_2-x_1|+|y_2-y_1|$。

我们将曼哈顿距离化开：

D = ∣ x 2 − x 1 ∣ + ∣ y 2 − y 1 ∣ = m a x ( x 2 − x 1 + y 2 − y 1 , x 1 − x 2 + y 2 − y 1 , x 1 − x 2 + y 1 − y 2 , x 2 − x 1 + y 1 − y 2 ) = m a x { ∣ ( x 2 + y 2 ) − ( x 1 + y 1 ) ∣ ,   ∣ ( x 2 − y 2 ) − ( x 1 − y 1 ) ∣ } D=|x_2-x_1|+|y_2-y_1|=max(x_2-x_1+y_2-y_1,x_1-x_2+y_2-y_1,x_1-x_2+y_1-y_2,x_2-x_1+y_1-y_2)=max\{|(x_2+y_2)-(x_1+y_1)|,\ |(x_2-y_2)-(x_1-y_1)|\} D=∣x2​−x1​∣+∣y2​−y1​∣=max(x2​−x1​+y2​−y1​,x1​−x2​+y2​−y1​,x1​−x2​+y1​−y2​,x2​−x1​+y1​−y2​)=max{∣(x2​+y2​)−(x1​+y1​)∣, ∣(x2​−y2​)−(x1​−y1​)∣}

上面这个形式其实和切比雪夫距离$D_{chess}=max(|x_2-x_1|,|y_2-y_1|)$形式上非常像，我们可以发现，曼哈顿意义下的$(x_1,y_1)(x_2,y_2)$两点的距离转化为切比雪夫意义下的$(x_2+y_2,x_2-y_2)(x_1+y_1,x_1-y_1)$两点的距离。

我们可以推广一下，给出曼哈顿意义下的坐标$(x,y)$，可以转化成切比雪夫意义下的坐标$(x+y,x-y)$。

反过来，给了切比雪夫意义下的坐标$(x,y)$，可以转化成曼哈顿意义下的坐标$(\frac{x+y}{2},\frac{x-y}{2})$。

有了上面这个转化怎么用呢？

我们现在假设$dis(i,j)$为曼哈顿意义下$i$到$j$的曼哈顿距离。如果当前的终点是$j$，那么有$ans={\sum }_{i=1}^{n}dis(i,j)$，这个要是求的话，复杂度是$O(n^2)$。我们化简一下：

${\sum }_{i=1}^{n}dis(i,j)=dis(1,j)+dis(2,j)+\cdots +dis(n,j)$

我们现在用$dis(i,j)$中的一部分$|x_i-x_j|$举例化简：

${\sum }_{i=1}^n\Delta x=|x_1-x_j|+|x_2-x_j|+\cdots +|x_j-x_j|+|x_{j+1}-x_j|+\cdots +|x_n-x_j|$

我们将横坐标排序成递增，容易发现$|x_j-x_j|$之前的式子都可以拆掉绝对值化简；$|x_{j+1}-x_j|$及以后的式子可以不用动，这样上面的式子就化简成：

$={\sum }_{i=1}^j(x_j-x_i)+{\sum }_{i=j+1}^n(x_i-x_j)$

这个式子可以发现就是前缀和。$dis(i,j)$中的$\Delta y$可以同理化简。

因此这道题目的整体思路就是：

读入切比雪夫意义下的坐标并化为曼哈顿坐标，坐标可以先不除以2，到最后将答案除以2也可以。

将$x,y$分别排序为递增，并求出前缀和。

然后枚举终点，使用二分将$x_i,y_i$在递增序列的位置找出，根据上面的式子，即可求出答案。

实现细节见代码：

#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
const int MAXN = 1e6 + 10;

int x[MAXN], y[MAXN], gx[MAXN], gy[MAXN], sumx[MAXN], sumy[MAXN], n;
int solve(int now) {
	int rx = lower_bound(gx + 1, gx + 1 + n, x[now]) - gx;
	int ry = lower_bound(gy + 1, gy + 1 + n, y[now]) - gy;
	int ans = rx * x[now] - sumx[rx] + sumx[n] - sumx[rx] - (n - rx) * x[now];
	ans += ry * y[now] - sumy[ry] + sumy[n] - sumy[ry] - (n - ry) * y[now];
	return ans;
} 
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0), cout.tie(0);
	cin >> n;
	for (int i = 1; i <= n; i++) {
		int x1, y1;
		cin >> x1 >> y1;
		x[i] = gx[i] = x1 + y1; // 这里不÷2，避免精度问题
		y[i] = gy[i] = x1 - y1;
	}
	sort(gx + 1, gx + 1 + n);
	sort(gy + 1, gy + 1 + n);
	for (int i = 1; i <= n; i++) {
		sumx[i] = sumx[i - 1] + gx[i];
		sumy[i] = sumy[i - 1] + gy[i];
	}
	int ans = 1e16;
	for (int i = 1; i <= n; i++) {
		ans = min(ans, solve(i));
	}
	cout << ans / 2 << endl;
	return 0;
}