【算法设计与分析】基于Go语言实现动态规划法解决TSP问题

news2024/11/26 2:54:03

 本文针对于最近正在学习的Go语言,以及算法课实验所需内容进行Coding,一举两得!     

一、前言 

        由于这个实验不要求向之前的实验一样做到那种连线的可视化,故可以用图形界面不那么好实现的语言进行编写,考虑到Go语言的方兴未艾,所以采用此种语言解决问题。


 二、问题

        TSP问题的大致解法,老师在课上已经说过了,清华大学出版社的《算法设计与分析》(第二版,然而书上伪代码存在一些疏漏)里面也有所阐述,这里不做细致解释。


三、代码分析

        主要可分为三个部分,输入、输出、计算。

1.输入

        输入部分需要一个整型变量存点(城市)的数量,一个矩阵存点到点的距离,另外增设一个矩阵存到某个点的最近的前驱。这里还有一个重要的问题是如何做出这些点的子集,也就是所要画的图(实验结果)的表头横向内容,代码如下:

 ·········
    var nums int
	// 读取二维数组的行数和列数
	fmt.Print("请输入(城市)点数: ")
	fmt.Scanln(&nums)
	// 初始化一个二维数组
	arc := make([][]int, nums)
	for i := range arc {
		arc[i] = make([]int, nums)
	}
	// 从控制台读取二维数组的值
	fmt.Println("请输入二维数组的元素,每行输入完毕后按回车键:")
	for i := 0; i < nums; i++ {
		for j := 0; j < nums; j++ {
			var putin int
			fmt.Scanf("%d", &putin)
			if putin == 0 {
				putin = 2004
			}
			arc[i][j] = putin
		}
		fmt.Scanln() // 跳过每行输入后的换行符
	}
	var LengthOfd int = int(math.Pow(2, float64(nums-1)))
	//下面的这个d就是那个动态规划法的表
	d := make([][]int, nums)
	for i := range d {
		d[i] = make([]int, LengthOfd)
	}
	//同样设置一个跟d一样的矩阵,来存最近的前驱
	front := make([][]int, nums)
	for i := range front {
		front[i] = make([]int, LengthOfd)
	}
	//初始化设置成很大的
	for i := range d {
		for j := range d[i] {
			d[i][j] = 2004
			front[i][j] = 2004
		}
	}
	for i := 0; i < nums; i++ {
		d[i][0] = arc[i][0]
		front[i][0] = 0
	}
	// 创建一维存所有要做子集的点。
	numName := make([]int, nums) //生成1,2,3
	for i := 1; i <= nums; i++ {
		numName[i-1] = i
	}
	numName = numName[:len(numName)-1]
	subset := subsets(numName) //生成子集
	sort.Slice(subset, func(i, j int) bool {
		return len(subset[i]) < len(subset[j])
	})
	fmt.Println(subset)
·······

        这里前面几个变量我不加以赘述,简单的创建和初始化(如果你用的其他语言写这道题,相信你能做到这个),这里说一下最小子集的寻找,我参考了LeetCode上一道题(力扣上的子集问题)以及CSDN上相关博主给出的解答(子集问题的GO语言其中一解,这里选择解题代码并未考虑时间及空间复杂度,大家可以试着采用leetcode上更快更好的代码),代码大意是采用了深度优先搜索的方式进行生成,下面是本题借用函数:

//来自于本站其他用户博文
func subsets(nums []int) [][]int {

    l := list.New()
    result := list.New()
    
    for i := 0; i <= len(nums); i++ {
         dfs(nums, 0, i, l, result)
    }
   
    arr := make([][]int, result.Len())
    k := 0
    for e := result.Front(); e != nil; e = e.Next(){
        curl := e.Value.(*list.List)
        arr[k] = make([]int, curl.Len())
        k++
    }

    i := 0;
    for e := result.Front(); e != nil; e = e.Next() {
        curl := e.Value.(*list.List)
        j := 0
        for p := curl.Front(); p != nil; p = p.Next() {
           arr[i][j] = p.Value.(int)
           j++
        }
        
        i++;
    }

    return arr
}

func dfs(nums []int, start int, len int, l *list.List, result *list.List) {

    if start == len {
        a := list.New()
        for e := l.Front(); e != nil; e = e.Next() {
            a.PushBack(e.Value)
        }
        result.PushBack(a)
        return
    }

    for i := start; i < len; i++ {
        l.PushBack(nums[i])
        dfs(nums, i+1, len, l, result)
        b := l.Back()
        l.Remove(b)
    }
}

         之后就是对得到的子集排序,因为上面代码跑出来的子集并非是有序的。

        这样便得到了计算所需的所有变量。

2.计算 

        计算这个方面可以参考书中的伪代码,值得注意的是在判断大小的一些地方需要改变传参。大家可通过书中样例矩阵以及表格进行相关推导,不难发现第0列被初始化,第1到3列依靠0列更新,3到5列依靠1到3列更新,第七列特别只有第0行需要更新,而我们的最终结果需要的表格的横向表头的二维数组长这样:

        所以,我们不难发现,它的长度似乎和计算有一些联系:在遍历{1}的时候,填充2和3,依靠0到1的数值,所以计算d[2][1] = arc[2][1] + d[1][0],d[3][1] = arc[3][1] + d[1][0]。依照此种规律,便可结合下文的代码内容分析:

	//下面才刚刚开始tsp
	for j := 1; j < len(subset); j++ { //以d中的列开始扫描,也就是上面的subset(V),表头内容
		for i := 1; i < nums; i++ { //挨个查找看看那个数字没在V的里面,真没在就去赋值,正在查找横表头有无
			judge := false
			for _, theNum := range subset[j] {
				if theNum == i {
					judge = true //在的在的,就不用操作了
					continue
				}
			}
			if judge == false { //wok,真没在

				var edge int = 1314

				for _, theNum := range subset[j] { //theNum,表头中含有的内容
					temp := arc[i][theNum] + d[theNum][j-theNum-len(subset[j])+1]
					if temp < edge {
						edge = temp
						d[i][j] = edge
						front[i][j] = theNum
					}
					//d[i][j] = int(math.Min(float64(d[i][j]), float64(arc[i][theNum]+d[theNum][j-1])))
				}
			}
		}
	}
	//求解最终结果
	TheLastEdge := 520
	for k := 1; k < nums; k++ {
		temp := arc[0][k] + d[k][LengthOfd-1-int(math.Pow(2, float64(k-1)))]
		if temp < TheLastEdge {
			TheLastEdge = temp
			d[0][LengthOfd-1] = TheLastEdge
			front[0][LengthOfd-1] = k
		}
	}

         最后求解最后的那个框,例如本题目就是{1,2,3}下面的第0行的内容。和上面写到的,平常的点求该问题的解法如出一辙,这里并不是很好理解,计算公式甚至可以理解成是我在纯粹的找规律凑数。值得注意的是,在计算的同时front矩阵也在记录他们的上一个点(前驱),这个在后面输出发挥着重要作用。

3.输出

        这里要注意以下路径是怎么打印出来的。i是直接指向最后更新的那个框框的纵坐标,j是横坐标,count用来计数确保不会在移动的途中迷路,根据书上表格,我们并不难发现下面的规律,j直接取front中的值会直接得到前驱节点的值,那么我们就应该在j行去定位这个前驱的下一个前驱,那么我们就只差寻找这个地方的纵坐标,相当巧合的是,现在的i减去现在的j的值再减去这次遍历计数器的值,正好到了我们要寻找的那个地方。(这里在草稿纸上列出我们需要加起来的点,恰好可以得出普适的规律)。

	fmt.Print("TSP最短的路径是:", "0")
	//打印路径
	for i, j, count := LengthOfd-1, 0, 0; ; {
		j = front[j][i]
		i = i - j - count
		fmt.Print("->", j)
		count++
		if i <= 0 {
			fmt.Println("->0")
			break
		}
	}
	//画表
	for i := 0; i < len(subset); i++ {
		str := fmt.Sprint(subset[i])
		fmt.Printf("%10s", str)
		//fmt.Print(subset[i], "\t\t")
	}
	fmt.Println()
	for i := 0; i < nums; i++ {
		for j := 0; j < LengthOfd; j++ {
			bye := "-"
			if d[i][j] == 2004 {
				fmt.Printf("%10s", bye)
			} else {
				fmt.Printf("%10d", d[i][j])
			}

		}
		fmt.Println()
	}

	fmt.Println("最短路径是:", d[0][LengthOfd-1])

        大致就是这样得到了最后结果,然后再对齐一下表格。


四、代码

// TSP Problem
// Created By DDD on 2024.5.25
//

package main

import (
	"container/list"
	"fmt"
	"math"
	"sort"
)

func subsets(nums []int) [][]int {

	l := list.New()
	result := list.New()

	for i := 0; i <= len(nums); i++ {
		dfs(nums, 0, i, l, result)
	}

	arr := make([][]int, result.Len())
	k := 0
	for e := result.Front(); e != nil; e = e.Next() {
		curl := e.Value.(*list.List)
		arr[k] = make([]int, curl.Len())
		k++
	}

	i := 0
	for e := result.Front(); e != nil; e = e.Next() {
		curl := e.Value.(*list.List)
		j := 0
		for p := curl.Front(); p != nil; p = p.Next() {
			arr[i][j] = p.Value.(int)
			j++
		}

		i++
	}

	return arr
}

func dfs(nums []int, start int, len int, l *list.List, result *list.List) {

	if start == len {
		a := list.New()
		for e := l.Front(); e != nil; e = e.Next() {
			a.PushBack(e.Value)
		}
		result.PushBack(a)
		return
	}

	for i := start; i < len; i++ {
		l.PushBack(nums[i])
		dfs(nums, i+1, len, l, result)
		b := l.Back()
		l.Remove(b)
	}
}

func main() {
	var nums int
	// 读取二维数组的行数和列数
	fmt.Print("请输入(城市)点数: ")
	fmt.Scanln(&nums)
	// 初始化一个二维数组
	arc := make([][]int, nums)
	for i := range arc {
		arc[i] = make([]int, nums)
	}
	// 从控制台读取二维数组的值
	fmt.Println("请输入二维数组的元素,每行输入完毕后按回车键:")
	for i := 0; i < nums; i++ {
		for j := 0; j < nums; j++ {
			var putin int
			fmt.Scanf("%d", &putin)
			if putin == 0 {
				putin = 2004
			}
			arc[i][j] = putin
		}
		fmt.Scanln() // 跳过每行输入后的换行符
	}
	var LengthOfd int = int(math.Pow(2, float64(nums-1)))
	//下面的这个d就是那个动态规划法的表
	d := make([][]int, nums)
	for i := range d {
		d[i] = make([]int, LengthOfd)
	}
	//同样设置一个跟d一样的矩阵,来存最近的前驱
	front := make([][]int, nums)
	for i := range front {
		front[i] = make([]int, LengthOfd)
	}
	//初始化设置成很大的
	for i := range d {
		for j := range d[i] {
			d[i][j] = 2004
			front[i][j] = 2004
		}
	}
	for i := 0; i < nums; i++ {
		d[i][0] = arc[i][0]
		front[i][0] = 0
	}
	// 创建一维存所有要做子集的点。
	numName := make([]int, nums) //生成1,2,3
	for i := 1; i <= nums; i++ {
		numName[i-1] = i
	}
	numName = numName[:len(numName)-1]
	subset := subsets(numName) //生成子集
	sort.Slice(subset, func(i, j int) bool {
		return len(subset[i]) < len(subset[j])
	})
	fmt.Println(subset)
	//下面才刚刚开始tsp
	for j := 1; j < len(subset); j++ { //以d中的列开始扫描,也就是上面的subset(V),表头内容
		for i := 1; i < nums; i++ { //挨个查找看看那个数字没在V的里面,真没在就去赋值,正在查找横表头有无
			judge := false
			for _, theNum := range subset[j] {
				if theNum == i {
					judge = true //在的在的,就不用操作了
					continue
				}
			}
			if judge == false { //wok,真没在

				var edge int = 1314

				for _, theNum := range subset[j] { //theNum,表头中含有的内容
					temp := arc[i][theNum] + d[theNum][j-theNum-len(subset[j])+1]
					if temp < edge {
						edge = temp
						d[i][j] = edge
						front[i][j] = theNum
					}
					//d[i][j] = int(math.Min(float64(d[i][j]), float64(arc[i][theNum]+d[theNum][j-1])))
				}
			}
		}
	}
	//求解最终结果
	TheLastEdge := 520
	for k := 1; k < nums; k++ {
		temp := arc[0][k] + d[k][LengthOfd-1-int(math.Pow(2, float64(k-1)))]
		if temp < TheLastEdge {
			TheLastEdge = temp
			d[0][LengthOfd-1] = TheLastEdge
			front[0][LengthOfd-1] = k
		}
	}

	fmt.Print("TSP最短的路径是:", "0")
	//打印路径
	for i, j, count := LengthOfd-1, 0, 0; ; {
		j = front[j][i]
		i = i - j - count
		fmt.Print("->", j)
		count++
		if i <= 0 {
			fmt.Println("->0")
			break
		}
	}
	//画表
	for i := 0; i < len(subset); i++ {
		str := fmt.Sprint(subset[i])
		fmt.Printf("%10s", str)
		//fmt.Print(subset[i], "\t\t")
	}
	fmt.Println()
	for i := 0; i < nums; i++ {
		for j := 0; j < LengthOfd; j++ {
			bye := "-"
			if d[i][j] == 2004 {
				fmt.Printf("%10s", bye)
			} else {
				fmt.Printf("%10d", d[i][j])
			}

		}
		fmt.Println()
	}

	fmt.Println("最短路径是:", d[0][LengthOfd-1])
}

/*
4
0 3 6 7
5 0 2 3
6 4 0 2
3 7 5 0
*/

五、参考文献

 算法分析与设计课程实验——TSP问题与01背包问题的动态规划算法实现-CSDN博客


Go的主函数不需要写return

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1693953.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

基于.net开发的博客系统

基于.net开发可以批量上传md文件生成文章的博客系统 .NET 个人博客 基于.net开发的博客系统 个人博客系统&#xff0c;采用.net core微服务技术搭建&#xff0c;采用传统的MVC模式&#xff0c;使用EF core来对mysql数据库(sqlite数据库)进行CRUD操作项目 为什么要自己开发博客…

uniapp微信小程序解决type=“nickname“获取昵称,v-model绑定值为空问题!

解决获取 type"nickname"值为空问题 文章目录 解决获取 type"nickname"值为空问题效果图Demo解决方式通过表单收集内容通过 uni.createSelectorQuery 效果图 开发工具效果图&#xff0c;真机上还会显示键盘输入框 Demo 如果通过 v-model 结合 blur 获取不…

使用梦畅闹钟,结合自定义bat、vbs脚本等实现定时功能

梦畅闹钟-每隔一段时间运行一次程序 休息五分钟bat脚本&#xff08;播放音乐视频&#xff0c;并锁屏&#xff09; chcp 65001 echo 回车开始休息5分钟 pause explorer "https://www.bilibili.com/video/BV1RT411S7Tk/?p47" timeout /t 3 /nobreak rundll32.exe use…

Font shape `U/rsfs/m/n‘ in size <29.86> not available size <24.88>

解决方法&#xff1a;mathrsfs 删除这个包 其他可以参考&#xff1a;koma script - Size substitution with fontsize14 - TeX - LaTeX Stack Exchange

若依框架对于后端返回异常后怎么处理?

1、后端返回自定义异常serviceException 2、触发该异常后返回json数据 因为若依对请求和响应都封装了&#xff0c;所以根据返回值response获取不到Code值但若依提供了一个catch方法用来捕获返回异常的数据 3、处理的方法

python给图片加上图片水印

python给图片加上图片水印 作用效果代码 作用 给图片加上图片水印图片水印的透明度&#xff0c;位置可自定义 效果 原始图片&#xff1a; 水印图片&#xff1a; 添加水印后的图片&#xff1a; 代码 from PIL import Image, ImageDraw, ImageFontdef add_watermark(in…

shelll 正则表达式

sort sort命令对行内容进行排序 sort语法&#xff1a; 1.sort &#xff08;选项&#xff09; 参数 2.cat file | sort 选项 选项&#xff1a; -n 按照数字进行排序 -r 反向排序 -k 指定排序 -f 忽略大小写 会将小写字母转化成大写字母来比较 -b 忽略每行前面的空格 .........…

cake (C# Make) cross-platform build automation system

1、安装Cake.Tool dotnet new tool-manifest dotnet tool install -g Cake.Tool --version 4.0.0 dotnet-cake --help C:\Users\用户名\.config\dotnet-tools.json 2、新建项目CakeDemo(net8.0) 编写脚本&#xff1a; build.cake var target Argument("target",…

BWVS 靶场测试

一、PHP弱类型 is_numeric() 输入&#xff1a;127.0.0.1/BWVS/bug/php/code.php # 1、源代码分析 如果num不是数字&#xff0c;那么就输出num&#xff0c;同时如果num1&#xff0c;就输出flag。即num要是字符串又要是数字 # 2、函数分析&#xff1a; is_numeric()函数&…

怎么查看项目中antd的版本

使用antd时&#xff0c;有在线参考资料&#xff0c;但是需要根据项目需要&#xff0c;选择对应版本的参考资料。 antd在线参考资料&#xff1a; 组件总览 - Ant Design 如何查看当前项目中antd的版本呢&#xff1f; 在项目的终端中输入&#xff1a; npm list antd antd官网选择…

【Linux】详解线程控制之线程创建线程终止线程等待线程分离

一、线程创建 thread&#xff1a;这是一个指向pthread_t类型的指针&#xff0c;用于获取新创建线程的线程ID。在调用pthread_create后&#xff0c;这个指针会被设置为新线程的ID。 attr&#xff1a;这是一个指向pthread_attr_t类型的指针&#xff0c;用于设置线程的属性&#x…

零基础的粉丝有福了:逐键提示盲打更轻松

盲打就是不看键盘去打字&#xff0c;对于零基础的粉丝而言&#xff0c;盲打入门通常都是很难的&#xff0c;今天就给大家放个福利&#xff1a;从今天开始就能盲打&#xff0c;3天之后盲打就入门了。 真的有这么简单吗&#xff1f;是的&#xff0c;跟着我做就可以了。 首先&am…

Leetcode刷题笔记2:数组基础2

导语 leetcode刷题笔记记录&#xff0c;本篇博客记录数组基础1部分的题目&#xff0c;主要题目包括&#xff1a; 977.有序数组的平方 &#xff0c;209.长度最小的子数组 &#xff0c;59.螺旋矩阵II 知识点 滑动窗口 所谓滑动窗口&#xff0c;就是不断的调节子序列的起始位…

sysbench压测mysql性能测试命令和报告

sysbench压测mysql性能测试命令和报告 一、安装sysbench工具二、创建测试数据库三、基于sysbench构造测试表和测试数据四、数据库性能测试1、数据库读写性能测试2、数据库读性能测试3、数据库删除性能测试4、数据库更新索引字段性能测5、数据库更新非索引字段性能测试6、数据库…

[激光原理与应用-93]:激光焊接检测传感器中常用的聚焦镜、分色镜、分光镜、滤波镜

目录 一、聚焦镜 1.1 原理及作用 1.2 性能指标 1.3 应用 1.4 类型 二、分色镜 2.1 原理及应用 2.2 种类 2.3 特点 2.4 注意事项 2.5 性能指标 三、分光镜 ​编辑 3.1 分光镜的类型 3.2 分光镜的工作原理 3.3 分光镜的应用 3.4 分光镜的参数 3.5 分光镜的优点…

LabVIEW控制Trio控制器

将LabVIEW与Trio控制器结合&#xff0c;可以实现对复杂运动系统的控制和监测。以下是详细的方法和注意事项&#xff1a; 一、准备工作 软件安装&#xff1a; 安装LabVIEW开发环境&#xff0c;确保版本兼容性。 安装Trio控制器的相关驱动程序和软件&#xff0c;如Trio Motion …

Springboot集成GRPC

Springboot集成GRPC 一、springboot版本二、GRPC的pom依赖2.1 服务端2.2 客户端3.构建依赖 三、配置文件服务端客户端 四、 demo4.1 编写proto文件4.2 生成文件4.3 服务端重写方法4.4 客户端调用该方法 五、测试 一、springboot版本 <groupId>org.springframework.boot&l…

数理逻辑:1、预备知识

17.1 命题和联结词 ​ 命题&#xff1a;可以判定真假的陈述句。&#xff08;则悖论&#xff0c;祈使句&#xff0c;疑问句都不是命题&#xff09; ​ 原子命题&#xff1a;不能被分割为更小的命题的命题 例如&#xff1a; 2既是素数又是偶数 可以由$p: 2 是素数&#xff0c;…

AI模型抉择:开源VS闭源,谁主沉浮?

AI模型抉择&#xff1a;开源VS闭源&#xff0c;谁主沉浮&#xff1f; &#x1f604;生命不息&#xff0c;写作不止 &#x1f525; 继续踏上学习之路&#xff0c;学之分享笔记 &#x1f44a; 总有一天我也能像各位大佬一样 &#x1f3c6; 博客首页 怒放吧德德 To记录领地 &am…

C++字符编码 cppp-reiconv库使用详解

经常写一些控制台小程序&#xff0c;常常会遇到输出中文乱码的问题&#xff0c;在windwos下可以使用MultiByteToWideChar转换字符编码&#xff0c;但跨平台就需要cppp-reiconv这样的第三方字符编码处理库&#xff0c;且开源。 一、下载cppp-reiconv库的源码和静/动态库 GitHu…