poj3301 Texas Trip（三分）

题目链接

分析：
题目大意：给出$n$个点的坐标，现在要求一个正方形，完全包围n个点，并且正方形面积最小，求最小的正方形面积

面积随着角度的变化是一个单峰的函数（下凸函数），那么自然就可以想到是三分
按照题目要求正方形最小的面积，如果正方形是平行于$x$轴的，那么正方形面积是x的最大距离$*y$的最大的距离
旋转正方形，在0到90度内总会找到一个正方形面积的最小值

但是旋转正方形比较麻烦，我们可以考虑旋转坐标系
将坐标系旋转0到90度，按旋转的角度重新计算各点的坐标，然后找出x的差和y的差，计算面积

角度旋转公式：$x = x*cos(α) - y*sin(α),y = x*sin(α) + y*cos(α)$

我们就三分一个旋转角度，计算正方形面积即可

tip

注意三分的精度
角度是弧度制

poj上A不了，找了一些网上的题解也A不掉，可能是网站的问题？

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>

using namespace std;

const double Pi=acos(-1.0);
const double INF=1e18;
const double eps=1e-12;
const int N=32;
struct node{
    double x,y;
};
node po[N];
int n;

double rotate(double &x,double &y,double a,int i)
{
    x=po[i].x*cos(a)-po[i].y*sin(a);
    y=po[i].x*sin(a)+po[i].y*cos(a);
}

double cal(double a)
{
    double mxx=-INF,mnx=INF,mxy=-INF,mny=INF;
    for (int i=1;i<=n;i++)
    {
        double x,y;
        rotate(x,y,a,i);
        mxx=max(mxx,x); mnx=min(mnx,x);
        mxy=max(mxy,y); mny=min(mny,y);
    }
    double X=max(mxx-mnx,mxy-mny);
    return X*X;
}

double sanfen()
{
    double l=0,r=Pi;
    double m1,m2;
    while (r-l>eps)
    {
        m1=l+(r-l)/3;
        m2=r-(r-l)/3;
        if (cal(m1)>cal(m2))
           l=m1;
        else r=m2;
    }
    return cal((m1+m2)/2);
}

int main()
{
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d",&n);
        for (int i=1;i<=n;i++)
            scanf("%lf%lf",&po[i].x,&po[i].y);
        printf("%0.2lf\n",sanfen());
    }
    return 0;
}