线型网络

题目描述 Description
有 N ( <=20 ) 台 PC 放在机房内，现在要求由你选定一台 PC，用共 N-1 条网线从这台机器开始一台接一台地依次连接他们，最后接到哪个以及连接的顺序也是由你选定的，为了节省材料，网线都拉直。求最少需要一次性购买多长的网线。（说白了，就是找出 N 的一个排列 P1 P2 P3 ..PN 然后 P1 -> P2 -> P3 -> … -> PN 找出 |P1P2|+|P2P3|+…+|PN-1PN| 长度的最小值)

输入描述 Input Description
第一行 N ，下面 N 行，每行分别为机器的坐标(x,y) ( 实数 -100<=x,y<=100 )

输出描述 Output Description
最小的长度，保留两位小数。

样例输入 Sample Input
3
0 0
1 1
1 -1

样例输出 Sample Output
2.83

这道题目确定的算法，非常容易超时，普通的搜索要经过很好的优化才可能
过。所以，我们采取随机化的算法。
题目中告诉你的，是各台机器的坐标。所以我们第一步工作，就是把每两台机器之间的距离都计算出来，方法是勾股定理（不会不知道吧）。
然后，我们可以先随机生成一条当前的最优路径。当然这条路径就是最优路径的情况是少见的，我们必须要对此进行修补。然后采用交换的方法试，具体方法是用Bubble Sort的方法，去枚举每两个点，如果第a-1个点到第a个点的距离与第b+1个点到第b个点的距离之和大于第a-1个点到第b个点的距离与第b+1个点到第a个点的距离之和，则把这一段路的所有节点用反向。（为什么？）然后就这样对一条路径不断优化，直到不能优化。

然而一次生成修补的路径并不能代表最优，我们还需要把上述过程重复一定的次数，以达到最优。

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
double x[22],y[22];

double f[22][22];
int n;
int index[22];
double dist(double x1,double y1,double x2,double y2){
       return sqrt( (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2) );
}
void init(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%lf%lf",&x[i],&y[i]);
    for(int i=1;i<=n;i++)
      for(int j=1;j<=n;j++)
        f[i][j]=dist(x[i],y[i],x[j],y[j]);
    for(int i=1;i<=n;i++){
        f[0][i]=0;
        f[i][n+1]=0;
    }
    for(int i=1;i<=n;i++) index[i]=i;
}
void sv(int a,int b)
{
    int i,j;
    for (i=a,j=b;i<=j;i++,j--)
    {
        swap(index[i],index[j]);
    }
}
double solve(){
       for(int i=1;i<=100;i++){
          int a=rand()%(n)+1;
          int b=rand()%(n)+1;
          swap(index[a],index[b]);
       }

       for(int i=1;i<=n-1;i++){
           for(int j=i+1;j<=n;j++)
           if(f[index[ i-1] ][index[i] ]+f[ index[j] ][index[j+1] ]>
              f[index[ i-1] ][index[j] ]+f[ index[i] ][index[j+1] ] ){
                 sv(i,j);
              }
       }
       double t=0;
       for(int i=1;i<=n-1;i++){
          t+=f[index[i] ][ index[i+1]];
       }
       return t;
}

int main(){
    freopen("linec.in","r",stdin);
    freopen("linec.out","w",stdout);
    init();
    double ans=100000000;
    for(int i=1;i<=2000;i++){
        double t=solve();
        if(t<ans) ans=t;
    }
    printf("%.2lf\n",ans);
    return 0;
}