Codeforces 429D Tricky Function（平面最近点对）

最新推荐文章于 2020-05-19 10:39:34 发布

转载最新推荐文章于 2020-05-19 10:39:34 发布 · 149 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/cxhscst2/p/8016415.html

本文介绍了一种解决平面最近点对问题的方法，并通过具体的代码实现了解决方案。将问题转化为平面上两点间的距离最小值问题，使用分治策略进行高效求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接 Tricky Function

$f(i, j) = (i - j)^{2} + (s[i] - s[j])^{2}$

把$(i, s[i])$塞到平面直角坐标系里，于是转化成了平面最近点对问题。

#include <bits/stdc++.h>

using namespace std;

#define rep(i, a, b)	for (int i(a); i <= (b); ++i)
#define dec(i, a, b)	for (int i(a); i >= (b); --i)

typedef long long LL;

const int N = 100010;

struct Point{
	LL x, y;
	void scan(){ scanf("%lld%lld", &x, &y);}
} p[N], q[N];

int n;
LL c[N];
LL s;

bool cmp_x(const Point &a, const Point &b){ 
	return a.x < b.x;
}
bool cmp_y(const Point &a, const Point &b){
	return a.y < b.y;
}

LL dist(const Point &a, const Point &b){ 
	return (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
}


LL work(int l, int r){
	if (r == l) return 9e18;
	if (r - l == 1) return dist(p[l], p[r]);
	if (r - l == 2) return min(min(dist(p[l], p[r]), dist(p[l + 1], p[r])), dist(p[l], p[l + 1]));
	int mid = (l + r) >> 1, cnt = 0;
	LL ret = min(work(l, mid), work(mid + 1, r));

	rep(i, l, r) if (p[i].x < p[mid].x + sqrt(ret) && p[i].x > p[mid].x - sqrt(ret)) q[++cnt] = p[i];
	sort(q + 1, q + cnt + 1, cmp_y);
	rep(i, 1, cnt - 1) rep(j, i + 1, cnt){
		if ((q[j].y - q[i].y) * (q[j].y - q[i].y) > ret) break;
		ret = min(ret, dist(q[i], q[j]));
	}

	return ret;
}

int main(){

	scanf("%d", &n);
	rep(i, 1, n) scanf("%lld", c + i);
	rep(i, 1, n){
		s += c[i];
		p[i] = {i, s};
	}
	sort(p + 1, p + n + 1, cmp_x);
	printf("%lld\n", work(1, n));
	return 0;

}

转载于:https://www.cnblogs.com/cxhscst2/p/8016415.html