계속해서 TLE가 나고 있어 질문 올립니다.

지금 제 방법은 임의의 2점을 잡고나서 나머지 2점은 이분검색으로 찾는데 시간복잡도는 2N^2logN입니다. N이 2만이므로 당연히 TLE가 나겠지만, 더 줄일 수 있는 방법을 몰라 이렇게 글을 올립니다.

줄일 수 있는 방법이나 다른 획기적인 방안이 없을까요?

// 밑에는 제 코드입니다.

#include
#include
#define M_abs(X, Y) (X-Y) > 0 ? (X-Y) : (Y-X)
using std::sort;

struct XY{

int x, y;

};
typedef struct XY xy;
xy sq[20010];
int sq_seq[20010][2];
bool compare ( const xy &fi, const xy &se ){

if( fi.x == se.x )
    return fi.y < se.y;
return fi.x < se.x;

}

bool find_point( int p, int x, int y )
{
int fp = p - sq_seq[p][0];
int rp = p + sq_seq[p][1];
int mp;

while( fp <= rp ){

    mp = (fp + rp) / 2;
    if( sq[mp].y == y )
        return true;
    else if( sq[mp].y > y ) 
        rp = mp-1;
    else
        fp = mp+1;
}
return false;

}
int main()
{
int T;
int N, i, j;
int length_X, length_Y;
int max, min;
freopen("input.txt","rt",stdin);
freopen("output.txt","wt",stdout);

scanf("%d", &T);
while( T-- ){

    max = -1; min = 0x7fffffff;
    scanf("%d", &N);
    for( i=0;i<N;i++ ){

        scanf("%d %d", &sq[i].x, &sq[i].y);
    }

    sort( sq, sq+N, compare );

    /*for( i=0;i<N;i++ ){

        printf("%d ", sq[i].x);
    }
    printf("\n");*/
    sq_seq[0][0] = 0;
    sq_seq[N-1][1] = 0;

    for( i=1;i<N;i++ ){

        if( sq[i].x != sq[i-1].x )
            sq_seq[i][0] = 0;
        else
            sq_seq[i][0] = sq_seq[i-1][0] + 1;

        if( sq[N-i-1].x != sq[N-i].x )
            sq_seq[N-i-1][1] = 0;
        else
            sq_seq[N-i-1][1] = sq_seq[N-i][1] + 1;
    }

    /*printf("----------------sq_seq--------------\n");
    for( i=0;i<N;i++ ){

        printf("%d ",sq_seq[i][0] );
    }
    printf("\n");
    for( i=0;i<N;i++ ){

        printf("%d ",sq_seq[i][1]);
    }
    printf("\n");
    printf("----------------sq_seq_end----------\n");*/

    for( i=0;i<N;i++ ){

        for( j=i+1;j<N;j++ ){

            length_X = M_abs(sq[i].x, sq[j].x);
            length_Y = M_abs(sq[i].y, sq[j].y);

            if( length_X != length_Y ) continue;
            // sq[j].x AND sq[i].y 
            // sq[i].x AND sq[j].y;

            if( find_point( j, sq[j].x , sq[i].y ) && find_point( i, sq[i].x, sq[j].y ) ){

                max = max < length_X ? length_X : max;
                min = min > length_X ? length_X : min;
            }
        }
    }
    printf("%d %d\n",min,max);
}

return 0;

}