【BZOJ1814】Ural 1519 Formula 1 (插头dp)

题面

BZOJ

Vjudge

题解

戳这里

上面那个链接里面写的非常好啦。

然后说几个点吧。

首先是关于为什么只需要考虑三进制状态,因为哈密顿回路是不可能出现自交的,因此对于当前的轮廓线一定直接分割了哈密顿回路的一部分,不可能出现只考虑分割出来的情况下,存在插头的连通性直接交叉,否则一定不合法(比如说四个连续位置,你不可能\(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)的更多相关文章

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

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

  2. BZOJ1814: Ural 1519 Formula 1(插头Dp)

    Description Regardless of the fact, that Vologda could not get rights to hold the Winter Olympic gam ...

  3. 【BZOJ1814】Ural 1519 Formula 1 插头DP

    [BZOJ1814]Ural 1519 Formula 1 题意:一个 m * n 的棋盘,有的格子存在障碍,求经过所有非障碍格子的哈密顿回路个数.(n,m<=12) 题解:插头DP板子题,刷板 ...

  4. 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 ...

  5. bzoj 1814 Ural 1519 Formula 1 ——插头DP

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1814 普通的插头 DP .但是调了很久.注意如果合并两个 1 的话,不是 “把向右第一个 2 ...

  6. Ural 1519 Formula 1 插头DP

    这是一道经典的插头DP单回路模板题. 用最小表示法来记录连通性,由于二进制的速度,考虑使用8进制. 1.当同时存在左.上插头的时候,需要判断两插头所在连通块是否相同,若相同,只能在最后一个非障碍点相连 ...

  7. bzoj 1814: Ural 1519 Formula 1 插头dp经典题

    用的括号序列,听说比较快. 然并不会预处理,只会每回暴力找匹配的括号. #include<iostream> #include<cstdio> #include<cstr ...

  8. bzoj1814 Ural 1519 Formula 1(插头DP)

    对插头DP的理解还不是很透彻. 先说一下肤浅的理解吧. 插头DP使用范围:指数级复杂度,且适用于解决网格图连通性问题,如哈密顿回路等问题.插头一般指每相邻2个网格的接口. 题目难度:一般不可做. 使用 ...

  9. 插头DP讲解+[BZOJ1814]:Ural 1519 Formula 1(插头DP)

    1.什么是插头$DP$? 插头$DP$是$CDQ$大佬在$2008$年的论文中提出的,是基于状压$D$P的一种更高级的$DP$多用于处理联通问题(路径问题,简单回路问题,多回路问题,广义回路问题,生成 ...

  10. 【Ural】1519. Formula 1 插头DP

    [题目]1519. Formula 1 [题意]给定n*m个方格图,有一些障碍格,求非障碍格的哈密顿回路数量.n,m<=12. [算法]插头DP [题解]<基于连通性状态压缩的动态规划问题 ...

随机推荐

  1. 格式化angularjs日期'/Date(-62135596800000)/'

    在实现在angularjs时,发现经序列化后的日期需要格式化显示. 翻看以前的博客,似乎有写过一篇有关js方面的解决办法<格式化json日期'/Date(-62135596800000)/'&g ...

  2. Java 面试题 == 和 equals 的区别

    int和Integer的区别 1.Integer是int的包装类,int则是java的一种基本数据类型 2.Integer变量必须实例化后才能使用,而int变量不需要 3.Integer实际是对象的引 ...

  3. redis 配置 架构 基础

    redis 官网  redis.io   io为某国家域名后缀 有redis各种版本. java 版本 又分各种工具 clients 下 RedisClient为图形化管理界面 Jedis 才是jav ...

  4. 两个非常好的bootstrap模板,外送大话设计模式!

    两个非常好的bootstrap模板,外送大话设计模式! 下载地址:http://download.csdn.net/download/wolongbb/10198756

  5. Ubuntu轻松编译openJDK

    花了三天在windows上搞openJDK,对bash本来就不熟,加上各种莫名依赖和脚本里的bug,身心俱疲.最后make all的时候产生一个莫名其妙的错误说什么有warning且-Werror置为 ...

  6. Js基础---红宝书读书日记(1)-------基本类型和引用类型

    JS的变量可能包含两种不同数据类型的值,基本类型和引用类型; 基本类型是指简单的数据段,引用类型是指可能由多个值构成的对象; JS高级程序设计第三章介绍了变量分为 5种简单数据类型(string/nu ...

  7. CSS 列表实例

    CSS 列表属性允许你放置.改变列表项标志,或者将图像作为列表项标志.CSS 列表属性(list)属性 描述list-style 简写属性.用于把所有用于列表的属性设置于一个声明中.list-styl ...

  8. mysql操作命令梳理(5)-执行sql语句查询即mysql状态说明

    在日常mysql运维中,经常要查询当前mysql下正在执行的sql语句及其他在跑的mysql相关线程,这就用到mysql processlist这个命令了.mysql> show process ...

  9. oracle系统化学习笔记

    CentOS 4.x上安装Oracle 9i(3讲) oracle9i非常成熟,刚学主要是学9i比较经典 学会安装9i具有现实意义,先学完9i再学11g等比较好 1.安装centos 2.安装orac ...

  10. linux内核分析程序破解实践报告