ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

计算几何 点积叉积 点类模板

2021-09-20 19:32:22  阅读:210  来源: 互联网

标签:AC AB const point double return 点类 点积叉积 模板


模板摘自:ACM国际大学生程序设计竞赛 算法与实现 计算几何篇

计算几何

3.1多边形

3.1.1几何误差

const double eps = 1e-8;
int cmp (double x) {
    if (fabs(x) < eps) return 0;
    if (x > 0) return 1;
    return -1;
}

3.1.2计算几何点类

函数用法
double sqr(double x)计算平方
double det(const point &a, const point &b)计算两个向量的叉积
double dot(const point &a, const point &b)计算两个向量的点积
double dist(const point &a, const point &b)计算两个点的距离
point rotate_point(const point &p, double A)向量op绕原点逆时针旋转A(弧度)

点积

几何意义:|AC|*|AD|
A B → ∗ A C → = ∣ A B ∣ ⋅ ∣ A C ∣ ⋅ c o s θ \overrightarrow{AB}*\overrightarrow{AC} = |AB|\cdot|AC|\cdot cos \theta AB ∗AC =∣AB∣⋅∣AC∣⋅cosθ
在这里插入图片描述

如果点积是正的,那就代表在1,4象限

负的就代表在2,3象限

叉积

几何意义:|AC|*|BD|,即三角形abc面积的两倍
A B → ∗ A C → = ∣ A B ∣ ⋅ ∣ A C ∣ ⋅ s i n θ \overrightarrow{AB}*\overrightarrow{AC} = |AB|\cdot|AC|\cdot sin \theta AB ∗AC =∣AB∣⋅∣AC∣⋅sinθ
如果是正的,就代表在12象限

负的就是34象限

const double pi = acos(-1.0);
inline double sqr (double x) {
    return x*x;
}
struct point {
    double x, y;
    point() {}
    point(double x, double y) : x(a), y(b) {}
    void input() {
        scanf("%lf%lf", &x, &y);
//        cin>>x>>y;
    }
    friend point operator + (const point &a, const point &b) {
		return point(a.x + b.x, a.y + b.y);
	}
	friend point operator - (const point &a, const point &b) {
		return point(a.x - b.x, a.y - b.y);
	}
	friend bool operator == (const point &a, const point &b) {
		return cmp(a.x - b.x) == 0 && cmp(a.y - b.y) == 0;
	}
	friend point operator * (const point &a, const double &b) {
		return point(a.x * b, a.y * b);
	}
	friend point operator * (const double &a, const point &b) {
		return point(a * b.x, a * b.y);
	}
	friend point operator / (const point &a, const double &b) {
		return point(a.x / b, a.y / b);
	}
	double norm() {//计算长度
		return sqrt(sqr(x) + sqr(y));
	}
}
//计算两个向量的叉积
double det (const point &a, const point &b) {
    return a.x*b.y - a.y*b.x;
}
//计算两个向量的点积
double dot (const point &a, const point &b) {
    return a.x*b.x + a.y*b.y;
}
double dist (const point &a, const point &b) {
    return (a - b).norm();
}
double rotate_point (const point &p, double A) {
    double tx = p.x, ty = p.y;
    return point (tx * cos(A) - ty * sin(A), tx * sin(A) + ty * cos(A));
}

标签:AC,AB,const,point,double,return,点类,点积叉积,模板
来源: https://blog.csdn.net/Dueser/article/details/120394410

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有