WIKIOI 1344 线型网络 题解与分析

【题目链接】：

           http://www.wikioi.com/problem/1344/

【分析】：

        随机化搜索，不加随机只能过6个点。首先，预处理若ij毗邻的距离和，然后设定循环执行次数<即人品够好设定一次循环就能过......> ，每次随机一个编号数组node，然后在node角标i，j满足i<j中任意数对判断：<令x1=node[i],y1=node[j]>如果x1前面那个数和y1毗邻,以及y1后面那个数与x1毗邻更优,就可以将x1-y1区间的数以对称轴对称调转区间内的所有数，并设定标记，表示进行过这一步，若没有执行过这一步，说明现在的排序方法为本序列最优，统计答案，最后与标准最优答案比较输出最小即可

 

【代码】：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<time.h>
using namespace std;
#define MAX 23
#define IMAX 21474836
#define MAXSTEP 5//卡出的 
struct POINT{double x,y;};
POINT a[MAX];
int N,node[MAX];
double ans=IMAX,dist[MAX][MAX],ans1;
double DIST(double x1,double y1,double x2,double y2){return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));}
void swap(int &x1,int &y1){int t=x1;x1=y1,y1=t;}
bool check(int x1,int y1)
{
      if(dist[node[x1-1]][node[x1]]+dist[node[y1]][node[y1+1]]>dist[node[x1-1]][node[y1]]+dist[node[x1]][node[y1+1]])//如果x1前面那个数和y1毗邻,以及y1后面那个数与x1毗邻更优,就可以将x1-y1区间的数以对称轴对称调转区间内的所有数 
            return true;
      return false;
}
void change(int x1,int y1)//调转区间数 
{
      while(x1<y1)
      {
            swap(node[x1],node[y1]);
            x1++;
            y1--;
      }
}
void work()
{
      bool sign=true;
      for(int i=1;i<=N;i++)
            node[i]=i;
      for(int i=1;i<=N;i++)//随机序列 
      {
            int num=rand()%(N-i+1)+i;
            swap(node[i],node[num]);
      }
      while(sign)
      {
            sign=false;
            for(int i=1;i<=N;i++)
                  for(int j=i+1;j<=N;j++)
                  {
                        if(check(i,j))
                        {
                              change(i,j); 
                              sign=true;
                        } 
                  }      
      }
}
int main()
{
	  scanf("%d",&N);
      for(int i=1;i<=N;i++)
            scanf("%lf%lf",&a[i].x,&a[i].y);
      srand(time(NULL));
      for(int i=1;i<=N;i++)//预处理ij毗邻的距离 
            for(int j=1;j<=N;j++)
                  dist[i][j]=DIST(a[i].x,a[i].y,a[j].x,a[j].y);
                  
      for(int i=1;i<=MAXSTEP;i++)//限定随机次数 
      {
            ans1=0;
            work();
            for(int j=1;j<N;j++)//统计答案 
                  ans1+=dist[node[j]][node[j+1]];
            ans=min(ans,ans1);
      }
      printf("%.2lf\n",ans);
      return 0;
}

 

转载注明出处： http://blog.csdn.net/u011400953