模拟退火解决01背包问题-GO语言实现

在网上找了一圈模拟退火解决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
(第一次发文章若有格式错误请轻喷)

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值