题目链接
本题是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;
}