求模拟退火算法解旅行商问题的C++代码

高分求模拟退火法解 TSP问题 C++或C语言程序~

tsproblem.h
class tsproblem
{
private:
int *parrayc;
int *parrayj;
int icost;
int icurpoint, isize;
unsigned long ls;
unsigned long power(int i, int j);
int g(int i, unsigned long s);
public:
tsproblem(int size, int point);
tsproblem(int size, int point, int *array);
void calculate();
void display();
~tsproblem();
};

#include "tsproblem.h"
#include
#include
#include
//采用动态规划方法,用c(i, j)表示边v(i, j)的长度,g(i, s)表
//示从结点i开始,经过s中的每一结点一次且仅一次到达i的最短路
//径长度,则有如下递推式:
// g(i, s)=min{c(i, j)+g(j, s-{j})} (j在s集合中)
//为实现集合s的表示,采用如下方法:
//把s定义为unsigned long型数,设其二进制表示为:bn....b2b1b0
//则当b0为1时表示0在集合s中,b0为0时表示0不在集合s中,余类推
tsproblem::tsproblem(int size, int point)
{
int i, j;
isize=size;
icurpoint=point;
if(icurpoint>=isize){
cout<<"the curpoint is out of range!
";
return;
}
parrayc=new int[isize*isize];
for(i=0; i<isize; i++)
for(j=0; j<isize; j++){
if(i==j)
parrayc[i*isize+j]=0;
else
parrayc[i*isize+j]=rand()*isize/18000+1;
}
parrayj=new int[isize*power(2, isize)];
}

void tsproblem::calculate()
{
ls=0;
ls=~ls;
ls=ls<<isize;
ls=~ls;
ls-=power(2, icurpoint);
icost=g(icurpoint, ls);
}

int tsproblem::g(int i, unsigned long s)
{
int j, k, t, temp=32767;
if(s==0)
return parrayc[i*isize+icurpoint];
for(j=0; j<isize; j++){
if((s>>j)%2==0) continue;
else{
t=parrayc[i*isize+j]+g(j, s-power(2, j));
if(temp>t){
temp=t;
k=j;
}
}
}
parrayj[i*power(2, isize)+s]=k;
return temp;
}

unsigned long tsproblem::power(int i, int j)
{
unsigned long t=1;
for(int k=0; k<j; k++)
t=t*i;
return t;
}

void tsproblem::display()
{
int temp, i, j, n=power(2, isize);
unsigned long s=ls;
temp=icurpoint;
cout<<" the array of cost is:
";
for(i=0; i<isize; i++){
for(j=0; j<isize; j++)
cout<<setw(4)<<parrayc[i*isize+j];
cout<<endl<<"";
}
cout<<"
the shortest path is:
";
cout<<" "<<icurpoint;
while(s){
temp=parrayj[temp*n+s];
s-=power(2, temp);
cout<<" "<<temp;
}
cout<<" "<<icurpoint<<endl<<" the min cost is:"
<<icost<<endl;
}

tsproblem::~tsproblem()
{
if(parrayc)delete parrayc;
if(parrayj)delete parrayj;
}

那个论文你怎么做的？

我自己写的：
#include<stdlib.h>
#include<cmath>
#include<algorithm>
using namespace std;
static double Tmax = 10, Tmin = 0.1, r = 0.999999;
static int k = 100;
inline void random_sele(int &fl, int &fp, int arr[], int n)
{
do
{
fl = rand() % n;
fp = rand() % n;
} while (fl == arr[fp] || arr[fl] == arr[fp]);
};
inline bool accept(double r1, double r2, double T)
{
const unsigned int MASK((1 << 30) - 1);
double p = 1.0 / (1.0 + exp((r2 - r1) / T));
unsigned int temp1=((rand()<<20)^(rand()<<10)^rand())&MASK;
unsigned int temp2=((rand()<<20)^(rand()<<10)^rand())&MASK;
double res=(1.0*temp1+1.0*temp2/MASK);
return res < p*MASK;
}
void set_SA()
{
printf(" Tmax Tmin r k\n");
printf(" %.8f %.8f %.8f %d \n",Tmax,Tmin,r,k);
printf("Input Tmax,Tmin,r,k:\n");
scanf("%lf%lf%lf%d",&Tmax,&Tmin,&r,&k);
}
void TSP_SA(int arr[], double &evl, int n, double map[][2000])
{
double r1, r2, T = Tmax;
int i, fl, sl, fp, sp;
int show_s=0;
srand(11827);
while (T >= Tmin)
{
for (i = 0; i < k; i++)
{
random_sele(fl, fp, arr, n);
sl = arr[fl], sp = arr[fp];
r1 = map[fl][sl] + map[fp][sp] + map[sp][arr[sp]];
r2 = map[fl][sp] + map[sp][sl] + map[fp][arr[sp]];
if (accept(r1, r2, T))
{
arr[fp] = arr[sp];
arr[fl] = sp;
arr[sp] = sl;
sl = sp;
evl = evl + r2 - r1;
}
}
if(++show_s==1000000/k)
{
printf("T=%f evl=%f\n",T,evl);
show_s=0;
}
T *= r;
}
}

---

高陵县15860562205： 求模拟退火算法解旅行商问题的C++代码？
成王眨血塞： 我自己写的: #include #include #include using namespace std; static double Tmax = 10, Tmin = 0.1, r = 0.999999; static int k = 100; inline void random_sele(int &fl, int &fp, int arr[], int n) { do { fl = rand() % n; fp = rand() % n; } while (fl == arr[fp] || arr[...

高陵县15860562205： 谁能帮我翻译一下我的毕业论文的题目:基于模拟退火算法解决旅行商问题 - ？
成王眨血塞： Based on simulated annealing algorithm to solve traveling salesman problem

高陵县15860562205： 高分求模拟退火法解 TSP问题 C++或C语言程序 - ？
成王眨血塞： tsproblem.hclass tsproblem{private:int *parrayc;int *parrayj;int icost;int icurpoint, isize;unsigned long ls;unsigned long power(int i, int j);int g(int i, unsigned long s);public:tsproblem(int size, int point);tsproblem(int size, int point, int *...

高陵县15860562205： 模拟退火算法 Pascal&C++？
成王眨血塞：模拟退火的基本思想: (1) 初始化:初始温度T(充分大),初始解状态S(是算法迭代的起点), 每个T值的迭代次数L (2) 对k=1,……,L做第(3)至第6步: (3) 产生新解S′ (4) 计算增量Δt′=C(S′)-C(S),其中C(S)为评价函数 (5) 若Δt′<0则接受S′作为新的当前解,否则以概率exp(-Δt′/T)接受S′作为新的当前解. (6) 如果满足终止条件则输出当前解作为最优解,结束程序. 终止条件通常取为连续若干个新解都没有被接受时终止算法. (7) T逐渐减少,且T->0,然后转第2步.

高陵县15860562205： 急求,模拟退火遗传算法的MATLAB程序!谢谢？
成王眨血塞： 你真幸福.我刚刚编了一个模拟退火算法,计算旅行商问题:注意:一共三个文件,第一个是主程序,下面两个是子函数.% for d=1:50 %循环10次发现最小路径为4.115,循环50次有3次出现4.115 T_max=80; %input('please input the start ...

高陵县15860562205： 蛮力法求解旅行商问题 求C++代码 - ？
成王眨血塞： #include<stdlib.h> #include<cmath> #include<algorithm> using namespace std; static double Tmax = 10, Tmin = 0.1, r = 0.999999; static int k = 100; inline void random_sele(int &fl, int &fp, int arr[], int n) { do { fl = rand() % n; fp = rand() % n; } while ...

高陵县15860562205： C语言实现红外图像自适应对比度增强？
成王眨血塞： c/c++的模拟退火算法程序 Tc下编译. 密码是:asd http://ds8.fileflyer.com/d/b02df8ec-82cf-48da-8fa6-9e1c28eb7ecd/vn1U/Lna5OAn/anneal.rar 可以打开,没问题 ,密码是解压缩密码 给你发了看下邮箱吧不会吧,还有邮箱吗?

高陵县15860562205： 模拟退火算法解决旅行商问题部分论文请哪位日语达人帮忙翻译一下？
成王眨血塞： 首先介绍了旅行商问题(最初に旅行行商人を质问导入した) 模拟退火算法原理及其算法实现(シミュレーションのアニーリングのアルゴリズムの主义およびアルゴリズムの认识) 应用模拟退火算法对TSP进行研究(TSPにシミュレーショ...

高陵县15860562205： 求货郎担问题的matlab算法 - ？
成王眨血塞： 货郎担问题有很多解法,模拟退火,遗传算法,动态规划等.基于matlab TSP问题遗传算法的实现%TSP问题(又名:旅行商问题,货郎担问题)遗传算法通用matlab程序%D是距离矩阵,n为种群个数,建议取为城市个数的1~2倍,%C为停止...

高陵县15860562205： 求用C++编程的优化设计中坐标轮换法的程序 - ？
成王眨血塞： 要是有点c基础的话可以看看c++ primer 看看前几章的容器和泛型算法 用这个实现遗传算法解决旅行商问题 代码100行左右