Задача на построение прямой

#include <iostream>
#include <utility>
#include <vector>
#include <cmath>
 
constexpr double eps = 1e-8;
 
struct vec2 { double x, y; };
 
double cross(const vec2 & lh, const vec2 & rh)
{ return lh.x * rh.y - rh.x * lh.y; }
 
vec2 operator-(const vec2 & lh, const vec2 & rh)
{ return { lh.x - rh.x, lh.y - rh.y }; }
 
int signImpl(double val)
{ return val < 0.0 ? -1 : 1; }
 
int sign(double val)
{ return std::abs(val) < eps ? 0 : signImpl(val); }
 
template <class It>
std::pair<It, It> findLine(It first, It last)
{
    for(auto prev = first; prev != last; ++prev)
    {
        for(auto next = std::next(prev); next != last; ++next)
        {
            int counter{};
            auto const a = *next - *prev;
            for (auto curr = first; curr != last; ++curr)
            {
                if (curr != prev && curr != next)
                {
                    auto const b = *curr - *prev;
                    counter += sign(cross(a, b));
                }
            }
            if (counter == 0)
                return {prev, next};
        }
    }
    return {last, last};
}
 
std::vector<vec2> readPoints()
{
    std::vector<vec2> result;
    // TO DO : сам
    return result;
}
 
int main()
{
    auto points = readPoints();
    auto result = findLine(points.begin(), points.end());
    if (result.first != points.end() && result.second != points.end())
    {
        // распечатать результат (result.first->x и т.д.)
    }
    else
    {
        // Точек не найдено
    }
}