标签:return Point CThirdConvexHull ret 凸包 算法 convex Total
ThirdConvexHull.h
#pragma once
class CThirdConvexHull
{
public:
CThirdConvexHull();
~CThirdConvexHull();
public:
static AcGePoint3dArray testConvexHullPoints(AcGePoint3dArray& pts);
};
ThirdConvexHull.cpp
#include "stdafx.h"
#include "ThirdConvexHull.h"
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>
using std::sort;
#define MAXN 1002
int N, L;
double sqr(double a)
{
return a*a;
}
struct Point
{
double x, y;
inline Point operator- (const Point &t)
{
Point ret;
ret.x = x - t.x;
ret.y = y - t.y;
return ret;
}
inline Point operator+ (const Point &t)
{
Point ret;
ret.x = x + t.x;
ret.y = y + t.y;
return ret;
}
inline int det(const Point &t)
{
return x*t.y - t.x*y;
}
inline double dist(Point &t)
{
return sqrt(sqr(x - t.x) + sqr(y - t.y));
}
}d[MAXN];
bool cmpPoint(const Point &a, const Point &b)
{
if (a.x != b.x) return a.x < b.x;
return a.y < b.y;
}
int convex[MAXN], cTotal;
void get_convex_hull()
{
sort(d, d + N, cmpPoint);
int Total = 0, tmp;
for (int i = 0; i < N; ++i)
{
while ((Total > 1) &&
((d[convex[Total - 1]] - d[convex[Total - 2]]).det(
d[i] - d[convex[Total - 1]]) <= 0)) Total--;
convex[Total++] = i;
}
tmp = Total;
for (int i = N - 2; i >= 0; --i)
{
while ((Total > tmp) &&
((d[convex[Total - 1]] - d[convex[Total - 2]]).det(
d[i] - d[convex[Total - 1]]) <= 0)) Total--;
convex[Total++] = i;
}
cTotal = Total;
}
CThirdConvexHull::CThirdConvexHull()
{
}
CThirdConvexHull::~CThirdConvexHull()
{
}
AcGePoint3dArray CThirdConvexHull::testConvexHullPoints(AcGePoint3dArray& pts)
{
AcGePoint3dArray outpts;
N = pts.length();
for (int i = 0; i < N; ++i)
{
d[i].x = pts[i].x;
d[i].y = pts[i].y;
}
get_convex_hull();
double Ans = 0;
for (int i = 0; i < cTotal; ++i)
{
AcGePoint3d tempt;
tempt.x = d[convex[i]].x;
tempt.y = d[convex[i]].y;
tempt.z = 0;
outpts.append(tempt);
//acutPrintf(_T("\n%2f,%2f,0"), d[convex[i]].x, d[convex[i]].y);
}
return outpts;
}
标签:return,Point,CThirdConvexHull,ret,凸包,算法,convex,Total 来源: https://www.cnblogs.com/Pond-ZZC/p/12446971.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。