DP模板—状压DP

一、状压DP

利用二进制表示当前状态

二、例题

（1）蒙德里安的梦想

当我们先尽可能摆放完$1\ast 2$的方格时，剩下的空格的摆放方式是唯一的——只能用$2\ast 1$补完或不能构成合法方案

所以我们考虑只摆放$1\ast 2$方格的方案并分析能否构成一个合法方案，

状态表示：

$f[i][j]代表前i-1列已经用1\ast 2的方格摆放好，且从第i-1列伸出到第i列为状态j的所有方案$

$j的第k位为1表示(i,k-1)的位置已被填上$

状态转移：

对于状态$f[i][j]$，当$f[i-1][k]$满足$(k\&j)==0$并且第$i$列摆放上1*2的方格合法时，$f[i][j]+=f[i-1][k]$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=12,M=1<<N;
int n,m;
ll f[N][M];
bool vis[M];
vector<int>ve[M];
int main(){
    while(cin>>n>>m,n||m){
        //预处理出与状态j能构成合法方案的所有方案
        for(int i=0;i<1<<n;i++){
            int cnt=0,f=1;
            for(int j=0;j<n;j++){
                if(i>>j&1){
                    if(cnt&1){
                        f=0;
                        break;
                    }
                    else cnt=0;
                }
                else cnt++;
            }
            if(cnt&1) f=0;
            vis[i]=f;
        }
        
        for(int i=0;i<1<<n;i++){
            ve[i].clear();
            for(int j=0;j<1<<n;j++){
                if((i&j)==0&&vis[i|j]) ve[i].emplace_back(j);
            }
        }
        
        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(auto z:ve[j]) f[i][j]+=f[i-1][z];
            }
        }
        cout<<f[m][0]<<endl;
    }
    return 0;
}

最短Hamilton路径

$f[i][j]$表示最后走到$j$点状态为$i$的最小花费

#include<bits/stdc++.h>
using namespace std;
const int N=21,M=1<<N;
int n,m;
int g[N][N];
int f[M][N];
int main(){
    cin>>n;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++) cin>>g[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>>j&1){
                for(int z=0;z<n;z++){
                    if(i>>z&1){
                        f[i][j]=min(f[i][j],f[i-(1<<j)][z]+g[z][j]);
                    }
                }
            }
        }
    }
    cout<<f[(1<<n)-1][n-1];
    return 0;
}