对使用状态压缩和动态规划求hamilton最短路径的理解

hamilton问题

即:

给定一张 n(n≤20) 个点的带权无向图,点从 0~n-1 标号,求起点 0 到终点 n-1 的最短Hamilton路径。 Hamilton路径的定义是从 0 到 n-1 不重不漏地经过每个点恰好一次。
通俗来说就是,图表示大概为这个样子
来源:http://home.ustc.edu.cn/~xxxujian/homework8.pdf

本文中,代码里的^ 代表xor, 文本中的^ 代表次方

最短路径

本文主要解释通过状态压缩和动态规划优化后的枚举求最短路径算法,本质上还是枚举算法

预设操作二进制的一些小算法

  • (n>>k)&1 取出n在二进制状态中第k位的值(即把n右移k位,用与1按位与取出右移后数字的第0位也就是原数的第k位),在本例中即为取是否经过第k个点
  • n^(1<<k) 给n在二进制状态中的第k位取反,在本例中即代表如果经过那个点之前的情况

状态压缩

本算法中,状态压缩主要体现在用一个int变量代替bool数组来优化时空
因为一个bool在C++中占1byte,而一个int变量在C++占4bytes (来源: https://docs.microsoft.com/zh-tw/cpp/cpp/fundamental-types-cpp?view=msvc-160#sizes-of-built-in-types)
所以当bool数组中数量超过4个时使用int类型代替必会节省空间,而本例中最多有20个需要表示的状态(状态指经过该点或没经过,n <= 20), 所以最多int变量需要有20为,也就是2^20(1<<20的值),为1,048,576,而int最大为2147483647, 所以在本例适用.
所以在代表每个节点的状态的int中,二进制状态的第x位的值即代表是否经过第x点

动态规划

动态规划在本例中就是把寻找最短路径分割成很多个小任务,比如找到起点到第一个节点的的最短距离,然后找到第一个节点到第二个节点的最短距离,再把这两个相加就得出了起点到第二个节点的最短距离,从而推到终点算出起点到终点的最短距离。
所以就需要一个int数字来储存最短距离数据
本例算法中以f为变量名,大小为(2 ^ 20)*20,为什么是2 ^ 20是因为以一个最大为2^20的数表示当前经过了哪些点即有2 ^ 20种行进方法,而有20列是代表在当前状态(前面的数表示的经过状态)下,到每个点的最短距离.
比如: f[12][4] 就代表从起点经过第2第3个点后(12的二进制是1100)到第4个点的最短距离
在下面给出的代码中空出了f[0]中任意内容是因为是没有意义的,起点即为第0位,无论什么情况都要经过起点,也就是表达状态的数字无论什么情况都要大于1

预设

本例中把无穷大(也就是动态规划里每个最短路径的初始值)设置成0x3F及以上(和0x7F一样达到了10^9量级)基本数据在本例中不会达到这种量级,而且无穷大和无穷大相加还是等于无穷大,无穷大和其他数相加也是等于无穷大,并且也不会造成数据溢出。

输入

输入为一个n*n二维数组,第[i, j]项,也就是二维数组中第i行,第j列的数值表示i到j的权重(也可以理解为距离吧),在下面用weight作为变量名,大概输入长这样:
来源: http://home.ustc.edu.cn/~xxxujian/homework8.pdf
每个值代表该行到该列的点的权重(在本例中也可以理解成距离)或该列到该行的距离,如果不可以直接到达即为无穷

输出

为int

算法代码

#include <iostream>
#include <memory.h>
using namespace std;

int f[1<<20][20]; //创建一个2^20*20的数组储存最短距离数据
// n为有几个点, weight为带权无向图的权重
int hamilton(int n, int weight[20][20]){
   
	memset(f, 0x3f, sizeof(f)); // 初始化f每一项为无穷大
	f[1][0] = 0; // 设置最开始的最短距离,即起点(因为第一个下标的1在二进制下第0位是1也就是代表经过第0点即起点的情况下)到起点(第二个下标是0即第0点)的距离为0
	for(int i = 1; i < (1 
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
哈密顿回路是一种经过图中每个节点一次且仅一次的回路。哈密顿回路问题是一个NP完全问题,因此没有已知的多项式时间算法可以解决这个问题。不过,可以使用启发式算法来解决近似的问题。 下面是一个使用Java实现的近似算法: ```java import java.util.*; public class HamiltonianPath { private static int[][] graph; // 图 private static int[] path; // 存储路径 private static boolean[] visited; // 标记是否访问过 private static int n; // 节点数 public static void main(String[] args) { Scanner sc = new Scanner(System.in); n = sc.nextInt(); graph = new int[n][n]; path = new int[n]; visited = new boolean[n]; // 构建图 for(int i = 0; i < n; i++) { for(int j = 0; j < n; j++) { graph[i][j] = sc.nextInt(); } } // 从第一个节点出发 path[0] = 0; visited[0] = true; if(findHamiltonianPath(1)) { // 打印路径 for(int i = 0; i < n; i++) { System.out.print(path[i] + " "); } } else { System.out.println("No Hamiltonian Path exists"); } } // 查找哈密顿路径 private static boolean findHamiltonianPath(int pos) { // 如果已经遍历完所有节点 if(pos == n) { // 判断最后一个节点是否与第一个节点相邻 if(graph[path[pos - 1]][path[0]] == 1) { return true; } else { return false; } } // 遍历其它节点 for(int i = 1; i < n; i++) { if(isValid(i, pos)) { path[pos] = i; visited[i] = true; if(findHamiltonianPath(pos + 1)) { return true; } // 回溯 visited[i] = false; } } return false; } // 判断节点是否可达 private static boolean isValid(int node, int pos) { // 如果节点已经被访问过,返回false if(visited[node]) { return false; } // 如果前一个节点与当前节点不相邻,返回false if(graph[path[pos - 1]][node] == 0) { return false; } return true; } } ``` 在这个算法中,我们使用了回溯的方法来查找哈密顿路径。我们从第一个节点开始,依次尝试访问其它节点,直到找到一条哈密顿路径或者遍历完所有节点。在查找过程中,我们使用visited数组来标记节点是否已经被访问过,使用path数组来存储路径。isValid方法用来判断节点是否可达。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值