2019.01.24 bzoj2310: ParkII(轮廓线dp)
传送门
题意简述:给一个m*n的矩阵,每个格子有权值V(i,j) (可能为负数),要求找一条路径,使得每个点最多经过一次且点权值之和最大。
思路:我们将求回路时的状态定义改进一下。
现在由于求的是路径说明有可能出现单插头的情况,于是我们用四进制表示,000对应无插头,111对应左括号插头,222对应右括号插头,333对应单插头。
然后多了一些很多转移,于是我把找左右插头的函数给封装了起来 。
细节较多各位慢慢写吧。
然后由于写的hashhashhash表之类的常数太大的主要是懒得写数组转移了导致bzojbzojbzoj速度垫底了。
注意细节。
代码:
#include<bits/stdc++.h>
#define ri register int
#define change (f[cur].insert(stat,mx+a[i][j]))
using namespace std;
inline int read(){
int ans=0,w=1;
char ch=getchar();
while(!isdigit(ch)){if(ch=='-')w=-1;ch=getchar();}
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
return ans*w;
}
const int mod=1e6+7;
struct Statement{
int tot,mx[mod],sta[mod],idx[mod];
inline void clear(){memset(idx,-1,sizeof(idx)),tot=0;}
inline void insert(int stat,int mxn){
int pos=stat%mod;
if(!pos)++pos;
while(~idx[pos]&&sta[idx[pos]]!=stat)pos=pos==mod-1?1:pos+1;
if(~idx[pos])mx[idx[pos]]=max(mx[idx[pos]],mxn);
else mx[idx[pos]=++tot]=mxn,sta[tot]=stat;
}
}f[2];
int n,m,mp[105][10],a[105][10],ans=-1e9;
bool cur;
inline int getbit(int x,int p){return (x>>((p-1)<<1))&3;}
inline void update(int&x,int p,int v){x^=(getbit(x,p)^v)<<((p-1)<<1);}
inline int findl(int stat,int pos){
for(ri bit,cnt=-1,i=pos-1;i;--i){
bit=getbit(stat,i);
if(bit==1)++cnt;
if(bit==2)--cnt;
if(!cnt)return i;
}
}
inline int findr(int stat,int pos){
for(ri bit,cnt=1,i=pos+1;i<=m+1;++i){
bit=getbit(stat,i);
if(bit==1)++cnt;
if(bit==2)--cnt;
if(!cnt)return i;
}
}
inline void solve(){
f[cur=0].clear(),f[cur].insert(0,0);
for(ri i=1;i<=n;++i){
for(ri j=1;j<=m;++j){
f[(cur^=1)].clear();
for(ri tt=1,stat,mx,tmp,stmp,p,q;tt<=f[cur^1].tot;++tt){
stat=f[cur^1].sta[tt],mx=f[cur^1].mx[tt],stmp=stat,update(stmp,j,0),update(stmp,j+1,0);
p=getbit(stat,j),q=getbit(stat,j+1);
if(p==1&&q==2)continue;
if(!(p+q)){
f[cur].insert(stat,mx);
if(mp[i][j+1])update((tmp=stat),j+1,3),f[cur].insert(tmp,mx+a[i][j]);
if(mp[i+1][j])update((tmp=stat),j,3),f[cur].insert(tmp,mx+a[i][j]);
if(mp[i][j+1]&&mp[i+1][j])update(stat,j,1),update(stat,j+1,2),change;
continue;
}
if(!p){
if(mp[i][j+1])change;
if(mp[i+1][j])update((tmp=stat),j,q),update(tmp,j+1,0),f[cur].insert(tmp,mx+a[i][j]);
if(q^3)update(stat,j+1,0),update(stat,(q==1?findr(stat,j+1):findl(stat,j+1)),3),change;
else if(!stmp)ans=max(ans,mx+a[i][j]);
continue;
}
if(!q){
if(mp[i+1][j])change;
if(mp[i][j+1])update((tmp=stat),j,0),update(tmp,j+1,p),f[cur].insert(tmp,mx+a[i][j]);
if(p^3)update(stat,j,0),update(stat,(p==1?findr(stat,j):findl(stat,j)),3),change;
else if(!stmp)ans=max(ans,mx+a[i][j]);
continue;
}
stat=stmp;
if(p==2&&q==1){change;continue;}
if(p==1&&q==1){update(stat,findr(stat,j+1),1),change;continue;}
if(p==2&&q==2){update(stat,findl(stat,j),2),change;continue;}
if(p==1&&q==3){update(stat,findr(stat,j),3),change;continue;}
if(p==2&&q==3){update(stat,findl(stat,j),3),change;continue;}
if(p==3&&q==1){update(stat,findr(stat,j+1),3),change;continue;}
if(p==3&&q==2){update(stat,findl(stat,j+1),3),change;continue;}
if(p==3&&q==3)if(!stat)ans=max(ans,mx+a[i][j]);
}
}
for(ri j=1;j<=f[cur].tot;++j)f[cur].sta[j]<<=2;
}
}
int main(){
n=read(),m=read();
for(ri i=1;i<=n;++i)for(ri j=1;j<=m;++j)ans=max(ans,a[i][j]=read()),mp[i][j]=1;
solve(),cout<<ans;
return 0;
}
2019.01.24 bzoj2310: ParkII(轮廓线dp)的更多相关文章
- 2019.01.24 bzoj3125: CITY(轮廓线dp)
传送门 题意简述:给一个n∗mn*mn∗m的网格图,有的格子不能走,有的格子只能竖着走,有的格子只能横着走,问用一条回路覆盖所有能走的格子的方案数. 思路: 就是简单的轮廓线dpdpdp加了一点限制而 ...
- 2019.01.24 NOIP训练 旅行(轮廓线dp)
传送门 题意简述: 给一个n∗mn*mn∗m的有障碍的网格图,问你从左上角走到左下角并覆盖所有可行格子的路径条数. 思路: 路径不是很好算. 将图改造一下,在最前面添两列,第一列全部能通过,第二列只有 ...
- 梦想MxWeb3D协同设计平台 2019.01.24更新
SDK开发包下载地址:http://www.mxdraw.com/ndetail_10124.html1. 编写快速入门教程2. 重构前端代码,支持一个页面多个三维控件同时加载,或二维和三维同时加 ...
- 【NOI2019模拟2019.7.1】三格骨牌(轮廓线dp转杨图上钩子定理)
Description \(n,m<=1e4,mod ~1e9+7\) 题解: 显然右边那个图形只有旋转90°和270°后才能放置. 先考虑一个暴力的轮廓线dp: 假设已经放了编号前i的骨牌,那 ...
- Tensorflow学习笔记2019.01.22
tensorflow学习笔记2 edit by Strangewx 2019.01.04 4.1 机器学习基础 4.1.1 一般结构: 初始化模型参数:通常随机赋值,简单模型赋值0 训练数据:一般打乱 ...
- HDU 4802 && HDU 4803 贪心,高精 && HDU 4804 轮廓线dp && HDU 4805 计算几何 && HDU 4811 (13南京区域赛现场赛 题目重演A,B,C,D,J)
A.GPA(HDU4802): 给你一些字符串对应的权重,求加权平均,如果是N,P不计入统计 GPA Time Limit: 2000/1000 MS (Java/Others) Memory ...
- HDU - 4804 Campus Design(状压+轮廓线dp)
Campus Design Nanjing University of Science and Technology is celebrating its 60th anniversary. In o ...
- POJ 3254 Corn Fields (状压DP,轮廓线DP)
题意: 有一个n*m的矩阵(0<n,m<=12),有部分的格子可种草,有部分不可种,问有多少种不同的种草方案(完全不种也可以算1种,对答案取模后输出)? 思路: 明显的状压DP啦,只是怎样 ...
- 轮廓线DP POJ3254 && BZOJ 1087
补了一发轮廓线DP,发现完全没有必要从右往左设置状态,自然一点: 5 6 7 8 9 1 2 3 4 如此设置轮廓线标号,转移的时候直接把当前j位改成0或者1就行了.注意多记录些信息对简化代码是很有帮 ...
随机推荐
- 支付宝小程序开发之与微信小程序不同的地方
前言: 本文仅汇总微信小程序移植支付宝小程序过程中遇到的一些不同的地方,详细请参考官方开发文档. 网络请求: 对于网络请求,基本上改动不大,也就支付宝小程序没有responseType属性及响应码字段 ...
- f5会话保持导致站点访问故障
1.环境: vs Service Port为443和80,vip为10.160.100.66,pool为10.160.45.42/43的8080.80的vs重新到443. 2.问题: 单独访问源站都没 ...
- shell脚本小集锦
1) 如何向脚本传递参数 ? ./script argument 例子: 显示文件名称脚本 ./show.sh file1.txt cat show.sh #!/bin/bash 2) 如何在脚本中使 ...
- http协议的学习
TPC/IP协议是传输层协议,主要解决数据如何在网络中传输,而HTTP是应用层协议,主要解决如何包装数据.
- 运行SVO
安装与运行的所有文档:https://github.com/uzh-rpg/rpg_svo/wiki 有两种安装方式: 有ros:https://github.com/uzh-rpg/rpg_svo/ ...
- Windows系统崩溃后快速恢复Oracle数据库的妙招
Windows系统崩溃后快速恢复Oracle数据库,以下是操作步骤 假设oracle数据安装在d:\\oracle文件夹中,数据库名称orcl 1>将崩溃的数据库安装目录"d:\\or ...
- javascript 重构alert()
javascript问题,关于重构window.alert()后,然后调用window原本的window.alert()的方法 大神们,问个问题,如果在script标签的第一行散写,重构了window ...
- 关于js动画简单理解;
1.CSS样式提供了运动 过度的属性:transition 过度的属性值:attr ,time , liner , delay: 值分别是:属性(css),花费的时间,变化的速度(默认匀速), ...
- thinkphp两表联查并且分页
ThinkPHP中关联查询(即多表联合查询)可以使用 table() 方法或和join方法,具体使用如下例所示: 1.原生查询示例: $Model = new Model(); $sql = 'sel ...
- Navicat连接MySQL,出现2059 - authentication plugin 'caching_sha2_password'的解决方案
昨天当我把MySQL的安装程序下载并安装好,然后又下载了另外一个工具来使用它,该工具的名称是Navicat Premium,当我通过该工具连接MySQL Workbench的时候,无法连接,提示“20 ...