POJ 3301Texas Trip

题意：给坐标系中几个点，然后求出将所有点都包围起来的最小正方形的面积

思路：三分加向量旋转角公式。。之前做POJ上一题接触了向量旋转公式，这次做起了就轻松多了

不是很懂三分的可以看看这个大牛的文章http://hi.baidu.com/czyuan_acm/item/81b21d1910ea729c99ce33db

#include <iostream>
#include <stdio.h>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iomanip>
#define INF 0xfffffffff
#define PI acos(-1.0)
#define eps 1e-15

using namespace std;

struct Node
{
    double x,y;
}f[50];

int T;
int n;
double cal(double t)
{
    double x,y;
    double xmin,xmax,ymin,ymax;
    xmin=ymin=INF;
    xmax=ymax=-INF;
    for(int i=1;i<=n;i++)
    {
        x=cos(t)*f[i].x-sin(t)*f[i].y;
        y=sin(t)*f[i].x+cos(t)*f[i].y;
         //找到最大最小x,y
        xmin=min(x,xmin);
        xmax=max(x,xmax);
        ymin=min(y,ymin);
        ymax=max(y,ymax);
    }
    return max(xmax-xmin,ymax-ymin);//返回最大差值，这样才能保证所有点都包围进去
}

int main()
{

    while(~scanf("%d",&T))
    {

        while(T--)
        {
            scanf("%d",&n);
            for(int i=1;i<=n;i++)
            scanf("%lf%lf",&f[i].x,&f[i].y);

            double mid,mmid,ri=PI/2,le=0;

    while(fabs(ri-le)>=eps)
    {
        mid=(ri+le)/2.0;
        mmid=(mid+ri)/2.0;
        if(cal(mid)<=cal(mmid))
          ri=mmid;
          else
          le=mid;
    }
     printf("%.2lf\n",cal(le)*cal(le));

        }
    }
    return 0;
}