BZOJ1209 [HNOI2004]最佳包裹 三维凸包 计算几何
欢迎访问~原文出处——博客园-zhouzhendong
去博客园看该题解
题目传送门 - BZOJ1209
题目概括
给出立体的n个点。求三维凸包面积。
题解
增量法,看了一天,还是没有完全懂。
上板子!
代码
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
using namespace std;
const double Eps=1e-8;
const int N=100+5;
int Dcmp(double x){
if (fabs(x)<Eps)
return 0;
return x<0?-1:1;
}
struct Point{
double x,y,z;
Point (){}
Point (double x_,double y_,double z_){
x=x_,y=y_,z=z_;
}
}P[N],p[N];
Point operator - (Point a,Point b){
return Point(a.x-b.x,a.y-b.y,a.z-b.z);
}
Point operator + (Point a,Point b){
return Point(a.x+b.x,a.y+b.y,a.z+b.z);
}
Point operator * (Point a,double x){
return Point(a.x*x,a.y*x,a.z*x);
}
Point operator / (Point a,double x){
return Point(a.x/x,a.y/x,a.z/x);
}
Point cross(Point a,Point b){
return Point(a.y*b.z-a.z*b.y,a.z*b.x-a.x*b.z,a.x*b.y-a.y*b.x);
}
Point cross(Point a,Point b,Point c){
return cross(b-a,c-a);
}
double Dot(Point a,Point b){
return a.x*b.x+a.y*b.y+a.z*b.z;
}
double Length(Point P){
return sqrt(Dot(P,P));
}
double rand01(){
return (double)rand()/(double)RAND_MAX;
}
double randEps(){
return (rand01()-0.5)*Eps;
}
Point add_noise(Point P){
return Point(P.x+randEps(),P.y+randEps(),P.z+randEps());
}
struct Face{
int v[3];
Face(){}
Face(int a,int b,int c){
v[0]=a,v[1]=b,v[2]=c;
}
Point Normal(Point *P){
return cross(P[v[0]],P[v[1]],P[v[2]]);
}
int CanSee(Point *P,int i){
return Dot(P[i]-P[v[0]],Normal(P))>0?1:0;
}
};
struct ConvexHull{
vector <Face> cur;
bool vis[N][N];
void Increment(Point *P,int n){
vector <Face> next;
memset(vis,0,sizeof vis);
for (int i=1;i<=n;i++)
P[i]=add_noise(P[i]);
cur.clear();
cur.push_back(Face(1,2,3));
cur.push_back(Face(3,2,1));
for (int i=4;i<=n;i++){
next.clear();
for (int j=0;j<cur.size();j++){
Face F=cur[j];
int res=F.CanSee(P,i);
if (!res)
next.push_back(F);
for (int k=0;k<3;k++)
vis[F.v[k]][F.v[(k+1)%3]]=res;
}
for (int j=0;j<cur.size();j++)
for (int k=0;k<3;k++){
int a=cur[j].v[k],b=cur[j].v[(k+1)%3];
if (vis[a][b]&&!vis[b][a])
next.push_back(Face(a,b,i));
}
cur=next;
}
}
double Area(Point a,Point b,Point c){
return 0.5*Length(cross(a,b,c));
}
double Area(Face F,Point *P){
return Area(P[F.v[0]],P[F.v[1]],P[F.v[2]]);
}
double Area(Point *P){
double Answer=0;
for (int i=0;i<cur.size();i++)
Answer+=Area(cur[i],P);
return Answer;
}
}Hull;
int n;
int main(){
scanf("%d",&n);
for (int i=1;i<=n;i++){
scanf("%lf%lf%lf",&P[i].x,&P[i].y,&P[i].z);
p[i]=P[i];
}
Hull.Increment(p,n);
printf("%.6lf",Hull.Area(P));
return 0;
}
BZOJ1209 [HNOI2004]最佳包裹 三维凸包 计算几何的更多相关文章
- bzoj 1209: [HNOI2004]最佳包裹 三维凸包
1209: [HNOI2004]最佳包裹 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 160 Solved: 58[Submit][Status] ...
- bzoj 1964: hull 三维凸包 计算几何
1964: hull 三维凸包 Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 54 Solved: 39[Submit][Status][Discuss ...
- 洛谷P2287 [HNOI2004]最佳包裹(三维凸包)
题面 传送门 题解 左转板子,调个精度就能\(A\)了 //minamoto #include<bits/stdc++.h> #define R register #define fp(i ...
- BZOJ1209 最佳包裹 (三维凸包 增量法)
题意 求三维凸包的表面积. N≤100N\le100N≤100 题解 暴力往当前的凸包里加点.O(n2)O(n^2)O(n2).题解详见大佬博客 扰动函数shakeshakeshake是为了避免四点共 ...
- 洛谷P4502 [ZJOI2018]保镖(计算几何+三维凸包)
题面 传送门 题解 我对计蒜几盒一无所知 顺便\(xzy\)巨巨好强 前置芝士 三维凸包 啥?你不会三维凸包?快去把板子写了->这里 欧拉公式 \[V-E+F=2\] \(V:vertex\)顶 ...
- Hnoi2004 金属包裹
传送门 三维凸包模板题……只是听了听计算几何的课之后心血来潮想写的…… 我的做法很无脑是吧……暴力枚举三个点组成的三角形,然后枚举剩下的点,判断其余点是否都在这个三角形的同一侧,是的话则说明这个三角形 ...
- POJ 2225 / ZOJ 1438 / UVA 1438 Asteroids --三维凸包,求多面体重心
题意: 两个凸多面体,可以任意摆放,最多贴着,问他们重心的最短距离. 解法: 由于给出的是凸多面体,先构出两个三维凸包,再求其重心,求重心仿照求三角形重心的方式,然后再求两个多面体的重心到每个多面体的 ...
- hdu4273Rescue(三维凸包重心)
链接 模板题已不叫题.. 三维凸包+凸包重心+点到平面距离(体积/点积) 体积-->混合积(先点乘再叉乘) #include <iostream> #include<cstd ...
- hdu4449Building Design(三维凸包+平面旋转)
链接 看了几小时也没看懂代码表示的何意..无奈下来问问考研舍友. 还是考研舍友比较靠谱,分分钟解决了我的疑问. 可能三维的东西在纸面上真的不好表示,网上没有形象的题解,只有简单"明了&quo ...
随机推荐
- Tomcat中配置URIEncoding="UTF-8"来处理中文的方法
http://www.cnblogs.com/seabird1979/p/4837237.htmlTomcat中配置URIEncoding="UTF-8"来处理中文的处理打开 se ...
- u-boot移植(十)---代码修改---支持nor flash
一.问题定位 开发板重启后打印了2个提醒和一个错误,caches的提醒先不看,看看flash和nand下面的提醒,bad CRC,Using default enviroment,我们可以定位Usin ...
- CF912E Prime Gift
传送门 看到\(n\)只有16,可以把这些质数分成两半,然后预处理出这些数相乘得出的小于\(10^{18}\)的所有数,排个序,然后二分最终答案,再用两个指针从前往后和从后往前扫,进行\(two-po ...
- PXC中的GTIDs
基本环境:PXC 5.7.19 Row+Gtid,3节点 一.Galera GTID vs MySQL GTID 1.1.Galera GTID vs MySQL GTID Both kinds of ...
- android 面试事件分发
开始的Activity和最好的view是没有onInterceptTouchEvent方法的 不能进行对事件的拦截 总论 在我们点击屏幕时,会有下列事件发生: Activity调用dispathTou ...
- Java注解之Retention、Documented、Target、Inherited介绍
先看代码,后面一个个来解析: @Retention(RetentionPolicy.RUNTIME) @Target(value = {ElementType.METHOD, ElementType. ...
- MR架构
MapReduce框架结构 Map/Reduce是一个用于大规模数据处理的分布式计算模型,它最初是由Google工程师设计并实现的,Google已经将它完整的MapReduce论文公开发布了.其中对它 ...
- GitHub:Awesome-Hacking(黑客技能列表-恶意代码)
0 初衷 GitHub这一份黑客技能列表很不错,包含了多个方向的安全.但目前我关注只有逆向工程与恶意代码,所以其他的被暂时略过. 虽然很感谢作者的辛勤付出,但并不打算复制粘贴全套转载.逐条整理是为了从 ...
- 在全志平台调试博通的wifi驱动(类似ap6212)【转】
转自:http://blog.csdn.net/fenzhi1988/article/details/44809779 调试驱动之前,首先先看看驱动代码,了解代码大致工作流程,再根据硬件配置驱动,比如 ...
- Python3学习笔记21-实例属性和类属性
由于Python是动态语言,根据类创建的实例可以任意绑定属性. 给实例绑定属性的方法是通过实例变量,或者通过self变量: class Student(object): def __init__(se ...