Найти количество точек в каждой полуплоскости

/* Verevkin, CodeBlocks, кодировка файла: UTF-8 BOM */
// d:\Soft\CodeBlocks\share\CodeBlocks\templates\wizard\console\cpp\main.cpp
 
#include <stdlib.h>
#include <stdio.h>
#include <conio.h>
#include <string.h>
#include <math.h>
#include <assert.h>
 
typedef struct
{
  double x, y;
} TPoint;
 
#define PRECISION   1e-6 // точность (число, меньшее по модулю, считается за 0)
#define ZERO(x)     (fabs(x) < PRECISION) // примерно равно нулю
#define EQUAL(a, b) ZERO((a) - (b)) // примерно равно
#define SIGN(x)     (ZERO(x) ? 0 : ((x) > 0 ? 1 : -1))  // знак числа
 
// уравнение прямой y = kx + c по двум точкам.
// возвращает false, если это вертикальная линия (k = ∞)
bool line_eq(TPoint a, TPoint b, double &k, double &c)
{
  if (EQUAL(a.x, b.x)) return false;
  k = (b.y - a.y) / (b.x - a.x);
  c = a.y + a.x * (a.y - b.y) / (b.x - a.x);
  return true;
}
 
// функция возвратит -1, если точка p находтися слева от прямой,
// проведённой через точки a и b, +1 - если справа, и 0, если
// точка лежит на прямой.
int foo(TPoint a, TPoint b, TPoint p)
{
  // защита от дурака: 2 точки для построения оси не должны совпадать.
  assert(!EQUAL(a.x, b.x) || !EQUAL(a.y, b.y));
 
  // коэффициенты уравнения прямой
  double k, c;
 
  if (!line_eq(a, b, k, c))
    return SIGN(p.x - a.x); // вырожденный случай - ось вертикальная
 
  // выше прямой (справа, +1), ниже прямой (слева -1) или лежит на ней (0)
  // в зависимости от знака тангенса угла наклона.
  return (k < 0) ? SIGN(p.y - k * p.x - c) : SIGN(k * p.x + c - p.y);
}
 
// строка уравнения прямой "y = kx + c" по двум точкам.
char* line2str(TPoint a, TPoint b)
{
  static char buf[128];
  double k, c;
 
    if (!line_eq(a, b, k, c)) sprintf(buf, "x = %g", a.x); // вертикальная линия
    else if (!ZERO(k))
         switch (SIGN(c))
         {
           case -1: sprintf(buf, "y = %gx - %g", k, -c); break;
           case +1: sprintf(buf, "y = %gx + %g", k, c); break;
           case  0: sprintf(buf, "y = %gx", k); break;
           default: assert(false); // не может быть
         }
       else sprintf(buf, "y = %g", c);  // горизонтальная линия
  return buf;
}
 
// ---------------------------------------------------------------
#define AX -2
#define AY  2
#define BX -2
#define BY -2
 
int main(int argc, char** argv)
{
    system("chcp 65001 & cls");
 
    const char* RESULT_STR[] = {"левее прямой", "лежит на прямой", "правее прямой"};
 
    // точки, заданные от балды
    TPoint pts[] =
    {
      {-6, -1},
      { 0,  1},
      { 1,  0},
      { 1,  1},
      {-2, -2},
      {4.5, 0.5},
      {-1,  1},
      {-1, -1},
      {-1, -2},
      {-4,  4},
      {-2, -5}
    };
 
    // две точки, через которые проходит ось.
    // По условию их вводит юзер, но я ленивая жопа,
    // поэтому задал как константы.
    TPoint a = {AX, AY}, b = {BX,  BY};
 
    // кол-во точек слева/снизу от оси, лежащих на оси, справа/сверху от оси
    unsigned counts[3] = {0, 0, 0};
 
    printf(" > Через точки A = (%g,%g) и B = (%g,%g) проходит прямая %s.\n\n",
           a.x, a.y, b.x, b.y, line2str(a, b));
 
    for (auto p: pts)
    {
      int position = foo(a, b, p); // положение точки (-1 | 0 | +1)
      counts[1 + position]++; // подсчёт количества
      printf(" # Точка (%g,%g) - %s.\n", p.x, p.y, RESULT_STR[1 + position]);
    }
 
    printf("\n $ Точек левее: %u; правее: %u; лежащих на прямой: %u.\n", counts[0], counts[2], counts[1]);
 
    printf("\n >> Нажми что-нибудь для выхода... ");
    getch();
    return 0;
}