7.26 简单计算几何

2024-04-10 18:38
文章标签 简单 计算 几何 7.26

本文主要是介绍7.26 简单计算几何,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

向量加减法,长度,数乘


代码实现

Pt operator - (Pt a,Pt b){
return Pt(a.x-b.x,a.y-b.y);}
Pt operator + (Pt a, Pt b){
return Pt(a.x + b.x, a.y + b.y); }
double Length ( Point q){
return sqrt (q.x*q.x+q.y*q.y); }
Pt operator * ( double A, Pt p){
return Pt(p.x*A, p.y*A); }

 

点积

double dot (Pt a, Pt b) {
return a.x * b.x + a.y * b.y; }
double Angle ( Vector A, Vector B) {
return acos ( dot(A,B)/ Length (A)/ Length (B ));}

叉积

double det (Pt a, Pt b) {
return a.x * b.y - a.y * b.x; }

 

点在直线上的判断,直线间的关系

点在线段上的判断

bool PtOnSegment (Pt s, Pt t, Pt a) {
return ! det (a-s,a-t)&& min (s.x,t.x)<=a.x &&
a.x <= max(s.x,t.x)&& min(s.y, t.y) <= a.y &&
a.y <= max(s.y, t.y);
}
bool PtOnSegment (Pt p, Pt a, Pt b) {
return ! sgn (det (p-a,b-a ))&& sgn (dot (p-a,p-b)) <=0;
}

线段间的关系

这里写图片描述

点到线段的距离

求点P到线段AB的最短距离
1.判断线段PA和AB所成的夹角,如果是钝角,那么jPAj是
点到线段的最短距离。
2.判断线段PB和AB所成的夹角,如果是钝角,那么jPBj是
点到线段的最短距离。
3.线段PA和线段PB与AB所成的夹角都不为钝角,那么
点P到线段AB的距离是点P到直线AB的距离,这个距离可
以用面积法直接算出来。

求三角形的面积

判断点是否在三角形内

重心

Pt triangleMassCenter (Pt a, Pt b, Pt c) {
return (a+b+c) / 3.0;}

求多边形的面积

//代码实现多边形的有向面积,这里返回正值
double PolygonArea (Pt* p,int n){
double area =0;
for (int i=1;i<n -1; i++)
area += det (p[i]-p[0] ,p[i+1] -p [0]);
return fabs ( area /2);
}

多边形内格点数

Pick公式

//代码实现:求多边形边界上的点数:可以当成伪代码来看
int polygon_border_point_cnt ( const Polygon &p) {
int ans = 0;
int n = p. size ();
for (int i = 0; i < n; ++i)
ans += gcd (Abs (int (p[ next (i)].x-p[i].x)),
Abs(int(p[ next (i)].y-p[i].y )));
return ans ;
}

凸包 

定义:凸包就是把给定点包围在内部的、面积最小的凸多边形

算法:Andrew算法

首先把所有点按照从小到大排序(如果x相同,按照y从小到大排
序),删除重复点后得到序列p1; p2; :::,然后把p1和p2放到凸包
中。从p3开始,当新点在凸包“前进”方向的左边时继续,否则
依次删除最近加入凸包的点,直到新点在左边。
注意:从左到右和从右到左各扫描一次

 

 

 

这篇关于7.26 简单计算几何的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/891836

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 10130 简单背包

题意: 背包和 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <queue>#include <map>

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

poj 3304 几何

题目大意:给出n条线段两个端点的坐标,问所有线段投影到一条直线上,如果这些所有投影至少相交于一点就输出Yes!,否则输出No!。 解题思路:如果存在这样的直线,过投影相交点(或投影相交区域中的点)作直线的垂线,该垂线(也是直线)必定与每条线段相交,问题转化为问是否存在一条直线和所有线段相交。 若存在一条直线与所有线段相交,此时该直线必定经过这些线段的某两个端点,所以枚举任意两个端点即可。