插头dp板子题??

搞了我一晚上,还tm全是抄的标程。。

还有高精,哈希混入,还是我比较弱,orz各种dalao

有不明白的可以去看原论文。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define base (int)1e9
#define maxh 2650
using namespace std;
int n,m;
int get(int s,int x){
return (s>>((x-1)<<1))&3;
}//查询对于状态s,第x位的状态. 三进制,所以左移一位.&3,得0,1,2
void change(int &s,int x,int y){
int w=(x-1)<<1;
s&=~(3<<w);
s|=y<<w;
}//将s状态的第x位改为y
int find(int s,int x){
for(int i=x,d=0,pos=get(s,x)==1?1:-1;;i+=pos){
int nn=get(s,i);
if(nn==1) d++; if(nn==2) d--;
if(!d) return i;
}
}//找与该位匹配的括号
struct bignum{//高精没啥好说的...
int a[5];
void clear(){memset(a,0,sizeof a);}
bignum(){clear();}
void set(int x){clear();while(x){a[++a[0]]=x%base;x/=base;}}
void print(){
printf("%d",a[a[0]]);
for(int i=a[0]-1;i>0;i--)printf("%09d",a[i]);
}
bignum operator + (bignum b){
static bignum c;c.clear();c.a[0]=max(a[0],b.a[0])+1;
for(int i=1;i<=c.a[0];i++){
c.a[i]+=a[i]+b.a[i];
c.a[i+1]+=c.a[i]/base;
c.a[i]%=base;
}
while(!c.a[c.a[0]])c.a[0]--;
return c;
}
void operator += (bignum b){*this=*this+b;}
void operator = (int x){set(x);}
}ans;
struct hashtable{//骚骚的哈希
int key[maxh];
int cnt,hash[maxh];
bignum val[maxh];
void clear(){
memset(key,-1,sizeof key);
memset(hash,0,sizeof hash);
memset(val,0,sizeof val);
cnt=0;
}
void newnode(int i,int x){
hash[i]=++cnt;
key[cnt]=x;
}
bignum & operator[](const int hh){
int k;
for(int i=hh%maxh;;i++){
if(i==maxh)i=0;
if(!hash[i]) newnode(i,hh);
if(key[hash[i]]==hh){k=hash[i]; break;}
}
return val[k];
}
}f[2];
int now,last;
void dp(int i,int j){//插头dp
now^=1;last^=1;
int tot=f[last].cnt;
f[now].clear();
for(int k=1,x,y;k<=tot;k++){
int s=f[last].key[k];
bignum vv=f[last].val[k];
x=get(s,j); y=get(s,j+1);
if(!x&&!y){if(i!=n&&j!=m){change(s,j,1);change(s,j+1,2);f[now][s]+=vv;}}
else if(x==1&&y==1){
change(s,find(s,j+1),1);
change(s,j,0); change(s,j+1,0);
f[now][s]+=vv;
}
else if(x==2&&y==2){
change(s,find(s,j),2);
change(s,j,0); change(s,j+1,0);
f[now][s]+=vv;
}
else if(x==1&&y==2){
if(i==n&&j==m) ans+=vv;
}
else if(x==2&&y==1){
change(s,j,0); change(s,j+1,0);
f[now][s]+=vv;
}
else if(x){
if(i!=n) f[now][s]+=vv;
if(j!=m){change(s,j,0);change(s,j+1,x);f[now][s]+=vv;}
}
else if(y){
if(j!=m) f[now][s]+=vv;
if(i!=n){change(s,j+1,0);change(s,j,y);f[now][s]+=vv;}
}
}
}
int main(){
freopen("postman.in","r",stdin);
freopen("postman.out","w",stdout);
scanf("%d%d",&n,&m);
if(n==1||m==1){printf("1\n");return 0;}
if(n<m)swap(n,m);
f[0][0]=1;
now=0;last=1;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++)dp(i,j);
if(i!=n)for(int j=1,k=f[now].cnt;j<=k;j++)
f[now].key[j]<<=2;
}
ans+=ans;
ans.print();
return 0;
}

bzoj 1210 [HNOI2004] 邮递员 插头dp的更多相关文章

  1. BZOJ.1210.[HNOI2004]邮递员(插头DP Hash 高精)

    BZOJ 洛谷 http://www.cnblogs.com/LadyLex/p/7326874.html 插头DP.\(m+1\)个插头的状态需要用三进制表示:\(0\)表示无插头,\(1\)表示是 ...

  2. 【BZOJ1210】[HNOI2004]邮递员 插头DP+高精度

    [BZOJ1210][HNOI2004]邮递员 Description Smith在P市的邮政局工作,他每天的工作是从邮局出发,到自己所管辖的所有邮筒取信件,然后带回邮局.他所管辖的邮筒非常巧地排成了 ...

  3. vijos 1110小胖邮递员;bzoj 1210: [HNOI2004]邮递员

    Description Smith在P市的邮政局工作,他每天的工作是从邮局出发,到自己所管辖的所有邮筒取信件,然后带回邮局.他所管辖的邮筒非常巧地排成了一个m*n的点阵(点阵中的间距都是相等的).左上 ...

  4. 无聊的 邮递员 插头dp

    邮递员想知道,如果他每天都用不同路线走过10×20个点阵邮筒,他必须活过多少个世纪才能走遍所有方案? 7:00 改完T1,开始肝插头dp 7:10 放弃,颓博客 7:20 学习插头dp 7:21 放弃 ...

  5. bzoj 2331: [SCOI2011]地板 插头DP

    2331: [SCOI2011]地板 Time Limit: 5 Sec  Memory Limit: 128 MBSubmit: 541  Solved: 239[Submit][Status] D ...

  6. 【BZOJ】2310: ParkII 插头DP

    [题意]给定m*n的整数矩阵,求经过所有点至多一次路径的最大数值和.n<=8,m<=100. [算法]插头DP [题解]最小表示法确实十分通用,处理简单路径问题只需要状态多加一位表示独立插 ...

  7. BZOJ 2331 [SCOI2011]地板 ——插头DP

    [题目分析] 经典题目,插头DP. switch 套 switch 代码瞬间清爽了. [代码] #include <cstdio> #include <cstring> #in ...

  8. BZOJ 1212: [HNOI2004]L语言( dp + trie )

    因为单词很短...用trie然后每次dp暴力查找...用哈希+dp应该也是可以的.... ------------------------------------------------------- ...

  9. BZOJ 1207: [HNOI2004]打鼹鼠( dp )

    dp.. dp[ i ] = max( dp[ j ] + 1 ) ------------------------------------------------------------------ ...

随机推荐

  1. 重定向和servlet生命周期

    重定向(1)什么是重定向服务器通知浏览器向一个新的地址发送请求.注:可以发送一个302状态码和一个Location消息头.(该消息头包含了一个地址,称之为重定向地址),浏览器收到之后,会立即向重定向地 ...

  2. 管理xcode插件

    1.打开终端 2.输入     open ~/Library/Application\ Support/Developer/Shared/Xcode/Plug-ins 3.出现 4.想删除那个就随意吧 ...

  3. Spring Cloud 入门教程 - Eureka服务注册与发现

    简介 在微服务中,服务注册与发现对管理各个微服务子系统起着关键作用.随着系统水平扩展的越来越多,系统拆分为微服务的数量也会相应增加,那么管理和获取这些微服务的URL就会变得十分棘手,如果我们每新加一个 ...

  4. 在 javascript 中,为什么 [1,2] + [3,4] 不等于 [1,2,3,4]?

    在 stackoverflow 上有人提问:arrays - Why does [1,2] + [3,4] = "1,23,4" in JavaScript? 问题 我想将一个数组 ...

  5. 如何卸载Centos自带jdk

    1.搜索安装的jdk: rpm -qa|grep jdk 结果如下: java-1.7.0-openjdk-1.7.0.45-2.4.3.3.el6.x86_64 java-1.6.0-openjdk ...

  6. python笔记--2--字符串、正则表达式

    字符串 ASCII码采用1个字节来对字符进行编码,最多只能表示256个符号. UTF-8以3个字节表示中文 GB2312是我国制定的中文编码,使用1个字节表示英语,2个字节表示中文:GBK是GB231 ...

  7. 经典Console案例

    /*下面的示例演示 WindowLeft.WindowTop.WindowWidth.WindowHeight.BufferWidth.BufferHeight 和 CursorVisible 属性以 ...

  8. 【转】H.264RTP封包原理

    原文地址:H.264RTP封包原理   作者:cnp11 1.  引言  随着信息产业的发展,人们对信息资源的要求已经逐渐由文字和图片过渡到音频和视频,并越来越强调获取资源的实时性和互动性.但人们又面 ...

  9. python argparse用法总结

    转:python argparse用法总结 1. argparse介绍 argparse是python的一个命令行解析包,非常适合用来编写可读性非常好的程序. 2. 基本用法 prog.py是我在li ...

  10. 架构之微服务(zookeeper)

    ZooKeeper是一个分布式的,开放源码的分布式应用程序协调服务,它包含一个简单的原语集,分布式应用程序可以基于它实现同步服务,配置维护和命名服务等.Zookeeper是hadoop的一个子项目,其 ...