编译器:C++ (g++) Particle类型用于表示处于三维空间中的质点，其包括如下成员：① 三个浮点数成员x、y和z表示质点在三维空间中的坐标；② 浮点数成员mass表示质点的质量；③ 接受空间坐标和质量参数的构造函数；④ 自定义operator-()操作符函数用于计算两个质点间的万有引力，即表达式p1 – p2调用执行p1的operator-()操作符函数，计算质点p1和p2间的万有引力。 万有引用公式F = GMm/r²中，万有引力常数G = 6.67×10-11N·m²/kg²。 请依据上述描述，设计Particle类，使得下述代码能正常执行。 注意：天文数字都很大，建议浮点数全部采用long double类型。 裁判测试程序样例： #include <iostream> #include <cmath> using namespace std; //在此处定义Particle类及其重载operator+()操作符函数 int main() { long double x, y, z, mass; cin >> x >> y >> z >> mass; Particle p1(x,y,z,mass); cin >> x >> y >> z >> mass; const Particle p2(x,y,z,mass); long double force = p1 - p2; printf("Gravity between two objects = %.1Lf N",force); return 0; } 输入样例： 0 0 0 5.965e15 0 0 38400000 7.342e14 输出样例： Gravity between two objects = 198101.7 N 说明：输入输出中的所有数据均为标准单位，坐标为米，质量为kg。 请注意：函数题只需要提交相关代码片段，不要提交完整程序。

代码如下：

#include <iostream>
#include <cmath>
using namespace std;

class Particle{
private:
    long double x, y, z, mass;
public:
    Particle(long double x=0, long double y=0, long double z=0, long double mass=0){
        this->x = x;
        this->y = y;
        this->z = z;
        this->mass = mass;
    }
    long double operator-(const Particle& p2){
        const long double G = 6.67e-11;
        long double dx = x - p2.x;
        long double dy = y - p2.y;
        long double dz = z - p2.z;
        long double r = sqrt(dx*dx + dy*dy + dz*dz);
        return G * mass * p2.mass / (r * r);
    }
};

int main(){
    long double x, y, z, mass;
    cin >> x >> y >> z >> mass;
    Particle p1(x,y,z,mass);
    cin >> x >> y >> z >> mass;
    const Particle p2(x,y,z,mass);
    long double force = p1 - p2;
    printf("Gravity between two objects = %.1Lf N",force);
    return 0;
}