分支限界法求解旅行商问题

这篇博客介绍了如何运用分支限界法来解决2016年华为软件精英挑战赛中的一种特殊旅行商问题。问题涉及到带有权重的有向图,寻找从特定起点s到终点t的不成环路径,路径需经过指定的顶点集合V'。文章强调了图的特性、路径要求以及权重计算,并指出与经典旅行商问题的区别。作者使用C++实现了初步的分支限界算法,以求近似解,并在后续对代码进行了优化,以节省内存。
摘要由CSDN通过智能技术生成

问题来来自2016华为软件精英挑战赛。

 问题定义

给定一个带权重的有向图G=(V,E),V为顶点集,E为有向边集,每一条有向边均有一个权重。对于给定的顶点s、t,以及V的子集V',寻找从s到t的不成环有向路径P,使得P经过V'中所有的顶点(对经过V'中节点的顺序不做要求)。

若不存在这样的有向路径P,则输出无解,程序运行时间越短,则视为结果越优;若存在这样的有向路径P,则输出所得到的路径,路径的权重越小,则视为结果越优,在输出路径权重一样的前提下,程序运行时间越短,则视为结果越优。

说明:

1)图中所有权重均为[1,20]内的整数;

2)任一有向边的起点不等于终点;

3)连接顶点A至顶点B的有向边可能超过一条,其权重可能一样,也可能不一样;

4)该有向图的顶点不会超过600个,每个顶点出度(以该点为起点的有向边的数量)不超过8;

5)V'中元素个数不超过50;

6)从s到t的不成环有向路径P是指,P为由一系列有向边组成的从s至t的有向连通路径,且不允许重复经过任一节点;

7)路径的权重是指所有组成该路径的所有有向边的权重之和。


与经典旅行商问题相比,有几点小区别:

①是有向图

②起点终点是不同的点

③大图中的指定点集,难度更大


由于大赛对程序运行时间有限制,所以需要一个求近似解的算法。


但还是现用分支限界试试,用C++做了初步的实现,暂做记录。


3月7日下午两点半修改了代码,改用链表维护path,节省内存


#pragma once
#include <iostream>
#include <queue>
#include <unordered_map>
#include <unordered_set>
#include <set>
#include <vector>
#include <stack>
using namespace std;

struct Node {
	int id;
	Node* path_father;	//路径中父节点
	int already;		//path中V'中节点数目(默认source和destination不在V'中)
	int path_cost;
	Node(int i) : id(i), path_cost(0), already(0), path_father(NULL) {};
	Node(int i
分支限界法求解旅行商问题的步骤如下: 1. 确定问题的目标函数和约束条件,即TSP问题的目标是找到一条通过所有城市的最短路径。 2. 采用邻接矩阵表示城市之间的距离。 3. 定义结点的数据结构。每个结点表示一个旅行商的旅行路径,包括已访问的城市、未访问的城市、当前路径长度等信息。 4. 定义结点的比较函数,用于根据当前路径长度进行结点的排序。 5. 使用优先队列维护结点的集合,每次取出当前路径长度最小的结点。 6. 对于每个结点,枚举下一个要访问的城市,计算扩展后的路径长度,并生成新的结点加入优先队列中。 7. 不断重复以上步骤,直到找到一条符合要的最短路径或者优先队列为空。 以下是一个基于C语言的分支限界法求解TSP问题的示例代码: ``` #include <stdio.h> #include <stdlib.h> #include <string.h> #include <limits.h> #define N 5 int map[N][N] = { {0, 2, 9, 10, 5}, {2, 0, 6, 4, 8}, {9, 6, 0, 7, 6}, {10, 4, 7, 0, 8}, {5, 8, 6, 8, 0} }; typedef struct Node { int path[N]; int len; int visited[N]; } Node; int cmp(const void *a, const void *b) { Node *nodeA = *(Node **)a; Node *nodeB = *(Node **)b; return nodeA->len - nodeB->len; } void printPath(Node *node) { printf("Path: "); for (int i = 0; i < N; i++) { printf("%d ", node->path[i]); } printf("%d\n", node->path[0]); printf("Length: %d\n", node->len); } void solveTSP() { Node start; memset(&start, 0, sizeof(start)); start.path[0] = 0; start.visited[0] = 1; start.len = 0; Node *current = &start; Node *queue[N*N]; int front = 0, rear = 0; queue[rear++] = current; while (front < rear) { current = queue[front++]; if (front >= N*N) { break; } if (current->len >= queue[rear-1]->len) { continue; } for (int i = 1; i < N; i++) { if (current->visited[i] == 0) { Node *node = (Node *)malloc(sizeof(Node)); memcpy(node, current, sizeof(Node)); node->visited[i] = 1; node->path[node->len+1] = i; node->len += map[current->path[node->len]][i]; if (node->len >= queue[rear-1]->len) { free(node); continue; } if (node->len < INT_MAX) { queue[rear++] = node; } if (node->len < queue[0]->len) { qsort(queue, rear, sizeof(Node *), cmp); } } } } printPath(queue[0]); for (int i = 0; i < rear; i++) { free(queue[i]); } } int main() { solveTSP(); return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值