1 条题解

  • 1
    @ 2020-04-25 13:10:48

    如果考虑DP的话,那么应该如何划分状态?
    首先需要确定我们要考虑的因素,因为Hamiton路径每个点只能经过一次,并且要经过所有的点,所以我们应该考虑当前的状态经过哪些点,但是如果只考虑当前经过哪些点的话,我们没有办法确定当前的最短路径是多少,即我们不知道我们是从哪个点转移过来的,那么我们就需要考虑我们当前所在的是哪个点,那么我们就可以知道我们从哪里转移过来,边权是多少,我们就可以更新我们的最短路。

    我们需要考虑的因素:

    1.当前经过哪些点

    2.当前在那哪个点

    那么我们就应该考虑状态:

    我们用\(f[i][j]\)表示当前的状态,其中i是当前所经过点的集合,j是当前所在的点。

    考虑完状态之后我们就需要考虑如何转移:

    因为每两个点之间都存在路径,所以我们可以从任意一个点转移过来,但是我们要确保向当前状态转移的状态的当前状态点没有访问且访问了向当前状态转移的状态的点。

    即:

    \(f[i][j] = min([state_k][k]+w[k][j],f[i][j])\)

    其中\(state_k\) \(j\)未被访问,且包括\(k\)这个点。

    换言之就是\(state_k\)为集合\(i\)中\(j\)未被访问,但\(k\)被访问的集合,这样就与当前状态建立了联系。

    那么我们首先可以枚举集合\(i\),即枚举通过的点来扩展状态。

    接下来我们可以通过枚举合法的\(j\)与合法的\(k\)来确定状态的最优解。

    \(j\)合法:即\(i\)集合中包括\(j\)。

    \(k\)合法:即\(i\)集合中\(k\)已被访问。

    那么我们用什么来确定每个点是否被访问?

    因为数据规模比较小,所以我们可以用一个32位二进制整数n的第i位来保存第i个点是否被访问。

    状态压缩

    二进制状态压缩,是将一个长度为m的bool数组用一个m位二进制整数表示存储的方法。可以利用如下的方法实现原bool数组对应下表的存取。

    操作 | 运算 |

    -|-|
    取出整数n在二进制表示下第k位| (n >> k) & 1 |
    取出整数n在二进制表示下第0~k-1位(后k位)|n & ((1<<k)-1) (即将后k位全都设为1并用n进行与运算)|
    把整数n在二进制表示状态下的第k位取反| n^(1<<k)|
    把整数n在二进制表示下的第k为赋值1| n or (1<<k)|
    对整数n在二进制表示下的第k位赋值为0 |n &(~(1<<k))


    那么我们就可以得到完整的状态转移方程:

    \(f[i][j] = min(f[i xor(1<<j)][k] + w[k][j],f[i][j])\)

    且满足于\((i >> j) \) & \( 1 = 1\)

    \((i xor (i << j)) >> k \) & \(1 = 1\)


    其中,\((i >> j) \) & \( 1 = 1\)(是确保当前状态下i的第\(j\)位一定为1,即\(j\)点被访问。

    而 \(i xor(1<<j)\)则代表将\(i\)的第\(j\)位取反 即在\(i\)中\(j\)未被访问(其余不变)的状态。

    其复杂度为\(O(2^n*n)\)

    #include <iostream>
    #include <cstring>
    using namespace std;
    int f[1 << 20][25],w[25][25],n;
    int main(){
        cin >> n;
        for (int i = 0; i < n; i++)
            for (int j = 0; j < n; j++)
                cin >> w[i][j];
        
        memset(f,0x3f,sizeof(f));
        f[1][0] = 0;
        for (int i = 0; i <= 1 << n; i++){
            for (int j = 0; j < n; j++){
                if (i&(1<<j) == 1 << j){
                    for (int k = 0; k < n; k++){
                        if ((i^1<<j) >> k & 1){
                            f[i][j] = min(f[i][j],f[i^1<<j][k] + w[k][j]);
                        }
                    }
                }
            }
        }
        cout << f[(1<<n)-1][n - 1];
        return 0;
    }
    
  • 1

信息

ID
1003
难度
9
分类
(无)
标签
(无)
递交数
11
已通过
1
通过率
9%
被复制
1
上传者