POJ 2187(求凸包內兩點之間最長距離)

推薦博客 :借鑑

/*

分別列舉每個凸包的每條邊,再列舉另一個凸包的點,求得一對對踵點,與當前記錄的值相取捨。
遍歷到下一條邊,從上次便利到的點開始尋找下一對踵點,再取捨......。直至取捨出最長、最短距離。
*/
#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int maxn=50001;
struct node
{
    double  x,y;
}st[maxn],num[maxn];
int n,k,e;
double dis(node p1,node p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
int multi(node p1,node p2,node p3)//叉積
{
    return (p1.x-p3.x)*(p2.y-p3.y)-(p2.x-p3.x)*(p1.y-p3.y);
}
bool cmp(node p1,node p2)
{
    if(multi(p1,p2,num[0])>0) return true;
    if(multi(p1,p2,num[0])==0&&dis(p1,num[0])<dis(p2,num[0]))return true;
    return false;
}
void Graham()//凸包
{
    e=2;
    sort(num+1,num+n,cmp);
    st[0]=num[0];
    st[1]=num[1];
    st[2]=num[2];
    for(int i=3;i<n;i++)
    {
        while(e>1&&multi(num[i],st[e],st[e-1])>=0)
            e--;
        st[++e]=num[i];
    }
}
double rotating_calipers()
{
    e++;
    int q=1;
    double ans=0;
    st[e]=st[0];
    for(int p=0;p<e;p++)
    {
        while(multi(st[p+1],st[q+1],st[p])>multi(st[p+1],st[q],st[p]))
            q=(q+1)%e;
        ans=max(ans,max(dis(st[p],st[q]),dis(st[p+1],st[q+1])));
    }
    return ans;
}
int main()
{
    scanf("%d",&n);
    k=0;//K表示起始點
    for(int i=0;i<n;i++){
        scanf("%lf %lf",&num[i].x,&num[i].y);
        if(num[i].x<num[k].x||(num[i].x==num[k].x&&num[i].y<num[k].y)) k=i;
    }
    swap(num[0],num[k]);
    //printf("%d\n",k);
    Graham();
    printf("%.0lf\n",rotating_calipers());
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章