洛谷P4724 【模板】三维凸包
题面
题解
先理一下关于立体几何的基本芝士好了……顺便全都是从\(xzy\)巨巨的博客上抄来的
加减
三维向量加减和二维向量一样
模长
\(|a|=\sqrt{x^2+y^2+z^2}\)
点积
两个向量\(a,b\)的点积还是代表\(a\)在\(b\)上的投影长\(\times b\)的模长,也有$a\cdot b=|a||b|\cos<a,b> \(,坐标上有\)(x_1,y_1,z_1)\cdot (x_2,y_2,z_2)=(x_1x_2,y_1y_2,z_1z_2)$
叉积
是一个三维向量,其模长表示这两个三维向量形成的平行四边形的面积,坐标表示为\((y_1z_2-z_1y_2,z_1x_2-x_1z_2,x_1y_2-y_1x_2)\)
关于叉积的方向,借(dao)用(yong)一下\(xzy\)巨巨的图
对于\(a\times b\),伸出右手,食指指向\(a\),中指指向\(b\),大拇指所对的方向为叉积后的向量方向,如图中\(AD\)即为\(AC\times AB\)
法向量
指和一个平面垂直的向量,记为\(n\),显然\(n\)有无数个
任取平面上两个向量做叉积即可
判断点\(P\)是否在平面上
任取平面上一点\(A\),如果\(n\)和\(AP\)点积为\(0\)说明\(P\)在该平面上
点\(P\)到平面的距离
任取平面上一点\(A\),那么这个距离实际上就等于\(AP\)在\(n\)上的投影的长度,那么\(AP\)点积\(n\)再除以\(|n|\)即可
求凸包
为了避免四点共面的情况先扰动一下,假装现在每一个面最多只有三个点
维护凸包的时候我们维护凸包的每一个面,把这个面的所有点按逆时针排序
假设我们现在要新增一个点\(P\),那么从\(P\)向之前的凸包做投影
然后把所有能看得到的面都删掉,加入新的面就行了
复杂度\(O(n^2)\)
ps:关于复杂度的计算,首先有欧拉公式\(V-E+F=2\),其中\(V\)为点数,\(E\)为边数,\(F\)为面数,凸包上最多有\(n\)个点,而且易发现有\(3F=2E\),解得\(F=2n-4,E=3n-6\)。增量法的复杂度是\(V\times E=O(n^2)\)
//minamoto
#include<bits/stdc++.h>
#define R register
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
double readdb()
{
R double x=0,y=0.1,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(x=ch-'0';(ch=getc())>='0'&&ch<='9';x=x*10+ch-'0');
for(ch=='.'&&(ch=getc());ch>='0'&&ch<='9';x+=(ch-'0')*y,y*=0.1,ch=getc());
return x*f;
}
const int N=2005;const double eps=1e-8;
inline double Rd(){return 1.0*rand()/RAND_MAX;}
inline double reps(){return (Rd()-0.5)*eps;}
struct point{
double x,y,z;
point(){}
point(R double xx,R double yy,R double zz):x(xx),y(yy),z(zz){}
inline void init(){x=readdb()+reps(),y=readdb()+reps(),z=readdb()+reps();}
inline point operator -(const point &b)const{return point(x-b.x,y-b.y,z-b.z);}
inline point operator *(const point &b)const{return point(y*b.z-z*b.y,z*b.x-x*b.z,x*b.y-y*b.x);}
inline double operator ^(const point &b)const{return x*b.x+y*b.y+z*b.z;}
inline double norm()const{return sqrt(x*x+y*y+z*z);}
}p[N];
struct node{
int id[3];
node(){}
node(R int x,R int y,R int z){id[0]=x,id[1]=y,id[2]=z;}
inline point normal(){return (p[id[1]]-p[id[0]])*(p[id[2]]-p[id[0]]);}
inline double area(){return normal().norm()/2;}
inline bool ck(const point &b){return ((b-p[id[0]])^normal())>0;}
}f[N*3],st[N*3];
int vis[N][N],n,top,cnt,v;double res;
void Convex_3D(){
f[++cnt]=node(1,2,3),f[++cnt]=node(3,2,1);
fp(i,4,n){
top=0;
fp(j,1,cnt){
v=f[j].ck(p[i]),!v?(st[++top]=f[j],0):0;
vis[f[j].id[0]][f[j].id[1]]=
vis[f[j].id[1]][f[j].id[2]]=
vis[f[j].id[2]][f[j].id[0]]=v;
}
fp(j,1,cnt){
if(vis[f[j].id[0]][f[j].id[1]]&&!vis[f[j].id[1]][f[j].id[0]])st[++top]=node(i,f[j].id[0],f[j].id[1]);
if(vis[f[j].id[1]][f[j].id[2]]&&!vis[f[j].id[2]][f[j].id[1]])st[++top]=node(i,f[j].id[1],f[j].id[2]);
if(vis[f[j].id[2]][f[j].id[0]]&&!vis[f[j].id[0]][f[j].id[2]])st[++top]=node(i,f[j].id[2],f[j].id[0]);
}
fp(j,1,top)f[j]=st[j];cnt=top;
}
fp(i,1,cnt)res+=f[i].area();
printf("%.3lf\n",res);
}
int main(){
// freopen("testdata.in","r",stdin);
n=read();
fp(i,1,n)p[i].init();
Convex_3D();
return 0;
}
洛谷P4724 【模板】三维凸包的更多相关文章
- luogu P4724 模板 三维凸包
LINK:三维凸包 一个非常古老的知识点.估计也没啥用. 大体上了解了过程 能背下来就背下来吧. 一个bf:暴力枚举三个点 此时只需要判断所有的点都在这个面的另外一侧就可以说明这个面是三维凸包上的面了 ...
- 题解-洛谷P4724 【模板】三维凸包
洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...
- 洛谷P3373 [模板]线段树 2(区间增减.乘 区间求和)
To 洛谷.3373 [模板]线段树2 题目描述 如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.将某区间每一个数乘上x 3.求出某区间每一个数的和 输入输出格式 输入格 ...
- 洛谷P3375 [模板]KMP字符串匹配
To 洛谷.3375 KMP字符串匹配 题目描述 如题,给出两个字符串s1和s2,其中s2为s1的子串,求出s2在s1中所有出现的位置. 为了减少骗分的情况,接下来还要输出子串的前缀数组next.如果 ...
- LCT总结——概念篇+洛谷P3690[模板]Link Cut Tree(动态树)(LCT,Splay)
为了优化体验(其实是强迫症),蒟蒻把总结拆成了两篇,方便不同学习阶段的Dalao们切换. LCT总结--应用篇戳这里 概念.性质简述 首先介绍一下链剖分的概念(感谢laofu的讲课) 链剖分,是指一类 ...
- 【AC自动机】洛谷三道模板题
[题目链接] https://www.luogu.org/problem/P3808 [题意] 给定n个模式串和1个文本串,求有多少个模式串在文本串里出现过. [题解] 不再介绍基础知识了,就是裸的模 ...
- 洛谷-P5357-【模板】AC自动机(二次加强版)
题目传送门 -------------------------------------- 过年在家无聊补一下这周做的几道AC自动机的模板题 sol:AC自动机,还是要解决跳fail边产生的重复访问,但 ...
- 洛谷.1919.[模板]A*B Problem升级版(FFT)
题目链接:洛谷.BZOJ2179 //将乘数拆成 a0*10^n + a1*10^(n-1) + ... + a_n-1的形式 //可以发现多项式乘法就模拟了竖式乘法 所以用FFT即可 注意处理进位 ...
- 洛谷.3803.[模板]多项式乘法(FFT)
题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...
- 洛谷.3803.[模板]多项式乘法(NTT)
题目链接:洛谷.LOJ. 为什么和那些差那么多啊.. 在这里记一下原根 Definition 阶 若\(a,p\)互质,且\(p>1\),我们称使\(a^n\equiv 1\ (mod\ p)\ ...
随机推荐
- 安装glibc-2.14
下载glibc-2.14.tar.gz 解压: [root@jrgc130 software]# tar xf glibc-2.14.tar.gz [root@jrgc130 software]# c ...
- Eclipse修改XML默认打开方式
用Eclipse开发Android的时候 默认的XML是采用Android xml editor 打开,这个工具不够直观,如果想直接看文本的XML的话,可以通过如下方式修改 1.菜单:Window ...
- MySQL备份还原之三使用xtrabackup
1 xtrabackup安装 1)解压源码包 tar -xzvf percona-xtrabackup-2.1.7.tar.gz 2)安装perl环境(DBI/DBD) yum install per ...
- ssh免密连接远程服务器
ssh免密连接远程服务器 借助ssky-keygen和ssh-copy-id工具,通过4个简单的步骤实现无需输入密码登录远程Linux主机 1 生成密钥 通过内置的工具生成RSA算法加密的密钥 ssh ...
- Pacemaker实现双机热备
在互联网高速发展的今天,尤其在电子商务的发展,要求服务器能够提供不间断服务.在电子商务中,如果服务器宕机,造成的损失是不可估量的.要保证服务器不间断服务,就需要对服务器实现冗余.在众多的实现服务器冗余 ...
- iOS静态库的制作与引用
[iOS静态库的制作与引用] 1.Configuring Exported Headers To configure which headers are exported to clients, se ...
- Java IO简介
-------------siwuxie095 Java IO简介: IO 也写作"I/O",可理解为 In 和 Out,即 输入 与 输出 ...
- win10 Kinect2 Visualstudio2015 opencv3环境搭建
1.下载kinect SDK ( Kinect for Windows SDK 2.0 ): https://www.microsoft.com/en-us/download/details.asp ...
- 409. Longest Palindrome 最长对称串
[抄题]: Given a string which consists of lowercase or uppercase letters, find the length of the longes ...
- Mpich编程
一.简介 通过安装MPICH构建MPI编程环境,从而进行并行程序的开发.MPICH是MPI(Message-Passing Interface)的一个应用实现,支持最新的MPI-2接口标准,是用于并行 ...