UVa12121/LA4062 You are around me ...

原创已于 2025-08-15 12:40:03 修改 · 135 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#icpc #Dhaka 2007 #UVa #图论 #最邻近点对 #抽屉原理

于 2022-08-12 15:34:08 首次发布

UVa部分题目解题报告同时被 2 个专栏收录

116 篇文章

订阅专栏

计算几何

18 篇文章

订阅专栏

探讨了在平面上放置多个不重叠椭圆的问题，通过旋转坐标系并应用分治算法求解最大面积。

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

题目链接

本题是2007年 ACM ICPC 亚洲区域赛 Dhaka 赛区的 H 题。

题意

给出平面上N个点 $(1<N<15001)$ ，要求在这些点上放置偏心率为e $(0.2<e\le1)$ ，长轴与x轴夹角为θ的椭圆的圆心，且所有的椭圆不重叠。问：椭圆最大的面积是多少？

分析

椭圆的面积等公式： $S = $$\pi$$ab, b = a \sqrt{1-e^{2}}$ ，b为短轴，a为长轴。

求出最大可能的短轴值，即可求出最大可能的椭圆面积。对于具体的某两个点(xi,yi)和(xj,yj)，其横纵坐标差值dx=xi-xj，dy=yi-yj，对应最大可能的短轴值为 $b = \sqrt{e^{2}dx^{2}+dy^{2}}$ ，此为优化目标。这其实是横向坐标收缩后的两点间距离，所以需要一个预处理：先把所有的坐标逆时针旋转θ度，使所有长轴与x轴平行，再对x轴进行伸缩。之后再求最邻近点对，进而可求面积。

最临近点对是算法导论中讲过的一个经典的分治问题，可在 $O(n\log{n})$ 复杂度内求解。这里面有一个非常棒的分析，可以得到一个很有趣的结论，用于把时间复杂度优化到 $O(n\log{n})$ ，具体参见隔壁知乎。这里简单给出结论：已知当前最短距离为d时，若某个点在下图所示的 $d\times2d$ 矩形边界上，则此矩形内距离它小于d的点最多有6个。

可以结合抽屉原理证明此结论，按上图6等分矩形，则每一个小矩形的长为 $\frac{d}{2}$ ，宽为 $\frac{2d}{3}$ ，对角线长度为 $\frac{5d}{6} < d$ ，如果矩形内有超过7个及以上的点，则必然有某个小矩形内至少包含2个点，而它们的距离小于d，矛盾。

AC代码

#include <iostream>
#include <iomanip>
#include <algorithm>
#include <cmath>
using namespace std;

#define N 15010
#define dis(a, b) sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y))
int a[N], n, kase = 0; double p = M_PI / 180, d, e, t;
struct point {
    double x, y;
    bool operator< (const point& rhs) const {
        return x < rhs.x;
    }
} s[N];

bool cmp(int i, int j) {
    return s[i].y < s[j].y;
}

void closest_dis(int l, int r) {
    if (r-l == 1) {
        d = min(d, dis(s[l], s[r]));
    } else if (r-l > 1) {
        int m = (l+r)>>1, c = 0;
        closest_dis(l, m);
        closest_dis(m+1, r);
        for (int i=m; i>=l; --i) {
            if (s[m].x-s[i].x >= d) break;
            a[c++] = i;
        }
        for (int i=m+1; i<=r; ++i) {
            if (s[i].x-s[m].x >= d) break;
            a[c++] = i;
        }
        sort(a, a+c, cmp);
        for (int i=0; i<c; ++i) for (int j=i+1, k=i+7; j<k; ++j) {
            if (j>=c || s[a[j]].y>=s[a[i]].y+d) break;
            d = min(d, dis(s[a[i]], s[a[j]]));
        }
    }
}

double solve() {
    t *= p; e = sqrt(1 - e*e);
    double co = cos(t), si = sin(t);
    for (int i=0; i<n; ++i) {
        double a, b; cin >> a >> b;
        s[i].x = e * (a*co + b*si);
        s[i].y = b*co - a*si;
    }
    if (e == 0) return 0;
    d = __DBL_MAX__;
    sort(s, s+n);
    closest_dis(0, n-1);
    return M_PI * d * d / 4 / e;
}

int main() {
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(6);
    while (cin>>n>>e>>t && n) cout << "Case " << ++kase << ':' << endl << solve() << endl;
    return 0;
}