main.cpp

#include "Image.h"
#include "Welzl.h"
#include <time.h>

// Ãåíåðàöèÿ òî÷åê â ðåêòàíãëå.
std::vector<Point> Generate(const Point &lb, const Point &rt, unsigned int count)
{
  std::vector<Point> points;
  points.reserve(count);

  Size size(static_cast<unsigned int>(rt.x - lb.x), static_cast<unsigned int>(rt.y - lb.y));

  for(unsigned int i = 0; i < count; ++i)
  {
    points.push_back(Point(lb.x + rand() % size.width, lb.y + rand() % size.height));
  }

  return points;
}


int main()
{
  srand(static_cast<unsigned int>(time(0)));
  Size size(1000, 1000);

  Image image;
  image.Resize(size);
  image.Fill(0xFFFFFFFF);

  std::vector<Point> points = Generate(Point(500, 300), Point(700, 700), 10000000);

  Welzl w;

  printf("begin\n");

  Circle c = w.Process(points);

  printf("end\n");

  image.DrawCircle(c.center, c.radius, 0xFF0000FF);

  for(auto it = points.begin(); it != points.end(); ++it)
  {
    image.DrawPoint(*it, 0x00FF00FF);
  }

  image.Save("img.png");

  return 0;
}