垃圾插头DP,照着打都调了我一下午,淦!!!
学这个玩意纯粹是因为模拟赛考了一道,要不然碰都不会碰......
我觉得插头DP的主要难度在于实现,而不是理解算法原理......
不说废话了,进入正题。
解决的问题
学一个新知识点,肯定先要了解这个知识点能够解决什么样的问题。(废话)
其实插头DP就是一类特殊的状压DP,他解决的就是在小数据规模下连通性的问题,数据太大插头DP也不好做,毕竟本质还是状压。
举几个简单的例子,就比如说格点图的哈密顿路径计数,求棋盘的黑白染色方案满足相同颜色之间形成一个连通块的方案数,以及特定图的生成树计数等等,这些你看起来连暴力都不是很好搞的东西,插头DP就可以做。(说实话,就这复杂度,说暴力都有人信)
接下来我们先来了解几个概念。
基本概念
插头
如果一个格子某个方向有插头,表示这个格子在这个方向与相邻格子是相连的。
简单点来说就是记录当前格子与相邻的格子是否连通。
轮廓线
就是已决策格子与未决策格子的分界线。(下图红色的线)
注意一下,这里分界线的边数是列数加一,别搞错了。
大概就这两个概念,我们现在来看看插头DP如何实现
实现
我们就借一道板题来说。
因为数据较小,但是搜索过不了,于是考虑状压。(即插头DP)
可以逐格DP(某些题目可以逐行),设 f i , j , s t a t e f_{i,j,state} fi,j,state 表示当前的方案数。
通常的 s t a t e state state 有括号表示和最小表示,这里着重介绍泛用性更好的最小表示。我们用长度为 m + 1 m+1 m+1 的整形数组,记录轮廓线上每个插头的状态, 0 0 0 表示没有插头,并约定连通的插头用相同的数字进行标记。
由于要形成最小回路,可知最多只有两个插头互相连通,并且跟两组连通的插头之间不会交叉,既不会出现 2 1 2 1 2\ 1\ 2\ 1 2 1 2 1 的情况,并且每一个格子最多只有两个插头。
考虑如何转移。
对于上图黄色的格子,轮廓线为红色, 蓝色使我们DP完当前节点后的轮廓线(包括第一、二、四列的红线),我们分以下几种情况考虑。
- 跟上面和左边的格子都有插头,此时还要分两种情况。
- 如果这两个插头没连通,则这个格子则不能有右插头和下插头,即新的两边轮廓线上没有任何插头。
- 如果这两个插头连通了,若当前枚举的格子为最后一个格子,则直接连通,否则不能连通。(若连通则成为多个闭合回路了)
- 跟上面和左边的格子有一个插头,则此时跟下面和右边任意一个格子也有一个插头。
- 上面和左边都没有插头,则此时右边和下面都必须有插头。
注意考虑一下边界和障碍的条件。
当我们DP完一行时,轮廓线会变为下图形式:
此时我们在DP下一行之前记得将整个状态左移即可。
优化
对于普通的插头DP,时空复杂度也并不是非常优,此时就需要一点优化。
我们会发现由于他是闭合回路,它的状态会非常稀疏,于是我们可以用 hash 来优化。
但这个 hash 要用挂表法,即将你所有模数为 x x x 的状态都挂到一起,可以用链式前向星维护。
你每次遇到一个新状态就将它存到哈希里面,然后查找状态则将那一条 hash 链枚举一遍即可。
代码
cpp
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 200005, P = 9973, offset = 4, mask = (1 << offset) - 1;
int n, m, mp[20][20], ex, ey;
LL Val;
int in() {
char ch = getchar();
while (ch != '.' && ch != '*')
ch = getchar();
return ch == '*';
}
struct HashTable {//hash的实现
int last[P], next[N], cnt;
LL ans[N], state[N];
void clear() {
cnt = 0;
memset(last, 0, sizeof(last));
}
void push(LL s) {
LL x = s % P;
for (int i = last[x]; i; i = next[i])
if (state[i] == s) {
ans[i] += Val;
return ;
}
cnt++, state[cnt] = s, ans[cnt] = Val;
next[cnt] = last[x], last[x] = cnt;
}
void roll() {
for (int i = 1; i <= cnt; i++)
state[i] = state[i] << offset;
}
} H[2], *H0, *H1;
int bb[25], b[25];
LL encode() {//这是将数组加码成状态
LL s = 0;
memset(bb, -1, sizeof(bb));
int bn = 1;
bb[0] = 0;
for (int i = m + 1; i; i--) {
if (!~bb[b[i]])
bb[b[i]] = bn++;
s <<= offset;
s |= bb[b[i]];
}
return s;
}
void decode(LL s) {//这是将状态解码成数组
for (int i = 1; i <= m + 1; i++) {
b[i] = s & mask;
s >>= offset;
}
}
void push(int j, int down, int right) {//将状态改变,然后往下DP
b[j] = down, b[j + 1] = right;
H1 -> push(encode());
}
void PlugDP() {//DP
H0 = H, H1 = H + 1;
H1 -> clear(), Val = 1, H1 -> push(0);
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
swap(H0, H1), H1 -> clear();
for (int l = 1; l <= (H0 -> cnt); l++) {
decode(H0 -> state[l]);
Val = H0 -> ans[l];
int left = b[j], up = b[j + 1];
//上面和左边的插头类型
bool down = i != n, right = j != m;;
//这两个变量即是判断是否处于边界位置,能否有右插头和下插头
if (mp[i][j] == 1) {//如果当前格子为障碍
if (!left && !up) {//如果上面和左边都没有插头,则可以转移成下面和右边也都没有插头
push(j, 0, 0);
}
}
else if (left && up) {//上面和左边都有插头
if (left == up) {//两个插头连通
if (i == ex && j == ey)
push(j, 0, 0);
}
else {//两个插头不连通
for (int t = 1; t <= m + 1; t++)
if (b[t] == left)
b[t] = up;
push(j, 0, 0);
}
}
else if (left || up) {//上面和左边只有一个插头
int t = left | up;
if (down)
push(j, t, 0);
if (right)
push(j, 0, t);
}
else if (down && right)//两边都没有插头并且下面和右边都可以有插头
push(j, m, m);
}
}
H1 -> roll();//即全部右移
}
}
int main() {
scanf("%d%d", &n, &m);
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++) {
mp[i][j] = in();
if (!mp[i][j])
ex = i, ey = j;
}
PlugDP();
cout << (H1 -> cnt == 1 ? H1 -> ans[1] : 0) << endl;
return 0;
}