Выберите четыре разные точки, которые являются вершинами квадрата наибольшего периметра

var
   a:array [1..1000] of integer;
   i,j,k,l,n,mi,mj,mk,ml:integer;
   d1,d2,d3,d4,p,mp:real;
{по х}
function
   x(i:integer):real;
begin
   x:=a[2*i-1];
end;
{по y}
function
   y(i:integer):real;
begin
   y:=a[2*i]
end;
{длина отрезка(используем 2 выше описанные функции)}
function
   dl(i,j:integer):real;
begin
   dl:=sqrt(sqr(x(i)-x(j))+sqr(y(i)-y(j)))
end;
{begin}
begin
   writeln('Количество точек: ');
   readln(n);
   for i:=1 to n do
   begin
      write('x[',i,']=');
      readln(a[2*i-1]);
      write('y[',i,']=');
      readln(a[2*i]);
   end;
   for i:=1 to n-3 do
      for j:=i+1 to n-2 do
         for l:=j+1 to n-1 do
            for k:=l+1 to n do
            begin
               d1:=dl(j,i);;
               d2:=dl(l,j);
               d3:=dl(k,l);
               d4:=dl(k,i);
               if (d1=d2) and (d2=d3) and (d3=d4) and (d4=d1) then
                  if (d1<>0) then
                  begin
                  p:=d1+d2+d3+d4;
                  if p>mp then
                     begin
                        mi:=i;
                        mj:=j;
                     mk:=k;
                        ml:=l;
                  mp:=p;
                  end;
                  end;
            end;
   if p<>0 then
   begin
      writeln('Квадрат с максимальным периметром: ');
      writeln('Точки :',mi,':',mj,':',ml,':',mk);
      writeln('Периметр: ',mp:4:2);
   end
   else
      writeln('НЕТУ КВАДРАТА!');
   ReadLn;
 end.

uses crt;
const t=0.000001;
function Dlina(x1,y1,x2,y2:real):real;
begin
Dlina:=sqrt(sqr(x1-x2)+sqr(y1-y2));
end;
var p:array[1..2,1..4] of real;
    d:array[1..6] of real;
    i,j,k:byte;
    x:real;
begin
clrscr;
writeln('Vvedite koordinaty 1 tochki:');
readln(p[1,1],p[2,1]);
writeln('Vvedite koordinaty 2 tochki:');
readln(p[1,2],p[2,2]);
writeln('Vvedite koordinaty 3 tochki:');
readln(p[1,3],p[2,3]);
writeln('Vvedite koordinaty 4 tochki:');
readln(p[1,4],p[2,4]);
k:=0;
for i:=1 to 3 do
for j:=i+1 to 4 do
 begin
  k:=k+1;
  d[k]:=Dlina(p[1,i],p[2,i],p[1,j],p[2,j]);
 end;
for i:=1 to 5 do
for j:=i+1 to 6 do
if d[i]>d[j] then
 begin
  x:=d[i];
  d[i]:=d[j];
  d[j]:=x;
 end;
{длины сторн и диагоналей}
for i:=1 to 6 do
write(d[i]:5:2);
writeln;
if (abs(d[1]-d[2])<t)and(abs(d[2]-d[3])<t)and(abs(d[3]-d[4])<t)
and(abs(d[5]-d[6])<t)
then write('Yes!') else write('No!');
readln
end.

if (abs(d[1]-d[2])<t)and(abs(d[2]-d[3])<t)and(abs(d[3]-d[4])<t)
and(abs(d[5]-d[6])<t)

uses crt;
const p=0.0001;{точность сравнения длин}
      nmax=16;{больше все равно вводить не будешь}
type Point=record {тип - точка}
           x,y:real;
           end;
function Dlina(a,b:Point):real;
begin
Dlina:=sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
end;
var t:array[1..nmax] of Point;{массив точек}
    d:array[1..6] of real;{массив длин сторон и диагоналей}
    n,i,j,k,l,ii,jj,imx,jmx,kmx,lmx:byte;
    x,pr,mx:real;
begin
clrscr;
repeat
writeln('Количество точек от 4 до ',nmax,' n=');;
readln(n);
until n in [4..nmax];
writeln('Введите координаты точек:');
for i:=1 to n do
  begin
   writeln('Точка ',i);
   write('x=');readln(t[i].x);
   write('y=');readln(t[i].y);
  end;
clrscr;
writeln('Координаты:');
write('№');
for i:=1 to n do
write(i:4);
writeln;
write('X:');
for i:=1 to n do
write(t[i].x:4:1);
writeln;
write('Y:');
for i:=1 to n do
write(t[i].y:4:1);
writeln;
writeln;
imx:=0;mx:=0;
for i:=1 to n-3 do
for j:=i+1 to n-2 do
for k:=j+1 to n-1 do
for l:=k+1 to n do
 begin {длины сторон и диагоналей в массив}
  d[1]:=Dlina(t[i],t[j]);
  d[2]:=Dlina(t[i],t[k]);
  d[3]:=Dlina(t[i],t[l]);
  d[4]:=Dlina(t[j],t[k]);
  d[5]:=Dlina(t[j],t[l]);
  d[6]:=Dlina(t[k],t[l]);
  for ii:=1 to 5 do {сортируем по возрастанию}
  for jj:=ii+1 to 6 do
  if d[ii]>d[jj] then
   begin
    x:=d[ii];
    d[ii]:=d[jj];
    d[jj]:=x;
   end; {если первые 4 равны и последние 2 равны и они больше первых}
  if (abs(d[1]-d[2])<p)and(abs(d[2]-d[3])<p)and(abs(d[3]-d[4])<p)
  and(abs(d[5]-d[6])<p)and(d[5]>d[1]) then
    begin
     pr:=4*d[1];{периметр}
     if pr>mx then
      begin
       mx:=pr;{макс. периметр}
       imx:=i;{номера веpшин}
       jmx:=j;
       kmx:=k;
       lmx:=l;
      end;
    end;
 end;
if imx=0 then write('Квадратов нет!')
else
 begin
  writeln('Квадрат с максимальным периметром образуют точки:');
  writeln(imx,'- x=',t[imx].x:0:1,' y=',t[imx].y:0:1);
  writeln(jmx,'- x=',t[jmx].x:0:1,' y=',t[jmx].y:0:1);
  writeln(kmx,'- x=',t[kmx].x:0:1,' y=',t[kmx].y:0:1);
  writeln(lmx,'- x=',t[lmx].x:0:1,' y=',t[lmx].y:0:1);
 end;
readln
end.