大连海事大学现代优化技术作业节约算法解tsp问题的c语言实现.docxVIP

作业2：节约算法解TSP问题的C语言实现 1120131233陈鑫交通运输工程 问题描述 节约算法(Saving Algorithm)是解决配送问题的一种常见的构筑算法。其 基本思路是： 首先，所有顾客都肓接从配送中心之间的肓接运输开始； 然后，计算上述每两个直接运输如果改成连续运输后的距离(费用)节约额 Sjj = cQi + cOj - Cjj, i, j为任意两个顾客； 最后，按照任意两个顾客i, j之间连续运输的节约额Sy从大到小顺序来组 织连续运输，知道配送车辆等制约条件不能满足为止。 本程序里不涉及任何吋间窗或其他约束，适合20个点以内 算法实现 #include stdio.hn #include math.h void main() { longx[21]={0},y[21]={0},d[191]={0},e[191]={0},g[41]={0},s[42]={0); int i,j,h=l,k,l,o,m,n=0,p,q=3,r=0; int flag; double a[21]={0},b[21][21]={0},c[21][21]={0},f[191]={0}; char ch; printfC* C-W 算法求解 TSP 问题 \n\n\nH); printfC*请输入坐标(20个以内),坐标之间用空格隔开，按回车键结束输入： \nH); re:scanf(H%d,%dH,x[h],y[h]); ch=getchar(); if(ch!=,\n,) { h++; goto re; } for(i=l;ih+l；i++) a[i]=sqrt(x[i]*x[i]+y[i]*y[i]); for(i=l;ih;i++) { for(j=i+1 ;j vh+1 ;j++) { n++; b[i][j]=sqrt((x[i]-x[j])*(x[i]-x[j])4-(y[i]-y[j])*(y[i]-y|j])); c[i]Ul=a[i]+a[j]-b[i]U]； d[n]=i； e[n]=j; f[n]=c[i]U]; } } fbr(i=l;in;i++) { flag=0; for(j=l;jvn+l?i;j++) if(f[j]fU+l]) { k=f[j];f[j]=fU+l];f[j+l]=k; l=d[j];d[j]二 d[j+l];d[j+l]=l; o=e[j];e[j]=e[j+1 ];e[j+l]=o; flag 二 1 ； } if(flag==O) break; } printf(n\n,r); printf(”节约值排序:\nn); for(i=n;i0;i-) { if(f[i]!=O) printf(,%lf(%ld,%ld)--(%ld,%ld)\n,\f[i],x[d[i]],y[d[i]],x[e[i]],y[e[i]]); 1 g[l]=d[n]; g[2]=e[n]; for(m= 1 ;mvh;m++) { for(j=n-l;j0;j-) { P=0; for(i=l;i41;i++) { if(dUl==g[i]) p 二 p+1; if(e[j]==g[i]) P=P+1； if(P=l) { g[ql=d[j]； g[q+l]=e[j]; if(p==l) break; q=q+2; } printf(M\n 选出连接点:\nu); for(i=l ;i41 ;i=i+2) { if(g[i]!=O) { printf(“(％ld,%ld)??(％ld,%ld)\n“,x[g[i]],y[g[i]],x[g[i+l]],y[g[i+l]]); r=r+2; } } s[21]=g[l]; s[22]二g ⑵； for(i=3;ir;i=i+2) for(j=l;j42;j++) { if(g[i]==s[j]) { if(s[j?l]=O) s[j?l]=g[i+l]； else sU+l]=g[i+l]； break; if(s[j?l]=O) s[j-l]=g[i]； Else s[j+l]二 g[i]； break; printf(\n旅彳亍商路线：\nH); printf(,,(O,O)-H); for(i=l;i42;i++) { if(s[i]!=O) printf(,,(%ld,%ld)-,,,x[s[i]],y[s[i]]); } printf(,,(O,O)n); getch(); 运行截图如下: ■ X:\USERS\ADMINISTRATOR\DESiaOP\hh\DebugV^e^ 请谕入坐标（20个以内）,坐标之间用空格隔开.按回车键结束输入: 12.23 26.56 89.45 75,66 45.32 ［节约值排序： 174.00000