三维凸包之增量算法

最新推荐文章于 2024-07-20 15:19:53 发布

theArcticOcean

最新推荐文章于 2024-07-20 15:19:53 发布

阅读量4.6k

点赞数 2

分类专栏： algorithm_计算几何文章标签：凸包

本文链接：https://blog.csdn.net/thearcticocean/article/details/50646428

版权

本文深入探讨了三维空间中凸包的增量算法，详细阐述了算法的工作流程，揭示了如何通过逐步添加点来构建凸包，适用于理解和实现高效的空间几何计算。

摘要由CSDN通过智能技术生成

算法工作过程：

选取点集中的四个点构成一个初始的四面体：

选两个点p1, p2，再选一个点p3，p3和p1, p2不在同一条直线上，这样形成了一个面。再选一个点，不在同一个面上，形成一个四面体。

然后继续研究剩余的点：

点在四面体内部，略过。

点在四面体外部，删除此点能“看见”的面，扩大凸包体积。

POJ 3528 Ultimate Weapon

http://poj.org/problem?id=3528

大意：求解多面体表面积

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const double eps=1e-7;
const int N=5e2+10;
struct point{
    double x,y,z;
    point(){ x=0; y=0; z=0; }
    point(double _x,double _y,double _z):x(_x),y(_y),z(_z) {}
    point operator -(const point other){
        return point(x-other.x,y-other.y,z-other.z);
    }
    point operator *(const point other){
        return point(y*other.z-z*other.y,-x*other.z+z*other.x,x*other.y-y*other.x);
    }
    double operator ^(const point other){ //点乘
        return x*other.x+y*other.y+z*other.z;
    }
};
struct ch3D{
     struct fac{
         int a,b,c;  // 面上三点的编号
         bool ok;   // 是否属于凸包
     };
     int n;  //初始点数
     point P[N];  //初始点

     int cnt;  //凸包上的三角形
     fac F[N*6];  // 被创建的面不超过6N

     int vis[N][N];
     double vlen(point v){
         return sqrt(v.x*v.x+v.y*v.y+v.z*v.z);
     }
     double area(point p1,point p2,point p3){
         return vlen((p2-p1)*(p3-p1));   // 三角形面积*2
     }
     double volumn(point p1,point p2,point p3,point p4){
         return fabs((p2-p1)*(p3-p1)^(p4-p1));  //V的6倍
     }
     double pToPlane(point p,fac ff){  // 点在面的同向: +
         point m=P[ff.b]-P[ff.a];
         point n=P[ff.c]-P[ff.a];
         point v=p-P[ff.a];
         return m*n^v;
     }
     void deal(int p,int a,int b){
         int f=vis[a][b];
         fac add;
         if(F[f].ok){
            if(pToPlane(P[p],F[f])>ep