【学习笔记】插头 DP

合集 - 学习笔记(42)

1.【学习笔记】Kruskal 重构树2022-11-25 2.【学习笔记】网络流2022-12-15 3.【学习笔记】高级数据结构2022-12-21 4.【学习笔记】线性基2022-12-27 5.【学习笔记】Link Cut Tree2022-12-28 6.【学习笔记】字符串后缀算法2022-12-30 7.【学习笔记】字符串回文算法2023-01-12 8.【学习笔记】组合数学2023-01-27 9.【学习笔记】多项式 1：基础操作2023-01-30 10.【学习笔记】多项式 2：集合幂级数2023-02-03 11.【学习笔记】多项式 3：多项式运算2023-02-11 12.【学习笔记】Prufer 序列2023-02-15 13.【学习笔记】多项式 4：生成函数2023-02-15 14.【学习笔记】DP 套 DP2023-02-19 15.【学习笔记】图的连通性2023-02-22 16.【学习笔记】差分约束2023-03-20 17.【学习笔记】长链剖分2023-04-22 18.【学习笔记】2-SAT2023-04-23 19.【学习笔记】根号算法2023-06-06 20.【学习笔记】Primal-Dual 原始对偶算法2023-06-15 21.【学习笔记】Bostan-Mori 算法2023-07-02 22.【学习笔记】狄利克雷卷积与高级筛法2023-07-03 23.【学习笔记】DP 优化 1：基础优化2023-07-03 24.【学习笔记】DP 优化 2：动态 DP2023-07-04 25.【学习笔记】李超线段树2023-07-09 26.【学习笔记】优化建图2023-07-11 27.【学习笔记】Segment Tree Beats2023-07-11

28.【学习笔记】插头 DP2023-07-13

29.【学习笔记】任意模数多项式乘法2023-07-16 30.【学习笔记】SG 函数与 SG 定理2023-08-05 31.【学习笔记】类欧几里得算法2023-08-06 32.【学习笔记】狄利克雷前/后缀和/差分2023-08-11 33.【学习笔记】DSU on Tree2023-08-22 34.【学习笔记】DP 优化 3：闵可夫斯基和优化 DP2023-09-01 35.【学习笔记】笛卡尔树2023-09-05 36.【学习笔记】Miller-Rabin 算法2023-10-24 37.【学习笔记】DP 优化 4：决策单调性2023-11-15 38.【学习笔记】DP 优化 5：wqs 二分优化 DP2023-12-04 39.【学习笔记】边分治2023-12-28 40.【学习笔记】KMP 相关算法2024-01-10 41.【学习笔记】概率生成函数2024-01-11 42.【学习笔记】离散对数和剩余2024-01-25

插头 DP，是一类解决网格图上连通性问题的状压 DP。

括号表示法#

对于连通性问题，需要将当前的插头的状态进行压缩用于转移。

在哈密顿回路问题中，插头一定两两配对，且不存在相交关系，因此可以用括号表示法，将当前路径的左端插头设为 ( 并用 $1$ 表示，右端插头设为 ) 并用 $2$ 表示，不含插头的位置则用 $0$ 表示。

在哈密顿路径问题中，插头可能独立存在，将这样的路径用 $3$ 表示。

由于合法的括号序列个数小于 $4^{m+1}$ ，因此可以使用手写哈希表来减少空间，转移只和上一个方格相关，因此类似滚动数组维护两组信息即可。

状态转移#

以模板题为例。

注意到左插头和上插头（如果存在）转移后分别变成了上插头和左插头（如果存在），且在括号序列的位置不变，转移时只需要考虑当前的两个插头即可。

设当前决策方格的左插头为 $p_1$ ，上插头为 $p_2$ 。

特判当前方格存在障碍的情况，此时 $p_1,p_2$ 必须均为 $0$ ，转移后同样均为 $0$ 。

以下转移均以当前方格以及与当前方格存在插头的方格均不存在障碍为前提。

若 $p_1,p_2$ 均为 $0$ ，转移相当于新建一个连通块，当前方格需要向下向右都连通，转移后 $p_1=1,p_2=2$ 。

若 $p_1$ 不为 $0$ ， $p_2$ 为 $0$ ，转移相当于拓展一个连通块，向下或向右有两种情况。

若 $p_1$ 为 $0$ ， $p_2$ 不为 $0$ ，同理也有两种情况。

若 $p_1=1,p_2=1$ ，转移相当于合并两个连通块，此时 $p_2$ 对应的右括号变成左括号。

若 $p_1=2,p_2=2$ ，同理，此时 $p_1$ 对应的左括号变成右括号。

若 $p_1=2,p_2=1$ ，二者直接合并。

若 $p_1=1,p_2=2$ ，回路闭合。

当切换一行时，注意到左插头从最后一位移至第一位，需要改变状态。

点击查看代码

int n,m;
char mp[15][15];
int ex,ey;
int head[mod+10],nxt[maxn];
int st[2][maxn],cnt[2];
ll dp[2][maxn],ans;
inline void update(int now,int s,ll val){
    int p=s%mod+1;
    for(int i=head[p];i;i=nxt[i]){
        if(st[now][i]==s) return dp[now][i]+=val,void();
    }
    st[now][++cnt[now]]=s,dp[now][cnt[now]]=val;
    nxt[cnt[now]]=head[p],head[p]=cnt[now];
}
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;++i){
        scanf("%s",mp[i]+1);
        for(int j=1;j<=m;++j){
            if(mp[i][j]=='.') ex=i,ey=j;
        }
    }
    int now=0;
    update(now,0,1);
    for(int i=1;i<=n;++i){
        for(int k=1;k<=cnt[now];++k) st[now][k]<<=2;
        for(int j=1;j<=m;++j){
            memset(head,0,sizeof(head));
            now^=1;
            cnt[now]=0;
            for(int k=1;k<=cnt[now^1];++k){
                int s=st[now^1][k];
                ll val=dp[now^1][k];
                int p1=s>>2*(j-1)&3,p2=s>>2*j&3;
                int t=s-(p1<<2*(j-1))-(p2<<2*j);
                if(mp[i][j]=='*'){
                    if(!p1&&!p2) update(now,s,val);
                }
                else{
                    if(!p1&&!p2){
                        if(mp[i+1][j]=='.'&&mp[i][j+1]=='.') update(now,t+(1<<2*(j-1))+(2<<2*j),val);
                    }
                    else if(p1&&!p2){
                        if(mp[i+1][j]=='.') update(now,s,val);
                        if(mp[i][j+1]=='.') update(now,t+(p1<<2*j),val);
                    }
                    else if(!p1&&p2){
                        if(mp[i+1][j]=='.') update(now,t+(p2<<2*(j-1)),val);
                        if(mp[i][j+1]=='.') update(now,s,val);
                    }
                    else{
                        if(p1==1&&p2==1){
                            int tmp=1;
                            for(int l=j+1;l<=m;++l){
                                if((s>>2*l&3)==1) ++tmp;
                                if((s>>2*l&3)==2) --tmp;
                                if(!tmp){
                                    update(now,t-(1<<2*l),val);
                                    break;
                                }
                            }
                        }
                        else if(p1==2&&p2==2){
                            int tmp=1;
                            for(int l=j-2;l>=0;--l){
                                if((s>>2*l&3)==2) ++tmp;
                                if((s>>2*l&3)==1) --tmp;
                                if(!tmp){
                                    update(now,t+(1<<2*l),val);
                                    break;
                                }
                            }
                        }
                        else if(p1==2&&p2==1) update(now,t,val);
                        else{
                            if(i==ex&&j==ey) ans+=val;
                        }
                    }
                }
            }
        }
    }
    printf("%lld\n",ans);
    return 0;
}