插头 DP

垃圾插头DP,照着打都调了我一下午,淦!!!

学这个玩意纯粹是因为模拟赛考了一道,要不然碰都不会碰......

我觉得插头DP的主要难度在于实现,而不是理解算法原理......

不说废话了,进入正题。

解决的问题

学一个新知识点,肯定先要了解这个知识点能够解决什么样的问题。(废话)

其实插头DP就是一类特殊的状压DP,他解决的就是在小数据规模下连通性的问题,数据太大插头DP也不好做,毕竟本质还是状压。

举几个简单的例子,就比如说格点图的哈密顿路径计数,求棋盘的黑白染色方案满足相同颜色之间形成一个连通块的方案数,以及特定图的生成树计数等等,这些你看起来连暴力都不是很好搞的东西,插头DP就可以做。(说实话,就这复杂度,说暴力都有人信)

接下来我们先来了解几个概念。

基本概念


插头

如果一个格子某个方向有插头,表示这个格子在这个方向与相邻格子是相连的。

简单点来说就是记录当前格子与相邻的格子是否连通。

轮廓线

就是已决策格子与未决策格子的分界线。(下图红色的线)

注意一下,这里分界线的边数是列数加一,别搞错了。


大概就这两个概念,我们现在来看看插头DP如何实现

实现

我们就借一道板题来说。

Luogu P5056 【模板】插头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;
}
相关推荐
lulu_gh_yu29 分钟前
数据结构之排序补充
c语言·开发语言·数据结构·c++·学习·算法·排序算法
ULTRA??1 小时前
C加加中的结构化绑定(解包,折叠展开)
开发语言·c++
凌云行者2 小时前
OpenGL入门005——使用Shader类管理着色器
c++·cmake·opengl
凌云行者2 小时前
OpenGL入门006——着色器在纹理混合中的应用
c++·cmake·opengl
~yY…s<#>2 小时前
【刷题17】最小栈、栈的压入弹出、逆波兰表达式
c语言·数据结构·c++·算法·leetcode
可均可可3 小时前
C++之OpenCV入门到提高004:Mat 对象的使用
c++·opencv·mat·imread·imwrite
白子寰3 小时前
【C++打怪之路Lv14】- “多态“篇
开发语言·c++
yannan201903133 小时前
【算法】(Python)动态规划
python·算法·动态规划
小芒果_013 小时前
P11229 [CSP-J 2024] 小木棍
c++·算法·信息学奥赛
gkdpjj3 小时前
C++优选算法十 哈希表
c++·算法·散列表