题目背景

感谢@throusea 贡献的两组数据

题目描述

回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*)。她发现,把大池子视为01矩阵(0表示对应位置无鱼,1表示对应位置有鱼)有助于决定吃鱼策略。

在代表池子的01矩阵中,有很多的正方形子矩阵,如果某个正方形子矩阵的某条对角线上都有鱼,且此正方形子矩阵的其他地方无鱼,猫猫就可以从这个正方形子矩阵“对角线的一端”下口,只一吸,就能把对角线上的那一队鲜鱼吸入口中。

猫猫是个贪婪的家伙,所以她想一口吃掉尽量多的鱼。请你帮猫猫计算一下,她一口下去,最多可以吃掉多少条鱼?

输入输出格式

输入格式:

有多组输入数据,每组数据:

第一行有两个整数n和m(n,m≥1),描述池塘规模。接下来的n行,每行有m个数字(非“0”即“1”)。每两个数字之间用空格隔开。

对于30%的数据,有n,m≤100

对于60%的数据,有n,m≤1000

对于100%的数据,有n,m≤2500

输出格式:

只有一个整数——猫猫一口下去可以吃掉的鱼的数量,占一行,行末有回车。

输入输出样例

输入样例#1:

4 6
0 1 0 1 0 0
0 0 1 0 1 0
1 1 0 0 0 1
0 1 1 0 1 0

输出样例#1:

3

说明

右上角的

1 0 0
0 1 0
0 0 1

解析:

这道题跟P1387 最大正方形 比较像,具体做法都是去检验形成新正方形的可行性。

想了很久不知道怎么处理从上面和左边的转移,原本就是想统计\(0\)的个数的,可惜没尝试。


两次\(dp\),分别处理两条对角线的最优解。

以从左上到右下的对角线为例,需要预处理出任意位置的前面的\(0\)的数量\(pre1[i][j],pre2[i][j]\)。

显然,如果某处是\(1\),那么它的下面的和右边的点前面的\(0\)的数量就是\(0\),如果是\(0\)那么它的下面的和右边的点的前面的\(0\)就是此处的\(0\)的个数加\(1\)。

设\(dp[i][j]\)表示以\((i,j)\)为右下角顶点的满足题目条件的正方形的边长,由定义得在\((1\sim i,1\sim j)\)这个区域的正方形除了对角线上是\(1\)外,其他地方都是\(0\)。

所以在考虑状态转移时,对于\(dp[i][j]\),我们已经得知\(dp[i-1][j-1]\)的最大正方形边长,那么也就得知在\((i-dp[i][j]\sim i,j-dp[i][j]\sim j )\)这个区域中,除了对角线上是\(1\)外其余都是\(0\)。如果要让\((i,j)\)为右下角顶点的正方形最大,我们就需要考虑上面提到的\(pre\),若当前位置这个\(pre\)还没\(dp[i-1][j-1]\)大,那我们只好让以\((i,j)\)为右下角顶点的正方形边长为\(pre[i][j]\)中的较小值了。

于是有状态转移方程:

\(dp[i][j]=min(dp[i-1][j-1],min(pre1[i-1][j],pre2[i][j-1]))+1\)

另一条对角线同理。

注意由于\(dp\)是以递推的形式进行的,所以我们可以边输入便处理,加快速度。

当然不这样也不会T的\(QWQ\)

参考代码:

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<ctime>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#define N 2505
#define INF 0x3f3f3f3f
#define ri register int
using namespace std;
int a[N][N],dp[N][N],pre1[N][N],pre2[N][N],n,m;
inline int read()
{
int f=1,x=0;char c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
int main()
{
n=read();m=read();
int ans=0;
for(ri i=1;i<=n;i++)
for(ri j=1;j<=m;j++){
a[i][j]=read();
if(!a[i][j]){
pre1[i][j]=pre1[i-1][j]+1;
pre2[i][j]=pre2[i][j-1]+1;
}
else dp[i][j]=min(dp[i-1][j-1],min(pre1[i-1][j],pre2[i][j-1]))+1;
ans=max(dp[i][j],ans);
}
memset(dp,0,sizeof(dp));
memset(pre1,0,sizeof(pre1));
memset(pre2,0,sizeof(pre2));
for(ri i=1;i<=n;i++){
for(ri j=m;j>=1;j--){
if(!a[i][j]){
pre1[i][j]=pre1[i-1][j]+1;
pre2[i][j]=pre2[i][j+1]+1;
}
else dp[i][j]=min(dp[i-1][j+1],min(pre1[i-1][j],pre2[i][j+1]))+1;
ans=max(dp[i][j],ans);
}
}
cout<<ans<<endl;
return 0;
}

P1736 创意吃鱼法[二维dp]的更多相关文章

  1. P1736 创意吃鱼法 图的DP

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

  2. P1387 最大正方形&&P1736 创意吃鱼法

    P1387 最大正方形 P1736 创意吃鱼法 两道类似的$DP$ 转移方程基本上类似于$f[i][j]=min(f[i-1][j-1],min(f[i][j-1],f[i-1][j]))$ 考虑构成 ...

  3. 洛谷 P1736 创意吃鱼法(多维DP)

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

  4. 洛谷 P1736 创意吃鱼法 Label:dp || 前缀和

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

  5. 洛谷P1736 创意吃鱼法 dp

    正解:dp 解题报告: 早就想写dp的题目辣!我发现我的dp好差啊QAQ所以看到列表的小朋友写dp的题目就跟着他们的步伐做下题好辣QwQ 这题的话没有那——么难,大概说下趴QwQ 首先说下题意 前面一 ...

  6. P1387 最大正方形 && P1736 创意吃鱼法(DP)

    题目描述 在一个n*m的只包含0和1的矩阵里找出一个不包含0的最大正方形,输出边长. 输入输出格式 输入格式: 输入文件第一行为两个整数n,m(1<=n,m<=100),接下来n行,每行m ...

  7. 洛谷 P1736 创意吃鱼法

    题目描述 题目链接:https://www.luogu.org/problemnew/show/P1736 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢( ...

  8. P1736 创意吃鱼法

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

  9. P1736 创意吃鱼法 (动态规划)

    题目描述 回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中,然后开始思考:到底要以何种方法吃鱼呢(猫猫就是这么可爱,吃鱼也要想好吃法 ^_*).她发现,把大池子视为01矩阵(0表示对应位置无鱼,1 ...

随机推荐

  1. http响应时长分析

    curl -o /dev/null -s -w "time_namelookup:%{time_namelookup}::time_connect:%{time_connect}::time ...

  2. 修改centos的源, 使其能访问阿里云服务器

    1. 从拷贝以下文件到/etc/yum.repos.d/ 下; Centos-7.repo # CentOS-Base.repo # # The mirror system uses the conn ...

  3. Appium移动自动化测试-----(四)安装 appium Server

    我们可以在Appium官方网站上下载操作系统相应的Appium版本. https://bitbucket.org/appium/appium.app/downloads/ 当前最新版本为 Appium ...

  4. NMF包的安装

    win10 操作系统,R3.4版本 NMF包的下载 source("https://bioconductor.org/biocLite.R") biocLite("Bio ...

  5. [转帖]Hive 快速入门(全面)

    Hive 快速入门(全面) 2018-07-30 16:11:56 琅琊山二当家 阅读数 4343更多 分类专栏: hadoop 大数据   转载: https://www.codercto.com/ ...

  6. Redis--list类型操作命令

    列表 list Redis列表是简单的字符串列表,按照插入顺序排序.你可以添加一个元素导列 表的头部(左边)或者尾部(右边) 列表 list——基本命令 lpush 语法:lpush key valu ...

  7. 20191031:GIL全局解释锁

    20191031:GIL全局解释锁 总结关于GIL全局解释锁的个人理解 GIl全局解释锁,本身不是Python语言的特性,而是Python语言底层的c Python解释器的一个特性.在其他解释器中是没 ...

  8. ideaui安装破解、相关配置、JRebel破解

    前言: Ideaui 2019(官网 https://www.jetbrains.com/idea/?fromMenu) 安装软件就简单啦,下载选择路径就完事了,注意文件名全英文:   但是按照咱们传 ...

  9. (转)基于FFPMEG2.0版本的ffplay代码分析

    ref:http://zzhhui.blog.sohu.com/304810230.html 背景说明 FFmpeg是一个开源,免费,跨平台的视频和音频流方案,它提供了一套完整的录制.转换以及流化音视 ...

  10. Vue $emit $event 传值(子to父)

    事件名 始终使用 kebab-case 的事件名. 通过事件向父组件发送信息 子组件中EnFontsize.vue中$emit <button @click="$emit('enlar ...