计算几何基础-2

文章目录

  • 直线:
  • 图形:
  • 求垂足
  • 求两圆交点
  • 直线与圆交点
  • 多边形问题
    • 判断一个点是否在任意多边形内部
    • Pick定理
  • 凸包
    • 求点集的凸包
      • 水平法:
      • 增量法:
  • 半平面
    • 半平面交
    • 求半平面交

直线:

struct Line{point p,v;Line(){}Line(point _p.point _v):p(_p),v(_v){}
}L[N];

图形:

在这里插入图片描述

求垂足

在这里插入图片描述

point vp(Line l,point p)//求p点向l做垂线得到的垂足 
{double Dis=dis(l,p);//求p到l的距离Dis Dis=Dis*Dis;double pDis=(p-l.p).length();//求AP距离pDis pDis=pDis*pDis;double vDis=sqrt(pDis-Dis);//勾股定理得到AH距离 point ans=l.p+l.v*(vDis/l.v.length());//通过方向单位向量乘长度得到路径,加起点得到答案 return ans;
}

求两圆交点

在这里插入图片描述

  1. 先求余弦定理 cosp1
  2. 反函数 角p1
  3. p1p2旋转到p11T方向
  4. 调整长度
  5. p1+p1T向量=T
    代码:
point get_intersection(circle A,circle B)
{point AB=B.o-A.o;double costheta=AB.length() * AB.length() + A.r *A.r-B.r*B.r;costheta=costheta/(2*AB.length()*A.r);//余弦定理AB=AB.rotatr(acos(costheta));//旋转theta角return A.o+AB*(A.r/AB.length()); 
} 

直线与圆交点

在这里插入图片描述
代码:

point get_intersection(circle A,Line B)
{point H=vp(B,A.o);//求垂线 double OH=(H-A.o).length();double HY=sqrt(A.r * A.r -OH * OH);return H + B.v *(HY/B.v.length());// return H + B.v *(-HY/B.v.length());
}

点+单位向量 * 长度

多边形问题

平面上n个点收尾顺次连接组成的平面图形
可能是凸多边形或者凹多边形

判断一个点是否在任意多边形内部

从这个点出发引一条射线,如果这个射线与多边形有奇数个交点则在内部,否则在外部在这里插入图片描述

double calc_S()
{double ans=0;int n=poly.size();for(int i=1;i<n;i++)ans+=poly[i-1]^poly[i];ans+=poly[n-1]^poly[0];return ans/2;
}

Pick定理

对于顶点都是整点的多边形,设其面积为S,多边形内部的点数为a,边上的点数为b,那么满足:
S = a + b/2 - 1

凸包

给出一个二维平面内的点集,如果任意两个点的连线都在点集内,则这个点集是个凸集
对于给定的散点集X,包含X的所有凸集的交集S叫做X的凸包
点就是墙上的钉子,用一个橡皮筋套在外面,收缩后形成的凸多边形就是凸包

求点集的凸包

水平法:

凸包一般使用水平法求解
将凸包分为上凸壳和下凸壳两部分,分别求解
每一部分按照x的坐标排序,用单调栈维护,利用叉积的符号判断凹凸性
在这里插入图片描述

增量法:

在这里插入图片描述
判断情况:
此时B合法
在这里插入图片描述
此时B不合法
在这里插入图片描述
不合法的点一定是最后添加入S的,且是连续的
在这里插入图片描述

//判断上凸包
sort(poly.begin(),poly.end(),cmp);
int top=0,n=poly.size();
Stack[++top]=poly[0];
for(int i=1;i<n;i++)
{point C=poly[i];while(top>1){point B=Stack[top];point A=Stack[top-1];point AB=B-A;point AC=C-A;if((AB^AC)>eps)top--;else break;}Stack[++top]=C;
}

半平面

• 顾名思义,就是平面的一半。一条直线会把平面分成两部分,就 是两个半平面。对于半平面,我们可以用直线方程式如: ax+by+c > 0 表示,更常用的是用直线表示

半平面交

• 顾名思义,就是多个半平面求交集。其结果可能是一个凸多边形、 无穷平面、直线、线段、点等。 • 什么时候需要半平面交?
• 1. 二维线性规划 (高中数学)
• 2. 求多边形的核 • 多边形的核:如果多边形中存在一个区域使得在区域中可以看到 多边形中任意位置(反之亦然),则这个区域就是多边形的核

求半平面交

• 理论上有很多很多种求法,但是在实际应用中效率最高而且最好 写的是 S&I 增量法,也就是一个一个插入半平面并且更新答案。
• 核心思想:
• 1. 选取逆时针方向为正方向,把所有的直线变成向量。
• 2. 所有的向量按照极角排序,角度相同的保留左边的。
• 3. 按照顺序每次插入一个平面,删掉右面的部分,保留左边的部 分。
###具体方法:
1.选方向,排序(如上文)
• 2.用双端队列保存构成当前核的所有向量。
• 3.按照顺序遍历所有向量,每次加入判断影响
• 假设队列中最后两条直线的交点是左图所示。如果新加入了直线 后原来的交点在直线右侧,说明最后一条直线没有用,把它从队 列中删除。

在这里插入图片描述
看遍边6要不要保留,就看5和6的交点在直线的哪一侧
在这里插入图片描述

由p指向v的向量
ans存的交点
q存的边
在这里插入图片描述

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318480.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

eShopOnContainers 看微服务 ②:配置 启动

一、什么是dockerDocker 是一个开源项目&#xff0c;通过把应用程序打包为可移植的、自给自足的容器&#xff08;可以运行在云端或本地&#xff09;的方式&#xff0c;实现应用程序的自动化部署。使用 Docker 的时候&#xff0c;需要创建一个应用或服务&#xff0c;然后把它和它…

判断整除(opj)(动态规划)

解析 与取模结合的动归&#xff0c;正常做即可 问题 眼瞎&#xff01;&#xff01;&#xff01; 这个序列的每个数都必须用到&#xff01;&#xff01;&#xff01; if(f[i-1][j]) f[i][j]1;上面这行就是不对的&#xff01;&#xff01;&#xff01; 头疼 仔细审题 opj的题…

[USACO09FEB]Revamping Trails G

题意&#xff1a; 约翰一共有 N 个牧场.由 MM 条布满尘埃的小径连接。小径可以双向通行。每天早上约翰从牧场 1 出发到牧场 N 去给奶牛检查身体。 通过每条小径都需要消耗一定的时间。约翰打算升级其中 K 条小径&#xff0c;使之成为高速公路。在高速公路上的通行几乎是瞬间完…

NET Core微服务之路:自己动手实现Rpc服务框架,基于DotEasy.Rpc服务框架的介绍和集成...

本篇内容属于非实用性&#xff08;拿来即用&#xff09;介绍&#xff0c;如对框架设计没兴趣的朋友&#xff0c;请略过。 快一个月没有写博文了&#xff0c;最近忙着两件事;一&#xff1a;阅读刘墉先生的《说话的魅力》&#xff0c;以一种微妙的&#xff0c;你我大家都会经常遇…

花店橱窗布置(洛谷P1854)(动态规划)

传送门 文章目录解析问题代码解析 一道很正常的动态规划 dp[i][j]表示到第j个花瓶放了第j朵花的dp最优值 注意&#xff1a;是严格使第i朵放在j瓶 找到最优解递归输出即可 问题 又是初始化的问题&#xff01;&#xff01;&#xff01; 一开始把dp赋值成负无穷时落掉了j0的一行…

P4009 汽车加油行驶问题

题目描述&#xff1a; 题解&#xff1a; 看了很多题解&#xff0c;无论什么解法都绕不开分层图 在本题中加满油的车每次可以移动K步&#xff0c;那么我们就可以建立一个K1层的分层图&#xff0c;表示汽车油量k的状态&#xff08;油量0…k&#xff09;&#xff0c;然后根据题目…

.net core i上 K8S(五).netcore程序的hostip模式

正文上一章讲了pod的管理&#xff0c;今天再分享一个pod的访问方式1.Pod的HostIP模式Pod的HostIP模式&#xff0c;可以通过宿主机访问pod内的服务&#xff0c;创建yaml文件如下apiVersion: v1 kind: Pod metadata: name: netcore-podlabels:app: netcorepod spec:containers: …

状态压缩:枚举子集(最优组队)(ybtoj)(动态规划)

解析 很裸的状压dp 但是直接暴力的话状态2n,枚举2n 乘在一起会T诶 怎么办呢&#xff1f; 使用下面这个循环&#xff0c;就可以保证只会枚举当前状态s的子集 for(int i(s-1)&s;i;i(i-1)&s){........ }证明 举举例子就挺明显了 为什么不重不漏呢&#xff1f; 首先i肯…

【活动(深圳)】告别2018之12.22 大湾区.NET Meet 大会 ,同时有网络直播

今年的 Connect(); 主题更加聚焦开发者工具生产力、开源&#xff0c;以及无服务器&#xff08;Serverless&#xff09;云服务。Visual Studio 2019 AI 智能加持的 IntelliCode、实时代码协作共享 Live Share&#xff0c;.NET Core 3.0的预览版本附带了大量更新&#xff0c;旨在…

最短路径(状压dp)(ybtoj)

解析 “bug总有de完的一天” 头疼 暴力写的话本题显然复杂度是假的 有一个很好的思路优化时间复杂度 先用dp[k][i]表示**从第k个关键点到任意i点的最短路 跑k遍 SPFA或迪杰斯特拉 即可 然后转移时可以只考虑关键点 使状态转移数大大降低 细节 头真疼 边界条件&#xff1a; …

P2403 [SDOI2010]所驼门王的宝藏(强连通分量)(拓扑排序)

文章目录题目描述解析代码洛谷传送门题目描述 解析 看题目要求很容易想到强连通分量缩点加拓扑dp 但是问题在于存图 第一感就是和暴力和每个点连边 但那样无论点数和边数都很爆炸 随后我们发现这个图非常稀疏 所以我们可以只连有宝藏的点 然而这样边数会被一行横门这样的数据…

Xamarin.Forms之UserDialogs 重制版本

在 forms 里面&#xff0c;目前使用比较多的弹出组件是 Acr.UserDialogs &#xff0c;但是这个组件有些小问题&#xff0c;比如 loading .hide 会同时把 toast 给一起关掉&#xff0c;android 下的 toast 希望是 安卓原生的toast 样子&#xff0c;而不是 底部弹出一个横条&…

H.Minimum-cost Flow

H.Minimum-cost Flow 题目&#xff1a; 其实就是给出每条边的单位费用&#xff0c;q次查询&#xff0c;每次查询改变所有边的容量&#xff08;所有边容量一样&#xff09;&#xff0c;问最后流出1流量的最小花费是多少&#xff1f; 题解&#xff1a; 暴力做法肯定是每次询问…

eShopOnContainers 看微服务③:Identity Service

引言通常&#xff0c;服务所公开的资源和 API 必须仅限受信任的特定用户和客户端访问。那进行 API 级别信任决策的第一步就是身份认证——确定用户身份是否可靠。在微服务场景中&#xff0c;身份认证通常统一处理。一般有两种实现形式&#xff1a;基于API 网关中心化认证&#…

P2168 [NOI2015] 荷马史诗(哈夫曼编码树)

传送门 文章目录题目描述解析细节代码题目描述 解析 其实就是构造一棵树 另所有带权点都为叶子节点 其代价为权值与深度的乘积 求最小代价及最小代价下的最小深度 可以看成一开始有n棵小树 每次把k棵合并在一起 最后合成一棵大树就行了 每次合并的代价是k棵树的权值和 看到这…

B-Suffix Array

B-Suffix Array 题意&#xff1a; 一个字符串只含有a和b&#xff0c;先给出b数组的构造方式&#xff1a; 对于每个位置i来说&#xff1a; 如果存在一个位置j&#xff0c;使得j<i,且s[j] s[i],则b[i]i-j否则b[i]0 现在对字符串每个后缀都构造B数组&#xff0c;并按照字典…

线段树合并、分裂

基本概念&#xff1a; 如果需要维护许多个大小为 \(10^5\) 级别的多重集&#xff0c;可以看做给每一个多重集建立一棵线段树。线段树的合并、分裂就是多重集的累加、分开。 这里使用动态开点的方式存储线段树树。 如果一个节点为空&#xff0c;那么它的编号为 \(0\) 。 变量释义…

(神奇的)虚树(初步了解)

参考文献&#xff1a; 博文1 博文2 博文3 引入 在一类树上动态规划问题中,题目给出的询问往往包含树上的很多各节点,并保证总的点数规模小于某个值. 如果我们直接在整颗树上进行dp的话,时间复杂度与询问的次数有关,这显然是不可接受的,如果我们可以找到一种动态规划的方法,使…

P2048 [NOI2010] 超级钢琴(RMQ 贪心)

文章目录题目描述解析代码传送门题目描述 解析 首先&#xff0c;如果只有一个和弦&#xff0c;那么问题显然简单了 用前缀和结合ST表随便做做即可 然而 这次要求前k大的 怎么办呢&#xff1f; 参照之前有一道序列合并的做法 我们想到&#xff0c;可以先建一个优先队列&#xf…

微服务架构基础之Service Mesh

ServiceMesh(服务网格) 概念在社区里头非常火&#xff0c;有人提出 2018 年是 ServiceMesh 年&#xff0c;还有人提出 ServiceMesh 是下一代的微服务架构基础。那么到底什么是 ServiceMesh&#xff1f;它的诞生是为了解决什么问题&#xff1f;企业是否适合引入 ServiceMesh&…