凸包 基础练习例子(一)

来源:互联网 发布:男人床上粗鲁 知乎 编辑:程序博客网 时间:2024/06/06 04:45
      有一个牧场,牧场上有很多个供水装置,现在牧场的主人想要用篱笆把这些供水装置圈起来,以防止不是自己的牲畜来喝水,各个水池都标有各自的坐标,现在要你写一个程序利用最短的篱笆将这些供水装置圈起来!(篱笆足够多,并且长度可变)
输入
      第一行输入的是N,代表用N组测试数据(1<=N<=10),第二行输入的是m,代表本组测试数据共有m个供水装置(3<=m<=100)接下来m行代表的是各个供水装置的横纵坐标
输出
     输出各个篱笆经过各个供水装置的坐标点,并且按照x轴坐标值从小到大输出,如果x轴坐标值相同,再安照y轴坐标值从小到大输出
样例输入
1
4
0 0
1 1
2 3
3 0
样例输出
0 0
2 3
3 0
 
 
program:
 
 

#include<iostream>

#include<algorithm>

#include<cmath>

using namespace std;

 

struct point

{  

    int x,y;

};

point p[105],res[105];

double Dist(const point &arg1, const point &arg2)

{

    return sqrt( (arg1.x - arg2.x)*(arg1.x - arg2.x) + (arg1.y - arg2.y)*(arg1.y - arg2.y) );

}

bool multi(point p0,point p1,point p2)

{  

    return (p1.x-p0.x)*(p2.y-p0.y)>(p2.x-p0.x)*(p1.y-p0.y);

}

int mysort1(point a,point b)

{   if(a.y!=b.y) return a.y<b.y;

    if(a.y==b.y&&a.x!=b.x) return a.x<b.x;

}

bool cmp(const point &a,const point &b)

{

    point temp=p[0];

    double xmt=(a.x-temp.x)*(b.y-temp.y)-(b.x-temp.x)*(a.y-temp.y);

    if(xmt)                             //向量不共线就按逆时针旋转

        return xmt>0;

    return Dist(a,temp)>Dist(b,temp);//向量共线取最长的。

}

 

int main()

{   int n,k,len;

    cin>>k;

    

    while(k--)

    {

      cin>>n;

      for(int i=0;i<n;i++)

         cin>>p[i].x>>p[i].y;

      sort(p,p+n,mysort1);//排序,找到最左下角的点

       

      res[0]=p[0];

 

      sort(p+1,p+n,cmp);//按照极角排序

      res[1]=p[1];

      res[2]=p[2];

      int top=2;

      for(int i=3;i<n;i++)

      { 

            while(multi(p[i],res[top],res[top-1]))

                top--;

            res[++top]=p[i];

      }

      

      for(int i=0;i<=top;i++)

         cout<<res[i].x<<" "<<res[i].y<<endl;

    }

    

    system("pause");

    return 0;

}

原创粉丝点击