洛谷 P2701 [USACO5.3]巨大的牛棚Big Barn Label:二维数组前缀和 你够了 这次我用DP

题目背景

(USACO 5.3.4)

题目描述

农夫约翰想要在他的正方形农场上建造一座正方形大牛棚。他讨厌在他的农场中砍树,想找一个能够让他在空旷无树的地方修建牛棚的地方。我们假定,他的农场划分成 N x N 的方格。输入数据中包括有树的方格的列表。你的任务是计算并输出,在他的农场中,不需要砍树却能够修建的最大正方形牛棚。牛棚的边必须和水平轴或者垂直轴平行。

EXAMPLE

考虑下面的方格,它表示农夫约翰的农场,‘.'表示没有树的方格,‘#'表示有树的方格

1 2 3 4 5 6 7 8

1 . . . . . . . .

2 . # . . . # . .

3 . . . . . . . .

4 . . . . . . . .

5 . . . . . . . .

6 . . # . . . . .

7 . . . . . . . .

8 . . . . . . . .

最大的牛棚是 5 x 5 的,可以建造在方格右下角的两个位置其中一个。

输入输出格式

输入格式:

Line 1: 两个整数: N (1 <= N <= 1000),农场的大小,和 T (1 <= T <= 10,000)有树的方格的数量

Lines 2..T+1: 两个整数(1 <= 整数 <= N), 有树格子的横纵坐标

输出格式:

只由一行组成,约翰的牛棚的最大边长。

输入输出样例

输入样例#1:

8 3
2 2
2 6
6 3

输出样例#1:

5

说明

题目翻译来自NOCOW。

USACO Training Section 5.3

代码

 #include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std; int a[][],N,M;
int f[][],ans; int main(){
// freopen("01.in","r",stdin);//freopen("01.out","w",stdout); scanf("%d%d",&N,&M);
for(int i=;i<=M;i++){int x,y;scanf("%d%d",&x,&y);a[x][y]=;} for(int i=;i<=N;i++){
for(int j=;j<=N;j++){
if(a[i][j]) continue;
else f[i][j]=min(min(f[i-][j],f[i][j-]),f[i-][j-])+;
ans=max(f[i][j],ans);
}
} printf("%d\n",ans); fclose(stdin);fclose(stdout);return ;
}

我就给个思路,两种解法:
1. DP,用f[i,j]表示以当前格子为右→下角能构成的最大的正方形大小,方程f[i][j]=min(min(f[i-1][j],f[i][j-1]),f[i-1][j-1])+1;代码如上,此处不再赘述;
2. 二维数组前缀和,可能会超时,可以参见 洛谷P2733 ,跟这道题的区别就是没有区别 ,核心代码如下,如有不理解可参见博客:

for(int i=1;i<=N;i++){//二维数组前缀和
  for(int j=1;j<=N;j++) sum[i][j]=sum[i][j-1]+int(a[i][j]-'0');
  for(int j=1;j<=N;j++) sum[i][j]+=sum[i-1][j];
}

for(int k=2;k<=N;k++)//对边长进行枚举
  for(int i=1;i+k-1<=N;i++)
    for(int j=1;j+k-1<=N;j++)
      if(cal(i,j,i+k-1,j+k-1)==k*k)
        ans=k;

以下是两道运用前缀和的题目
[洛谷 P1736 创意吃鱼法](http://www.cnblogs.com/radiumlrb/p/5808527.html)
[洛谷 P1387 最大正方形 ](http://www.cnblogs.com/radiumlrb/p/5808285.html)

洛谷 P2701 [USACO5.3]巨大的牛棚Big Barn Label:二维数组前缀和 你够了 这次我用DP的相关教程结束。

相关推荐:

Python中如何实现NumPy数组的符号函数运算_利用sign函数判断正负

np.sign()返回整数1、-1、0分别表示正、负、零;复数返回单位复数;不区分±0.0;用作布尔掩码前须处理零值,性能依赖数据类型,不支持非数值类型。 NumPy的np.sign()到底返回什么值 np.sign()对数组元素逐个判断符号,返回整数:正数→1,负数→-1,零→0。注意它不区分+0...

如何在Python中对NumPy数组进行转置与轴换位操作?

NumPy数组转置与轴换位本质相同但接口语义不同:.T仅适用于二维数组的简单转置,三维及以上会完全反转轴序;精确控制需用np.transpose()、np.swapaxes()或np.moveaxis()。 NumPy数组的转置和轴换位本质是同一类操作,但接口不同、语义侧重不同:用.T或np.tra...

如何在 NumPy 中将二维数组沿新轴扩展为三维数组

本文介绍多种高效、内存友好的方法,将二维NumPy数组沿新增的第一维(axis=0)重复扩展为三维数组,涵盖np.expand_dims+np.repeat、np.tile、广播机制等主流方案,并对比其性能与适用场景。 本文介绍多种高效、内存友好的方法,将二维numpy数组沿新增的第一维(axis=...

如何使用 Bash 的 join 和 sed 实现基于前缀匹配的多文件合并

本文介绍如何通过Bash命令行工具(join+sed+进程替换)高效合并两个制表符分隔文件,使第一文件的首列(如a)与第二文件首列前缀(如a_rna_1)自动对齐并横向拼接,无需临时文件或编程语言介入。 本文介绍如何通过bash命令行工具(`join`+`sed`+进程替换)高效合并两个制表符分隔文...

如何用 Bash 的 join 和 sed 实现基于前缀匹配的多文件合并

本文介绍一种无需编程语言、仅用标准Linux工具(join+sed+sort)即可完成两列不完全相等但存在前缀关联的TSV文件合并的方法,适用于基因组学等场景中如a与a_rna_1的映射合并。 本文介绍一种无需编程语言、仅用标准linux工具(`join`+`sed`+`sort`)即可完成两列不完...

Python中如何计算NumPy数组的特征值与特征向量_使用linalg.eig

np.linalg.eig返回的特征向量是列向量,即v[:,i]对应第i个特征值w[i],但print(v)屏幕显示按行展示,易误读为“每行一个向量”;实际布局符合线性代数惯例,验证方式为检查A@v[:,0]≈w[0]*v[:,0]。 为什么linalg.eig返回的特征向量是列向量,但结果却是按行...