POJ 2079(计算几何+凸包+旋转卡壳法)
来源:互联网 发布:侧铣头编程讲解 编辑:程序博客网 时间:2024/05/22 05:22
问题描述:
Given n distinct points on a plane, your task is to find the triangle that have the maximum area, whose vertices are from the given points.
Input
The input consists of several test cases. The first line of each test case contains an integer n, indicating the number of points on the plane. Each of the following n lines contains two integer xi and yi, indicating the ith points. The last line of the input is an integer −1, indicating the end of input, which should not be processed. You may assume that 1 <= n <= 50000 and −10 4 <= xi, yi <= 10 4 for all i = 1 . . . n.
Output
For each test case, print a line containing the maximum area, which contains two digits after the decimal point. You may assume that there is always an answer which is greater than zero.
Sample Input
33 42 62 752 63 92 08 06 5-1Sample Output
0.5027.00题目题意:题目给了我们n个点,问这n个点中三个点能构成三角形最大面积是多少?
题目分析:我们以这n个点构成一个二维凸包,那么这三个点一定在上面,因为它们是最外面的点,现在就是在这个凸包上找到三个点求最大面积。
挑战设计竞赛上有种方法叫旋转卡壳法。
假定有 p q,r 三个点,我们先固定p,q,去求出使得三角形面积最大的点r,然后固定p,r,求出最大面积时的q,p就是遍历凸包所有的点
代码如下:
#include<iostream>#include<cstdio>#include<cmath>#include<algorithm>using namespace std;const int MAXN=50000+500;struct point{ int x,y;};point list[MAXN];//凸包的点int stack[MAXN],top;//stack保存的是最后凸包的点的序号double cross(point p0,point p1,point p2) //计算叉积 p0p1 X p0p2{ return (double ) (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);}double dis(point p1,point p2) //计算 p1p2的 距离{ return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));}bool cmp(point p1,point p2) //极角排序函数 , 角度相同则距离小的在前面{ int tmp=cross(list[0],p1,p2); if(tmp>0) return true; else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true; else return false;}void init(int n) //输入,并把 最左下方的点放在 list[0] 。并且进行极角排序{ int i,k; point p0; scanf("%d%d",&list[0].x,&list[0].y); p0.x=list[0].x; p0.y=list[0].y; k=0; for(i=1;i<n;i++) { scanf("%d%d",&list[i].x,&list[i].y); if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) ) { p0.x=list[i].x; p0.y=list[i].y; k=i; } } list[k]=list[0]; list[0]=p0; sort(list+1,list+n,cmp);}void graham(int n){ int i; if(n==1) {top=0;stack[0]=0;} if(n==2) { top=1; stack[0]=0; stack[1]=1; } if(n>2) { for(i=0;i<=1;i++) stack[i]=i; top=1; for(i=2;i<n;i++) { while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<=0) top--; top++; stack[top]=i; } }}int main(){ int n; while (scanf("%d",&n)!=EOF) { if (n==-1) break; init(n); graham(n); double ans=0; int p=1,q=2; stack[++top]=stack[0]; for (int i=0;i<top;i++) { while (cross(list[stack[i]],list[stack[p]],list[stack[q+1]])>cross(list[stack[i]],list[stack[p]],list[stack[q]])) q=(q+1)%top; ans=max(ans,fabs(cross(list[stack[i]],list[stack[p]],list[stack[q]]))/2.0); while (cross(list[stack[i]],list[stack[p+1]],list[stack[q]])>cross(list[stack[i]],list[stack[p]],list[stack[q]])) p=(p+1)%top; ans=max(ans,fabs(cross(list[stack[i]],list[stack[p]],list[stack[q]]))/2.0); } printf("%.2f\n",ans); } return 0;}
- POJ 2079(计算几何+凸包+旋转卡壳法)
- poj 2187 Beauty Contest 计算几何(凸包+旋转卡壳法)/最远点对
- 【转载】[Poj 2187]计算几何之凸包(三) {旋转卡壳初步}
- 【BZOJ1185】最小矩形覆盖 计算几何 凸包 旋转卡壳
- 计算几何之凸包(三) {旋转卡壳初步}
- [计算几何][凸包][旋转卡壳] 最远距离点对
- POJ 2187(计算几何+旋转卡壳法或暴力)
- POJ 2187 计算几何之旋转卡壳
- POJ 2079 Triangle 凸包+旋转卡壳
- poj 2079 Triangle 凸包+旋转卡壳
- <模板><计算几何>凸包-旋转卡壳最小面积外接矩形
- ZOJ 2419-- Triangle-凸包+旋转卡壳求最大面积三角形(计算几何)
- 计算几何之二维凸包:卷包裹算法、Graham Scan Algorithm、旋转卡壳算法
- HDU 5251 矩形面积 (计算几何+旋转卡壳求覆盖凸包的最小矩形面积)
- poj 2079 Triangle(凸包+旋转卡壳)
- POJ 2079 Triangle (凸包中的最大三角形&旋转卡壳)
- POJ 2079 Triangle(凸包-旋转卡壳)
- poj 2079 Triangle (二维凸包旋转卡壳)
- Spark Streaming缓存、Checkpoint机制
- CVE-2012-0158
- 开始-最初的梦想
- 【笔记】《WebGL编程指南》学习-第2章WebGL入门(6-改变点的颜色))
- Immutable.js与React和Redux实践
- POJ 2079(计算几何+凸包+旋转卡壳法)
- HDU 6215 Brute Force Sorting【双端链表】
- 详谈PHP垃圾回收机制
- 地球是平的?!
- < Android 源码分析(0)> Android 8.0 源码下载编译
- CVE-2012-4792
- linux下的yum环境的搭建问题CR
- 文章标题
- linux网络编程