Найти наименьший периметр треугольника

#include<iostream.h>
#include<stdlib.h>
#include<math.h>
#include<conio.h>
main()
{
const N=10;
int X[N],Y[N];
int i,j,k;
randomize();
//zapolnenie massivov sluchainym obrazom
for (i=0;i<N;i++)
  {
   X[i]=random(41)-20;
   Y[i]=random(41)-20;
  }
//Vyvod pervonachal'nyh dannyh
for (i=0;i<N;i++)
   cout<<i+1<<"Koordinaty tochki:"<<endl
       <<"("<<X[i]<<";"<<Y[i]<<")"<<endl;
double a,b,c;
for (i=0;i<N-2;i++)
  for (j=i+1;j<N-1;j++)
    for (k=j+1;k<N;k++)
 
a=sqrt((X[i]-X[j])*(X[i]-X[j])+(Y[i]-Y[j])*(Y[i]-Y[j]));
b=sqrt((X[j]-X[k])*(X[j]-X[k])+(Y[j]-Y[k])*(Y[j]-Y[k]));
c=sqrt((X[i]-X[k])*(X[i]-X[k])+(Y[i]-Y[k])*(Y[i]-Y[k]));
 
 
double min=0;
int i_min,j_min,k_min;
for (i=0;i<N-2;i++)
  for (j=i+1;j<N-1;j++)
    for (k=j+1;k<N;j++)
         a=sqrt((X[i_min]-X[j_min])*(X[i_min]-X[j_min])+(Y[i_min]-Y[j_min])*(Y[i_min]-Y[j_min]));
         b=sqrt((X[j_min]-X[k_min])*(X[j_min]-X[k_min])+(Y[j_min]-Y[k_min])*(Y[j_min]-Y[k_min]));
         c=sqrt((X[i_min]-X[k_min])*(X[i_min]-X[k_min])+(Y[i_min]-Y[k_min])*(Y[i_min]-Y[k_min]));
 
      if(a+b+c<min)
                 {min=a+b+c;
                  i_min=i;
                  j_min=j;
                  k_min=k;}
 
cout <<a+b+c<<endl;
getch();
}