传送门

三维凸包模板题……只是听了听计算几何的课之后心血来潮想写的……

我的做法很无脑是吧……暴力枚举三个点组成的三角形,然后枚举剩下的点,判断其余点是否都在这个三角形的同一侧,是的话则说明这个三角形是凸包的一个面。

理论复杂度应该是$O(n^4)$,不过看上去跑得飞快?人帅自带小常数哈哈

这个故事告诉我们:大力出奇迹……

 #include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<fstream>
#include<iomanip>
using namespace std;
const int maxn=;
const long double eps=1e-;
struct Point{
long double x,y,z;
Point(long double x=0.0,long double y=0.0,long double z=0.0):x(x),y(y),z(z){}
Point operator-(const Point &a)const{return Point(a.x-x,a.y-y,a.z-z);}//A - B = B到A的位移
Point operator/(const long double &a)const{return Point(x/a,y/a,z/a);}
}a[maxn];
typedef Point Vector;
long double noise();
long double Dot(const Vector&,const Vector&);
Vector Cross(const Vector&,const Vector&);
long double Length(const Vector&);
long double Area2(const Point&,const Point&,const Point&);
Vector LawVector(const Vector&,const Vector&);
long double Distance(const Point&,const Point&,const Vector&);
ifstream fin("enwrap.in");
ofstream fout("enwrap.out");
long double ans=0.0;
int n,t;
int main(){
fin>>n;
for(int i=;i<=n;i++){
fin>>a[i].x>>a[i].y>>a[i].z;
a[i].x+=noise();
a[i].y+=noise();
a[i].z+=noise();
}
Vector A;
long double d;
bool bad;
for(int i=;i<=n;i++)for(int j=;j<i;j++)for(int k=;k<j;k++){
bad=false;
t=;
A=LawVector(a[j]-a[i],a[k]-a[i]);
for(int l=;l<=n;l++)if(i!=l&&j!=l&&k!=l){
d=Distance(a[l],a[i],A);
if(!t)t=(d<?-:);
else if(t*d<){
bad=true;
break;
}
}
if(!bad)ans+=Area2(a[i],a[j],a[k]);
}
ans/=2.0;
fout<<fixed<<ans;
return ;
}
long double noise(){
static int a=,b=,p=,x=;
x=a*x+b;x%=p;
if(x<)x+=p;
return (long double)(x-p)/p*5e-;
}
inline long double Dot(const Vector &A,const Vector &B){return A.x*B.x+A.y*B.y+A.z*B.z;}
inline Vector Cross(const Vector &A,const Vector &B){return Vector(A.y*B.z-B.y*A.z,A.z*B.x-B.z*A.x,A.x*B.y-B.x*A.y);}
inline long double Length(const Vector &A){return sqrt(Dot(A,A));}
inline long double Area2(const Point &A,const Point &B,const Point &C){return Length(Cross(B-A,C-A));}
inline Vector LawVector(const Vector &A,const Vector &B){
Vector n=Cross(A,B);
return n/Length(n);
}//两个向量的叉积一定同时垂直于这两个向量
inline long double Distance(const Point &A,const Point &P,const Vector &n){return Dot(A-P,n);}//(P,n)是平面的点法式,n是单位向量
/*
三维凸包——暴力法
暴力枚举三个点组成的三角形,
判断其他点是否都在这个三角形的同侧,
是则说明这个三角形一定是凸包的一个面。
*/

话说一开始忘了怎么求平面的法向量,手推了个公式然后发现搞出NaN了……因为我推的公式默认法向量在z维的长度不为0,但其实有很多平面的法向量是垂直于z轴的,然后就除零爆炸了……又脑补了很久,后来才想起来两个向量的叉积必定同时垂直于这两个向量……这人没救了

还有一件事,这题卡精度,随机扰乱搞得太大会炸精度……计算几何毁我青春

Hnoi2004 金属包裹的更多相关文章

  1. bzoj 1209: [HNOI2004]最佳包裹 三维凸包

    1209: [HNOI2004]最佳包裹 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 160  Solved: 58[Submit][Status] ...

  2. BZOJ1209 [HNOI2004]最佳包裹 三维凸包 计算几何

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ1209 题目概括 给出立体的n个点.求三维凸包面积. 题解 增量法,看了一天,还是没有完全懂. 上板 ...

  3. 洛谷P2287 [HNOI2004]最佳包裹(三维凸包)

    题面 传送门 题解 左转板子,调个精度就能\(A\)了 //minamoto #include<bits/stdc++.h> #define R register #define fp(i ...

  4. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  5. BZOJ 1212: [HNOI2004]L语言 [AC自动机 DP]

    1212: [HNOI2004]L语言 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1367  Solved: 598[Submit][Status ...

  6. 如何让include标签包裹的布局置于屏幕最下方?

    如何让一个Layout 始终在屏幕的下方 我想让<include layout="@layout/bottom" />一直在屏幕下,怎么做? 1.相对布局中用属性  a ...

  7. float包裹性与破坏性及清除浮动几种方法

    1. “清除浮动” ??准确的描述应该是“清除浮动造成的影响”  学习视频录制推荐的教程<CSS深入理解之float浮动> 2.如何清除浮动造成的影响??? 栗子 块级div元素包含一个内 ...

  8. Vijos1450 包裹快递[二分答案]

    背景 小K成功地破解了密文.但是乘车到X国的时候,发现钱包被偷了,于是无奈之下只好作快递员来攒足路费去Orz教主…… 描述 一个快递公司要将n个包裹分别送到n个地方,并分配给邮递员小K一个事先设定好的 ...

  9. python--基础学习(五)参数位置传递、关键字传递、包裹传递及解包裹

    python系列均基于python3.4环境 1.位置传递和关键字传递 代码示例 #位置传递 def fun(a,b,c): print("a: {0}, b: {1}, c: {2}&qu ...

随机推荐

  1. P4097 [HEOI2013]Segment 李超线段树

    $ \color{#0066ff}{ 题目描述 }$ 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 i 条被插入的线段的标号为 i 给定一个数 k,询问与直线 x = k 相交的线 ...

  2. P3267 [JLOI2016/SHOI2016]侦察守卫

    $ \color{#0066ff}{ 题目描述 }$ 小R和B神正在玩一款游戏.这款游戏的地图由N个点和N-1条无向边组成,每条无向边连接两个点,且地图是连通的.换句话说,游戏的地图是一棵有N个节点的 ...

  3. 用 gdb 和 qemu 调试 grub

    因为qemu内置了gdbserver,所以我们可以用gdb调试qemu虚拟机上执行的代码,而且不受客户机系统限制. 以下内容是我调试 grub 0.97 时的一份笔记. 准备 qemu, gdb,以及 ...

  4. PHP类库生成pdf代码实例

    require_once('./tcpdf/tcpdf.php'); //引入库文件        $pdf = new TCPDF('P', 'mm', 'A4', true, 'UTF-8', f ...

  5. springAOP实现方法运行时间统计

    aop的before和after,寻思分别在两个方法里获得当前时间,最后一减就可以了. 因此,今天就探讨了一下这个问题,和大家分享一下. 创建maven工程,引入spring的依赖包,配置好appli ...

  6. 1016 部分A+B (15 分)

    // 题目不难,感觉变量定义的有点儿多了#include <iostream> #include <cmath> using namespace std; int main() ...

  7. python 小练习('''产生了多少个对象''')('''模拟两个英雄,昵称,攻击力,生命值''')

    '''小练习 计数器,产生了多少个对象''' # class Student: # school = 'luffycity' # count = 0 # def __init__(self, name ...

  8. vue-cli 启动过项目步骤

    一. 安装 node.js 安装完成后,可以命令行工具中输入 node -v 和 npm -v,如果能显示出版本号,就说明安装成功. 二.安装webpack npm install webpack - ...

  9. es第三篇:Search APIs

    大多数search API都是可以操作多个索引的,除了explain API. 当执行一个search API时,可以指定routing参数,去搜索特定的主分片及其副本分片.routing参数值可以是 ...

  10. windows 下创建 sqlite 数据库

    说明:windows 下执行创建 sqlite 数据库命令后数据库文件不会马上生成,需要创建表以后才会生成. 1.将 sqlite3.exe 文件放在任何位置(如放在 d:\tools )2.在 CM ...