百度之星初赛1006(计算几何:能包含凸包的最小矩形面积)

矩形面积

 

 Accepts: 717

 

 Submissions: 1619

 Time Limit: 2000/1000 MS (Java/Others)

 

 Memory Limit: 32768/32768 K (Java/Others)

Problem Description

小度熊有一个桌面，小度熊剪了很多矩形放在桌面上，小度熊想知道能把这些矩形包围起来的面积最小的矩形的面积是多少。

Input

第一行一个正整数 T，代表测试数据组数( 1≤T≤20 )，接下来 T 组测试数据。

每组测试数据占若干行，第一行一个正整数  N(1≤N<≤1000) ，代表矩形的数量。接下来 N 行，每行 8 个整数 x1,y1,x2,y2,x3,y3,x4,y4 ，代表矩形的四个点坐标，坐标绝对值不会超过10000。

Output

对于每组测试数据，输出两行：

第一行输出"Case #i:"，i 代表第 i 组测试数据。 第二行包含1 个数字，代表面积最小的矩形的面积，结果保留到整数位。

Sample Input

2
2
5 10 5 8 3 10 3 8
8 8 8 6 7 8 7 6
1
0 0 2 2 2 0 0 2

Sample Output

Case #1:
17
Case #2:
4

问题可以转化为平面上有4*n个点，求能够包含这些点的最小矩形面积，直接上模板（蓝儿弱没怎么做计算几何，都不知道，先mark一下）

#include<stdio.h>
#include<cmath>
#include<algorithm>
#define eps 1e-8
#define N 50010
using namespace std;
struct Point
{
    double x,y;
    Point() {}
    Point(double x0,double y0):x(x0),y(y0) {}
};
Point p[N];
int con[N];
int cn;
int n;
struct Line
{
    Point a,b;
    Line() {}
    Line(Point a0,Point b0):a(a0),b(b0) {}
};
double Xmult(Point o,Point a,Point b)
{
    return (a.x-o.x)*(b.y-o.y)-(b.x-o.x)*(a.y-o.y);
}
double Dmult(Point o,Point a,Point b)
{
    return (a.x-o.x)*(b.x-o.x)+(a.y-o.y)*(b.y-o.y);
}
int Sig(double a)
{
    return a<-eps?-1:a>eps;
}
double Dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int cmp(Point a,Point b)
{
    double d=Xmult(p[0],a,b);
    if(d>0)
        return 1;
    if(d==0 && Dis(p[0],a)<Dis(p[0],b))
        return 1;
    return 0;
}
double min(double a,double b)
{
    return a<b?a:b;
}
void Graham()
{
    int i,ind=0;
    for(i=1; i<n; i++)
        if(p[ind].y>p[i].y || (p[ind].y==p[i].y) && p[ind].x>p[i].x)
            ind=i;
    swap(p[ind],p[0]);
    sort(p+1,p+n,cmp);
    con[0]=0;
    con[1]=1;
    cn=1;
    for(i=2; i<n; i++)
    {
        while(cn>0 && Sig(Xmult(p[con[cn-1]],p[con[cn]],p[i]))<=0)
            cn--;
        con[++cn]=i;
    }
    int tmp=cn;
    for(i=n-2; i>=0; i--)
    {
        while(cn>tmp && Sig(Xmult(p[con[cn-1]],p[con[cn]],p[i]))<=0)
            cn--;
        con[++cn]=i;
    }
}
double Solve()
{
    int t,r,l;
    double ans=999999999;
    t=r=1;
    if(cn<3)
        return 0;
    for(int i=0; i<cn; i++)
    {
        while(Sig( Xmult(p[con[i]],p[con[i+1]],p[con[t+1]])-
                   Xmult(p[con[i]],p[con[i+1]],p[con[t]])   )>0)
            t=(t+1)%cn;
        while(Sig( Dmult(p[con[i]],p[con[i+1]],p[con[r+1]])-
                   Dmult(p[con[i]],p[con[i+1]],p[con[r]])   )>0)
            r=(r+1)%cn;
        if(!i) l=r;
        while(Sig( Dmult(p[con[i]],p[con[i+1]],p[con[l+1]])-
                   Dmult(p[con[i]],p[con[i+1]],p[con[l]])   )<=0)
            l=(l+1)%cn;
        double d=Dis(p[con[i]],p[con[i+1]]);
        double tmp=Xmult(p[con[i]],p[con[i+1]],p[con[t]])*
                   ( Dmult(p[con[i]],p[con[i+1]],p[con[r]])-
                     Dmult(p[con[i]],p[con[i+1]],p[con[l]]) )/d/d;
        ans=min(ans,tmp);
    }
    return ans;
}
int main()
{
    int i; int T, cas = 1;
    scanf("%d", &T);
    while(T --)
    {
        scanf("%d",&n);
        n *= 4;
        for(i=0; i<n; i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        Graham();
        printf("Case #%d:\n%.0f\n",cas ++, Solve());
    }
    return 0;
}