HDU1007
Posted qq-1615160629
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU1007相关的知识,希望对你有一定的参考价值。
layout: post
title: "HDU 1007"
date: 2017-04-23 14:27:18 +0800
categories: ACM
tags: Partition
author: SteveDevin
mathjax: true
- content
{:toc}
思路:
分治裸题, 求平面最接近点对: 对于要求最接近点对的区间S[l, r], 取 m = (l + r) >> 1, 则最短距离为以x或y排序后的区间 S[l, m], S[m + 1, r], 以及左右两个区间之间某两个点距离的最小值.
设 min(S[l, m], S[m + 1, r]) = ret; 则两个区间之间的最短距离只需考虑 x 坐标 在 [m - ret, m + ret]的点即可, 进而可以证明只需考虑[m - ret, m + ret]中的6个点即可.
我先是枚举[m - ret, m + ret] 所有的点求最短距离, TLE. 换成区间内最近的6个点后AC (感觉这样不严谨, 数据太水?)
代码
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn = 100000 + 2;
struct Point{
double x;
double y;
}points[maxn];
int n;
bool comp(Point &a, Point &b)
{
if(a.x == b.x) return a.y < b.y;
return a.x < b.x;
}
double calc(int a, int b)
{
return sqrt(fabs(pow(points[a].x - points[b].x, 2) + pow(points[a].y - points[b].y, 2)));
}
double CPair(int l, int r)
{
if(r - l + 1 < 2) return 1 << 28;
int m = (l + r) >> 1;
double ret = 1 << 28;
ret = CPair(l, m);
ret = min(ret, CPair(m + 1, r));
int L = m, R = 0;
//while(L && points[L].x > points[m].x - ret) L--;
// while(R < r && points[R].x < points[m].x + ret) R++;
L = max(m - 3, l);
R = min(m + 3, r);
for(int j = L; j <= m; j++)
for(int k = m + 1; k <= r; k++)
ret = min(ret, calc(j, k));
return ret;
}
int main()
{
while(~scanf("%d", &n) && n)
{
for(int i = 0; i < n; i++)
scanf("%lf%lf", &points[i].x, &points[i].y);
sort(points, points + n, comp);
printf("%.2f
", CPair(0, n - 1) / 2);
}
}
Minds:
感觉自己越来越水了, 之前刷了那么长时间的题, 才不过两个月左右没刷题就废了.
距离省赛还有两个周, 好慌啊……
以上是关于HDU1007的主要内容,如果未能解决你的问题,请参考以下文章