标签: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. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。