EOJ 1019. 着弹点

Time limit per test: 2.0 seconds

Memory limit: 256 megabytes

炮兵某部进行实弹射击 , 对一个平面区域里连续开炮 , 得到了很多的弹坑 . 当射击完成后 , 作为技术人员的你 ,

想要得到一个重要的参数 , 就是相隔距离最大的炮弹着弹点的距离 .

Input

多组数据，每组第一行 n 代表点数，接着 n 行为着弹点的坐标，坐标为正整数，不超过 1000000。

2<=n<=30000。

Output

每组一行，求出相隔距离最大的两个着弹点之间的距离，保留 2 位小数

Examples

input

3
1 2
3 4
5 6

output

5.66

 

 

裸的旋转卡壳，找了一个模版结果是错的。。。。

后来换了另一个才过。。。

其实这个题也可以直接暴力，双向求3000对也能出答案。。。

 1 #include <cstdio>
 2 #include <algorithm>
 3 #include <cmath>
 4 #include <iostream>
 5 #include <cstdlib>
 6 #define N 30005
 7 using namespace std;
 8 struct point {
 9     double x, y;
10 
11     point(double x = 0, double y = 0) : x(x), y(y) {}
12 
13     bool operator<(const point &c) const { return x < c.x || (x == c.x && y < c.y); }
14 
15     point operator-(const point &c) const { return point(x - c.x, y - c.y); }
16 
17     double operator*(const point &c) const { return x * c.y - y * c.x; }
18 
19     double operator|(const point &c) const { return (x - c.x) * (x - c.x) + (y - c.y) * (y - c.y); }
20 } a[N], p[N];
21 int n;
22 int convex_hull(point a[], point p[], int tot) {
23     sort(a + 1, a + 1 + tot);
24     int top = 0;
25     for (int i = 1; i <= tot; i++) {
26         while (top > 1 && (p[top] - p[top - 1]) * (a[i] - p[top - 1]) <= 0) top--;
27         p[++top] = a[i];
28     }
29     int k = top;
30     for (int i = tot - 1; i >= 1; i--) {
31         while (top > k && (p[top] - p[top - 1]) * (a[i] - p[top - 1]) <= 0) top--;
32         p[++top] = a[i];
33     }
34     return top - 1;
35 }
36 double rotating_calipers(point p[], int top) {
37     double ans = 0;
38     p[top + 1] = p[1];
39     int now = 2;
40     for (int i = 1; i <= top; i++) {
41         while ((p[i + 1] - p[i]) * (p[now] - p[i]) < (p[i + 1] - p[i]) * (p[now + 1] - p[i])) {
42             now++;
43             if (now == top + 1) now = 1;
44         }
45         ans = max(ans, (p[now] | p[i]));
46     }
47     return sqrt(ans);
48 }
49 
50 int main() {
51     while (scanf("%d", &n)!=EOF) {
52         for (int i = 1; i <= n; i++)
53             scanf("%lf%lf", &a[i].x, &a[i].y);
54         int top = convex_hull(a, p, n);
55         printf("%.2lf\n", rotating_calipers(p, top));
56     }
57     return 0;
58 }

View Code

 

转载于:https://www.cnblogs.com/htwx/articles/8029821.html