计算面积 - 多边形相关（面积，点的位置）

三角形：海伦公式

$S=\frac{1}{2}(x1y2+x2y3+x3y1-x1y3-x2y1-x3y2);$

//三角形面积
double sTriangulation1(vector<double>x, vector<double>y) {
	//海伦公式
	double S =0.5* ((x[0]*y[1] + x[1]*y[2] + x[2]*y[0]) - (x[1]*y[0] + x[2]*y[1] + x[0]*y[2]));
	if (S < 0)
		return -S;
	return S;
}

double sTriangulation2(vector<double>x, vector<double>y) {
     //向量  与原点
	int n = x.size();
	double sum = 0;
	for (int i = 0; i < n; i++) {
		sum += 0.5 * (x[i]*y[(i+1)%n]-x[(i+1)%n]*y[i]);
	}
	if (sum < 0)
		sum = -sum;
	return sum;
}

多边形面积

考虑有凹边形，统一使用 向量 计算面积的方式，有正负
思路 ： 向量法 （点数组需要时顺/逆时针排序的）

• $S_{ABC}=\frac{1}{2}\overrightarrow{AB}\ast \overrightarrow{BC}$
• 以 原点(0,0) 为 第一个点，那么三角学的面积可以写成
  $S_{ABC}=S_{OBC}-S_{OAB}-S_{OAC}$
  这里默认 BC是最远的边，但一般不知是哪个边，无法提前确认+ -，故都使用 向量去做
• $S_{ABC}=S_{OAB}+S_{OBC}+S_{OCA}$
  (上一个结束点 是下一个开始点 循环到第一个开始点)
• 同理 五边形ABCDE有：
  $S_{ABCDE}=S_{OAB}+S_{OBC}+S_{OCD}+S_{ODE}+S_{OEA}$
  对于一个$S_{OAB}$的计算公式； 其中$A(x_1,y_1),B(x_2,y_2)$：
  $S=\frac{1}{2}\overrightarrow{0A}\ast \overrightarrow{0B}=\frac{1}{2}（x_1\ast y_2-x_2\ast y_1）$

//x={1,0,-1,0},y={0,1,0,-1}; S=2.0
//x={0,1,0},y={0,0,1}; S=0.5


double SumTriangulation(vector<double>x, vector<double>y) {

	//看成任意两边间与原点的面积（有正负）//N个 最后一个是OEA
	int n = x.size();//n个点   有n-1次
	vector<double> sum;
	double s = 0.0;
	for (int i = 0; i < n-1; i++) {
		s= 0.5 * (x[i] * y[i + 1] - x[i + 1]*y[i]);
		sum.push_back(s);
	}
	//最后一个三角形 OEA
	s= 0.5 * (x[n-1] * y[0] - x[0] * y[n-1]);
	sum.push_back(s);
	s = 0;
	//输出看看
	for (int i = 0; i < sum.size(); i++) {
		cout << sum[i] << " ";
		s += sum[i];
	}
	cout <<endl;
	return s;
}

P点是否在多边形内

思路：射线法 （默认时顺/逆时针排序的）

• 通过p点引出一条 水平右射线
• 判断射线与 多边形的边的交点和数量，奇数就是在内
• 特殊边界：
  • 交点刚好是多边形的点，不算
  • 交点在左边，不算
• 加速判断 P(x,y) 与 边 $(x_1,y_1),(x_2,y_2)$
  • $y>=max(y_1,y_2)$， 没有
  • $y<=min(y_1,y_2)$，没有 (这两包含了 交点刚好是边的情况)
  • $x>=max(x_1,x_2)$，没有
  • $y_1==y_2$ ，边水平，没有
• 利用两点式，求出 交点 X
  两点式：$\frac{y-y_1}{X-x_1}=\frac{y_2-y_1}{x_2-x_1}=$
  转化为：$X=\frac{(y-y_1)\ast (x_2-x_1)}{(y-y_1)}+x_1$

//点在多边面内
//map<pair<int,int>>(x,y)
bool PtInPolygon(vector<double>x, vector<double>y, double px, double py) {

	int n = x.size();
	int count = 0;
	for (int i = 0; i < n; i++) {
		//判断 是否跳过
		if(y[i]==y[(i + 1) % n]) continue;
		if(py>= max(y[i], y[(i + 1) % n])) continue;
		if(py<= min(y[i], y[(i + 1) % n])) continue;
		if(px>= max(x[i], x[(i + 1) % n])) continue;
		
		//两点式
		double X = ((y[(i + 1) % n] - y[i]) * (x[(i + 1) % n])/(py-y[i]) - x[i]) + x[i];

		if (X > px) //还是有出现在左边的情况 去掉
			count++;
	}
	if (count % 2 == 1)
		return true;
	return false;
}