题目大意:
在二维坐标系中给出一些点,求能覆盖他们的最小正方形的面积(正方形的边不一定平行坐标轴)
解题思路:
对于一个点,若坐标轴旋转a度(弧度制),那么X‘=X*cos(a)-Y*sin(a);Y‘=Y*cos(a)+X*sin(a);
对于角度三分,对于正方形面积是个单峰函数在【0,pi】。有最小值。
下面是代码:
#include <set> #include <map> #include <queue> #include <math.h> #include <vector> #include <string> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <iostream> #include <algorithm> #define eps 1e-8 #define pi acos(-1.0) #define inf 107374182 #define inf64 1152921504606846976 #define lc l,m,tr<<1 #define rc m + 1,r,tr<<1|1 #define iabs(x) ((x) > 0 ? (x) : -(x)) #define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE)) #define clearall(A, X) memset(A, X, sizeof(A)) #define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE)) #define memcopyall(A, X) memcpy(A , X ,sizeof(X)) #define max( x, y ) ( ((x) > (y)) ? (x) : (y) ) #define min( x, y ) ( ((x) < (y)) ? (x) : (y) ) using namespace std; struct node { double x,y; }point[305]; int n; double does(double a) { double maxx=-10000000,maxy=-100000000,minx=100000000,miny=100000000,tx,ty; for(int i=0;i<n;i++) { tx=point[i].x*cos(a)-point[i].y*sin(a); ty=point[i].y*cos(a)+point[i].x*sin(a); maxx=max(maxx,tx); maxy=max(maxy,ty); minx=min(minx,tx); miny=min(miny,ty); } return max(maxx-minx,maxy-miny); } int main() { int T; double l,r,lmid,rmid,ans1,ans2; scanf("%d",&T); while(T--) { scanf("%d",&n); for(int i=0;i<n;i++) { scanf("%lf%lf",&point[i].x,&point[i].y); } l=0.0; r=pi; while(r-l>eps) { lmid=(l+r)/2; rmid=(lmid+r)/2; ans1=does(lmid); ans2=does(rmid); if(ans1<=ans2)r=rmid; else l=lmid; } printf("%.2lf\n",ans1*ans1); } return 0; }
POJ 3301 Texas Trip,布布扣,bubuko.com
原文:http://blog.csdn.net/lin375691011/article/details/38511577