在网上找了一圈模拟退火解决01背包问题的算法,结果只有一个matlab代码,于是自己改成了golang代码
由于本人学golang不到一年,此代码也只是使用了少量的样例,对于原matlab代码中的操作和golang中的操作可能有点不同(比如一个数组直接等于另一个数组)可能导致其实现方法只是凑巧命中,欢迎各位大佬改良,以下是改的代码:
func SA(d, k []int, restriction int) {
//温度衰减速度,马氏链长度,物品数量
tem, mar_length, num := 0.95, 1000, len(k)
sol_new := randnum(num)
sol_current, sol_best := make([]int, num), make([]int, num)
sol_best[0] = 0
E_current, E_best := math.MaxInt64, math.MaxInt64
tf, t := 3.0, 97.0
for t > tf {
for i := 1; i < mar_length; i++ {
//产生随机扰动
temp1 := rand.Intn(num)
sol_new[temp1] = sol_new[temp1] ^ 1
for {
s := dot(sol_new, d, 1) > restriction
if s {
//如果不满足约束条件随机放弃一个物品
find(sol_new[:], 1)
} else {
break
}
}
//计算背包中物品的价值模拟退火算法只能求最小值,所以价值取负号
E_new := dot(sol_new, k, 0)
if E_new < E_current {
E_current = E_new
sol_current = sol_new
if E_new < E_best {
E_best = E_new
sol_best = sol_new
}
} else {
if rand.Float64() < math.Exp(float64(-(E_new-E_best)/int(t))) {
E_current = E_new
sol_current = sol_new
} else {
sol_new = sol_current
}
}
}
t = t * tem
}
fmt.Println(-E_best)
// return -E_best
}
// */
// 随机生成一个固定长度只有0和1的数组
func randnum(n int) []int {
num := make([]int, n)
for key, _ := range num {
num[key] = rand.Intn(2)
}
return num
}
//模拟点乘
func dot(n, m []int, flag int) int {
dot := 0
for i := 0; i < len(n); i++ {
dot += n[i] * m[i]
}
if flag == 1 {
return dot
} else {
return -dot
}
}
//用来随机丢弃一个物品
func find(n []int, m int) {
for {
i := rand.Intn(len(n))
if n[i] == m {
n[i] = m ^ 1
return
}
}
}
测试代码
w, v := make([]int, 3), make([]int, 3)
w[0], w[1], w[2] = 10, 8, 5
v[0], v[1], v[2] = 5, 4, 1
SA(w, v, 16)
原代码地址:
https://www.cnblogs.com/haimishasha/p/9785567.html
(第一次发文章若有格式错误请轻喷)