POJ 3301 Texas Trip 计算几何 三分法求极值点

题目链接：http://poj.org/problem?id=3301

题目大意：平面上有n个点，求一个可以覆盖所有点的面积最小的正方形。

由于正方形的边与x轴的夹角未知，但是如果已知正方形的两条平行边与x轴的夹角，那么正方形的面积是好求的；正方形的面积是边与x轴夹角的函数。问题在于怎么求出alpha使得正方形的面积最大。这是一道典型的求函数极值的问题，前提是函数是单峰的。组合数学科上老师讲过用优选法求函数的极值，即0.618法。。改天再用这种方法做些优化

总结一下，求函数零点用二分，求函数极值用三分（即0.618法）

已知两个定点的坐标和分别过这两个点的两条平行线与x轴正半轴的夹角，求两平行线之间的距离的推导并不难，高中知识就够。

/*AC*/
#include <cstdio>
#include <iostream>
#include <cmath>
using namespace std;

struct Point{
    int x, y;
};

Point arr[30];
int n;
const double eps = 0.00000005;

inline double dis1(Point p1, Point p2, double a){ //a是两条平行线与x轴正半轴的夹角
    return fabs((p2.y - p1.y)*cos(a) - (p2.x - p1.x)*sin(a));
}

inline double dis2(Point p1, Point p2, double a){ //a' = pi - (pi/2 - a);
    return fabs((p2.y - p1.y)*sin(a) + (p2.x - p1.x)*cos(a));
}

double cal(double a){
    double maxd1 = 0.0, maxd2 = 0.0;
    for(int i=0; i<n-1; i++){
        for(int j=i+1; j<n; j++){
            double d1 = dis1(arr[i], arr[j], a);
            double d2 = dis2(arr[i], arr[j], a);
            if(d1 > maxd1) maxd1 = d1;
            if(d2 > maxd2) maxd2 = d2;
        }
    }
    return max(maxd1, maxd2);
}

int main(){
    int T;
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        for(int i=0; i<n; i++)
            scanf("%d %d", &arr[i].x, &arr[i].y);
        double ans = 1000.0, lb = 0.0, ub = asin(1.0); //0到pi/2的范围内搜索，找出极值点
        while(ub - lb > eps){
            double mid1 = lb + (ub-lb)/3;
            double mid2 = ub - (ub-lb)/3;
            double d1 = cal(mid1);
            double d2 = cal(mid2);
            if(d1 < d2){
                ans = d1;
                ub = mid2;
            }
            else{
                ans = d2;
                lb = mid1;
            }
        }
        printf("%.2lf\n", ans*ans);
    }
    return 0;
}