деление множество точек на две равные части

#include <iostream>
#include <vector>
 
const double circleRadius = 4;
std::pair<double, double> circleCenter(5,5);
 
void initVector(std::vector< std::pair<double, double> >& v);
inline bool comparePoint(const std::pair<double, double>& p1, const std::pair<double, double>& p2, const std::pair<double, double>& p3);
inline double length(const std::pair<double, double>& p1, const std::pair<double, double>& p2);
void print(std::vector< std::pair<double, double> >& v);
void sort(std::vector< std::pair<double, double> >& v);
 
int main()
{
    std::vector< std::pair<double, double> > v;
    std::vector< std::pair<double, double> > copy;
    std::vector< std::pair<double, double> >::iterator it;
    
    //заливаем множество точек в вектор
    initVector(v);
 
    //оставляем только те точки, которые в окружность попадают. Окружность я определяю сам
    for(it = v.begin(); it != v.end(); ++it)
    {
        if(circleRadius >= length(circleCenter, *it))
        {
            copy.push_back(*it);
        }
    }
 
    //сортирую все точки которые попали в окружность по Х-координате
    sort(copy);
    //print(copy);
 
    //поделив вектор, пополам получим точку, которая лежит по середине остальных. Ну и вторую получим прибавив + 1 к середине вектора, можно и - 1 сделать :)
    //т. е. я поделил по вертикали множество точек!
    //в условии сказано что приблизительно должно быть равное количество.
    std::cout << copy[copy.size()/2 - 1].first << copy[copy.size()/2 - 1].second << std::endl;
    std::cout << copy[copy.size()/2].first << copy[copy.size()/2].second << std::endl;
 
    return 0;
}
 
void initVector(std::vector< std::pair<double, double> >& v)
{
    v.push_back(std::pair<double, double>(2,2));
    v.push_back(std::pair<double, double>(5,8));
    v.push_back(std::pair<double, double>(7,3));
    v.push_back(std::pair<double, double>(6,3));
    v.push_back(std::pair<double, double>(7,8));
    v.push_back(std::pair<double, double>(3,5));
    v.push_back(std::pair<double, double>(6,6));
    v.push_back(std::pair<double, double>(10,10));
}
 
bool comparePoint(const std::pair<double, double>& p1, const std::pair<double, double>& c, const std::pair<double, double>& p2)
{
    return ((c.first - p1.first)*(p2.second - c.second) - (p2.first - c.first)*(c.second - p1.second)) == 0;
}
 
double length(const std::pair<double, double>& p1, const std::pair<double, double>& p2)
{
    return sqrt((p2.first - p1.first)*(p2.first - p1.first) + (p2.second - p1.second)*(p2.second - p1.second));
}
 
void print(std::vector< std::pair<double, double> >& v)
{
    std::vector< std::pair<double, double> >::iterator it;
    for(it = v.begin(); it != v.end(); ++it)
    {
        std::cout << '[' << it->first << ',' << it->second << ']' << std::endl;
    }
}
 
void sort(std::vector< std::pair<double, double> >& v)
{
    for ( int i = 0; i < v.size(); i++ )
        for ( int j = 1; j < v.size(); j++ )
            if ( v[j].first < v[j-1].first )
                std::swap(v[j], v[j - 1]);
}

#include <iostream>
#include <cmath>
#define M 1000
#define R 100
#include <cstdlib>
#include <ctime>
#include <vector>
using namespace std;
struct coordinati{
    float a;
    float b;
    void set_coord(float& x, float& y)
    {
        a = x;
        b = y;
    }
};
int main()
{
        coordinati coord;
        vector<coordinati> vec;
    srand(time(NULL));
    int i = 0;
 
    float x, y;
    while(i != M)
    {
        x = (rand()%(R *2)) - R;
        y = (rand()%(R *2)) - R;
 
        if(sqrt(x*x + y*y) <= R)
        {
            coord.set_coord(x,y);
            vec.push_back(coord);
        }
        i++;
    }
    for(vector<coordinati>::size_type i = 0; i != vec.size(); i++)
    {
                for(vector<coordinati>::size_type j = 0; j != vec.size(); j++)
                {
                    if(vec[i].a == (- vec[j].a)&& vec[i].b == -(vec[j].b))
                    cout << "Cherez tochku a[" << vec[i].a <<"][" << vec[i].b << "]" <<
                    " i tuchku b[" << vec[j].a <<"][" << vec[j].b <<"]" <<
                    "mozhno razdelitj krug popolam" << endl;
 
                }
    }
return 0;
}