洛谷P1736 创意吃鱼法【DP】

时空限制 1000ms / 128MB

题目描述

回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中，然后开始思考：到底要以何种方法吃鱼呢（猫猫就是这么可爱，吃鱼也要想好吃法 ^_*）。她发现，把大池子视为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

输出格式：

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

---

题目分析

用两个dp数组$dp[i][j]$表示以$(i,j)$为右/左下角的符合要求的正方形的最大边长
不难想到两个$dp[i][j]$需要分别从$dp[i-1][j-1]/dp[i-1][j+1]$转移得来

但有一个比较容易犯的思维漏洞

若不能全部继承$dp[i-1][j-1]$，则$dp[i][j]=1$
这是错的！错的！！错的！！！

比如像这样的矩阵
1 0 0
0 1 0
1 0 1
假设已有$dp[2][2]=2$，发现不能有$dp[3][3]=dp[2][2]+1=3$
但是却可以有$dp[3][3]=2$

从这里我们可以发现
$dp[i][j]$是从$1$~$dp[i-1][j-1]$中那个最大的符合条件的数继承来的
而这个数显然具有可二分性

于是正确的思路应该是在区间$[1,dp[i-1][j-1]]$二分出最大的符合条件的数x
令$dp[i][j]=x+1$

然而这题数据弱，直接枚举区间1~dp[i-1][j-1]中的数也过了

---

#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=2510;
int n,m;
int a[maxn][maxn];
int dp1[maxn][maxn],dp2[maxn][maxn],ans;

int main()
{
    n=read();m=read();
    for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j)
    a[i][j]=read();
    
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=m;++j)
        {
    	    if(!a[i][j])continue; int k;
            
            for(k=1;k<=dp1[i-1][j-1];++k)
            if(a[i-k][j]||a[i][j-k]) break;
            dp1[i][j]=k;
            
            for(k=1;k<=dp2[i-1][j+1];++k)
            if(a[i-k][j]||a[i][j+k]) break;
            dp2[i][j]=k;
            
            ans=max(ans,max(dp1[i][j],dp2[i][j]));
        }
    }
    printf("%d",ans);
    return 0;
}