#include<fstream.h>
#include<iomanip.h>// setw
#include<iostream.h>
#include<stdio.h>
#include<stdlib.h>//randomize();
#include<math.h>

void main(void)
{
int x[20],y[20],Z0[20],Z[20],Zmin[20];
float W[20][20];

ifstream ff("mami.txt");//Считываем координаты городов
int k=0;//Число городов
while(!ff.eof())
    {
    ff>>x[++k];
    ff>>y[k];
    };

int i,j;
for (i=1;i<k;i++)
cout<<i<<": "<<x[i]<<", "<<y[i]<<endl;

for (i=1;i<k;i++)
 {
 for (j=1;j<k;j++)
 W[i][j]=sqrt(pow(x[i]-x[j],2)+pow(y[i]-y[j],2));
 cout<<endl;
 }
for (j=1;j<k;j++) W[j][j]=999;//Бесконечное расстояние
for (j=1;j<k;j++) Z0[j]=j;//Первая последовательность-случайная

float L0=999,T=100,alpha=0.9,L1,P,P0;
int z1,i1,i2;

randomize();
for (int m=1;m<400;m++)//Основной цикл
{
for (j=1;j<k;j++) Z[j]=Z0[j];
   i1=random(k-1)+1; i2=random(k-1)+1;
   z1=Z[i1]; Z[i1]=Z[i2]; Z[i2]=z1;//Меняем местами 2 случайных города
   L1=W[Z[k-1]][Z[1]];
   for (i=1;i<k-1;i++) L1+=W[Z[i]][Z[i+1]];//Длина гамильтонова цикла
   if (L1<=L0){L0=L1;//Если  длина маршрута уменьшилась  - выбираем этот маршрут
   for (j=1;j<k;j++)Z0[j]=Z[j];
   for (j=1;j<k;j++)Zmin[j]=Z[j];//Запоминаем лучший маршрут
   };
   if (L1>L0){//Если больше, то в зависимости от Р, выбираем или нет
       P0=random(100)/100.;
       P=exp(-(L1-L0)/T);//Вероятность выбора худшего решения
       if (P0<P) for (j=1;j<k;j++)Z0[j]=Z[j];
   };
T=alpha*T;//Снижение температуры
}
cout<<"Lmin="<<L0<<endl;
cout<<"T="<<T<<endl;
for (j=1;j<k;j++) cout<<setw(8)<<Zmin[j];
cout<<endl;
getchar();
ff.close();
}