dd爱探险（状压DP）

题目描述:

链接：dd爱探险
来源：牛客网

星际中有n个空间站，任意两个空间站间可以相互跳跃，由空间站x跳跃到空间站y所需要的代价为P[x][y],注意不保证p[x][y]=p[y][x]，dd可以任意选择出发的空间站，并通过恰好n−1次跳跃把所有空间站跳完，并且dd必须选择2次跳跃，其中一次跳跃中进行重力加速，另一次跳跃中进行反重力加速，重力加速会导致当前跳跃代价变为0，反重力加速会导致当前跳跃代价翻倍(乘2)，问跳完所有空间站所需要最小代价

输入描述:

第一行一个数n(3≤n≤16) 接下来n行每行n个数，第x行第y个数表示p[x][y]
(0≤p[x][y]≤100000)

输出描述:

一个数,表示最小代价

示例1 :

输入：

3
0 2 1
2 0 1
3 1 0

输出：

2

说明：

1->2重力加速
2->3反重力加速
代价0+1*2=2

思路：

借鉴牛客竞赛 @sudo_Yiii
状压DP
最短Hamilton路径，加拓展
mp存一下两点间距离，dp[i][j][k]表示在第i中情况下到达点j的状态k。
k = 0表示未经过加速
k = 1表示经过一次反重力
k = 2表示经过一次重力
k = 3表示两次加速都已经经历过

dp[i][j][1] = min(dp[i][j][1], min(dp[i - (1 << j)][k][0] + mp[k][j] * 2, dp[i - (1 << j)][k][1] + mp[k][j]));
表示到这个点 经历一次反重力 与 已经经过一次反重力的经历一次正常 相比较
dp[i][j][2] = min(dp[i][j][2], min(dp[i - (1 << j)][k][0], dp[i - (1 << j)][k][2]) + mp[k][j]);
同理表示这个点 经历一次重力 与 已经经过一次重力的经历一次正常 相比较

dp[i][j][3] = min(dp[i][j][3], min(dp[i - (1 << j)][k][3] + mp[k][j], min(dp[i - (1 << j)][k][1], dp[i - (1 << j)][k][2] + mp[k][j] * 2)));
最终状态：经历过反重力的经历一次重力 与 经历过重力的经历一次反重力 与 都经历过的经历一次正常 三者相比较
最后取最小值即可

代码：

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;

int mp[16][16], dp[(1 << 16) + 10][16][4];

int main()
{
    int n;
    cin >> n;
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++) cin >> mp[i][j];
    memset(dp, 127,sizeof dp);
    for(int i = 0;i < n;i++) dp[1 << i][i][0] = 0;
    for(int i = 1;i < (1 << n);i++)
        for(int j = 0;j < n;j++){
            if(i >> j & 1){
                for(int k = 0;k < n;k++){
                    if((i >> k) & 1){
                        dp[i][j][0] = min(dp[i][j][0], dp[i - (1 << j)][k][0] + mp[k][j]);
                        dp[i][j][1] = min(dp[i][j][1], min(dp[i - (1 << j)][k][0] + mp[k][j] * 2, dp[i - (1 << j)][k][1] + mp[k][j]));
                        dp[i][j][2] = min(dp[i][j][2], min(dp[i - (1 << j)][k][0], dp[i - (1 << j)][k][2]) + mp[k][j]);
                        dp[i][j][3] = min(dp[i][j][3], min(dp[i - (1 << j)][k][3] + mp[k][j], min(dp[i - (1 << j)][k][1], dp[i - (1 << j)][k][2] + mp[k][j] * 2)));
                    }
                }
            }
        }
    int ans = 0xffffff;
    for(int i = 0;i < n;i++) ans = min(ans, dp[(1 << n) - 1][i][3]);
    cout << ans << endl;
    return 0;
}