计算几何专项：UVa 12307

应用旋转卡壳求最小面积外接矩阵和最小周长外接矩阵，在网上找到了实现策略，可是不会自己实现。参考了别人的代码，搞了一上午把细节搞懂了，奉上AC代码。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;
const double eps=1e-6;
const double inf=1e10;
int dcmp(double x)
{if(fabs(x)<eps) return 0;else return x<0?-1:1;
}
struct point
{double x,y;point(double x=0,double y=0):x(x),y(y){}
};
point operator-(point a,point b){return point(a.x-b.x,a.y-b.y);}
point operator+(point a,point b){return point(a.x+b.x,a.y+b.y);}
point operator*(point a,double p){return point(a.x*p,a.y*p);}
bool operator<(const point& a,const point& b)
{return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
bool operator==(const point& a,const point& b)
{return dcmp(a.x-b.x)==0&&dcmp(a.y-b.y)==0;
}
double cross(point a,point b){return a.x*b.y-a.y*b.x;}
double dot(point a,point b){return a.x*b.x+a.y*b.y;}
double length(point a){return sqrt(dot(a,a));}
point normal(point a)
{double l=length(a);return point(a.x/l,a.y/l);
}
double distoline(point p,point a,point b)
{point v1=b-a,v2=p-a;return fabs(cross(v1,v2))/length(v1);
}
vector<point> p;
double mina,minp;
vector<point> convex(vector<point>& p)
{sort(p.begin(),p.end());int n=p.size();vector<point> ch(n+1);int m=0;for(int i=0;i<n;i++){while(m>1&&cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;ch[m++]=p[i];}int k=m;for(int i=n-2;i>=0;i--){while(m>k&&cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0) m--;ch[m++]=p[i];}if(n>1) m--;ch.resize(m);return ch;
}
void rotating_calipers(vector<point>& points)
{vector<point> p=convex(points);int n=p.size();p.push_back(p[0]);mina=inf;minp=inf;int l=1,r=1,u=1;for(int i=0;i<n;i++){point edge=normal(p[(i+1)%n]-p[i]);while(dot(edge,p[r%n]-p[i])<dot(edge,p[(r+1)%n]-p[i])) r++;while(u<r||cross(edge,p[u%n]-p[i])<cross(edge,p[(u+1)%n]-p[i])) u++;while(l<u||dot(edge,p[l%n]-p[i])>dot(edge,p[(l+1)%n]-p[i])) l++;double w=dot(edge,p[r%n]-p[i])-dot(edge,p[l%n]-p[i]);double h=distoline(p[u%n],p[i],p[(i+1)%n]);mina=min(mina,w*h);minp=min(minp,2*(w+h));}
}
int main()
{int n;while(~scanf("%d",&n)){if(!n) break;p.clear();for(int i=0;i<n;i++){double x,y;scanf("%lf%lf",&x,&y);p.push_back(point(x,y));}rotating_calipers(p);printf("%.2f %.2f\n",mina,minp);}return 0;
}

计算几何专项：UVa 12307相关推荐

ACM计算几何专项练习题目总结
//第一期计算几何题的特点与做题要领: 1.大部分不会很难,少部分题目思路很巧妙 2.做计算几何题目,模板很重要,模板必须高度可靠. 3.要注意代码的组织,因为计算几何的题目很容易上两百行代码,里面 ...
POJ 计算几何专项训练（1）【2318】【2398】【3304】【2653】【1556】【1066】...
POJ 2318 TOYS 题意是在一个大矩形里有n条分割线把矩形分割成n+1部分.再给出一些玩具的坐标.要求统计每个部分内有多少个玩具. 具体做法就是二分求解出当前玩具右方的第一条线.这可以用叉积判 ...
计算几何(二分) - Crossed Ladders - UVA 10566
计算几何(二分) - Crossed Ladders - UVA 10566 题意: 如上图,给定三个浮点数x,y,c,计算两个房子之间的间距.如上图,给定三个浮点数x,y,c,计算两个房子之间的间距 ...
计算几何(三角形的外接圆+圆内接正多边形) - Temple of Dune - UVA 10439
计算几何(三角形的外接圆+圆内接正多边形) - Temple of Dune - UVA 10439 题意: 二维平面内 , 给定 3 个点的坐标 , 二维平面内,给定3个点的坐标, ...
计算几何(梅涅劳斯定理) - Triangle Hazard - UVA 12165
计算几何(梅涅劳斯定理) - Triangle Hazard - UVA 12165 题意: 如上图,给定△PQR的三个顶点坐标以及6个比例m1,m2,...,m6,如上图,给定\triangle P ...
计算几何(三角形内切圆 + 海伦公式) - InCircle - UVA 11524
计算几何(三角形内切圆 + 海伦公式) - InCircle - UVA 11524 题意: 对于△ABC,对于\triangle ABC,对于△ABC, 假设其内切圆与其相切于P.Q.R三点,假设其 ...
UVA 11178 Morley’s Theorem（莫雷定理计算几何）
Morley's Theorem Input: Standard Input Output: Standard Output Morley's theorem states that that the ...
uva 12426 Counting Triangles 计算几何
题意:给一个凸多边形,求以凸多边形的顶点为顶点的三角形中,有多少个面积小于给定数值K? 思路:旋转卡壳+二分答案时间复杂度:n^2logn 1 #include<iostream> 2 ...
UVA 11178-Morley's Theorem(计算几何_莫雷定理)
Problem D Morley's Theorem Input: Standard Input Output: Standard Output Morley's theorem states tha ...
也许是史上最不良心的低阶计算几何讲解和习题集？？
-3.在此声明: 笔者极端厌恶计算几何,以至于直到今天之前都只打过几个计算几何的模板~~~~~ 不过鉴于今年是18年,所以感觉SD很有可能考计算几何(18年是什么理由啊喂) 于是滚过来整理计算几何的资 ...

计算几何专项：UVa 12307

计算几何专项：UVa 12307相关推荐

最新文章

热门文章