【ACWing】2135. 马步距离（配数学证明）

题目地址：

https://www.acwing.com/problem/content/description/2137/

在国际象棋和中国象棋中，马的移动规则相同，都是走“日”字，我们将这种移动方式称为马步移动。如下图所示，从标号为$0$的点出发，可以经过一步马步移动达到标号为$1$的点，经过两步马步移动达到标号为$2$的点。

任给平面上的两点$p$和$s$，它们的坐标分别为$(x_p,y_p)$和$(x_s,y_s)$，其中，$x_p,y_p,x_s,y_s$均为整数。从$(x_p,y_p)$出发经过一步马步移动可以达到$(x_p+1,y_p+2),(x_p+2,y_p+1),(x_p+1,y_p-2),(x_p+2,y_p-1),(x_p-1,y_p+2),(x_p-2,y_p+1),(x_p-1,y_p-2),(x_p-2,y_p-1)$。假设棋盘充分大，并且坐标可以为负数。现在请你求出从点$p$到点$s$至少需要经过多少次马步移动？

输入格式：
只包含$4$个整数，它们彼此用空格隔开，分别为$x_p,y_p,x_s,y_s$。

输出格式：
输出一个整数，表示从点$p$到点$s$至少需要经过的马步移动次数。

数据范围：
$-10^7<x_p,y_p,x_s,y_s<10^7$

由对称性，步数与从$(0,0)$到$(x,y)=(|x_p-x_s|,|y_p-y_s|)$的步数相同，进一步地，若$x<y$，步数与到$(y,x)$也相同。所以问题转化为从$(0,0)$到$(x,y),x\ge y$的最小步数是多少，设步数为$f(x,y)$。我们可以先考虑步数比较少的情形，试图找出规律。我们只关注第一象限中$x\ge y$的区域：

这个区域，又可以进一步细化为$x\ge 2y$的区域和$2y>x\ge y$的区域，前者我们称为$1$号区域，而对于后者，在图中我们可以发现，当$2y>x\ge y$的时候，由于在$x\ge y$的区域内，当$x\ge 5$的时候，有$f(x,y)=f(x+1,y-1)$，所以我们可以从$(x,y)$这个点向右下方走，直到走到$1$号区域为止，所以本质上我们只需要对$1$号区域求解即可。

接下来只考虑$1$号区域。容易发现，当$x=2y$的时候，$f(x,y)=f(2y,y)=y$，并且$f(2y,y-d)=y+d\%2$，即$1$号区域内，当$x$是偶数的时候，有$f(x,y)=x/2+(x/2-y)\%2$；当$x$是奇数的时候，由于$f(x,y)=f(x+1,y+1)$，所以有$f(x,y)=(x+1)/2+((x+1)/2-(y+1))\%2$，可以验证$y=0$的时候公式也成立。此外还有几个特殊情况需要特判，分别是$f(0,0)=0,f(1,0)=3,f(2,2)=4$。后两个特例的产生原因是路径会经过$x\ge y$这个区域之外。

简要证明一下。我们只考虑$x\ge y$的区域。首先暴力解出$f(x,y)\le 4$的区域，这一步可以BFS来做。
先证明不在$1$号区域的情形，这个可以用数学归纳法来考虑，对$f(x,y)$进行归纳，并且假设最短路径可以完全含在$x\ge y$的区域里。对于某个$f(x,y)\ge 5$的格子，设$k=f(x,y)$，那么由归纳假设，对于最少步数是$k-1$的点是成立上面的命题的，那么对于$(x,y)$这个点，其必然是从某个最少步数是$k-1$的点扩展而来，比如是从$(a,b)$扩展而来，找到最右下角的那个$(a,b)$，这样这个$(a,b)$能走到$(x+1,y-1)$，所以$f(x+1,y-1)=f(x,y)$。
再证明$1$号区域的情形，容易证明$x=2y$这条直线上，$f(x,y)=y$。先用数学归纳法证明$x$为偶数的情形，可以类似上面的方法证明出$f(x,y)=f(x,y-2)=f(x,y-4)=...$。接下来我们两列两列考虑，在图中是可以发现规律的，对于$x$是偶数的那一列，从上到下是$f(2y,y),f(2y,y)+1,f(2y,y),f(2y,y)+1,...$，对于$x$是奇数的那一列，最上面的数字可以由$f(x,y)=f(x+1,y+1)$定出，然后向下依次是$f(2y,y),f(2y,y)+1,f(2y,y),f(2y,y)+1,...$，规律一样。接下来可以用数学归纳法证明这个两列的规律。先证明两列两列所需步数是单调上升的，这一点比较容易，接下来可以按照马的走法模拟一遍，即可证明。

代码如下：

#include <iostream>
using namespace std;

inline long f(long x, long y) {
    if (!x && !y) return 0;
    if (x == 1 && y == 0) return 3;
    if (x == 2 && y == 2) return 4;
    if (y <= x - y) {
        if (x % 2 == 0) return x / 2 + (x / 2 - y) % 2;
        else return (x + 1) / 2 + ((x + 1) / 2 - (y + 1)) % 2;
    }
    
    long d = (2 * y - x) / 3;
    if (!d) d = 1;
    x += d, y -= d;
    return f(x, y);
}

int main() {
    long x, y, a, b;
    scanf("%ld%ld%ld%ld", &x, &y, &a, &b);
    long c = abs(x - a), d = abs(y - b);
    if (c < d) swap(c, d);
    printf("%ld\n", abs(f(c, d)));
    return 0;
}

时空复杂度$O(1)$。