程序分析:
给定平面上一点p(x0,y0),判断该点是否在三角形ABC中,三角形顶点坐标分别为A(xa,xb)、B(xb,yb)、C(xc,yc)。可以使用面积法来判断,方法如下:其中S(A,B,C)表示三角形ABC的面积。
(1)若abs( S(A,B,C) ) = abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ,则P在三角形ABC的内部或边上;如果还有abs( S(P,B,C) )、abs( S(A,P,C) ) 和abs( S(A,B,P) )全都大于0,则说明P在三角形ABC的内部,否则P在三角形ABC的边上,具体为:S(P,B,C)为0,则说明P在BC边上,S(A,P,C)为0,则说明P在AC边上,S(A,B,P)为0,则说明P在AB边上;
(2)若abs( S(A,B,C) ) < abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) ,则P在三角形ABC的外部;
(3)对abs( S(A,B,C) ) > abs( S(P,B,C) ) + abs( S(A,P,C) ) + abs( S(A,B,P) ) 情况在理论上是不存在的。
对于如何求平面中三角形的面积?可以使用叉乘法来实现,即S(A,B,C)为向量AB叉乘AC所得向量模的1/2,再对该值求绝对值就是三角形ABC的面积。
C++代码:
#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cmath>
using namespace std;
#define EPS 1e-6
class Point {
public:
void setXY(int x, int y) {
this->x = x;
this->y = y;
}
int getX() {
return this->x;
}
int getY() {
return this->y;
}
private:
int x;
int y;
};
double Area(Point a, Point b, Point c) {
int x1 = b.getX() - a.getX();
int y1 = b.getY() - a.getY();
int x2 = c.getX() - a.getX();
int y2 = c.getY() - a.getY();
return abs((x1*y2 - x2*y1) / 2.0);
}
void IsInsideTriangle(Point a, Point b, Point c, Point p) {
double area = Area(a, b, c);
double area1 = Area(a, b, p);
double area2 = Area(a, c, p);
double area3 = Area(b, c, p);
double sumarea = area1 + area2 + area3;
if (abs(area - sumarea) < EPS) {
if (abs(area1) < EPS || abs(area2) < EPS || abs(area3) < EPS) {
cout << "The point is on the triangle." << endl;
}
else {
cout << "The point is inside the triangle." << endl;
}
}
else {
cout << "The point is outside the triangle." << endl;
}
}
int main() {
Point a, b, c, p;
int x, y;
cout << "The first point of triangle:";
cin >> x >> y;
a.setXY(x, y);
cout << "The second point of triangle:";
cin >> x >> y;
b.setXY(x, y);
cout << "The third point of triangle:";
cin >> x >> y;
c.setXY(x, y);
cout << "The point:";
cin >> x >> y;
p.setXY(x, y);
IsInsideTriangle(a, b, c, p);
return 0;
}
C++运行结果: