Решение задачи Array 138

#include <stdio.h>
 
int main(void)
{
    float a[10][2];
    int n;
 
    printf("N: ");
    scanf("%i",&n);
 
    int i;
    for (i=0; i<n; ++i){
        printf("a[%i]:\n",i+1);
        printf("  x : ");
        scanf("%f",&a[i][0]);
        printf("  y : ");
        scanf("%f",&a[i][1]);
    }
 
    int p1=0,p2=1,p3=2,i2,i3;
    float p,pmin=0;
    pmin+= sqrt(pow(a[0][0]-a[1][0],2)+pow(a[0][1]-a[1][1],2));
    pmin+= sqrt(pow(a[0][0]-a[2][0],2)+pow(a[0][1]-a[2][1],2));
    pmin+= sqrt(pow(a[1][0]-a[2][0],2)+pow(a[1][1]-a[2][1],2));
 
 
    for (i=0; i<n; ++i)
        for (i2=i+1; i2<n; ++i2)
            for (i3=i2+1; i3<n; ++i3){
                p=0;
                p+= sqrt(pow(a[i][0]-a[i2][0],2)+pow(a[i][1]-a[i2][1],2));
                p+= sqrt(pow(a[i][0]-a[i3][0],2)+pow(a[i][1]-a[i3][1],2));
                p+= sqrt(pow(a[i2][0]-a[i3][0],2)+pow(a[i2][1]-a[i3][1],2));
                if (p < pmin) {
                    p1=i;
                    p2=i2;
                    p3=i3;
                    pmin=p;
                }
            }
 
 
    printf("P: %f\n",pmin);
    printf("A %i :\n x: %f\n y: %f\n",p1+1,a[p1][0],a[p1][1]);
    printf("A %i :\n x: %f\n y: %f\n",p2+1,a[p2][0],a[p2][1]);
    printf("A %i :\n x: %f\n y: %f\n",p3+1,a[p3][0],a[p3][1]);
    return 0;
}

program Array138;
var
 A: array [1..10,1..2] of real;
 P,Pmin:real;
 N,i,i2,i3,Point1,Point2,Point3: integer;
 
begin
 Write('N: ');
 Readln(N);
 
 Writeln('A: ');
 for i:=1 to N do
  begin
   write(i,'.x :');
   readln(a[i,1]);
   write(i,'.y :');
   readln(a[i,2]);
  end;
 
 Pmin:=Pmin+sqrt(sqr(A[1,1]-A[2,1])+sqr(A[1,2]-A[2,2]));
 Pmin:=Pmin+sqrt(sqr(A[2,1]-A[3,1])+sqr(A[2,2]-A[3,2]));
 Pmin:=Pmin+sqrt(sqr(A[1,1]-A[3,1])+sqr(A[1,2]-A[3,2]));
 Point1:=1;
 Point2:=2;
 Point3:=3;
 
 for i:=1 to N do
  for i2:=i+1 to N do
   for i3:=i2+1 to N do
    begin
     P:=0;
     P:=P+sqrt(sqr(A[i,1]-A[i2,1])+sqr(A[i,2]-A[i2,2]));
     P:=P+sqrt(sqr(A[i2,1]-A[i3,1])+sqr(A[i2,2]-A[i3,2]));
     P:=P+sqrt(sqr(A[i,1]-A[i3,1])+sqr(A[i,2]-A[i3,2]));
     if Pmin>P then
      begin
       Pmin:=P;
       Point1:=i;
       Point2:=i2;
       Point3:=i3;
      end;
 
    end;
 
  Writeln ('Pmax:',Pmin,' (',A[Point1,1],':',A[Point1,2],')(',A[Point2,1],':',A[Point2,2],')(',A[Point3,1],':',A[Point3,2],')');
end.