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

#include <iostream>
#include <vector>
struct point {
    double x;
    double y;
};
struct polygon {
    int i1;
    int i2;
    int i3;
    int count;
};
int FindCountVertexIn(const point& a, const point& b, const point& c, const point& pt) {
    if (
        ((b.x - a.x) * (pt.y - a.y) - (b.y - a.y) * (pt.x - a.x) > 0) &&
        ((c.x - b.x) * (pt.y - b.y) - (c.y - b.y) * (pt.x - b.x) > 0) &&
        ((a.x - c.x) * (pt.y - c.y) - (a.y - c.y) * (pt.x - c.x) > 0)
    )
        return 1;
    return 0;
}
int main() {
    const int SIZE = 10;
    point mass[SIZE] = {
                        {0, 0}, {0, 1}, {1, 1}, {1, 0}, {2, 2},
                        {-1, -1}, {-1, 0}, {0, -1}, {-1, 1}, {1, -1}
                        };
    polygon pl_temp;
    std::vector<polygon> coll;
    for(int i = 0; i < SIZE - 2; ++i)
        for(int k = i + 1; k < SIZE - 1; ++k)
            for(int m = k + 1; m < SIZE; ++m) {
                pl_temp.i1 = i;
                pl_temp.i2 = k;
                pl_temp.i3 = m;
                pl_temp.count = 0;
                for(int s = 0; s < SIZE; ++s) {
                    if (
                        (s == i) ||
                        (s == k) ||
                        (s == m)
                        )
                        continue;
                    pl_temp.count += FindCountVertexIn(
                                                        mass[i],
                                                        mass[k],
                                                        mass[m],
                                                        mass[s]
                                                        );
                }
                coll.push_back(pl_temp);
            }
    std::cout << "====================== debug ======================\n";
    for(int n = 0; n < coll.size(); ++n)
        std::cout << coll[n].i1 << ' ' << coll[n].i2 << ' ' << coll[n].i3 << ' ' << coll[n].count << '\n';
    std::cout << "==================== end debug ====================\n";
    int id_max = 0;
    for(int n = 0; n < coll.size(); ++n)
        if (coll[n].count > coll[id_max].count)
            id_max = n;
    std::cout << "polygon(x, y, z, count) : ";
    std::cout << coll[id_max].i1 << ' '
              << coll[id_max].i2 << ' '
              << coll[id_max].i3 << ' '
              << coll[id_max].count << '\n';
 
    return 0;
}