洛谷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)\ ...
随机推荐
- SSD知识
不管什么接口的SSD,一般都由以下部分组成:主控,Flash,板,壳,品牌.下面本佬就这些部分一一发帖,仅供娱乐参考,不作任何推荐和偏向,有不同见解请直接发表,有任何错误,请直接指正,不为吵架,只为娱 ...
- Java之匿名内部类和包装类
匿名内部类 作用: 假如某个类只使用一次,则可以使用匿名内部类,无需再新建该类 我们上下代码: package com.learn.chap03.sec16; /** * 定义接口 */ public ...
- linux下编译wpa_supplicant …
linux下编译wpa_supplicant 收藏 前一段时间只在vs2005下编译成功过.经过近一段时间的琢磨,今天终于在linux下成功编译了wpa_supplicant. 挺简单的事情折腾了这么 ...
- 文件流的处理 C C++
C: 借鉴C 文件操作库函数总结 eg: #include <fstream> #include <iostream> #include <cstdio> usin ...
- sudo apt install libreadline-dev Reading package lists... Error!
luo@luo-ThinkPad-W540:~$ luo@luo-ThinkPad-W540:~$ luo@luo-ThinkPad-W540:~$ luo@luo-ThinkPad-W540:~$ ...
- SQL Server误区30日谈 第26天 SQL Server中存在真正的“事务嵌套”
误区 #26: SQL Server中存在真正的“事务嵌套”错误 嵌套事务可不会像其语法表现的那样看起来允许事务嵌套.我真不知道为什么有人会这样写代码,我唯一能够想到的就是某个哥们对SQL Serve ...
- 更改windows服务的配置文件(app.config)必须重启服务才能生效吗?
这个问题是前一阶段写windows服务碰到的.本来在写获取配置文件的某个配置的值的时候,通常我都是写类似下面的这么一个静态方法来获取: 1: /// <summary> 2: /// 获取 ...
- 29.MAX() 函数
MAX() 函数 MAX 函数返回一列中的最大值.NULL 值不包括在计算中. SQL MAX() 语法 SELECT MAX(column_name) FROM table_name 注释:MIN ...
- fiddler抓包时显示Tunnel to......443
打开手机浏览器,输入http://192.168.0.65:8888/FiddlerRoot.cer
- jqgrid扩展 获取表单数据
$.fn.GetPostData = function () { var data = {}; var k = false; $(this).find(".datacontrol" ...