dfs
回溯搜索。(哈密顿回路 / TSP问题)
给一个无向图,给定起点,再给几个商店(最多10
个),要求找一条从起点出发、经过所有商店、再回到起点的最短路。
因为商店数量很少,肯定是搜索了。可以用dfs
,也可以用状压dp
?(不太懂,没试,像这个一样?)
需要先最短路预处理一下,因为只会用到起点和这些商店之间的最短路,所以先把这些关键点的下标映射到一个数组,然后以每个关键点为源点来一遍最短路。
然后就是写这个dfs
了,从起点出发,用一个vis[]
来避免单次寻路中重复访问商店,最后都商店都访问完了再和起点连一下、记录结果。这个dfs
的参数有三个,这种题dfs
的写法务必掌握(真的很基础了。。)。其实很像n皇后问题(参考这个),递归函数自带了回溯的过程,让dfs
夹在两个标记设置中间。这样就枚举出了所有可能的路线,时间复杂度O(n!)(瞎猜的 )
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>
#include <string>
#include <queue>
using namespace std;
const int INF = 1e9;
const int MAXN = 1e5;
int N, M, T, X;
int S[11]; // 包括起点,最多11个
struct Edge
{
int n, w;
};
vector<Edge> ve;
vector<int> v[MAXN];
int d[11][MAXN]; // 以(起点以及所有商店)为源点执行最短路
bool inq[MAXN];
bool vis[11]; // dfs判断是否已访问,一次连续的向下探索过程不能重复访问
int ans;
void init()
{
ve.clear();
for (int i = 0; i < N; i++) v[i].clear();
memset(vis, 0, sizeof vis);
ans = INF;
}
void spfa(int s)
{
queue<int> q;
memset(inq, 0, sizeof inq);
fill(d[s], d[s] + N, INF);
q.push(S[s]);
inq[S[s]] = true;
d[s][S[s]] = 0;
for (; !q.empty();)
{
int t = q.front();
q.pop();
inq[t] = false;
for (int i = 0; i < v[t].size(); i++)
{
int n = ve[v[t][i]].n;
int w = ve[v[t][i]].w;
if (d[s][t] + w < d[s][n])
{
d[s][n] = d[s][t] + w;
if (!inq[n])
{
q.push(n);
inq[n] = true;
}
}
}
}
}
void dfs(int n, int sum, int step)
{
if (step == X)
{
ans = min(ans, sum + d[n][0]);
return;
}
for (int i = 1; i <= X; i++) // 中间没必要访问起点,全访问完了最后回起点就行了
{
if (!vis[i])
{
vis[i] = true;
//for (int j = step; j--;) cout << " ";
//cout << "(" << n << "," << sum << "," << step << ")调用(" << i << "," << sum + d[n][S[i]] << "," << step + 1 << ")" << endl;
dfs(i, sum + d[n][S[i]], step + 1);
//for (int j = step; j--;) cout << " ";
//cout << "(" << n << "," << sum << "," << step << ")调用结束(" << i << "," << sum + d[n][S[i]] << "," << step + 1 << ")" << endl;
vis[i] = false;
}
}
}
int main()
{
int a, b, c;
scanf("%d", &T);
for (; T--;)
{
scanf("%d%d", &N, &M);
init();
for (int i = 0; i < M; i++)
{
scanf("%d%d%d", &a, &b, &c);
ve.push_back(Edge{ b,c });
ve.push_back(Edge{ a,c });
v[a].push_back(i << 1);
v[b].push_back(i << 1 | 1);
}
S[0] = 0;
spfa(0);
scanf("%d", &X);
for (int i = 1; i <= X; i++)
{
scanf("%d", &S[i]);
spfa(i);
}
dfs(0, 0, 0); // 总共最多就10个点
printf("%d\n", ans);
}
return 0;
}