状态压缩dp

状态压缩dp一般用来解决n<=20的情况。

状态压缩dp主要有两种类型的题目——棋盘式和集合式。

棋盘式

对于棋盘式的状态压缩dp，我们的做法是将棋盘抽象为一个状态（通常用二进制进行抽象），用最后一次摆放的位置来划分状态，实现转移。

291. 蒙德里安的梦想 - AcWing题库

棋盘式dp的经典题目。

通用的做法是把每一行抽象为一个二进制数。

对于本题，我们要先明白一个重要的性质：所有矩形只有横向和纵向两种摆放方式，在摆放完某一个方向的所有矩形后，我们的摆放方式就被唯一地确定下来了。因此，我们只需要考虑某个方向的矩形如何摆放。

以横向为例。

在本题中，对于每一列，我们用“该列是否为一个矩形的末尾”来表示状态。f(i,j)表示，当前已经摆放完前i列，第i列状态为j。

在状态转移时，我们以“i-1列的状态”为划分依据进行转移，我们枚举所有i-1列的合法状态，转移到f(i,j)中。

什么样的状态是合法的呢？

在本题中，可能存在的不合法情况有：

ac代码如下：

#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
typedef long long ll;
const int N=15,M=1<<11;
int n,m;
ll f[N][M];
vector<int> state[M];
int main()
{
    while(cin>>n>>m&&n!=0&&m!=0)
    {
        memset(f,0,sizeof(f));
        for(int i=0;i<1<<n;i++)
        {
            state[i].clear();
            
            for(int j=0;j<1<<n;j++)
            {
                if((i&j)!=0) continue;
                int t=i|j;
                int cnt=0;
                bool is_valid=1;
                for(int k=0;k<n;k++)
                {
                    if((t>>k&1)==0)
                        cnt++;
                    else
                    {
                        if(cnt%2) 
                            is_valid=0;
                        cnt=0;
                    }
                }
                if(cnt%2) is_valid=0;
                if(is_valid) state[i].push_back(j);
            }
        }
        f[0][0]=1;
        f[1][0]=1;
        for(int i=2;i<=m+1;i++)
            for(int j=0;j<1<<n;j++)
            {
                for(int k=0;k<state[j].size();k++)
                {
                    f[i][j]+=f[i-1][state[j][k]];
                }
            }
        printf("%lld\n",f[m+1][0]);
    }
}

集合式

91. 最短Hamilton路径 - AcWing题库

集合式dp的经典题目。同样地，我们将当前走过了哪些点作为状态，状态用一个二进制数表示。

f(i,j)表示当前位于第j个点，走过状态是i。

我们以“上一个点是什么”作为状态划分依据。遍历状态i，当某个点状态为1时，可以通过f(i,j)=min(f(i-上一个点,上一个点)+a[上一个点][j],f(i,j))来进行状态转移。

ac代码如下：

#include<iostream>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=25,M=1<<20;
int n;
int a[N][N];
ll f[M][N];
int main()
{
    memset(f,0x3f,sizeof(f));
    scanf("%d",&n);
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {
            scanf("%d",&a[i][j]);
        }
    f[1][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(k!=j&&i>>k&1)
                    {
                        f[i][j]=min(f[i-(1<<j)][k]+a[k][j],f[i][j]);
                    }
                }
        }
    }
    
    printf("%lld",f[(1<<n)-1][n-1]);
}