题意:用最小的正方形覆盖所有的点,输出此正方形的面积。
思路:一条边长与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;
}
该博客探讨了如何用最小的正方形覆盖所有点的问题,通过分析正方形边长与x轴平行时,旋转角度对其投影长度的影响,得出这是一个单峰函数。作者提出采用三分搜索的方法来解决,但由于实际操作中旋转正方形复杂,因此选择旋转坐标系以简化计算。文章介绍了坐标旋转的数学表达式。
720

被折叠的 条评论
为什么被折叠?



