UVA10173(求凸包的面积最小外接矩形)

题目：10173 - Smallest Bounding Rectangle

求凸包的最小外接矩形的面积。

思路：

旋转卡壳
给定点集S,求S的最小覆盖矩形
最小覆盖矩形的四条边上，其中一条边有至少两个点，其他边上至少有一个点。
然后沿着凸包的边旋转，维护矩形另外三条边上的点。

#include<stdio.h>
#include<cmath>
#include<algorithm>
#define eps 1e-8
#define N 50010
using namespace std;
struct Point
{
double x,y;
Point(){}
Point(double x0,double y0):x(x0),y(y0){}
};
Point p[N];
int con[N];
int cn;
int n;
struct Line
{
Point a,b;
Line(){}
Line(Point a0,Point b0):a(a0),b(b0){}
};
double Xmult(Point o,Point a,Point b)
{
return (a.x-o.x)*(b.y-o.y)-(b.x-o.x)*(a.y-o.y);
}
double Dmult(Point o,Point a,Point b)
{
return (a.x-o.x)*(b.x-o.x)+(a.y-o.y)*(b.y-o.y);
}
int Sig(double a)
{
return a<-eps?-1:a>eps;
}
double Dis(Point a,Point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int cmp(Point a,Point b)
{
double d=Xmult(p[0],a,b);
if(d>0)
return 1;
if(d==0 && Dis(p[0],a)<Dis(p[0],b))
return 1;
return 0;
}
double min(double a,double b)
{
return a<b?a:b;
}
void Graham()
{
int i,ind=0;
for(i=1;i<n;i++)
if(p[ind].y>p[i].y || (p[ind].y==p[i].y) && p[ind].x>p[i].x)
ind=i;
swap(p[ind],p[0]);
sort(p+1,p+n,cmp);
con[0]=0;
con[1]=1;
cn=1;
for(i=2;i<n;i++)
{
while(cn>0 && Sig(Xmult(p[con[cn-1]],p[con[cn]],p[i]))<=0)
cn--;
con[++cn]=i;
}
int tmp=cn;
for(i=n-2;i>=0;i--)
{
while(cn>tmp && Sig(Xmult(p[con[cn-1]],p[con[cn]],p[i]))<=0)
cn--;
con[++cn]=i;
}
}
double Solve()
{
int t,r,l;
double ans=999999999;
t=r=1;
if(cn<3)
return 0;
for(int i=0;i<cn;i++)
{
while(Sig( Xmult(p[con[i]],p[con[i+1]],p[con[t+1]])-
Xmult(p[con[i]],p[con[i+1]],p[con[t]])   )>0)
t=(t+1)%cn;
while(Sig( Dmult(p[con[i]],p[con[i+1]],p[con[r+1]])-
Dmult(p[con[i]],p[con[i+1]],p[con[r]])   )>0)
r=(r+1)%cn;
if(!i) l=r;
while(Sig( Dmult(p[con[i]],p[con[i+1]],p[con[l+1]])-
Dmult(p[con[i]],p[con[i+1]],p[con[l]])   )<=0)
l=(l+1)%cn;
double d=Dis(p[con[i]],p[con[i+1]]);
double tmp=Xmult(p[con[i]],p[con[i+1]],p[con[t]])*
( Dmult(p[con[i]],p[con[i+1]],p[con[r]])-
Dmult(p[con[i]],p[con[i+1]],p[con[l]]) )/d/d;
ans=min(ans,tmp);
}
return ans;
}
int main()
{
int i;
while(scanf("%d",&n) && n)
{
for(i=0;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
Graham();
printf("%.4f\n",Solve());
}
return 0;
}

UVA10173(求凸包的面积最小外接矩形)相关推荐

opencv findContours()轮廓特征分析大全(求面积、周长、几何矩、质心、凸包、最小外接矩形、最小外接三角形、最小外接椭圆等)
文章目录前言一.效果二.opencv对应函数介绍 1.轮廓面积 2.轮廓周长 3.轮廓几何矩 4.轮廓的最大外接矩形 5.轮廓的凸包 6.轮廓的最小外接矩形 7.轮廓的最小外接三角形 8.轮廓的 ...
C++/Qt：最小外接矩形
以凸包的某条边作为矩形的一条边,求包围所有点的面积最小外接矩形.(感觉是对的)下图中蓝色为凸包,红色为MBR. 分享给有需要的人,代码质量勿喷. 一.求凸包,返回 xjList ...
三维空间点集的最小外接矩形
文章目录引言说明二维空间内的方案方案选定新坐标轴变换矩阵 G G G 后续引言前段时间毕设卡在求一系列点的最小外接矩形, 卡了好久在导师的帮助下终于解决了. 来分享下思路说明本篇所 ...
旋转卡壳算法求最小外接矩形代码
旋转卡壳原理:旋转卡壳详解_大学要有梦想的博客-CSDN博客_旋转卡壳思路: 1.选择卡壳算法用于求凸多边形的最小外接矩形 1.多边形最小的外接矩形一定是以多边形的的一条边为底的一部分 2.通过这条 ...
opencv 图像轮廓特征图像面积，轮廓周长，外接矩形、最小外接矩形、最小外接圆、拟合椭圆
找出图像轮廓 contours, hierarchy = cv.findContours(thresh, 3, 2) 画出图像轮廓 cnt = contours[1] cv.drawContours( ...
求取SHP文件的最小外接矩形并裁剪图像
目的: 求取shp文件中每一个形状的最小外接矩形. 根据每一个形状的最小外接矩形裁剪图像. 已知数据: 一个shp文件,包含若干个形状. 2.shp文件对应的影像. 工具 ARCGIS10.4 pyt ...
opencv 轮廓查找，凸包，最小外接矩形，最小外接圆，最小外接椭圆
本章内容: 1. 轮廓查找 2. 绘制轮廓 3. 凸包 4.最小外接矩形 5.最小外接圆 6.最小外接椭圆 1.搜索轮廓 2.绘制轮廓输出结果 3.凸包输出结果 4.最小外接矩形输出结果 5.最 ...
【千律】OpenCV基础：图像外接矩形、最小外接矩形、凸包、外接圆、拟合椭圆的绘制
环境:Python3.8 和 OpenCV 内容:图像外接矩形.最小外接矩形.凸包.外接圆.拟合椭圆的绘制 import cv2 as cv import numpy as np import mat ...
matlab外接矩形,matlab求二值圖像最小外接矩形
全部來自ilovematlab中文論壇,其中matlab較新的版本不再支持其中的有些函數,或者說有些函數的參數做了改變,因此適當修改后貼出來.注:無意侵權,記錄在此自我學習,順帶服務群眾. 首先是尋找 ...

UVA10173(求凸包的面积最小外接矩形)

UVA10173(求凸包的面积最小外接矩形)相关推荐

最新文章

热门文章