2018 Wannafly summer camp Day3--Knight

Knight
题目描述：
有一张无限大的棋盘，你要将马从((0,0))移到((n,m))。
每一步中，如果马在((x,y)(x,y)),你可以将它移动到 ((x+1,y+2)(x+1,y+2)),
((x+1,y-2)(x+1,y−2)),((x-1,y+2)(x−1,y+2)),((x-1,y-2)(x−1,y−2)),
((x+2,y+1)(x+2,y+1)),((x+2,y-1(x+2,y−1)),((x-2,y+1)(x−2,y+1)或(x-2,y-1)(x−2,y−1))。
你需要最小化移动步数。
输入：
第一行一个整数tt表示数据组数 ((1leq tleq 1000))。

每组数据一行两个整数(n,m (|n|,|m| leq 10^9))。

输出：
每组数据输出一行一个整数表示最小步数。
样例输入
2
0 4
4 2
样例输出
2
2

• 由于数据有(10^9)，所以BFS被毙了(～￣▽￣)～，没想到什么好的做法，所以BFS打表找规律=￣ω￣=。

• 打表结果及代码
  

#include<iostream>
#include <queue>
using namespace std;
int dir[8][2] = {
	{1,2},{1,-2},{-1,2},{-1,-2},
	{2,1},{2,-1},{-2,1},{-2,-1}
};
int n, m;
int maze[1100][1100];
bool vis[1100][1100];
struct Point {
	int x, y, step;
	Point(int _x, int _y, int _step) :
		x(_x), y(_y), step(_step) {}
};
void bfs(int sx, int sy)
{
	queue<Point>q;
	q.push(Point(sx, sy, 0));
	vis[sx][sy] = 1;
	maze[sx][sy] = 0;
	while (!q.empty())
	{
		int x = q.front().x;
		int y = q.front().y;
		int step = q.front().step;
		maze[x][y] = step;
		q.pop();
		for (int i = 0; i < 8; i++)
		{
			int tx = x + dir[i][0];
			int ty = y + dir[i][1];
			if (!vis[tx][ty]&&tx<61&&ty<61&&tx>=0&&ty>=0)
			{
				q.push(Point(tx, ty, step + 1));
				vis[tx][ty] = 1;
			}
		}
	}
}
int main() {
	//freopen("1.txt", "w", stdout);
	bfs(30, 30);
	for (int i = 0; i < 60; i++) {
		for (int j = 0; j <60; j++) {
			cout << maze[i][j] << " ";
		}
		cout << endl;
	}
	return 0;
}

• 从上面看，很明显是有规律的，据说大佬能一眼就看出来，以前我是不信的，直到现场有dalao花了4分钟拿了一血……<@_@>蒟蒻只能慢慢推了。首先先把上面的数据放到Excel里面，先预处理一下，将每个答案作为点，以起点为原点建立平面直角坐标系，结果如下：
  

之前我犯了一个错误，BFS起点放到数组边界上去了，应该放到偏中心的位置，把表打出来。将答案统一起来看，从2开始，所有相同的答案围成了一个八边形，这个八边形与坐标轴平行的边都是4层，不平行的都是3层，同时答案基本是向外递增的这样看的时候会发现两个特殊的地方，一个是((0,1),(1,0),(-1,0),(0,-1))这四个点为3，((2,0),(0,2),(0,-2),(-2,0))着四个点4，所以将这些点加入特判。
不难看出，这个表关于坐标轴对称(图中蓝色线)，同时也关于(y=+-x)对称(图中橙色线)，所以(x)轴正半轴为起点，逆时针划分为8个区域，每个区域都一样，只需要考虑1号区域就行了。
现在考虑的为1号区域，希望找到递增的答案之间存在的关系，这个关系为(y=x/2),可以发现这条直线上的整点正好是答案的递增：(0，0）->(2,1)->(4,2).....->(x,floor(x/2)))。将这条直线画出来。（floor()是对一个数值向下取整）
现在看(y=x/2)下方的点，满足关系(y<x/2),也就是(y<x-y)(精度问题，计算时应该用double)，而且下方的点都是在刚才所说的八边形的4层边上，所以可以发现将这些点作如下变换后可以将横坐标和(y=x/2)对应:

double(x-y-y)/4.0*2;

最后将上面这个值取反(+x-y)就是答案。同理可以推出(y=x/2)上方的点，满足关系(y>x/2),在刚才所说的八边形的3层边上，最后推出

double(x-y-y)/3.0*2;

• Code

#include <cstdio>
#include <cstring>
#include <cmath>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
typedef long long ll;
ll fun(ll x, ll y) {

	if (x == 1 && y == 0) {
		return 3;
	}
	if (x == 2 && y == 2) {
		return 4;
	}
	ll delta = x - y;
	if (y>delta) {
		return delta - 2 * floor(((double)(delta-y)) / 3.0);
	}
	else {
		return delta - 2 * floor(((double)(delta-y)) / 4.0);
	}
}

int main()
{

	int t;
	cin >> t;
	while (t--)
	{
		ll x, y;
		cin >> x >> y;
		x = abs(x);
		y = abs(y);
		if (x < y) {
			swap(x, y);
		}
		cout << fun(x, y) << endl;
	}

	return 0;
}

• 最后，为正经题解
  Knight：
  不妨假设(x>=y>=0)。
  当(x<=2y) 时，定义每一步的冗余值(w_i=3-dx-dy)，那么(Σw_i=Σ(2-dx)=3*步数-x-y)，显然我们只需要最小化冗余值。我们先只用(+2,+1)(若x 为奇数则加一步(+1,+2))走到(x,y’)，然后通过将(+2,+1)替换为2 个(+1,+2)使得(0<=y-y’<3)。
  若(y-y’=0)，则冗余值为0，显然最小。
  若(y-y’=1)，则将(+1,+2)替换为(+2,+1)和(-1,+2)或将2 个(+2,+1)替换为(+1,+2),(+1,+2),(+2,-1)，冗余值为2，显然最小。（此处需要特判(2,2)）
  若(y-y’=2)，则加上((+2,+1)和(-2,+1))，冗余值为4，由于不存在冗余值为1的步，所以最小。
  当(x>2y) 时，定义每一步的冗余值(w_i=2-dx)，那么(Σw_i=Σ(2-dx)=2*步数-x)，显然我们只需要最小化冗余值。我们先只使用(+2,+1)走到(2y,y)，然后用
  (+2,+1)和(+2,-1)走到((x’,y)使得0<=x-x’<4)。
  若(x-x’=0)则冗余值为0，显然最小。
  若(x-x’=1) 则将之前的(+2,+1)改为(+1,+2)和(+2,-1)，冗余值为1，显然最
  小。（此处需要特判(1,0)）若(x-x’=2) 则加上(+1,+2)和(+1,-2)，冗余值为2，由x/2+y 的奇偶性可知
  最小。
  若(x-x’=3) 则加上(+2,+1),(+2,+1),(-1,-2)，冗余值为3，由x/2+y 的奇偶性可知最小。
  时间复杂度O(t)