凸包模板
来源:互联网 发布:php网盘系统源码 编辑:程序博客网 时间:2024/04/29 18:29
#include<iostream>
#include<cmath>
#include<stack>
using namespace std;
#define Max 1000
int multi(int p1,int p2,int p3);
int Point[Max][2] , a[Max][2] , n;
int multi(int p1,int p2,int p3)
{
int v;
int v1x=Point[p2][0]-Point[p1][0],
v1y=Point[p2][1]-Point[p1][1],
v2x=Point[p3][0]-Point[p1][0],
v2y=Point[p3][1]-Point[p1][1];
if(p1==p2) return p3;
if(p1==p3) return p2;
else{
v=v1x*v2y-v1y*v2x;
if(v>0) return p2;
if(v<0) return p3;
if(v==0){
if((v1x*v1x+v1y*v1y)>(v2x*v2x+v2y*v2y))
return p2;
else return p3;
}
}
}
int qiutubao(int point[Max][2],int &n,int a[Max][2])//point是输入的点集 求出的凸包存到a中
{ //返回 的值表示凸包上顶点个数
int i;
int m=Point[0][1],t=0;
for(int k=1;k<n;k++){
if (Point[k][1]<m)
{
m=Point[k][1];t=k;
}
else if(Point[k][1]==m)
{
if (Point[k][0]<Point[t][0])
{m=Point[k][1];t=k;}
}
}
Point[n][0] = Point[t][0];
Point[n][1] = Point[t][1];
int q=n;
stack<int> p;
p.push(t);
int kk=t,g=0;
while(true)
{
for(int i=1;i<=n;i++)
g= multi(kk,g,i);
if(g==q) break;
p.push(g); kk=g; g=0;
}
int get; i=0;
while(!p.empty()) //出栈并存到a中
{
get = p.top();
a[i][0]=Point[get][0];a[i][1]=Point[get][1];
i++; p.pop();
}
return i;
}
int main()
{
int i=0,l;
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%d%d",&Point[i][0],&Point[i][1]);
l=qiutubao(Point,n,a);
cout<<endl;
for(i=0;i<l;i++)
cout<<a[i][0]<<' '<<a[i][1]<<endl;
return 0;
}
- 凸包模板
- 1392 凸包模板
- 凸包模板
- 凸包模板题
- 更新凸包模板
- 三维凸包模板
- 凸包模板
- 三维凸包模板
- 凸包模板
- 凸包模板
- 三维凸包模板
- 三维凸包模板
- 凸包模板
- 凸包 andrew模板
- 凸包面积模板
- 模板,凸包
- 凸包模板
- 凸包模板!
- Struts2与Struts1的对比
- 来自百度MM
- 不要使用矿泉水瓶接热水喝——致癌
- Linux下的调试工具《转》
- Virtual Memory II: the return of objrmap
- 凸包模板
- JAVA的环境变量的配置
- 利用LUCENE求相似文档
- BS下自己~~
- 遍历二叉树的递归和非递归算法
- 动易0day--Region.asp注入
- sql server 2005安装安装不上问题,
- 深入了解CLR的加载过程
- 会议记录(7月26日——8月25日)