二进制状压dp CH - 0103 最短Hamilton路径

今天初识状态压缩dp，这是一个很有代表性的题目

题意如下

Click_here_to_visit

描述

给定一张 n(n≤20) 个点的带权无向图，点从 0~n-1 标号，求起点 0 到终点 n-1 的最短Hamilton路径。
Hamilton路径的定义是从 0 到 n-1 不重不漏地经过每个点恰好一次。

输入格式

第一行一个整数n。

接下来n行每行n个整数，其中第i行第j个整数表示点i到j的距离（一个不超过10^7的正整数，记为a[i,j]）。

对于任意的x,y,z，数据保证 a[x,x]=0，a[x,y]=a[y,x] 并且 a[x,y]+a[y,z]>=a[x,z]。

输出格式

一个整数，表示最短Hamilton路径的长度。

有注意到n <= 20还是蛮小的 （一点都不小） ，所以我们可以开一种很奇怪的数组 （哪里奇怪了） 。为了方便的表示每个时刻哪些点已经走过，哪些没有，我们引进一个n位的二进制数，若其第i位上是1，则表示第i位已经被走过，反之亦然。这个二进制i是算法的精髓，每种状态都对应着一个i,每种状态又能延伸出新的状态，所以我们可以构造一个 诡异的 二维数数组f[i][j], 其具体意义是在状态‘i’上处于j点的最短Hamilton路径长度

状态转移方程：f[i][j] = min{f[i ^ 1 << j][k] + map[j][k]}

详情请看有注释又友好的代码：

#pragma G++ optimize(2)
#pragma GCC optimize(2)
#include<iostream>
#include<algorithm>
#include<stack>
#include<vector>
#include<iostream>
#include<climits>
#include<queue>
#include<cassert>
#include<iomanip>
#include<cmath>
#include<string> 
#include<cstdio>
#include<cstring>
#define _rep(i, a, b) for(ll i = (a); i <= (b); ++i)
#define _rev(i, a, b) for(ll i = (a); i >= (b); --i)
#define _for(i, a, b) for(ll i = (a); i <(b); ++i)
#define _rof(i, a, b) for(ll i = (a); i >(b); --i)
#define maxn 22
#define maxm 109
#define oo 0x3f
#define ll long long
#define met(a,b) memset((a),(b), sizeof(a))
#define db double
#define eps 1e-8
#define lowbit(x) x & (-x)
//#define DEBUG
using namespace std;
int map[maxn][maxn], n;
int f[1 << 20][maxn];
int main() {
 ios::sync_with_stdio(0);
 cin >> n;
 _for(i, 0, n) {
  _for(j, 0, n) {
   cin >> map[i][j];
  }
 }
 met(f, oo);
 f[1][0] = 0;
 _for(i, 1, 1 << n)//第一维,枚举所有“被经过的点的状态” 
 _for(j, 0, n)//枚举当前所在的位置
  if (i >> j & 1) // 如果被经过的点的状态里面有j
  _for(k, 0, n)//枚举上一个点
   if ((i ^ 1 << j) >> k & 1)// 没过j点的状态过了k点 
    f[i][j] = min(f[i][j], f[i ^ 1 << j][k] + map[k][j]);
     cout << f[(1 << n) - 1][n - 1];  
}