\(\color{#0066ff}{ 题目描述 }\)

给出n*m的方格,有些格子不能铺线,其它格子必须铺,形成一个闭合回路。问有多少种铺法?

\(\color{#0066ff}{输入格式}\)

第1行,n,m(2<=n,m<=12)

从第2行到第n+1行,每行一段字符串(m个字符),"*"表不能铺线,"."表必须铺

\(\color{#0066ff}{输出格式}\)

输出一个整数,表示总方案数

\(\color{#0066ff}{输入样例}\)

4 4
**..
....
....
....

\(\color{#0066ff}{输出样例}\)

2

\(\color{#0066ff}{数据范围与提示}\)

none

\(\color{#0066ff}{ 题解 }\)

插头DP本来以为多niubility的算法原来本质还是个DP,就是情况多了点qwq

状压分割线,有三种情况,无,左插头,右插头(详见洛谷题解懒得写了)

只要明白状态分析出所有情况即可

#include<bits/stdc++.h>
#define LL long long
LL in() {
char ch; LL x = 0, f = 1;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
return x * f;
}
std::unordered_map<LL, LL> f[2];
//我TM就不写hash表
int n, m, s, t;
bool mp[100][100];
char getch() {
char ch = getchar();
while(ch != '*' && ch != '.') ch = getchar();
return ch;
}
void init() {
n = in(), m = in();
for(int i = 1; i <= n; i++)
for(int j = 1; j <= m; j++) {
mp[i][j] = getch() == '.';
if(mp[i][j]) s = i, t = j;
}
}
LL pos(int v, int x) { return (v << (x << 1)); }
//返回第x大块的v状态(两个二进制来状压)
LL work() {
int now = 0, nxt = 1;
f[0][0] = 1;
LL U = (1LL << ((m + 1) << 1)) - 1;
//全集
for(int i = 1; i <= n; i++) {
for(int j = 1; j <= m; j++) {
f[nxt].clear();
for(auto &k:f[now]) {
LL S = k.first, val = k.second;
LL L = (S >> ((j - 1) << 1)) & 3, R = (S >> (j << 1)) & 3;
//分割线(L是当前竖着的那个,R是紧接着横着的那个)
if(!mp[i][j]) {
if(!L && !R) f[nxt][S] += val;
continue;
}
// 0 0
if(!L && !R) {
if(mp[i][j + 1] && mp[i + 1][j]) f[nxt][S ^ pos(1, j - 1) ^ pos(2, j)] += val;
//0 0 -> 1 2
}
//2 1
else if(L == 2 && R == 1) {
//2 1 -> 0 0
f[nxt][S ^ pos(L, j - 1) ^ pos(R, j)] += val;
}
// 0 1 // 1 0 // 0 2 // 2 0
else if(!L || !R) {
//0 1 // 0 2
if(!L) {
//拐弯
if(mp[i][j + 1]) f[nxt][S] += val;
//不拐弯
if(mp[i + 1][j]) f[nxt][S ^ pos(L, j - 1) ^ pos(L, j) ^ pos(R, j - 1) ^ pos(R, j)] += val;
}
//同上
else if(!R) {
if(mp[i][j + 1]) f[nxt][S ^ pos(L, j - 1) ^ pos(L, j) ^ pos(R, j - 1) ^ pos(R, j)] += val;
if(mp[i + 1][j]) f[nxt][S] += val;
}
}
//1 1 // 2 2
else if(L == R) {
// 1 1
if(L == 1) {
int du = 0;
//1 1 -> 0 0 但是源头接口处右插头变成左插头
for(int p = j; ; p++) {
LL o = (S >> (p << 1)) & 3;
if(o == 1) du++;
if(o == 2) du--;
if(!du) {
//原来状态消去,弄上新状态
f[nxt][S ^ pos(L, j - 1) ^ pos(R, j) ^ pos(2, p) ^ pos(1, p)] += val;
break;
}
}
}
//根上面差不多
else if(L == 2) {
int du = 0;
for(int p = j - 1; ; p--) {
LL o = (S >> (p << 1)) & 3;
if(o == 1) du--;
if(o == 2) du++;
if(!du) {
f[nxt][S ^ pos(L, j - 1) ^ pos(R, j) ^ pos(2, p) ^ pos(1, p)] += val;
break;
}
}
}
}
//本状态当且仅当是终点,用于封口
else if(L == 1 && R == 2 && i == s && j == t) return val;
}
std::swap(now, nxt);
}
f[nxt].clear();
//末尾的竖分割线到了下一行就变成了行首的分割线,把状态《《给分割线腾出地方
for(auto &k:f[now]) f[nxt][(k.first << 2) & U] += k.second;
std::swap(now, nxt);
}
return 0;
}
int main() {
init();
printf("%lld\n", work());
return 0;
}

P5056 【模板】插头dp的更多相关文章

  1. 模板—插头dp(Ural 1519 Formula 1)

    括号表示法: 据说比下一个要快而且灵活. #include<iostream> #include<cstring> #include<cstdio> #define ...

  2. 插头DP模板

    /* 插头dp模板 抄的GNAQ 的 括号表示法 */ #include<cstdio> #include<algorithm> #include<cstring> ...

  3. 模板:插头dp

    前言: 严格来讲有关dp的都不应该叫做模板,因为dp太活了,但是一是为了整理插头dp的知识,二是插头dp有良好的套路性,所以姑且还叫做模板吧. 这里先推荐一波CDQ的论文和这篇博客http://www ...

  4. bzoj1814 Ural 1519 Formula 1(插头dp模板题)

    1814: Ural 1519 Formula 1 Time Limit: 1 Sec  Memory Limit: 64 MBSubmit: 924  Solved: 351[Submit][Sta ...

  5. LG5056 【模板】插头dp

    题意 题目背景 ural 1519 陈丹琦<基于连通性状态压缩的动态规划问题>中的例题 题目描述 给出n*m的方格,有些格子不能铺线,其它格子必须铺,形成一个闭合回路.问有多少种铺法? 输 ...

  6. 【模板】插头dp

    题目描述 题解: 插头$dp$中经典的回路问题. 首先了解一下插头. 一个格子,上下左右四条边对应四个插头.就像这样: 四个插头. 一个完整的哈密顿回路,经过的格子一定用且仅用了两个插头. 所以所有被 ...

  7. [学习笔记]插头dp

    基于连通性的状压dp 巧妙之处:插头已经可以表示内部所有状态了. 就是讨论麻烦一些. 简介 转移方法:逐格转移,分类讨论 记录状态方法:最小表示法(每次要重新编号,对于一类没用“回路路径”之类的题,可 ...

  8. 插头dp小结

    插头dp: \(A:\)插头dp是什么? \(B:\)一种基于连通性状态压缩的动态规划问题 \(A:\)请问有什么应用呢? \(B:\)各种网格覆盖问题,范围允许状压解决,常用于计算方案数与联通块权值 ...

  9. 插头dp

    插头dp 感受: 我觉得重点是理解,算法并不是直接想出怎样由一种方案变成另一种方案.而是方案本来就在那里,我们只是枚举状态统计了答案. 看看cdq的讲义什么的,一开始可能觉得状态很多,但其实灰常简单 ...

随机推荐

  1. Java 中的关键字和保留字

    关键字: Java 语言中已经事先定义好了的,有着特殊含义和用途 访问控制 类.方法和变量修饰符 程序控制 异常处理 包相关 基本类型 变量引用 public abstract break try i ...

  2. vBulletin 5.x 版本通杀远程代码执行漏洞复现

    漏洞介绍 vBulletin中存在一个文件包含问题,可使恶意访问者包含来自 vBulletin 服务器的文件并且执行任意 PHP 代码.未经验证的恶意访问者可通过向index.php发出包含 rout ...

  3. Dev控件GridControl实现CheckBox列和ComBox列

    1.在sql语句中添加空白行,如select c1,c2 null c3 from xxx; 2.将sql语句查询结果与gdc绑定CmmFrm.BestFitGridViewColumnsWidth( ...

  4. Windows安装mysql 5.7.*.zip步骤

    1.去官网上下载.zip格式的文件. 2.解压到一个文件夹,这里我用D:\MySql表示 3.在D:\MySql\mysql-5.7.17-winx64下新建my.ini配置文件 黄色背景色的地方需要 ...

  5. 10-20C#基础---一维、二维数组&&冒泡排序

    一.一维数组 1.定义:是某一种数据类型的数据的组合,数组用来分组基本类型或相同类型的对象.数组中的实体叫做数组的元素或成员. 2. 格式:int[ ] shuzu=new int[ 6];存放int ...

  6. DAY10-MYSQL库操作

    一 系统数据库 information_schema: 虚拟库,不占用磁盘空间,存储的是数据库启动后的一些参数,如用户表信息.列信息.权限信息.字符信息等performance_schema: MyS ...

  7. C#WinForm如何调整控件的Tab按键顺序

    在日常生活中,很多用户都会有使用Tab键的习惯.而在C#的WinForm开发中,Tab按键的顺序默认是你拖拽进窗体的顺序.那么我们如何修改这个顺序呢?答案如下(以VS2010为例). 只需要点击[视图 ...

  8. [hdu1251]统计难题(trie模板题)

    题意:返回字典中所有以测试串为前缀的字符串总数. 解题关键:trie模板题,由AC自动机的板子稍加改造而来. #include<cstdio> #include<cstring> ...

  9. ArcGIS10拓扑规则-面规则(转)

    ArcGIS10拓扑规则-面规则 原创 2013年12月27日 10:20:44 标签: ArcGIS 1879 ARCGIS 10 里提供的拓扑规则共32种,下面一一介绍: 首先介绍的对于面图层的拓 ...

  10. String/StringBuilder 类 判断QQ号码

    1.1. 训练描述:[方法.String类]  一.需求说明:请用户输入一个“QQ号码”,我们来判断这个QQ号码是否正确. 要求:使用方法来完成判断功能. 1.2. 操作步骤描述 建立MainApp类 ...