动态规划——状态压缩DP

我所理解的状态压缩DP就是将所有状态用二进制数表示，然后根据前一个合法的二进制数表示的状态进行状态转移，得到最终答案。

蒙德里安的梦想

• 题目描述

求把 N×M 的棋盘分割成若干个 1×2 的的长方形，有多少种方案。
例如当 N=2，M=4 时，共有 5 种方案。当 N=2，M=3 时，共有 3 种方案。

• 分析

这道题的核心是要发现所有的情况之和横着放的小矩形的情况有关。当所有横着放的矩形确定以后，竖着放的小矩形只能从剩余的空间中塞进去。
所以可以考虑每一列，在每一列中有横着放或者前一列中横着放插过来的小矩形的行为1，没有的为0，则每一列可以表示成一个二进制的数，每个二进制数都表示一个状态。
然后需要考虑所有合法的情况，首先，前一列插过来的和本列横着放的小矩形不能在同一行，即j&k=0，此外，每一列连续的0的个数应该是偶数个，这样才能塞下剩余的竖着放的小矩形。

• 代码

#include<bits/stdc++.h>
using namespace std;
const int N = 12,M = 1<<N;
long long f[N][M];
bool st[M];
int n,m;

int main()
{
    while(scanf("%d%d",&n,&m)&&(n||m)){
    //遍历找到所有符合第二种的合法情况。
        for(int i=0;i<1<<n;i++){
            int cnt = 0;
            st[i] = true;
            for(int j=0;j<n;j++){
                if(i>>j&1){
                    if(cnt&1)st[i] = false;
                    cnt = 0;
                }
                else cnt++;
            }
            if(cnt&1)st[i] = false;
        }
        //初始化，每次状态弄为0，f[0][0] = 1,表示第0列放0个矩形有1种情况。
        memset(f,0,sizeof(f));
        f[0][0] = 1;
        for(int i=1;i<=m;i++){
            for(int j=0;j<1<<n;j++){
                for(int k=0;k<1<<n;k++){
                    if((j&k)==0&&(st[j|k]))f[i][j] += f[i-1][k];
                }
            }
        }
        //最终答案为第m列没有外伸的情况。
        printf("%lld\n",f[m][0]);
    }
    
    return 0;
}

最短Hamilton路径

• 题目描述

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

• DP分析

由于此题要走过所有点的最短路径，所以不能用前边学过的图论的算法，需要另辟蹊径。用f[i,j]表示从0走到j的所有路径的集合，其中j表示的是点的标号，i是一个二进制数，其从低到高某位上为1就表示取第该位对应的点。
因此，f[i,j] = min(f[i-1<<j,k]+w[k][j]。
最终得答案是所有点都取到即i=1<<n-1;

• 代码

#include<bits/stdc++.h>
using namespace std;
const int N = 20,M = 1<<N;
int f[M][N],w[N][N];
int n;


int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)scanf("%d",&w[i][j]);
    //初始化为最大值   
    memset(f,0x3f,sizeof(f)); 
    //从0出发，0到0路径为0
    f[1][0] = 0;
    for(int i=0;i<1<<n;i++)
        for(int j=0;j<n;j++)
            if(i>>j&1)
                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]);
    
    printf("%d",f[(1<<n)-1][n-1]);
    return 0;
}