closest pair 알고리즘 문제 관련 질문입니다.

/*
dev cpp 4.9.9.2 에서 작성했습니다.
편의상 전체 평면을 P, 분할된 좌측 평면을 PL, 우측평면을 PR이라고 하겠습니다.
*/ 
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
//두점간의 거릴를 구하는 함수를 매크로함수로 구현했습니다. 
#define DIST(X,Y,P,Q) ((((X)-(P))*((X)-(P)))+(((Y)-(Q))*((Y)-(Q)))) 
//절대값을 반환하는 함수 입니다. MINOR(PL의 최소거리,PR의 최소거리) 안에 있는 square을 만들때 사용합니다. 
#define ABSOLUTE(A) ((A)<(0))?(-1*(A)):(A)  
//두값중 작은 값을 반환하는 함수입니다.
#define MINOR(A,B) ((A)<(B))?(A):(B) 
using namespace std;
struct POINT{long long xpos;long long ypos;};
 
long long CNT=0;
 
//점이 저장될 변수입니다. 
vector<POINT> SET;
//PL에서 1점, PR에서 한점을 뽑을때 그 점들이 위치할 수 있는 집합을 저장할 변수입니다. 
vector<POINT> SQUARE;
long long MIN=1000000000;
//SET을 x좌표 오름 차순으로 정렬하기 위한 함수 입니다.
bool comp(const POINT P1,const POINT P2){return P1.xpos<P2.xpos;} 
 
//PL,PR을 가로지르며 생성되는 직선을 조사해볼 집합을 생성하는 함수입니다.
//현재 점의 인덱스(CDEX)가 주어졌을때, CDEX로 부터 상하좌우로 MINOR(PL의 최소거리,PR의 최소거리)내의 점들을 조사하여 square에 삽입합니다. 
void mkSQUARE(int current,double min){
     int i,len=SET.size();
     for(i=current-1;i>=0;i--){
                               if((SET[current].xpos-SET[i].xpos)>min){break;}
                               if(ABSOLUTE(SET[current].ypos-SET[i].ypos)<=min){SQUARE.push_back(SET[i]);}
     }
     for(i=current+1;i<len;i++){
                                if((SET[i].xpos-SET[current].xpos)>min){break;}
                                if(ABSOLUTE(SET[i].ypos-SET[current].ypos)<=min){SQUARE.push_back(SET[i]);}
     }
     return;
}
//구해진 square에 대해 최소거리를 brute force로 계산합니다. 
void combi(int N){
     int i,j;
     for(i=0;i<N-1;i++){
                        for(j=i+1;j<N;j++){
                                           MIN=MINOR(MIN,DIST(SQUARE[i].xpos,SQUARE[i].ypos,SQUARE[j].xpos,SQUARE[j].ypos));
     }}
     return;
} 
 
//SET의 원소가 3개 이하일경우 BRUTEFORCE를 사용합니다. 
long long small(int LDEX,int RDEX){
     long long ret=1000000000;
     if((RDEX-LDEX)==1){ret=MINOR(ret,DIST(SET[LDEX].xpos,SET[LDEX].ypos,SET[LDEX+1].xpos,SET[LDEX+1].ypos));}
     if((RDEX-LDEX)==2){
                        ret=MINOR(ret,DIST(SET[LDEX].xpos,SET[LDEX].ypos,SET[LDEX+1].xpos,SET[LDEX+1].ypos));
                        ret=MINOR(ret,DIST(SET[LDEX+2].xpos,SET[LDEX+2].ypos,SET[LDEX+1].xpos,SET[LDEX+1].ypos));
     }
     return ret;
}
//답을 출력하는 함수입니다.
//구간의 시작 인덱스 LDEX와 끝 인덱스 RDEX를 매개변수로 취합니다. 
long long mkANS(int LDEX,int RDEX){
     long long ret=1000000000;
     //SET의 원소가 3개 이하인 경우 small을 호출합니다. 
     if((RDEX-LDEX)<3){return small(LDEX,RDEX);}
     //SET의 원소가 3개보다 많을때 다음을 호출합니다. 
     else if(LDEX<RDEX){
                        //ret에 MINOR(PL의 최소거리,PR의 최소거리)을 저장합니다. 
                        ret=MINOR(mkANS(LDEX,(LDEX+RDEX)/2),mkANS((LDEX+RDEX)/2+1,RDEX));
                        MIN=1000000000;SQUARE.clear();
                        //가운데 점의 인덱스 (LDEX+RDEX)/2와 MINOR(PL의 최소거리,PR의 최소거리)인 sqrt(ret)을 이용하여 SQUARE을 구성합니다. 
                        mkSQUARE((LDEX+RDEX)/2,sqrt(ret));
                        //구해진 square에 대해 brute force로 최소값을 구합니다. 
                        combi(SQUARE.size());
                        ret=MINOR(ret,MIN);
     }
     return ret;
}
int main(void){
    int i,N;
    POINT temp;
    cin>>N;
    for(i=0;i<N;i++){
                     cin>>temp.xpos>>temp.ypos;
                     SET.push_back(temp);
    }
    //입력받은 SET을 x오름 차순으로 정렬합니다. 
    sort(SET.begin(),SET.end(),comp);
    cout<<mkANS(0,SET.size()-1)<<endl;
    system("PAUSE");
    return 0;
}