#include <iostream>
#include <math.h>

using namespace std;
 
int main()
{
    setlocale(LC_ALL, "Russian");
double x[5]={0.5,1,1.5,10,12.5}, y[5]={1.5,0.5,1,7,9};

cout<<"          X    Y"<<endl;
for(int i=0;i<5;i++)
{
cout<<"Объект"<<i<<"  "<<x[i]<<"    "<<y[i]<<endl;
}

double maxX,maxY;
double minX=maxX=x[0], minY=maxY=y[0];
for (int i=1;i<5;i++)
{
   if (minX>x[i])//минимум х
      minX=x[i];
   if (maxX<x[i])//максимум х
      maxX=x[i];
   if (minY>y[i])//минимум y
      minY=y[i];
   if (maxY<y[i])//максимум y
      maxY=y[i];
}

double R0,R1,k;
cout<<"Введите множитель(от 0 до 1) k=";
cin>>k;
R0=pow((maxX-minX)*(maxX-minX)+(maxY-minY)*(maxY-minY),0.5)/2;//общий радиус
R1=R0*k;//радиус группы вместо 0.3 можно использовать переменную k запрет ввода чисел больше 1 и отрицательных
cout<<"Нахождение общего радиуса R0"<<endl;
cout<<"R0="<<R0<<endl;
cout<<"Радиуса для классов R1=R0*k"<<endl;
cout<<"R1="<<R1<<endl;
cout<<endl;
cout<<"Нахождене длины L от Объект2 ("<<x[1]<<":"<<y[1]<<") обозначим номер точки Ti";
//Рандомный выбор точки из массива
double d[5];
for (int i=0;i<5;i++)
{
d[i]=abs(pow((x[1]-x[i])*(x[1]-x[i])+(y[1]-y[i])*(y[1]-y[i]),0.5));
}
cout<<endl;


for(int i=0;i<5;i++)
{
cout<<"L(Т"<<i<<")="<<d[i]<<endl;
}
cout<<endl;


int x1[5];
for(int i=0;i<5;i++)
{
if (R1>=d[i]) x1[i]=1,cout<<"Класс1: "<<x[i]<<" : "<<y[i]<<endl;
else if (R1<=d[i]) x1[i]=0;
}



cout<<"Находим центр тяжести для найденного класса:"<<endl;
//Среднее арифметическое X и Y и получем координаты центра тяжести
int z;
z=0;
for(int i=0;i<5;i++)
{
   z+=x1[i];
}
double Cx=0,Cy=0;
for(int i=0;i<5;i++)
{
Cx+=x[i]*x1[i];
Cy+=y[i]*x1[i];
}
Cx=Cx/z;
Cy=Cy/z;
cout<<"Центр тяжести(координаты)= "<<Cx<<":"<<Cy<<endl;
double d1[4];
d1[0]=abs(pow((Cx-x[0])*(Cx-x[0])+(Cy-y[0])*(Cy-y[0]),0.5));
d1[1]=abs(pow((Cx-x[2])*(Cx-x[2])+(Cy-y[2])*(Cy-y[2]),0.5));
d1[2]=abs(pow((Cx-x[3])*(Cx-x[3])+(Cy-y[3])*(Cy-y[3]),0.5));
d1[3]=abs(pow((Cx-x[4])*(Cx-x[4])+(Cy-y[4])*(Cy-y[4]),0.5));
cout<<endl;
system("pause");
return 0;
}

procedure TForm1.Button4Click(Sender: TObject);
     var i,j,d,n,m,f,l: integer;
         r: array of array of real;
         s,e,min, max: real;
         tk,p: array [1..10] of real;
     begin
       n:= StrToInt(Edit1.Text);
       m:= StrToInt(Edit2.Text);
       e:= StrToFloat(Edit3.Text);
       SetLength(r,n);
       for i:=0 to n-1 do
       SetLength(r[i],m);
       for i:=0 to n-1 do begin
       for j:=0 to m-1 do begin
       r[i][j]:= StrToFloat(StringGrid3.Cells[j+1,i+1]);
       end;end;
       min:=500;
       for i:=0 to n-1 do begin
         max:=-100;
           for j:=0 to m-1 do begin
             if(r[i][j] > max) then max:=r[i][j];
           end;
         if(max < min) then min:=max;
       end;
       s:=min-e;
       for i:= 0 to n-1 do tk[i]:=0;
       f:=1;
       repeat
         max:=-1;
         for i:=0 to n-1 do begin
           p[i]:=0;
           for j:=0 to m-1 do begin
             if (r[i,j] < s) and(tk[j]=0) and(tk[i]=0)then p[i]:=p[i]+1;
           end;
           if (p[i]> max) then begin
             max:=p[i];
             l:=i;
           end;
         end;
         tk[l]:=f;
         for i:=0 to n-1 do if (r[l][i]<s) then tk[i]:=f;
         f:=f+1;
         d:=0;
         for i:=0 to n-1 do if (tk[i]=0) then d:=1;
       until  d=0;
       for i:=0 to f-1 do begin
         Memo1.Lines.Add('в таксон с центром '+IntToStr(i)+' вошли:');
           for j:=0 to n-1 do begin
             if (tk[j]=i) then Memo1.Lines.Add(FloatToStr(tk[j]));;
           end;
       end;
     end;