状压dp

蒙德里安的梦想

求把NM的棋盘分割成若干个12的的长方形，有多少种方案。

例如当N=2，M=4时，共有5种方案。当N=2，M=3时，共有3种方案。

输入格式
输入包含多组测试用例。

每组测试用例占一行，包含两个整数N和M。

当输入用例N=0，M=0时，表示输入终止，且该用例无需处理。

输出格式
每个测试用例输出一个结果，每个结果占一行。

数据范围 1≤N,M≤11

输入样例：
1 2
1 3
1 4
2 2
2 3
2 4
2 11
4 11
0 0
输出样例：
1
0
1
2
3
5
144
51205

代码：
 

#include <iostream>
#include <cstring>

using namespace std;

const int N = 12, M = 1 << N;
int n, m;
long long f[N][M]; // f[i][j] : 第i行，j状态的方案数， j由为一个二进制数，用来表示选择竖长方形的状态，
bool st[M];

int main()
{
    while(cin >> n >> m  && (m || n))
    {
        memset(f, 0, sizeof f);
        memset(st, false, sizeof st);
        
        for(int i = 0; i < (1 << m); i++) // 用数组储存每个状态是否存在落单的列
        {
            st[i] = true;
            int cnt = 0;
            
            for(int j = 0; j < m; j++)
            {
                if((i >> j & 1) == 0) cnt ++;
                else
                {
                    if(cnt & 1)
                    {
                        st[i] = false;
                        break;
                    }
                    cnt = 0;
                }
            }
            if(cnt & 1) st[i] = false;
        }
        
        for(int i = 0; i <= (1 << m); i++) // 初始第一列
            if(st[i]) f[1][i] = 1;
        
        for(int i = 2; i <= n ; i++)
            for(int j = 0; j < (1 << m); j++)
                for(int k = 0; k < (1 << m); k++)
                    if((k & j) == 0 && st[k | j]) // 欲使j能由k转移而来，需满足k，j不存在冲突的选择，
                        f[i][j] += f[i-1][k];     // 且选完后不能存在单独的列，否则将无法放置横长方形。

        cout << f[n][0] << endl; // 最后一列已无法选择任何竖长方形，否则将超出。
    }
    
    return 0;
}

 
 

最短Hamilton路径

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

输入格式
第一行输入整数n。

接下来n行每行n个整数，其中第i行第j个整数表示点i到j的距离（记为a[i,j]）。

对于任意的x,y,z，数据保证 a[x,x]=0，a[x,y]=a[y,x] 并且 a[x,y]+a[y,z]>=a[x,z]。

输出格式
输出一个整数，表示最短Hamilton路径的长度。

数据范围 1 ≤ n ≤ 20 ， 0 ≤ a[i,j] ≤ 107

输入样例：
5
0 2 4 5 1
2 0 6 5 3
4 6 0 8 3
5 5 8 0 5
1 3 3 5 0
输出样例：
18

代码
 

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 20, M = 1 << N;
int a[N][N], f[M][N]; // f[i][j] : 从0->j的且状态为i的路径的最小值，i表示所经过的点。

int main()
{
    int n;
    cin >> n;
    
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            cin >> a[i][j];
  
    memset(f, 0x3f, sizeof f);
    f[1][0] = 0; 
    
    for(int i = 1; i < (1 << n); i++) // 枚举每个状态
        for(int j = 0; j < n; j++)
            if((i >> j) & 1) // 0->j的路径必然经过第j个点
                for(int k = 0; k < n; k++)
                    if((i >> k) & 1 && j != k) //状态转移:0->j由0->k, k->j转移而来
                         f[i][j] = min(f[i][j], f[i - (1 << j)][k] + a[k][j]);
    
    cout << f[(1 << n) - 1][n-1];
    return 0;
}