【BZOJ1814】Ural 1519 Formula 1 (插头dp)
【BZOJ1814】Ural 1519 Formula 1 (插头dp)
题面
题解
戳这里
上面那个链接里面写的非常好啦。
然后说几个点吧。
首先是关于为什么只需要考虑三进制状态,因为哈密顿回路是不可能出现自交的,因此对于当前的轮廓线一定直接分割了哈密顿回路的一部分,不可能出现只考虑分割出来的情况下,存在插头的连通性直接交叉,否则一定不合法(比如说四个连续位置,你不可能\(1,3\)匹配,因为这样子画一条路径出来,无论如何都会和\(2,4\)的路径相交)。因此,我们把一组匹配看成左右括号,再加上有的位置可能不存在插头,所以是左括号、右括号、无插头三种情况。
然后是关于转移中的一部分转移,当左侧和上方的边都有插头连进来的时候,此时当前格子内的插头就固定了,同时转移的时候就得到了两个无插头的位置。注意一下,如果当前位置分别是\(1,1\)或者\(2,2\)的时候,你还需要修改其他位置的插头使他们满足左右匹配的关系,举个例子,本来是\(1122\)的匹配关系,你现在把\(22\)给连在一起了,于是这个状态就变成了\(1100\),所以你需要把第二个\(1\)变成\(2\),也就是这个状态应该是\(1200\)。
我的程序直接刚的三进制的,自己本机测没有问题,但是常数比较大,交到\(bzoj\)上就直接\(TLE\)了。可能直接改成四进制用位运算反而会快些。懒得改了,就直接贴代码吧。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
#define ll long long
#define MAX 15
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int n,m,ln,lm,g[MAX][MAX];
char ch[MAX];
ll f[13][13][42000];
int bin[MAX];
int zt[42000],q[1600000],tot;
ll ans;
void dfs(int x,int c,int S)
{
if(c<0||c>m-x+1)return;
if(x>m){zt[++tot]=S;q[S]=tot;return;}
dfs(x+1,c,S);
dfs(x+1,c+1,S+bin[x]);
dfs(x+1,c-1,S+bin[x]+bin[x]);
}
int nxt(int x,int p)
{
for(int i=p,c=0;i<=m;++i)
{
if(x/bin[i]%3==2)--c;
if(x/bin[i]%3==1)++c;
if(!c)return i;
}
return m+1;
}
int pre(int x,int p)
{
for(int i=p,c=0;~i;--i)
{
if(x/bin[i]%3==1)--c;
if(x/bin[i]%3==2)++c;
if(!c)return i;
}
return m+1;
}
int main()
{
n=read();m=read();
for(int i=1;i<=n;++i)
{
scanf("%s",ch+1);
for(int j=1;j<=m;++j)
if(ch[j]=='.')g[i][j]=1,ln=i,lm=j;
}
bin[0]=1;for(int i=1;i<=m;++i)bin[i]=bin[i-1]*3;
dfs(0,0,0);f[1][0][q[0]]=1;
for(int i=1;i<=n;++i)
{
for(int j=1;j<=m;++j)
for(int k=1;k<=tot;++k)
{
int left=zt[k]/bin[j-1]%3,up=zt[k]/bin[j]%3;
if(!g[i][j]){if(!left&&!up)f[i][j][k]+=f[i][j-1][k];}
else
{
if(!left&&!up&&g[i][j+1]&&g[i+1][j])
f[i][j][q[zt[k]+bin[j-1]+2*bin[j]]]+=f[i][j-1][k];
if(!left&&up)
{
if(g[i][j+1])f[i][j][k]+=f[i][j-1][k];
if(g[i+1][j])f[i][j][q[zt[k]+up*(bin[j-1]-bin[j])]]+=f[i][j-1][k];
}
if(left&&!up)
{
if(g[i+1][j])f[i][j][k]+=f[i][j-1][k];
if(g[i][j+1])f[i][j][q[zt[k]+left*(bin[j]-bin[j-1])]]+=f[i][j-1][k];
}
if(left==1&&up==1)f[i][j][q[zt[k]-bin[j-1]-bin[j]-bin[nxt(zt[k],j)]]]+=f[i][j-1][k];
if(left==2&&up==1)f[i][j][q[zt[k]-2*bin[j-1]-bin[j]]]+=f[i][j-1][k];
if(left==2&&up==2)f[i][j][q[zt[k]-2*bin[j-1]-2*bin[j]+bin[pre(zt[k],j-1)]]]+=f[i][j-1][k];
if(left==1&&up==2&&i==ln&&j==lm)ans+=f[i][j-1][k];
}
}
for(int j=1;j<=tot;++j)
if(zt[j]%3==0)f[i+1][0][j]+=f[i][m][q[zt[j]/3]];
}
printf("%lld\n",ans);
return 0;
}
【BZOJ1814】Ural 1519 Formula 1 (插头dp)的更多相关文章
- bzoj1814 Ural 1519 Formula 1(插头dp模板题)
1814: Ural 1519 Formula 1 Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 924 Solved: 351[Submit][Sta ...
- BZOJ1814: Ural 1519 Formula 1(插头Dp)
Description Regardless of the fact, that Vologda could not get rights to hold the Winter Olympic gam ...
- 【BZOJ1814】Ural 1519 Formula 1 插头DP
[BZOJ1814]Ural 1519 Formula 1 题意:一个 m * n 的棋盘,有的格子存在障碍,求经过所有非障碍格子的哈密顿回路个数.(n,m<=12) 题解:插头DP板子题,刷板 ...
- bzoj 1814 Ural 1519 Formula 1 插头DP
1814: Ural 1519 Formula 1 Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 942 Solved: 356[Submit][Sta ...
- bzoj 1814 Ural 1519 Formula 1 ——插头DP
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1814 普通的插头 DP .但是调了很久.注意如果合并两个 1 的话,不是 “把向右第一个 2 ...
- Ural 1519 Formula 1 插头DP
这是一道经典的插头DP单回路模板题. 用最小表示法来记录连通性,由于二进制的速度,考虑使用8进制. 1.当同时存在左.上插头的时候,需要判断两插头所在连通块是否相同,若相同,只能在最后一个非障碍点相连 ...
- bzoj 1814: Ural 1519 Formula 1 插头dp经典题
用的括号序列,听说比较快. 然并不会预处理,只会每回暴力找匹配的括号. #include<iostream> #include<cstdio> #include<cstr ...
- bzoj1814 Ural 1519 Formula 1(插头DP)
对插头DP的理解还不是很透彻. 先说一下肤浅的理解吧. 插头DP使用范围:指数级复杂度,且适用于解决网格图连通性问题,如哈密顿回路等问题.插头一般指每相邻2个网格的接口. 题目难度:一般不可做. 使用 ...
- 插头DP讲解+[BZOJ1814]:Ural 1519 Formula 1(插头DP)
1.什么是插头$DP$? 插头$DP$是$CDQ$大佬在$2008$年的论文中提出的,是基于状压$D$P的一种更高级的$DP$多用于处理联通问题(路径问题,简单回路问题,多回路问题,广义回路问题,生成 ...
- 【Ural】1519. Formula 1 插头DP
[题目]1519. Formula 1 [题意]给定n*m个方格图,有一些障碍格,求非障碍格的哈密顿回路数量.n,m<=12. [算法]插头DP [题解]<基于连通性状态压缩的动态规划问题 ...
随机推荐
- (一)在 Blend 中绘制形状和路径
原文:(一)在 Blend 中绘制形状和路径 https://docs.microsoft.com/zh-cn/previous-versions/jj170881(v=vs.120) 在 Blend ...
- 手机端@media的屏幕适配
@media only screen and (width: 320px) { html { font-size: 16px; }} @media only screen and (width: 36 ...
- (7)学习笔记 ) ASP.NET CORE微服务 Micro-Service ---- 利用Polly+AOP+依赖注入封装的降级框架
创建简单的熔断降级框架 要达到的目标是: 参与降级的方法参数要一样,当HelloAsync执行出错的时候执行HelloFallBackAsync方法. public class Person { [H ...
- Shell编程基础篇-上
1.1 前言 1.1.1 为什么学Shell Shell脚本语言是实现Linux/UNIX系统管理及自动化运维所必备的重要工具, Linux/UNIX系统的底层及基础应用软件的核心大都涉及Shell脚 ...
- [UWP 自定义控件]了解模板化控件(10):原则与技巧
1. 原则 推荐以符合以下原则的方式编写模板化控件: 选择合适的父类:选择合适的父类可以节省大量的工作,从UWP自带的控件中选择父类是最安全的做法,通常的选择是Control.ContentContr ...
- ant+Jacoco 统计tomcat远程部署后项目接口自动化测试或者功能测试代码覆盖率
1.安装ant 环境,https://ant.apache.org/bindownload.cgi 2.下载jacoco包 https://www.eclemma.org/jacoco/ ,解压后, ...
- json中获取key值
<script type="text/javascript"> getJson('age'); function getJson(key){ var jsonObj={ ...
- Final 个人最终作业。
1.对软件工程M1/M2做一个总结 在M1阶段,我在C705组.M1阶段我与黄漠源同学结对,一起完成提取关键词算法的优化.最初我们一起测试提取关键词算法功能的实现效果,随后我主要负责从网络上搜寻并整理 ...
- 个人博客作业Week7(阅读文章,心得体会)
Alpha阶段结束了,内心可以说是五味杂陈.不是说我们的产品拿不上台面那般差劲,复杂的心绪主要来源于和别的队的比较,别的队才刚刚发布没多久访问量和注册量就破百了,并且还发起了找bug送红包的活动.可能 ...
- Linux内核分析——进程的描述和进程的创建
进程的描述和进程的创建 一. 进程的描述 (一)进程控制块PCB——task_struct 1.操作系统的三大管理功能包括: (1)进程管理 (2)内存管理 (3)文件系统 2.PCB task_st ...