POJ3301最小正方形覆盖（三分）

关注

原创最新推荐文章于 2024-10-26 16:15:10 发布 1.5k 阅读 AI 写同款· GEO 优化 ›

poj 3301 Texas Trip(旋转+三分) 题目：http://poj.org/problem?id=3301 大意：给出一些点求出最小的正方形面积把它们都覆盖掉。要解决它我们需要会一些基础知识：坐标旋转公式，三分搜索算法。逆时针旋转：x=xcosθ+ysinθ; y=ycosθ-xsinθ 顺时针旋转：x=xcosθ-ysinθ; y=ycosθ+xsinθ 阅读详情

题意：用最小的正方形覆盖所有的点，输出此正方形的面积。

思路：一条边长与x轴平行的正方形，旋转的角度与投影到x轴的长度的变化关系是一个单峰函数，所以用三分解题，但旋转正方形不方便，所以旋转坐标系。

x=a[i].x*cos(degree)-a[i].y*sin(degree);
y=a[i].x*sin(degree)+a[i].y*cos(degree);

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
using namespace std;
const double eps =1e-12;
const int maxn = 50;
struct Node{
	int x;
	int y;
}a[maxn];
int n;
double C(double degree){
	double x,y,maxx,maxy,minx,miny;
	maxx = maxy = -INF;
	minx = miny = INF;
	for(int i=0;i<n;i++){
		x=a[i].x*cos(degree)-a[i].y*sin(degree);
		y=a[i].x*sin(degree)+a[i].y*cos(degree);
		maxx=max(maxx,x);
		maxy=max(maxy,y);
		minx=min(minx,x);
		miny=min(miny,y);
	}
	return max(maxx-minx,maxy-miny);
}
double bsearch(){
	double left=0,right=PI/2.0;
	double midl,midr;
	while(right-left>eps){
		midl=(left+right)/2;
		midr=(midl+right)/2;
		if(C(midl)<C(midr)){
			right=midr;
		}
		else{
			left=midl;
		}
	}
	return C(left);
}
int main()
{
	int T;
	cin>>T;
	while(T--){
		cin>>n;
		for(int i=0;i<n;i++){
			cin>>a[i].x>>a[i].y;
		}
		double low = bsearch();
		printf("%.2lf\n",low*low);
	}
	return 0;
}

收录于