poj3301

参考三分讲解。

http://hi.baidu.com/czyuan_acm/item/81b21d1910ea729c99ce33db

#include <iostream>
#include <algorithm>
#include <math.h>
#include <iomanip>

using namespace std;

#define pi acos(-1.0)
#define eps 1e-8
#define INF 1<<30

struct point{ int x, y; }p[35];
int n, t;

double calc(double angle)
{
    double xmax, xmin, ymax, ymin, y, x;
    xmax = ymax = -INF;
    xmin = ymin = INF;
    for (int i = 0; i < n; i++)//rotate
    {
        x = cos(angle)*p[i].x - sin(angle)*p[i].y;
        y = sin(angle)*p[i].x + cos(angle)*p[i].y;
        xmax = max(xmax, x);
        xmin = min(xmin, x);
        ymax = max(ymax, y);
        ymin = min(ymin, y);
    }
    return max(xmax - xmin, ymax - ymin);
}
void solve()
{
    double low = 0, high = pi/4, mid, midmid;
    while (fabs(high - low) > eps)// the process
    {
        mid = (high + low) / 2;
        midmid = (mid + high) / 2;
        if (calc(mid) <= calc(midmid))
            high = midmid;
        else low = mid;
    }
    double ans = calc(mid);
    cout << fixed << setprecision(2) << ans * ans << endl;
}
int main()
{
    cin >> t;
    while (t--)
    {
        memset(p, 0, sizeof(p));
        cin >> n;
        for (int i = 0; i < n; i++)
        {
            cin >> p[i].x >> p[i].y;
        }
        solve();
    }
    return 0;
}