d_1.r

这是07年的大学生建模比赛最为完整的解答 使用了lingo 和R软件

### 读数据
rc<-read.csv("D_data.csv")

### 将数据写成矩阵
str<-1:16
source("read_matrix.R")
dist<-read_matrix(data=rc, str=str)

### Floyd 算法求任意两点的最短路
source("floyd.R")
dist<-floyd(dist)
dist

#### 重新开始计算 ####################################
D<-dist[c(1, str+1), c(1, str+1)]
D

### 写文件, 为LINGO运算使用
source("write_file.R")
write_file(str, D)

### 此时用LINGO软件计算最短路
 
#####################################################
ru<-c(10,15,6,9,13,6,11,4,13,17,11,2,11,21,13,14); sum(ru)
chu<-c(9,14,5,10,9,10,13,9,15,9,6,7,13,15,10,16); sum(chu)
dif<-chu-ru; dif

#####################################################
### 计算小圈的TSP
str<-2:7             ### 转重新开始计算
str<-c(3,2,6,7,5,4)
str<-c(4,5,7,6,2,3)

str<-c(8:10, 15, 16) ### 转重新开始计算
str<-c(10,9,8,16,15)

str<-c(1, 11:14)     ### 转重新开始计算
str<-c(13,1, 14, 12, 11)
# str<-c(11,12,14,1,13)

sum(ru[str])
sum(chu[str])

#####################################################
### 计算小圈的TSP
str<-c(4,5,2,3,1,13)

str<-6:10
### str<-10:6

str<-c(14, 15, 16, 11, 12)

#####################################################
### 计算车辆载货数
w0<-sum(ru[str]);w0
n<-length(str); w<-numeric(0)
w[1]<-w0+dif[str[1]]
for (i in 2:n)
  w[i]<-w[i-1]+dif[str[i]]
w

### 计算空车率
Wt<-65
K0<-(Wt-w0)/Wt;K0
K<-numeric(0)
for (i in 1:n)
   K[i]<-(Wt-w[i])/Wt
K

### 计算空车损失
D<-dist[c(1, str+1), c(1, str+1)]; D
L<-K0*D[1,2]
for (i in 1:(n-1))
   L<-L+K[i]*D[i+1, i+2]
L<-L+K[n]*D[(n+1),1]
L*2

###########################################
### 空车率矩阵
Km<-dist* c(8, dif+8);Km
D<-Km[c(1, str+1), c(1, str+1)]
D

### 写文件, 为LINGO运算使用
source("write_file.R")
write_file(str, D)

### 此时用LINGO软件D1_TSP计算最短路


####################################
### 计算结果
str<-c(4,5,2,3,1,13);    19.32308
str<-6:10;               28.89231
str<-c(14,16,15, 11,12); 33.75385

####################################
### 计算结果
str<-c(3,2,7,6,5,4)
str<-c(12,11,14,1,13)

str<-c(1,2,3,12,13,14)
str<-c(12,14,13,1,2,3)