Задача на поиск самой близкой пары точек располагающейся по разные стороны от оси ОХ

#include<iostream>
#include<math.h>
 
class Point
{
  private:
      int x;
      int y;
  public:
    Point(int _x, int _y) : x(_x), y(_y) {};
    int getx() { return x; };
    int gety() { return y; };
    friend double distance(Point*, Point*);
};
 
double distance(Point* r, Point* q)
{
  double dist=(r->x-q->x)*(r->x-q->x)+(r->y-q->y)*(r->y-q->y);
  return(sqrt(dist));
}
 
double first_distance(Point** p)
{
  double d=0.0;
  int i=0; int j=0;
  Point* temp;
  for(i=0; p[i]!=NULL; i++)
  {
    temp=p[i];
    for(j=0; p[j]!=NULL; j++)
    {
      if((temp->gety()*p[j]->gety())<0)
        {
          d=distance(temp, p[j]);
          return d;
        }
    }
  }
  if (d==0.0)
    exit(-1);
}
 
int main(int argc, char* argv[])
{
  double dist_temp;
  double d;
  Point** p;
  Point* temp;
  int x,y;
  int i=0; int j=0; int j1;
  int k=0;
 
  if(argc<2)
    return(-1);
 
  p=new Point*[argc];
  while(++i<argc)
  {
    sscanf(argv[i], "%d%*c%d", &x, &y);
    p[j]=new Point(x,y);
    j++;
  }
  p[j]=NULL;
 
  d=first_distance(p);
  for(i=0; p[i]!=NULL; i++)
  {
    temp=p[i];
    for(j=0; p[j]!=NULL; j++)
    {
      if((temp->gety()*p[j]->gety())<0)
      {
        dist_temp=distance(temp,p[j]);
        if(dist_temp<d)
          d=dist_temp;
      }
    }
  }
  printf("The minimal distance is %f\n", d);
 
  j=0;
  for(i=0; p[i]!=NULL; i++)
  {
    temp=p[i];
    for(j1=j; p[j]!=NULL; j1++)
    {
      if((temp->gety()*p[j1]->gety())<0)
      {
        dist_temp=distance(temp,p[j1]);
        if(dist_temp==d)
        {
        ++k;
        printf("The %d pair is (%d;%d), (%d;%d)\n", k, temp->getx(), temp->gety(), p[j1]->getx(), p[j1]->gety());
        }
      }
    }
    j++;
  }
 
  for(j=0; p[j]!=NULL; j++)
    delete p[j];
  delete [ ]p;
 
  return 0;
}