【洛谷 P8668】[蓝桥杯 2018 省 B] 螺旋折线 题解(数学+平面几何)

107 阅读2分钟

[蓝桥杯 2018 省 B] 螺旋折线

题目描述

如图所示的螺旋折线经过平面上所有整点恰好一次。

对于整点 (X,Y)(X, Y),我们定义它到原点的距离 dis(X,Y)\text{dis}(X, Y) 是从原点到 (X,Y)(X, Y) 的螺旋折线段的长度。

例如 dis(0,1)=3\text{dis}(0, 1)=3dis(2,1)=9\text{dis}(-2, -1)=9

给出整点坐标 (X,Y)(X, Y),你能计算出 dis(X,Y)\text{dis}(X, Y) 吗?

输入格式

XXYY

输出格式

输出 dis(X,Y)\text{dis}(X, Y)

样例 #1

样例输入 #1

0 1

样例输出 #1

3

提示

对于 40%40\%的数据,1000X,Y1000-1000\le X,Y\le 1000

对于 70%70\% 的数据,105X,Y105-10^5\le X,Y \le 10^5

对于 100%100\% 的数据,109X,Y109-10^9\le X,Y \le 10^9


思路

暴力模拟会超时,不妨想想其他办法。

函数d(int t)计算从原点到点(-t, t)的螺旋折线距离,使用等差数列求和公式2t * (2t + 1) * 2 / 2

函数dis(int x, int y)是主要的计算函数,根据输入的坐标点(x, y),计算出螺旋折线距离。

观察到有三条分界线:y=xy = xy=xy=-xy=x+1y=x+1 将图形分为四个部分。可以对每个区域分类讨论,通过从原点到点(-t, t)的螺旋折线距离,换算出从原点到任意点(x, y)的螺旋折线距离。

  • 当点位于左侧区域,并且满足条件x0x \geq 0xy-x \leq yy<xy < x时,此时的螺旋折线距离可以通过公式d(t)+3tyd(t) + 3 * t - y计算,其中t=xt = x

  • 当点位于上侧区域,并且满足条件y>0y > 0y<x-y < xxyx \leq y时,此时的螺旋折线距离可以通过公式d(t)+t+xd(t) + t + x计算,其中t=yt = y

  • 当点位于右侧区域,并且满足条件x<0x < 0x+1<yx + 1 < yyxy \leq -x时,此时的螺旋折线距离可以通过公式d(t)t+yd(t) - t + y计算,其中t=xt = -x

  • 当点位于下侧区域,此时的螺旋折线距离可以通过公式d(t)3tx+1d(t) - 3 * t - x + 1计算,其中t=y+1t = -y + 1

main函数中,使用cin从输入中读取坐标点(x, y),然后调用dis(int x, int y)计算螺旋折线距离,并使用cout输出结果。

注意

对于 100%100\% 的数据,109X,Y109-10^9\le X,Y \le 10^9。dis 可能会超出 int 的范围。需要开 long long,否则无法通过部分测试点。


AC代码

#include <algorithm>
#include <iostream>
#define mp make_pair
#define AUTHOR "HEX9CF"
using namespace std;
using ll = long long;

const int N = 1e6 + 7;
const int INF = 0x3f3f3f3f;
const ll MOD = 1e9 + 7;

// 从原点到(-t, t)
ll d(int t) {
	// 2t * (2t + 1) * 2 / 2
	return (1ll * 2 * t * (2 * t - 1));
}

ll dis(int x, int y) {
	if (x >= 0 && -x <= y && y < x) {
		// 右
		int t = x;
		// d + 2t + (t - y)
		return (1ll * d(t) + 3 * t - y);
	} else if (y > 0 && -y < x && x <= y) {
		// 上
		int t = y;
		// d + (t - x)
		return (1ll * d(t) + t + x);
	} else if (x < 0 && x + 1 < y && y <= -x) {
		// 左
		int t = -x;
		// d - (t - y)
		return (1ll * d(t) - t + y);
	} else {
		// 下
		int t = -y + 1;
		// d - (2t - 1) - (x + t)
		return (1ll * d(t) - 3 * t - x + 1);
	}
}

int main() {
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);

	int x, y;
	cin >> x >> y;
	cout << dis(x, y) << "\n";

	return 0;
}