状态压缩dp、轮廓线、插头dp——从入门到不会

题目清单

• POJ1185 炮兵布阵（状态压缩dp）
• HDU1693 闭合线路统计（插头dp）
• POJ2411 平面骨牌密铺（状态压缩dp/轮廓线更新）
• HDU1565 矩阵选数（轮廓线更新）
• HDU2167 矩阵选数（比上一题更严格的条件）（轮廓线更新）

POJ1185 该题是状态压缩dp的入门题，通常的做法是用一个01串表示行上某个位置是否安排炮兵。在状态转移时遍历前两行的状态并安排炮兵向下一行转移：对于每一行进行二重循环，分别遍历两行的状态，并进行检测转移合法性。然而考虑到状态数过多（2^10），实际循环遍历轮数是它的平方。所以我们对状态进行预处理，事先将不合法的状态（行内安排间隔小于2）剔除，重新编码，能显著减少状态数（在列数<10的情况下，状态数小于65）。

#include 
#include 
#include 
using namespace std;
const int maxn=1024;
int pro[maxn], rev[65], tot;
int dp[2][65][65];
int mp[110];
int n, m;

inline bool judge(int x) { return (x&(x>>1))||(x&(x>>2)); }
// 检查合法状态
void init() {
    tot=0;
    for (int i=0;i<(1<=0;--j) {
            mp[i]<<=1;
            mp[i]+=(tmp[j]=='H');
        }
    }

    memset(dp, 0xff, sizeof(dp));
    dp[0][0][0]=0;
    for (int i=1;i<=n;++i) {
        if (i>1) memset(dp[i%2], 0xff, sizeof(dp[i%2]));
        for (int u=0;u2&&(last&mp[i-2])) continue;
            if (judge(last)) continue;
            for (int v=0;v1&&(cur&mp[i-1])) continue;
                if (judge(cur)) continue;
                for (int w=0;w

POJ1693 在一个相对较小的矩阵中有些位置不可以走，现在只能上、下、左、右访问合法位置，走过的路不可以再走。现在在遵守以上约定的情况下，在矩阵中安排回路，求总的安排种类数。虽然题目中要求每个回路中必须包含一个或多个禁止通行的方块，但是完全没考虑这样的情况提交之后也ac了……

我们考虑这样的状态：对于某个格子u，它所处的行在它左侧的所有格子和在它上方的所有格子已经被统计过，在统计过与未被统计过的格子中间划一条“边缘线”（轮廓线），沿着这个轮廓线，被统计过的格子是否与未被统计的格子之间存在“插头”，将这个记作state，它可以用一个01串表示。所谓“插头”就是是否有跨过轮廓线的路径。dp[u][state]表示统计到u时轮廓线上插头状态为state的时候，总的路径安排种类数。

插头dp01.png

#include 
#include 
#include 
using namespace std;
typedef long long ll;

const int N=13;
int mp[N][N], n, m;
ll dp[2][1<11
                                dp[cur][st|(3<>j)&3)!=3) {              // 01/10->01/10
                                dp[cur][st^(3<>j)&3)==1||((st>>j)&3)==2)     //01 10->10
                                dp[cur][((st|(1<>j)&3)==3)                //11->00
                                dp[cur][(st^(3<

POJ2411 平面骨牌密铺，如果选用和插头dp类似的轮廓线转移方式做，状态转移会比较方便，因为状态转移到另一种状态的种类数有限。此次记录轮廓线下方（与上题不同，轮廓线上有段竖直的线状态不计）的方块是否已经被占用，进行状态转移。大体与上一题类似。

#include 
#include 
//#include 
#include 
using namespace std;
const int N=12;
typedef long long ll;
ll dp[2][1<

POJ1565 给定的矩阵中有一些数，我们可以选择一些数出来。现在规定不可以选择相邻的数（上下左右方向有公共边），问一次可以选的数的和最大为多少。

这题的转移方式和上一题“骨牌密铺”有点类似，同样使用轮廓线转移来写，记录的是轮廓线下方（注意是下方而不是下方和右方）方格是否可选的状态。

#include 
#include 
#include 
using namespace std;
const int N=21;
int a[N][N], dp[2][1<

POJ2167 与上题类似，不过这道题将“相邻”定义为周围八个格子。这时我们就不得不记录更多信息：对于宽度为n的图的轮廓线，我们将轮廓线下方和右方共n+1个方格是否可以选择的状态记录下来，并进行状态转移和计数。

#include 
#include 
#include 
using namespace std;
const int N=16;
int a[N][N], n, dp[2][1<0?nst|(1<<(j-1)):nst;
                            nst&=((2<=0){
                        ans=max(ans,dp[cur][st]);
                    }
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}

未完待续

最小表示法