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 ...
随机推荐
- modbus-vcr介绍
相关链接:modbus-vcr modbus-vcr是一个Ettercap的插件,被使用在缺少数据完整性的工业控制系统协议方面. 这个Ettercap插件执行一个MITM攻击在使用Modbus协议的系 ...
- C# 简单线程实例
1.简单线程实例 以及委托(同步委托.异步委托) using System; using System.Collections.Generic; using System.Linq; using Sy ...
- 利用 python requests完成接口文件上传
最近在准备一个公开课,主题就是利用不同的语言和不同的工具去实现文件的上传和下载. 在利用Jmeter去实现功能的时候,以及利用loadrunner去写脚本的时候,都很顺利,没有任何问题,当我尝试用Py ...
- js 执行上下文理解
前端基础进阶(三):变量对象详解http://www.jianshu.com/p/330b1505e41d 1.创建阶段 a.生成变量对象 1.创建arguments对象 2.functio ...
- Shiro+Spring+SpringMVC+Mybatis整合
Demo源码地址,v1.0分支 https://github.com/jxjy/hr
- VMware12虚拟机中Ubuntu16.04安装CPU版本Caffe
首先,可以自行下载VMware12进行安装,基本上都是直接点击‘下一步’直到安装完成,这里重点讲一下Ubuntu16及Caffe的安装步骤 第一步: 下载Ubuntu16.04版本的文件,这里给出链接 ...
- springboot项目发布到独立的tomcat中运行&打成jar包运行
springboot的打包方式依赖于插件:(下面插件打出的包与普通的包目录结构有区别) <plugin> <groupId>org.springframework.boot&l ...
- Maven安装配置操作
1)下载maven安装包并解压: 2)环境变量配置: 3)编辑环境变量Path,追加%MAVEN_HOME%\bin; 4)maven安装配置后进行dos命令检查:在cmd中输入 mvn -v 5)配 ...
- C++学习8-面向对象编程基础(模板)
模板 模板是一种工具,模板可以使程序员能建立具有通用类型的函数库与类库: 模板具有两种不同的形式: 函数模板 类模板 函数模板 当一个add()函数接收两个参数,因为某种特定情况,所传入的实参数据类型 ...
- 【转】Python之日期与时间处理模块(date和datetime)
[转]Python之日期与时间处理模块(date和datetime) 本节内容 前言 相关术语的解释 时间的表现形式 time模块 datetime模块 时间格式码 总结 前言 在开发工作中,我们经常 ...