【计算几何03】距离计算
感谢OI wiki对本文的帮助
欧氏距离
欧氏距离，一般也称作欧几里得距离。在平面直角坐标系中，设点 \(A,B\) 的坐标分别为 \(A(x_1,y_1),B(x_2,y_2)\) ，则两点间的欧氏距离为：
举个例子，若在平面直角坐标系中，有两点 \(A(6,5),B(2,2)\) ，通过公式，我们很容易得到 \(A,B\) 两点间的欧氏距离：
除此之外， \(P(x,y)\) 到原点的欧氏距离可以用公式表示为：
那么，三维空间中两点的欧氏距离公式呢？我们来观察下图。
我们很容易发现，在 \(\triangle ADC\) 中， \(\angle ADC = 90^\circ\) ；在 \(\triangle ACB\) 中， \(\angle ACB = 90^\circ\) 。
由此可得，三维空间中欧氏距离的距离公式为：
NOIP2017 提高组 奶酪 就运用了这一知识，可以作为欧氏距离的例题。
以此类推，我们就得到了 \(n\) 维空间中欧氏距离的距离公式：对于 \(\vec A(x_{11}, x_{12}, \cdots,x_{1n}) ,~ \vec B(x_{21}, x_{22}, \cdots,x_{2n})\) ，有
欧氏距离虽然很有用，但也有明显的缺点。两个整点计算其欧氏距离时，往往答案是浮点型，会存在一定误差。
曼哈顿距离
在二维空间内，两个点之间的曼哈顿距离（Manhattan distance）为它们横坐标之差的绝对值与纵坐标之差的绝对值之和。设点 \(A(x_1,y_1),B(x_2,y_2)\) ，则 \(A,B\) 之间的曼哈顿距离用公式可以表示为：
观察下图：
在 \(A,B\) 间，黄线、橙线都表示曼哈顿距离，而红线、蓝线表示等价的曼哈顿距离，绿线表示欧氏距离。
同样的栗子，在下图中 \(A,B\) 的坐标分别为 \(A(25,20),B(10,10)\) 。
通过公式，我们很容易得到 \(A,B\) 两点间的曼哈顿距离：
经过推导，我们得到 \(n\) 维空间 的曼哈顿距离公式为：
除了公式之外，曼哈顿距离还具有以下数学性质：
-
非负性
曼哈顿距离是一个非负数。
\(d(i,j)\geq 0\)
-
统一性
点到自身的曼哈顿距离为 \(0\) 。
\(d(i,i) = 0\)
-
对称性
\(A\) 到 \(B\) 与 \(B\) 到 \(A\) 的曼哈顿距离相等，且是对称函数。
\(d(i,j) = d(j,i)\)
-
三角不等式
从点 \(i\) 到 \(j\) 的直接距离不会大于途经的任何其它点 \(k\) 的距离。
\(d(i,j)\leq d(i,k)+d(k,j)\)
了解基本的以后来看一道例题 ↓
例题 1 洛谷P5098
根据题意，对于式子 \(|x_1-x_2|+|y_1-y_2|\) ，我们可以假设 \(x_1 - x_2 \geq 0\) ，根据 \(y_1 - y_2\) 的符号分成两种情况：
-
\((y_1 - y_2 \geq 0)\rightarrow |x_1-x_2|+|y_1-y_2|=x_1 + y_1 - (x_2 + y_2)\)
-
\((y_1 - y_2 \lt 0)\rightarrow |x_1-x_2|+|y_1-y_2|=x_1 - y_1 - (x_2 - y_2)\)
只要分别求出 \(x+y, x-y\) 的最大值和最小值即能得出答案。
参考代码1
#include <bits/stdc++.h>
using namespace std;
template
inline void read(T &x) {
x = 0;
char c = getchar();
bool f = 0;
for (; !isdigit(c); c = getchar()) f ^= c == '-';
for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + (c ^ 48);
x = f ? -x : x;
}
int n, x, y, minx = 0x7fffffff, maxx, miny = 0x7fffffff, maxy;
int main() {
read(n);
for (int i = 1; i <= n; i++) {
read(x), read(y);
minx = min(minx, x + y), maxx = max(maxx, x + y);
miny = min(miny, x - y), maxy = max(maxy, x - y);
}cpp
printf("%d\n", max(maxx - minx, maxy - miny));
return 0;
}
其实还有第二种做法，那就是把曼哈顿距离转化为切比雪夫距离求解，最后部分会讲到。
切比雪夫距离
切比雪夫距离（Chebyshev distance）是向量空间中的一种度量，二个点之间的距离定义是其各坐标数值差绝对值的最大值。——来源： 切比雪夫距离 - 维基百科
在二维空间内，两个点之间的切比雪夫距离为它们横坐标之差的绝对值与纵坐标之差的绝对值的最大值。设点 \(A(x_1,y_1),B(x_2,y_2)\) ，则 \(A,B\) 之间的切比雪夫距离用公式可以表示为：
仍然是这个栗子，下图中 \(A,B\) 的坐标分别为 \(A(6,5),B(2,2)\) 。
\(n\) 维空间中切比雪夫距离的距离公式：
曼哈顿距离与切比雪夫距离的相互转化
首先，我们考虑画出平面直角坐标系上所有到原点的曼哈顿距离为 \(1\) 的点。
通过公式，我们很容易得到方程 \(|x| + |y| = 1\) 。
将绝对值展开，得到 \(4\) 个 一次函数，分别是：
将这 \(4\) 个函数画到平面直角坐标系上，得到一个边长为 \(\sqrt{2}\) 的正方形，如下图所示：
正方形边界上所有的点到原点的 曼哈顿距离 都是 \(1\) 。
同理，我们再考虑画出平面直角坐标系上所有到原点的 切比雪夫距离 为 \(1\) 的点。
通过公式，我们知道 \(\max(|x|,|y|)=1\) 。
我们将式子展开，也同样可以得到可以得到 \(4\) 条 线段，分别是：
画到平面直角坐标系上，可以得到一个边长为 \(2\) 的正方形，如下图所示：
正方形边界上所有的点到原点的切比雪夫距离都是 \(1\) 。
将这两幅图对比，我们会神奇地发现：
这 \(2\) 个正方形是相似图形。
所以，曼哈顿距离与切比雪夫距离之间会不会有联系呢？
接下来我们简略证明一下：
假设 \(A(x_1,y_1),B(x_2,y_2)\) ，
我们把曼哈顿距离中的绝对值拆开，能够得到四个值，这四个值中的最大值是两个非负数之和，即曼哈顿距离。则 \(A,B\) 两点的曼哈顿距离为：
我们很容易发现，这就是 \((x_1 + y_1,x_1 - y_1), (x_2 + y_2,x_2 - y_2)\) 两点之间的切比雪夫距离。
所以将每一个点 \((x,y)\) 转化为 \((x + y, x - y)\) ，新坐标系下的切比雪夫距离即为原坐标系下的曼哈顿距离。
同理， \(A,B\) 两点的切比雪夫距离为：
而这就是 \((\dfrac{x_1 + y_1}{2},\dfrac{x_1 - y_1}{2}), (\dfrac{x_2 + y_2}{2},\dfrac{x_2 - y_2}{2})\) 两点之间的曼哈顿距离。
所以将每一个点 \((x,y)\) 转化为 \((\dfrac{x + y}{2},\dfrac{x - y}{2})\) ，新坐标系下的曼哈顿距离即为原坐标系下的切比雪夫距离。
结论
- 曼哈顿坐标系是通过切比雪夫坐标系旋转 \(45^\circ\) 后，再缩小到原来的一半得到的。
- 将一个点 \((x,y)\) 的坐标变为 \((x + y, x - y)\) 后，
- 原坐标系中的曼哈顿距离等于新坐标系中的切比雪夫距离
- 将一个点 \((x,y)\) 的坐标变为 \((\dfrac{x + y}{2},\dfrac{x - y}{2})\) 后，
原坐标系中的切比雪夫距离等于新坐标系中的曼哈顿距离
碰到求切比雪夫距离或曼哈顿距离的题目时，我们往往可以相互转化来求解。两种距离在不同的题目中有不同的优缺点，应该灵活运用。
例题 2
P4648「IOI2007」pairs 动物对数 （曼哈顿距离转切比雪夫距离）
P3964「TJOI2013」松鼠聚会 （切比雪夫距离转曼哈顿距离）
最后给出 P5098「USACO04OPEN」Cave Cows 3 的第二种解法：
我们考虑将题目所求的曼哈顿距离转化为切比雪夫距离，即把每个点的坐标 \((x,y)\) 变为 \((x + y, x - y)\) 。
所求的答案就变为 \(\max_{i,j\in n}\begin{Bmatrix} \max\begin{Bmatrix} |x_i - x_j|,|y_i - y_j|\end{Bmatrix}\end{Bmatrix}\)
现要使得横坐标之差和纵坐标之差最大，只需要预处理出 \(x,y\) 的最大值和最小值即可。
Code
#include <bits/stdc++.h>
using namespace std;
template
inline void read(T &x) {
x = 0;
char c = getchar();
bool f = 0;
for (; !isdigit(c); c = getchar()) f ^= c == '-';
for (; isdigit(c); c = getchar()) x = (x << 3) + (x << 1) + (c ^ 48);
x = f ? -x : x;
}
int n, x, y, a, b, minx = 0x7fffffff, maxx, miny = 0x7fffffff, maxy;
int main() {
read(n);
for (int i = 1; i <= n; i++) {
read(a), read(b);
x = a + b, y = a - b;
minx = min(minx, x), maxx = max(maxx, x);
miny = min(miny, y), maxy = max(maxy, y);
}
printf("%d\n", max(maxx - minx, maxy - miny));
return 0;
}
对比两份代码，我们又能够发现，两种不同的思路，写出来的代码却是完全等价的，是不是很神奇呢？当然，更高深的东西需要大家另行研究。
\(L_m\) 距离
一般地，我们定义平面上两点 \(A(x_1, y_1)\) , \(B(x_2, y_2)\) 之间的 \(L_m\) 距离为
\(d(L_m) = (|x_1-x_2|^m+|y1-y2|^m)^{\frac{1}{m}}\)
特殊的， \(L_2\) 距离就是欧几里得距离， \(L_1\) 距离就是曼哈顿距离。
汉明距离
汉明距离是两个字符串之间的距离，它表示两个长度相同的字符串对应位字符不同的数量
我们可以简单的认为对两个串进行异或运算，结果为 1 的数量就是两个串的汉明距离。
部分内容搬运自 浅谈三种常见的距离算法 ，感谢作者 xuxing 的授权。
