问题描述

对于三个点p0、p1、p2，请按照下列情况进行输出：
p0、p1、p2成逆时针方向 COUNTER_CLOCKWISE
p0、p1、p2成顺时针方向 CLOCKWISE
p2、p0、p1依次排列在同一直线上 ONLINE_BACK
p0、p1、p2依次排列在同一直线上 ONLINE_FRONT
p2在线段p0p1上 ON_SEGMENT

输入：
x p 0 x_{p0} xp0 y p 0 y_{p0} yp0 x p 1 x_{p1} xp1 y p 1 y_{p1} yp1
q q q
x p 2 0 x_{p2_0} xp20 y p 2 0 y_{p2_0} yp20
x p 2 1 x_{p2_1} xp21 y p 2 1 y_{p2_1} yp21
…
x p 2 q − 1 x_{p2_{q-1}} xp2q−1 y p 2 q − 1 y_{p2_{q-1}} yp2q−1
第1行输入p0、p1的坐标。接下来给出q个p2的坐标用作问题。
输出：
根据各问题输出上述状态之一，每个问题占1行。
限制：
1 ≤ q ≤ 1000
-10000 ≤ x i , y i x_i,y_i xi,yi ≤ 10000
p0、p1不是同一个点。

输入示例

输出示例

COUNTER_CLOCKWISE
CLOCKWISE
ONLINE_BACK
ON_SEGMENT
ONLINE_FRONT

讲解

向量a、b外积的方向与a、b所在平面垂直且满足右手螺旋定则，因此我们可以根据外积的方向判断向量a、b的位置关系。

(1)p2位于p0→p1的逆时针方向
(2)p2位于p0→p1的顺时针方向
(3)p2位于直线p0p1上，且顺序为p2→p0→p1
(4)p2位于直线p0p1上，且顺序为p0→p1→p2
(5)p2位于线段p0p1上，顺序为p0→p2→p1

设p0到p1的向量为a，p0到p2的向量为b，则各种情况的判断方法如下：

(1)外积的大小cross(a, b)为正时，可确定b在a的逆时针位置
(2)外积的大小cross(a, b)为负时，可确定b在a的顺时针位置
(3)(1)(2)都不符合时，表示p2位于直线p0p1上（不一定在线段p0p1上）。cosθ在θ大于90度或小于-90度时为负，因此a与b的内积dot(a, b)为负时，可确定p2位于线段p0p1后方，即顺序为p2→p0→p1。
(4)不符合(3)时，p2的位置为p0→p1→p2或p0→p2→p1。因此。如果b的大小大于a的，则可确定p2位于线段p0p1上。
(5)不符合(4)时，可确定p2位于线段p0p1上

下面的程序将三个点p0、p1、p2作为参数，返回点p2与向量p0→p1的位置关系。

逆时针方向CCW：

static const int COUNTER_CLOCKWISE = 1;
static const int CLOCKWISE = -1;
static const int ONLINE_BACK = 2;
static const int ONLINE_FRONT = -2;
static const int ON_SEGMENT = 0;int ccw(Point p0, Point p1, Point p2) {Vector a = p1 - p0;Vector b = p2 - p0;if( cross(a, b) > EPS ) return COUNTER_CLOCKWISE;if( cross(a, b) < -EPS ) return CLOCKWISE;if( dot(a, b) < -EPS ) return ONLINE_BACK;if( a.norm() < b.norm() ) return ONLINE_FRONT;return ON_SEGMENT;
}

AC代码如下

#include<stdio.h>
#include<iostream>
#include<cmath>
using namespace std;#define EPS (1e-10)
#define equals(a, b) (fabs((a) - (b)) < EPS)class Point {//Point类，点 public:double x, y;Point(double x = 0, double y = 0): x(x), y(y) {}Point operator + (Point p) { return Point(x + p.x, y + p.y); }Point operator - (Point p) { return Point(x - p.x, y - p.y); }Point operator * (double a) { return Point(a * x, a * y); }Point operator / (double a) { return Point(x / a, y / a); }double abs() { return sqrt(norm()); }double norm() { return x * x + y * y; }bool operator < (const Point &p) const {return x != p.x ? x < p.x : y < p.y;}bool operator == (const Point &p) const {return fabs(x - p.x) < EPS && fabs(y - p.y) < EPS;}
};typedef Point Vector;//Vector类，向量 struct Segment{//Segment 线段 Point p1, p2;
};
double dot(Vector a, Vector b) {//内积 return a.x * b.x + a.y * b.y;
}double cross(Vector a, Vector b) {//外积 return a.x*b.y - a.y*b.x;
}static const int COUNTER_CLOCKWISE = 1;
static const int CLOCKWISE = -1;
static const int ONLINE_BACK = 2;
static const int ONLINE_FRONT = -2;
static const int ON_SEGMENT = 0;int ccw(Point p0, Point p1, Point p2) {Vector a = p1 - p0;Vector b = p2 - p0;if( cross(a, b) > EPS ) return COUNTER_CLOCKWISE;if( cross(a, b) < -EPS ) return CLOCKWISE;if( dot(a, b) < -EPS ) return ONLINE_BACK;if( a.norm() < b.norm() ) return ONLINE_FRONT;return ON_SEGMENT;
}int main(){Point p0, p1, p2;cin>>p0.x>>p0.y>>p1.x>>p1.y;int q;cin>>q;while(q--){cin>>p2.x>>p2.y;switch(ccw(p0, p1, p2)){case 1: cout<<"COUNTER_CLOCKWISE"<<endl; break;case -1: cout<<"CLOCKWISE"<<endl; break;case 2: cout<<"ONLINE_BACK"<<endl; break;case -2: cout<<"ONLINE_FRONT"<<endl; break;case 0: cout<<"SEGMENT"<<endl; break;}}
}

注：以上本文未涉及代码的详细解释参见：计算几何学

计算几何学 | 逆时针方向 | Counter-Clockwise | C/C++实现相关推荐

计算几何学 | 知识点及C++代码实现汇总
博客上看到的非常优质的计算几何学c++代码实现. 代码可直接访问github仓库计算几何学计算几何学 | 几何对象的基本元素与表现 | Basic 计算几何基础知识整理大全 | Vector | ...
算法导论之计算几何学
计算几何学是计算机科学的一个分支,专门研究集合问题的解决的算法.计算几何学的问题一般输入关于一组集合对象的描述,如一组点.一组线段:输出是对问题的回答,如直线是否相交.三维空间和高维空间很难视觉化,这 ...
计算几何学(Computational Geometry)
计算几何学(Computational Geometry) 分类:算法学习 | 作者:酷~行天下 | 发表于2011/11/27 2条评论 3,299 views 看名字就知道是研究几何问题的算法. ...
计算几何学 | 圆与直线的交点 | Cross Points of a Circle and a Line | C/C++实现
问题描述求圆c与直线 l l l的交点. 输入: 输入按照下述格式给出: c x cx cx c y cy cy r r r q q q L i n e 1 Line_1 Line1 L i n ...
Counter-Clockwise 三点位置关系计算几何学
题目链接:https://vjudge.net/problem/Aizu-CGL_1_C #include<bits/stdc++.h> using namespace std; #def ...
机器学习 - [源码实现决策树小专题]决策树中混杂度数值度量的Python编程实现（信息熵和基尼系数的计算）
混杂度数值度量的Python编程实现李俊才的 CSDN 博客:https://blog.csdn.net/qq_28550263?spm=1001.2101.3001.5343 邮箱 :29114 ...
华院计算｜维托里斯—最长寿的数学家
百岁长寿的数学家并不多见,大家熟识的两位中国人是108岁的浙江大学朱良璧教授(1913年11月–2021年11月20日)和101岁的复旦大学苏步青教授(1902年9月23日–2003年3月17日).朱 ...
量子计算机可以计算葛立恒数,葛立恒数
葛立恒数曾经被视为在正式数学证明中出现过最大的数,后来则被TREE(3)取代. 中文名葛立恒数外文名 Graham's Number性质拉姆齐理论问题发现人葛立恒现知下限 12 葛 ...
数学术语的英汉对照(权威,全面)
abbreviation 简写符号:简写 abscissa 横坐标 absolute complement 绝对补集 absolute error 绝对误差 absolute inequality 绝 ...

计算几何学 | 逆时针方向 | Counter-Clockwise | C/C++实现

问题描述

输入示例

输出示例

讲解

AC代码如下

计算几何学 | 逆时针方向 | Counter-Clockwise | C/C++实现相关推荐

最新文章

热门文章