ABC 180 E - Traveling Salesman among Aerial Cities (状态压缩DP)

经典TSP问题，状压一下就可以了，代码是xxy大佬的。

// #pragma GCC optimize(2)
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <random>
#include <cctype>
#include <bitset>
#include <string>
#include <vector>
#include <queue>
#include <cmath>
#include <stack>
#include <set>
#include <map>
#define IO                       \
    ios::sync_with_stdio(false); \
    // cout.tie(0);
using namespace std;
// int dis[8][2] = {0, 1, 1, 0, 0, -1, -1, 0, 1, -1, 1, 1, -1, 1, -1, -1};
typedef unsigned long long ULL;
typedef long long LL;
typedef pair<int, int> P;
const int maxn = 2e5 + 10;
const int maxm = 1e6 + 10;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int inf = 0x3f3f3f3f;
const LL mod = 1e9 + 7;
const double eps = 1e-8;
const double pi = acos(-1);
// int dis[4][2] = {1, 0, 0, -1, 0, 1, -1, 0};
// int m[13] = {0, 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
struct Node
{
    LL a, b, c;
} a[maxn];
LL dis[20][20];
LL bit[20]; // 记录 1 的位置
LL dp[1 << 20][20];
LL cal(int i, int j)
{
    return abs(a[i].a - a[j].a) + abs(a[i].b - a[j].b) + max(0LL, a[j].c - a[i].c);
}
int main()
{
#ifdef WXY
    freopen("in.txt", "r", stdin);
//	 freopen("out.txt", "w", stdout);
#endif
    IO;
    int n;
    cin >> n;
    for (int i = 1; i <= n; i++)
        cin >> a[i].a >> a[i].b >> a[i].c;

    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            dis[i][j] = cal(i, j);
    // 预处理 两点间的距离
    bit[0] = 1;
    for (int i = 1; i <= n; i++)
        bit[i] = bit[i - 1] << 1;

    int all = bit[n] - 1; // 所有种情况 all 为全一的情况
    for (int i = 1; i <= all; i++)
        for (int j = 1; j <= n; j++)
            dp[i][j] = 1e9; // 初始化

    dp[1][1] = 0;                 // 当这条路上只有一个点且这个点为1的时候，显然还没有发生转移，花费为 0
    for (int i = 1; i < all; i++) // 所有状态 用 01 串来枚举
    {
        for (int j = 1; j <= n; j++)
        {
            if (i & bit[j - 1]) // 假如点 j 在状态集合中，那么可以通过这个点转移
            {
                // 用 bit[j-1] 找到对应位置的 1
                for (int k = 1; k <= n; k++)
                {
                    if (!(i & bit[k - 1])) // 当前点还没有在路径中，那么可以转移到这个点
                        dp[i | bit[k - 1]][k] = min(dp[i][j] + dis[j][k], dp[i | bit[k - 1]][k]);
                    // 更新当前最终点为 k 的时候的最优解
                }
            }
        }
    }
    LL ans = INF;
    for (int i = 2; i <= n; i++)
        ans = min(ans, dp[all][i] + dis[i][1]); //枚举一个点 使 1号点经过所有点后到这个点加上从这个点回到 1 号点最小
    cout << ans;
    return 0;
}