空间中任意三个点组成三角形面积

问题1：遍历所有可能的3个点

1 2 3

for(int i = 0; i < n; i++)    for(int j = i + 1; j < n; j++)      for(int k = j + 1; k < n; k++)

问题2：判断3个点是否能组成三角形

1 2 3 4 5 6 7 8

double a = dis(i, j);     //计算两点距离 double b = dis(i, k); double c = dis(k, j); if(a < (b + c) && b < (a + c) && c < (a + b)) {   return true; } return false;

问题3：判断颜色是否相同或全不同

1 2 3 4 5 6 7 8 9

if(V[i].c == V[j].c && V[j].c == V[k].c) {    return true; } else if(V[i].c != V[j].c &&V[i].c != V[k].c&&V[k].c != V[j].c) {   return true; } return false;

问题4：计算三角形面积 (海伦公式 百度~)

double a = L3(i, j);//求三边的长度

double b = L3(i, k);

double c = L3(k, j);

double p = (a + b + c) / 2;

return sqrt(p * (p - a) * (p - b) * (p - c));

#include<bits/stdc++.h> 
using namespace std;
 
/*
3
R 0 0 0
R 0 4 0
B 0 0 3
*/
 
struct node{
    char c;
    int x, y, z;
};
 
double MAX = 0;
int n;
vector<node> V;
 
//计算两点之间的举例
double L3(int i, int j){
    return sqrt(double((V[i].x - V[j].x)*(V[i].x - V[j].x) + (V[i].y - V[j].y)*(V[i].y - V[j].y) + (V[i].z - V[j].z)*(V[i].z - V[j].z)));
}
 
//判断是否是三角形
bool isSan(int i, int j, int k){
    double a = L3(i, j);
    double b = L3(i, k);
    double c = L3(k, j);
    if (a < (b + c) && b < (a + c) && c < (a + b))
    {
        return true;
    }
    return false;
}
 
//判断颜色
bool isColour(int i, int j, int k){
    if (V[i].c == V[j].c && V[j].c == V[k].c)
    {
        return true;
    }
    else if (V[i].c != V[j].c &&V[i].c != V[k].c&&V[k].c != V[j].c)
    {
        return true;
    }
    return false;
}
 
//计算三角形面积
double CmputeArea(int i, int j, int k){
    double a = L3(i, j);
    double b = L3(i, k);
    double c = L3(k, j);
    double p = (a + b + c) / 2;
    return  sqrt(p * (p - a) * (p - b) * (p - c));
}
 
//遍历所有可能的三个点
void run(){
    for (int i = 0; i < n; i++)
    {
        for (int j = i + 1; j < n; j++)
        {
            for (int k = j + 1; k < n; k++)
            {
                if (isSan(i, j, k) && isColour(i, j, k))
                {
                    double tArea = CmputeArea(i, j, k);
                    if (tArea > MAX)
                    {
                        MAX = tArea;
                    }
                }
            }
        }
    }
}
 
int main(){
    cin >> n;
 
    for (int i = 0; i < n; i++)
    {
        node t;
        cin >> t.c >> t.x >> t.y >> t.z;
        V.push_back(t);
    }
    run();
 
    printf("%.5lf", MAX);
 
}