#include <iostream>
#include <vector>

#include "CGAL/Exact_predicates_exact_constructions_kernel.h"
typedef CGAL::Exact_predicates_exact_constructions_kernel Kernel;
typedef Kernel::Point_2 Point_2;

#include "CGAL/Polygon_2.h"
typedef CGAL::Polygon_2<Kernel> Polygon_2;

#include "CGAL/Polygon_set_2.h"
typedef CGAL::Polygon_set_2<Kernel> Polygon_set_2;

int main(int argc, char* argv[])
{
    auto startTime = std::chrono::high_resolution_clock::now();
    const size_t polygonSize = 2000;
    std::vector<Polygon_2> vecPolygons(polygonSize);
    for (size_t i = 0; i < polygonSize; ++i)
    {
        vecPolygons[i].push_back(Point_2(i - 0.5, i - 0.5));
        vecPolygons[i].push_back(Point_2(i + 0.5, i - 0.5));
        vecPolygons[i].push_back(Point_2(i + 0.5, i + 0.5));
        vecPolygons[i].push_back(Point_2(i - 0.5, i + 0.5));
    }
    Polygon_set_2 polygonSet;
    polygonSet.join(vecPolygons.begin(), vecPolygons.end());

    for (size_t i = 0; i < polygonSize; ++i)
    {
        polygonSet.do_intersect(vecPolygons[i]);
    }
    
    auto endTime = std::chrono::high_resolution_clock::now();
    auto diffTime = endTime - startTime;
    std::cout << "Total cost: " << diffTime.count() * 1e-9 << std::endl;
    return 0;
}