传送门

题意简述:给一个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)的更多相关文章

  1. 2019.01.24 bzoj3125: CITY(轮廓线dp)

    传送门 题意简述:给一个n∗mn*mn∗m的网格图,有的格子不能走,有的格子只能竖着走,有的格子只能横着走,问用一条回路覆盖所有能走的格子的方案数. 思路: 就是简单的轮廓线dpdpdp加了一点限制而 ...

  2. 2019.01.24 NOIP训练 旅行(轮廓线dp)

    传送门 题意简述: 给一个n∗mn*mn∗m的有障碍的网格图,问你从左上角走到左下角并覆盖所有可行格子的路径条数. 思路: 路径不是很好算. 将图改造一下,在最前面添两列,第一列全部能通过,第二列只有 ...

  3. 梦想MxWeb3D协同设计平台 2019.01.24更新

    SDK开发包下载地址:http://www.mxdraw.com/ndetail_10124.html1.  编写快速入门教程2.  重构前端代码,支持一个页面多个三维控件同时加载,或二维和三维同时加 ...

  4. 【NOI2019模拟2019.7.1】三格骨牌(轮廓线dp转杨图上钩子定理)

    Description \(n,m<=1e4,mod ~1e9+7\) 题解: 显然右边那个图形只有旋转90°和270°后才能放置. 先考虑一个暴力的轮廓线dp: 假设已经放了编号前i的骨牌,那 ...

  5. Tensorflow学习笔记2019.01.22

    tensorflow学习笔记2 edit by Strangewx 2019.01.04 4.1 机器学习基础 4.1.1 一般结构: 初始化模型参数:通常随机赋值,简单模型赋值0 训练数据:一般打乱 ...

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

  7. HDU - 4804 Campus Design(状压+轮廓线dp)

    Campus Design Nanjing University of Science and Technology is celebrating its 60th anniversary. In o ...

  8. POJ 3254 Corn Fields (状压DP,轮廓线DP)

    题意: 有一个n*m的矩阵(0<n,m<=12),有部分的格子可种草,有部分不可种,问有多少种不同的种草方案(完全不种也可以算1种,对答案取模后输出)? 思路: 明显的状压DP啦,只是怎样 ...

  9. 轮廓线DP POJ3254 && BZOJ 1087

    补了一发轮廓线DP,发现完全没有必要从右往左设置状态,自然一点: 5 6 7 8 9 1 2 3 4 如此设置轮廓线标号,转移的时候直接把当前j位改成0或者1就行了.注意多记录些信息对简化代码是很有帮 ...

随机推荐

  1. JS判断变量类型

    typeof v 只能用于识别基础类型,不能识别对象 v instanceof MyClass 判断类型 Object.prototype.toString.call(v.p) === "[ ...

  2. JS参数转发

    在没有装饰器之前不方便. 可以用Reflect.apply. cls = function f() { let obj = {}; obj.show = function(a, b) { consol ...

  3. leetcode输入输出加速

    C++兼容C的输入输出,即cin与scanf混用文件指针不会出错,cout亦同,导致cin有额外开销. 可以用std::ios::sync_with_stdio(false);手动关闭. cin.ti ...

  4. 【nginx】配置详解

    nginx作为web server,是最常用的网络服务器.来看下它的配置文件. 实验1:最基础的配置(可访问静态资源) 包括listen,server_name,root,index,access_l ...

  5. [剑指Offer]59-队列的最大值(题目二待补)

    题目一:滑动窗口的最大值 题目链接 https://www.nowcoder.com/practice/1624bc35a45c42c0bc17d17fa0cba788?tpId=13&tqI ...

  6. iOS - 常用本机URL跳转设置

    import UIKit class ViewController: UIViewController { override func touchesBegan(_ touches: Set<U ...

  7. JS 获取屏幕的宽度和高度,各种方式

      Javascript: 网页可见区域宽: document.body.clientWidth网页可见区域高: document.body.clientHeight网页可见区域宽: document ...

  8. java.security.InvalidKeyException: Illegal key size 解决办法

    下载对应的文件并替换到指定目录 Java Cryptography Extension (JCE) Unlimited Strength Jurisdiction Policy Files 6 Jav ...

  9. Oauth2.0 认证的Web api例子

    Oauth2.0的解释 OAuth(开放授权)是一个开放标准,允许用户授权第三方移动应用访问他们存储在另外的服务提供者上的信息,而不需要将用户名和密码提供给第三方移动应用或分享他们数据的所有内容.OA ...

  10. c语言使用指针交换数值

    练习题:将两个int类型数值交换 #include <stdio.h> void swap(int*,int*); int main(void){ , hex = 0x5f1043; sw ...