Составить треугольник наибольшей площади, используя в качестве сторон три отрезка из заданных

#include <bits/stdc++.h>
using namespace std;
int main (){
    double n,p,q=0,s=0,a[2000];
    int k,i,j;
    cin>>n;
    for(int i=1; i<=n; i++)cin>>a[i];
 
    for(i=1; i<=n-2; i++){
        for(j=i+1; j<=n-1; j++){
            for(k=j+1; k<=n; k++){
        p=(a[i]+a[j]+a[k])/2.0;
        q=p*(p-a[i])*(p-a[j])*(p-a[k]);
            if(q>s){s=q;}
            }
        }
    }
    if (s>0) {cout<<fixed<<setprecision(3)<<(sqrt(s));} else {cout<<0;}
    return 0;
}

const int N = 1E3;
int a[N];
 
int comp (const void * a, const void * b) { return ( *(int*)b - *(int*)a ); }
double sq(double x, double y, double z) { if (x > (y+z)) return 0; double p = (x + y + z) / 2; return sqrt(p * (p-x) * (p-y) * (p-z)); }
 
int main() {
    int n; std::cin >> n; for (int i=0; i<n; i++) std::cin >> a[i];
    qsort (a, n, sizeof(int), comp);
    double r = 0;
    for (int i=0; i<n-2; i++) { double s = sq(a[i], a[i+1], a[i+2]); r = std::max(r, s); }
    if (r>0) std::cout << std::fixed << std::setprecision(3) << r; else std::cout << 0;
}