1. 程式人生 > >空間中任意三個點組成三角形面積

空間中任意三個點組成三角形面積

問題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);
 
}