生活随笔
收集整理的這篇文章主要介紹了
poj 2079(旋转卡壳)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題意:求出平面內的點集所組成的面積最大的三角形。
解題思路:考慮凸包+旋轉卡殼。面積最大的三角形的三點必定在凸包的頂點上,只不過這里要注意,三角形的邊不一定就是凸包的邊,有可能三角形相鄰兩點是橫跨凸包的。
關鍵是如何找三個頂點。這里采用的類似于尺取法,先固定一個點,剩下的兩個指針依次逆時針方向旋轉,找到最大的面積。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;const int maxn = 50000;
const double eps = 1e-8;
struct Point
{double x,y;
}p[maxn],Stack[maxn];
int n,top;double Cross(Point a,Point b,Point c)
{ return (a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y);
} double dis(Point a,Point b)
{ return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
} int cmp(Point a,Point b)
{ if(Cross(a,b,p[0]) > eps) return 1; if(fabs(Cross(a,b,p[0])) < eps && dis(b,p[0]) - dis(a,p[0]) < eps) return 1; return 0;
}int Graham()
{sort(p+1,p+n,cmp);top = 2; Stack[0] = p[0]; Stack[1] = p[1]; Stack[2] = p[2];for(int i = 3; i < n; i++){ while(top >= 1 && Cross(p[i],Stack[top],Stack[top-1]) > eps){ top--; } Stack[++top]=p[i]; } return top;
}double rotating_calipers()//旋轉卡殼
{int i,j=1,k=2;double ans=0;Stack[++top]=Stack[0];for(i=0; i<top; i++){/*邊上的兩點可以不是凸包上的點,WA無數次*/while(fabs(Cross(Stack[(k+1)%top],Stack[i],Stack[j])) > fabs(Cross(Stack[k],Stack[i],Stack[j])))k=(k+1)%top;while(fabs(Cross(Stack[k],Stack[i],Stack[(j+1)%top])) > fabs(Cross(Stack[k],Stack[i],Stack[j])))j=(j+1)%top;ans=max(ans,fabs(Cross(Stack[k],Stack[i],Stack[j])));}return ans;
}int main()
{while(scanf("%d",&n),n != -1){for(int i = 0; i < n; i++)scanf("%lf%lf",&p[i].x,&p[i].y);int k = 0; for(int i = 1;i < n; i++){ if(p[k].y > p[i].y || (p[k].y == p[i].y) && (p[k].x > p[i].x)){ k=i; } } swap(p[0],p[k]);Graham();printf("%.2lf\n",rotating_calipers()/2);}return 0;
}
總結
以上是生活随笔為你收集整理的poj 2079(旋转卡壳)的全部內容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。