poj3301 三分

传送门

题意：给定n个点，求能覆盖住这些点的最小的正方形的面积

思路：这个题的想法我觉得还是很巧妙的，一般来说我们考虑点不动，来找这个正方形，然而就算计算个凸包好像也不好实现。这个题通过旋转点，而正方形都是水平竖直的，旋转点寻找最大最小的横坐标之差与最大最小的纵坐标之差的最大值就是当前旋转角度下得到的正方形的边长（表达力为负啊……），每个角度都对应一个正方形的边长，并且是一个单峰函数，所以用三分来求得这个最小边长就可以了

旋转角度的转化可以用和差公式，如果原来的点为(x1,y1)，转化后的点为(x2,y2)，角度为x，那么

x2=x1*cos(x)-y1*sin(x)

y2=y1*cos(x)+x1*sin(x)

完整代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
typedef long long LL;
const double PI=acos(-1.0);
const double eps=1e-9;
int n;
struct point
{
    double x,y;
}data[35];
double Function(double a)
{
    double maxx=-600,maxy=-600;
    double minx=600,miny=600;
    double nx,ny;
    for(int i=0;i<n;i++)
    {
        nx=data[i].x*cos(a)-data[i].y*sin(a);
        ny=data[i].y*cos(a)+data[i].x*sin(a);
        maxx=max(maxx,nx);
        maxy=max(maxy,ny);
        minx=min(minx,nx);
        miny=min(miny,ny);
    }
    return max((maxx-minx),(maxy-miny));
}
double trisection(double l,double r)
{
    while(r-l>eps)
    {
        double ll=l+(r-l)/3;
        double rr=l+2*(r-l)/3;
        double xl=Function(ll);
        double xr=Function(rr);
        if(xl<xr)
            r=rr;
        else l=ll;
    }
    return l;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&data[i].x,&data[i].y);
        double left=0;
        double right=PI/2;
        double res=trisection(left,right);
        double ans=Function(res)*Function(res);
        printf("%.2f\n",ans);
    }
    return 0;
}