计算两个多边形的重叠面积（C++）

计算两个多边形的重叠面积（C++）

法一

#include <iostream>
#include <cmath>
#include <cstring>

using namespace std;

const int maxn = 300;//maxn表示多边形的最大顶点数量，设为300，即多边形最多有300个顶点
const double eps = 1e-6;//eps表示一个极小的数值，设为1e-6，即0.000001，通常用于判断两个浮点数是否相等。
//位置标识
int dcmp(double x)
{
   
    if(x > eps) return 1;
    return x < -eps ? -1 : 0;
}
// 定义二维点的结构体
struct Point
{
   
    double x, y;
};
// 计算向量(p1, p2)和(p1, p3)的叉积
double cross(Point a, Point b, Point c)  
{
   
    return (a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y);  //叉积公式
}

//计算线段ab和cd的交点坐标
Point intersection(Point a, Point b, Point c, Point d)
{
   
    Point p = a;
    // 计算交点坐标，公式为二维向量叉积的结果
    double t =((a.x-c.x)*(c.y-d.y)-(a.y-c.y)*(c.x-d.x))/((a.x-b.x)*(c.y-d.y)-(a.y-b.y)*(c.x-d.x)); 
    p.x +=(b.x-a.x)*t; 
    p.y +=(b.y-a.y)*t;
    cout << "intersection p.x=" << p.x << ", p.y=" << p.y << endl;
    return p;
}

//计算多边形面积，将多边形拆解成连续三个顶点组合成的多个三角形进行计算，这个循环计算一次其实是计算两次多边形的面积。
double PolygonArea(Point p[], int n)
{
   
    if(n < 3) return 0.0;
    double s = p[0].y * (p[n - 1].x - p[1].x);
    for(int i = 1; i < n - 1; ++ i) {
   
        s += p[i].y * (p[i - 1].x - p[i + 1].x);
        // cout << "p[i-1].x =" << p[i-1].x << ", p[i-1].y=" << p[i-1].y << endl;
        // cout << "p[i].x =" << p[i].x << ", p[i].y=" << p[i].y << endl;
        // cout << "p[i+1].x =" << p[i+1].x << ", p[i+1].y=" << p[i+1].y << endl;
    }
    s += p[n - 1].y * (p[n - 2].x - p[0].x);
    cout << "s =" << s << endl;
    return fabs(s * 0.5);
}
// 计算两个凸多边形的交集面积
double CPIA(Point a[], Point b[], int na, int nb)  //ConvexPolygonIntersectArea
{
   
    Point p[20], tmp[20];//定义两个点集数组，p存储b数组中的所有点，tmp存储在计算交集时得到的新点集
    int tn, sflag, eflag;//定义临时变量：新点集中的点数，起点和终点在连续点之间的位置标志
    memcpy(p,b,sizeof(Point)*(nb));// 把点集b的所有点复制到p中，初始化p数组
    for(int i = 0; i < na && nb > 2; i++)// 遍历a数组中的所有点，如果b中的点数小于等于2，则交集面积为0
    {
   
    	// 计算连续两个点组成的矢量线段在p数组中的位置
        if (i == na - 1) {
   //对于a中的最后一个点，计算a[0]和a[na-1]组成的矢量线段在p数组中的位置
    		sflag = dcmp(cross(a[0], p[0],a[i]));//计算a[0], a[na-1]所在矢量线段和p[0]的位置关系
    	} else {
   //计算a[i], a[i+1]所在矢量线段和p[0]的位置关系
    		sflag = dcmp(cross(a[i + 1], p[0],a[i]));
    	}
        for(int j = tn = 0; j < nb; j++, sflag = eflag)//遍历p数组中的所有点，计算在两个矢量线段之间的点，并存储到新的tmp数组中
        {
   
            if(sflag>=0) {
   
            	tmp[tn++] = p[j];
            }
            if (i == na - 1) {
   //计算下一个连续点在矢量线段的位置
            	if (j == nb -1) {
   //如果已经到达p数组中的最后一个点，则计算a[0], a[na-1]所在矢量线段和p[0]的位置关系
            		eflag = dcmp(cross(a[0], p[0], a[i]));
				} else {
   //计算a[0], a[na-1]所在矢量线段和p[j], p[j+1]的位置关系
					eflag = dcmp(cross(a[0], p[j + 1], a[i])); //计算下一个连续点在矢量线段的位置
				}
			} else {
   
				if (j == nb -1) {
   //计算a[i], a[i+1]所在矢量线段和p[0]的位置关系
					eflag = dcmp(cross(a[i + 1], p[0], a[i]));
				} else {
   //计算a[i], a[i+1]所在矢量线段和p[j], p[j+1]的位置关系
					eflag = dcmp(cross(a[i + 1], p[j + 1], a[i]));
				}
			}
            if((sflag ^ eflag) == -2) {
     //1和-1的异或为-2，也就是两个点分别在矢量线段的两侧
            	if (i == na - 1) {
   
            		if (j == nb -1) {
   
            			tmp[tn++] = intersection(a[i], a[0], p[j], p[0]); //求交点
            		} else {
   
            			tmp[tn++] = intersection(a[i], a[0], p[j], p[j + 1]);
            		}
				} else {
   
					if (j == nb -1) {
   
						tmp[tn++] = intersection(a[i], a[i + 1], p[j], p[0]);
					} else {
   
						tmp[tn++] = intersection(a[i], a[i + 1], p[j], p[j + 1]);
					}
				}
            }
        }
        memcpy(p, tmp, sizeof(Point) * tn);
        nb = tn, p[nb] = p[0];
    }
    if(nb < 3) return 0.0;
    return PolygonArea(p, nb);
}

double SPIA(Point a[], Point b[], int na, int nb)    //SimplePolygonIntersectArea