3707: 圈地

3707: 圈地

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 229  Solved: 85
[Submit][Status][Discuss]

Description

2维平面上有n个木桩，黄学长有一次圈地的机会并得到圈到的土地，为了体现他的高风亮节，他要使他圈到的土地面积尽量小。圈地需要圈一个至少3个点的多边形，多边形的顶点就是一个木桩，圈得的土地就是这个多边形内部的土地。（因为黄学长非常的神，所以他允许圈出的第n点共线，那样面积算0）

Input

第一行一个整数n，表示木桩个数。
接下来n行，每行2个整数表示一个木桩的坐标，坐标两两不同。

Output

仅一行，表示最小圈得的土地面积，保留2位小数。

Sample Input

3
0 0
0 1
1 0

Sample Output

0.50

HINT

对于100%的数据，n<=1000。

Source

[Submit][Status][Discuss]
网上正解一大堆啦。。。但是一个都看不懂不过此题有一个有趣的暴力方法，将所有点按照x坐标第一关键字，y坐标第二关键字排序每大概50个点分成一个小块，块内暴力枚举尝试更新答案但是这样的做法肯定是挂的，因为两块之间的点对并没有枚举不过，可以随机旋转坐标轴，再重复上面的操作，执行100次一般就能出正解了

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>using namespace std; const int maxn = 1010;typedef double DB; struct Point{    DB x,y; Point(){}    Point (DB x,DB y): x(x),y(y){}    bool operator < (const Point &B) const    {        if (x < B.x) return 1;        if (x > B.x) return 0;        return y < B.y;    }    Point operator - (const Point &B) {return Point(x - B.x,y - B.y);}}p[maxn],q[maxn];typedef Point Vector; int n;DB Ans = 1E16; int max(const int &x,const int &y) {return x > y ? x : y;}DB Cross(const Vector &v1,const Vector &v2) {return v1.x * v2.y - v2.x * v1.y;} void Calc(int l,int r){    for (int i = l; i <= r; i++)        for (int j = i + 1; j <= r; j++)            for (int k = j + 1; k <= r; k++)                Ans = min(Ans,fabs(Cross(q[j] - q[i],q[k] - q[i])));} int main(){    #ifdef DMC        freopen("DMC.txt","r",stdin);    #endif         cin >> n; int siz = 50;    for (int i = 1; i <= n; i++) scanf("%lf%lf",&p[i].x,&p[i].y);    for (int I = 0; I < 100; I++)    {        DB alpha = (DB)(rand() - rand()) / 233.00;        DB Cos = cos(alpha),Sin = sin(alpha); //printf("%.8lf %.8lf\n",Cos,Sin);        for (int i = 1; i <= n; i++)            q[i] = Point(p[i].x * Cos - p[i].y * Sin,p[i].x * Sin + p[i].y * Cos);        sort(q + 1,q + n + 1); int last = 1,tot = 0;        for (int i = 1; i <= n; i++)        {            ++tot;            if (tot == siz || i == n) Calc(last,i),tot = 0,last = i + 1;        }    }    printf("%.2f\n",Ans / 2.00);    return 0;}