PAT1150 Travelling Salesman Problem

题目

The “travelling salesman problem” asks the following question: “Given a list of cities and the distances between each pair of cities, what is the shortest possible route that visits each city and returns to the origin city?” It is an NP-hard problem in combinatorial optimization, important in operations research and theoretical computer science. (Quoted from “https://en.wikipedia.org/wiki/Travelling_salesman_problem”.)

In this problem, you are supposed to find, from a given list of cycles, the one that is the closest to the solution of a travelling salesman problem.

Input Specification:

Each input file contains one test case. For each case, the first line contains 2 positive integers N (2<N≤200), the number of cities, and M, the number of edges in an undirected graph. Then M lines follow, each describes an edge in the format City1 City2 Dist, where the cities are numbered from 1 to N and the distance Dist is positive and is no more than 100. The next line gives a positive integer K which is the number of paths, followed by K lines of paths, each in the format:

n C1 C2 ……Cn​​

where n is the number of cities in the list, and C​i‘s are the cities on a path.

Output Specification:

For each path, print in a line Path X: TotalDist (Description) where X is the index (starting from 1) of that path, TotalDist its total distance (if this distance does not exist, output NA instead), and Description is one of the following:

TS simple cycle if it is a simple cycle that visits every city;
TS cycle if it is a cycle that visits every city, but not a simple cycle;
Not a TS cycle if it is NOT a cycle that visits every city.
Finally print in a line Shortest Dist(X) = TotalDist where X is the index of the cycle that is the closest to the solution of a travelling salesman problem, and TotalDist is its total distance. It is guaranteed that such a solution is unique.

Sample Input:

6 10
6 2 1
3 4 1
1 5 1
2 5 1
3 1 8
4 1 6
1 6 1
6 3 1
1 2 1
4 5 1
7
7 5 1 4 3 6 2 5
7 6 1 3 4 5 2 6
6 5 1 4 3 6 2
9 6 2 1 6 3 4 5 2 6
4 1 2 5 1
7 6 1 2 5 4 3 1
7 6 3 2 5 4 1 6

Sample Output:

Path 1: 11 (TS simple cycle)
Path 2: 13 (TS simple cycle)
Path 3: 10 (Not a TS cycle)
Path 4: 8 (TS cycle)
Path 5: 3 (Not a TS cycle)
Path 6: 13 (Not a TS cycle)
Path 7: NA (Not a TS cycle)
Shortest Dist(4) = 8

代码

#include <iostream>
#include <vector>
#include <algorithm>
#include <set> 
using namespace std;
const int INF = 1e9;
const int MAXV = 205;
int G[MAXV][MAXV];

int travel(vector<int> path, int &flag, int n){
	set<int> travelPath;
	int res = 0;
	int len = path.size();
	int begin = path[0];
	int end = path[len-1];
	if(begin == end)  flag = 0;// 简单圈 
	else flag = 2;
	int pre = path[0];
	travelPath.insert(pre);
	int cur;
	for(int i=1; i<len; i++){
		cur = path[i];
		travelPath.insert(cur);
		if(cur == begin && i < len - 1)  flag = 1;// 不止一个圈 
		if(G[pre][cur] != INF){
			res += G[pre][cur];
			pre = cur;
		}
		else{
			flag = 3;// 不可达 
			return INF;
		}
	}
	
	// 如果没有访问完 
	if(travelPath.size() < n)  flag = 2;
	return res; 	
}
int main(){
	int n;// 城市个数 
	int m;// 无向图的边
	int v1, v2, dis;
	scanf("%d %d", &n, &m);
	fill(G[0], G[0] + MAXV * MAXV, INF);
	for(int i=0; i<m; i++){
		scanf("%d %d %d", &v1, &v2, &dis);
		G[v1][v2] = dis;
		G[v2][v1] = dis;
	}
	int k, num, flag, minDis = INF, minNum = 0;
	scanf("%d", &k);
	vector<int> path;
	set<int> cmPath;
	for(int i=1; i<=k; i++){
		path.clear();
		flag = -1;
		scanf("%d", &num);
		path.resize(num);
		for(int j=0; j<num; j++){
			scanf("%d", &path[j]);
		}
		
		int totalDist = travel(path, flag, n);
		
		if(flag == 0){
			printf("Path %d: %d (TS simple cycle)\n", i, totalDist);
			// 更新最短路径 
			if(totalDist < minDis){
				minDis = totalDist;
				minNum = i;
			}
		}
		else if(flag == 1){
			printf("Path %d: %d (TS cycle)\n", i, totalDist);
			// 更新最短路径 
			if(totalDist < minDis){
				minDis = totalDist;
				minNum = i;
			}
		}
		else if(flag == 2){
			printf("Path %d: %d (Not a TS cycle)\n", i, totalDist);
		}
		else if(flag == 3){
			printf("Path %d: NA (Not a TS cycle)\n", i);
		}
		
	}
	printf("Shortest Dist(%d) = %d\n", minNum, minDis);
	return 0;
}