UVa12121/LA4062 You are around me ...

题目链接

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

题意

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

分析

        椭圆的面积等公式：，b为短轴，a为长轴。

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

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

        可以结合抽屉原理证明此结论，按上图6等分矩形，则每一个小矩形的长为，宽为，对角线长度为，如果矩形内有超过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;
}