POJ 1228 Grandpa’s Estate(給定凸包上的一部分頂點 判斷表示的凸包是否唯一)

來源:互聯網
上載者:User

題目連結:http://poj.org/problem?id=1228

這個題目不好想,開始老是wa,後來不知道改了哪裡就AC了

這個題目從範例可以看出給的凸包山的點不是按照順序來的,是隨便的一些點

那麼我們就在這些點上建立凸包,那麼在建立凸包的時候我們不能隨便去掉邊

也就是即使是可以去除的點也不能去掉,也就是傳統的graphm方法不能預先選擇

結束之後要想這個凸包表示唯一那麼就不能在外部加任何一個點使其還能長成

一個更大的凸包,這個就要求一個頂點與其前面兩個頂點或者後面兩個頂點至少

有邊能共線,否則你就可以在外面插入一個頂點增加凸包,滿足則不能!

#include <iostream>#include <stdio.h>#include <string.h>#include <algorithm>#include <cmath>using namespace std;#define eps 1e-6#define PI 3.14159265struct point{double x;double y;}po[1500],temp;int n,pos;bool zero(double a){return fabs(a) < eps;}double dis(point &a,point &b)//返回兩點之間距離的平方{return (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y);}double across(point &a,point &b,point &c)//求a b and a c 的X積{return (b.x-a.x)*(c.y-a.y) - (b.y-a.y)*(c.x-a.x);}int cmp(const void *a,const void *b){    int x;    x=across(po[0],*(point*)a,*(point*)b);    if(zero(x))    return dis(po[0],*(point*)a) < dis(po[0],*(point*)b) ? -1 : 1;return across(po[0],*(point*)a,*(point*)b) > 1e-8 ? -1 : 1;}int graham(int num){    bool flag=true;int i,j,k=2;po[num]=po[0];//fangbiannum++;for(i=3;i<num;i++){while(across(po[k-1],po[k],po[i]) < -eps){k--;}po[++k]=po[i];//就這個迴圈結束,不需要了!}n=k;k=0;for(i=0;i<n;i++)if(zero(across(po[i],po[i+1],po[i+2]))){    k=1;    continue;}else{    flag=false;    if(k==0)    {printf("NO\n");k=10;return 0;}    k=0;}if(flag){    printf("NO\n");  return 0;}if(k!=10)printf("YES\n");return 0;}int main(){int i,j,k,t;point my_temp,temp1;scanf("%d",&t);while(t--){    scanf("%d",&n);scanf("%lf%lf",&po[0].x,&po[0].y);temp=po[0];pos=0;for(i=1;i<n;i++){scanf("%lf%lf",&po[i].x,&po[i].y);if(po[i].y < temp.y)temp=po[i],pos=i;if(po[i].y==temp.y && po[i].x >temp.x){temp1=temp; temp=po[i],pos=i;}}temp=temp1;if(n<6){    printf("NO\n");    continue;}my_temp=po[0];po[0]=po[pos];po[pos]=my_temp;qsort(po+1,n-1,sizeof(po[0]),cmp);graham(n);}return 0;}

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.