Codevs 1159 最大全0子矩阵 悬线法!!!!

1159 最大全0子矩阵

时间限制: 1 s    空间限制: 128000 KB    题目等级 : 黄金 Gold

题目描述 Description

在一个0,1方阵中找出其中最大的全0子矩阵,所谓最大是指O的个数最多。

输入描述 Input Description

输入文件第一行为整数N,其中1<=N<=2000,为方阵的大小,紧接着N行每行均有N个0或1,相邻两数间严格用一个空格隔开。

输出描述 Output Description

输出文件仅一行包含一个整数表示要求的最大的全零子矩阵中零的个数。

样例输入 Sample Input

5
0 1 0 1 0
0 0 0 0 0
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0

样例输出 Sample Output

9

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 #include<iostream>
 5 #include<cstdio>
 6 using namespace std;
 7 const int maxn=2333;
 8 int h[maxn][maxn],l[maxn][maxn],r[maxn][maxn];
 9 int n,maxl,maxr,ans,martrix[maxn][maxn];
10 int main() {
11     scanf("%d",&n);
12     for(int i=1;i<=n;++i)
13         for(int j=1;j<=n;++j)
14             scanf("%d",&martrix[i][j]);
15     for(int i=1;i<=n;++i)
16         r[0][i]=n;
17     for(int i=1;i<=n;++i) {
18         maxr=n;
19         maxl=1;
20         for(int j=1;j<=n;++j)
21             if(martrix[i][j]) {
22                 maxl=j+1;
23                 h[i][j]=l[i][j]=0;
24             }
25             else {
26                 h[i][j]=h[i-1][j]+1;
27                 l[i][j]=max(maxl,l[i-1][j]);
28             }
29         for(int j=n;j>0;--j)
30             if(martrix[i][j]) {
31                 maxr=j-1;
32                 r[i][j]=n;
33             }
34             else {
35                 r[i][j]=min(maxr,r[i-1][j]);
36                 ans=max(ans,(r[i][j]-l[i][j]+1)*h[i][j]);
37             }
38     }
39     printf("%d\n",ans);
40     return 0;
41 }

看到这个题目的时候妥妥的蒙蔽了~~,按照正方形敲出来就A了一个点,题解居然是悬线法!!!!(跑去学了悬线法)

就是记录一下能扩展到的最靠左、右、上的边界,搞一下就出来了。。。

悬线法

时间复杂度O(NM) 空间复杂度O(NM)

定义

有效竖线:除了两个端点外,不覆盖任何一个障碍点的竖直线段。 
悬线:上端覆盖了一个障碍点或者到达整个矩形上边界的有效线段。 
每个悬线都与它底部的点一一对应,矩形中的每一个点(矩形顶部的点除外)都对应了一个悬线。 
悬线的个数=(N-1)*M; 
如果把一个极大子矩形按照横坐标的不同切割成多个与y轴平行的线段,那么其中至少有一个悬线。 
如果把一个悬线向左右两个方向尽可能的移动,那么就得到了一个矩形,我们称它为悬线对应的矩形。 
悬线对应的矩形不一定是极大子矩形,因为下边界可能还可以向下扩展。 
设计算法: 
原理:所有悬线对应矩形的集合一定包含了极大子矩形的集合。 
通过枚举所有的悬线,找出所有的极大子矩形。 
算法规模: 
悬线个数=(N-1)×M 
极大子矩形个数≤悬线个数 
具体方法: 
设 H[i,j]为点(i,j)对应的悬线的长度。 
L[i,j]为点(i,j)对应的悬线向左最多能够移动到的位置。 
R[i,j]为点(i,j)对应的悬线向右最多能够移动到的位置。

!考虑点(i,j)对应的悬线与点(i-1,j)对应的悬线的关系(递推思想):

如果(i-1,j)为障碍点,那么,(i,j)对应的悬线长度1,左右能移动到的位置是整个矩形的左右边界。 
即 H[i,j]=1, 
L[i,j]=0,R[i,j]=m

如果(i-1,j)不是障碍点,那么,(i,j)对应的悬线长度为(i-1,j)对应的悬线长度+1。 
即 H[i,j]=H[i-1,j]+1

•如果(i-1,j)不是障碍点,那么(i,j)对应的悬线左右能移动的位置要在(i-1,j)的基础上变化。 
L[i-1,j] 
L[i,j]=max 
(i-1,j)左边第一个障碍点的位置 
•同理,也可以得到R[i,j]的递推式 
R[i-1,j] 
R[i,j]=min 
(i-1,j)右边第一个障碍点的位置

时间: 2024-08-05 06:34:56

Codevs 1159 最大全0子矩阵 悬线法!!!!的相关文章

[codevs1159]最大全0子矩阵(悬线法)

解题关键:悬线法模板题.注意此模板用到了滚动数组. #include<cstdio> #include<cstring> #include<algorithm> #include<cstdlib> #include<iostream> #include<cmath> #define maxn 2002 using namespace std; typedef long long ll; int map[maxn][maxn],l[ma

1159 最大全0子矩阵

/*f(i,j)表示以(i,j)为右下角的最大全0子矩阵的边长若a[i][j]==1,f(i,j)=0否则:f(i,j)=min{ f(i-1,j),f(i,j-1),f(i-1,j-1) }+1 这样求得的是最大全0正方形子矩阵要求长方形矩阵,上述思路行不通假设以(i,j)为右下角的最大矩阵=12它可能是3*4.4*3.2*6.6*2.1*12.12*1按上述思路进行状态转移的话,取得最优值的方案不唯一时,所有的方案需要都记下,用于后续的状态转移. 在长方形全0子矩阵中,考察某个位置(i,j)

最大子段和||最大子矩阵和||最大全0子矩阵||最大全0子正方形

最大子段和(略) 定义一个最大值dp[i]表示以i结尾的最大子段和: 初始化: dp[0]=A[0]: dp[i]=max(dp[i-1]+A[i],A[i]) 即dp[i-1]+A[i]<0时 dp[i] = A[i]; 否则 dp[i]=dp[i-1]+A[i] 最大全0子矩阵 https://blog.csdn.net/Flere825/article/details/54605662 1159 最大全0子矩阵 /* #define N 2010 int a[N][N],ans,n; in

BZOJ 1057: [ZJOI2007]棋盘制作 悬线法求最大子矩阵+dp

1057: [ZJOI2007]棋盘制作 Description 国际象棋是世界上最古老的博弈游戏之一,和中国的围棋.象棋以及日本的将棋同享盛名.据说国际象棋起源于易经的思想,棋盘是一个8*8大小的黑白相间的方阵,对应八八六十四卦,黑白对应阴阳.而我们的主人公小Q,正是国际象棋的狂热爱好者.作为一个顶尖高手,他已不满足于普通的棋盘与规则,于是他跟他的好朋友小W决定将棋盘扩大以适应他们的新规则.小Q找到了一张由N*M个正方形的格子组成的矩形纸片,每个格子被涂有黑白两种颜色之一.小Q想在这种纸中裁减

bzoj1057: [ZJOI2007]棋盘制作(悬线法)

题目要求纵横坐标和奇偶性不同的点取值不同,于是我们把纵横坐标和奇偶性为1的点和0的点分别取反,就变成经典的最大全1子矩阵问题了,用悬线法解决. #include<iostream> #include<cstring> #include<cstdlib> #include<cstdio> using namespace std; const int maxn=2010,inf=1e9; int n,m,ans1,ans2; int h[maxn],mp[max

重(zhong)新学习悬线法

2017.9.28今天模拟赛T1就是求最大子矩阵的经典题 然而我已经好久没有写悬线了,以前悬线也是拉的,于是现在就系统的学习一波吧 给定一个N * M的01矩阵,求最大全0矩阵的大小 N^3的做法其实有很多种,前缀和乱搞什么的都可以 考虑N^2的做法 对于任意一个非1的点i,j,记其向上能到达的最长长度为up[i][j] 显然up[i][j] = up[i-1][j] + 1 (a[i][j] == 0) up[i][j] = 0 (a[i][j] = 1)我们称up[i][j]为一条悬线 显然

BZOJ 3039 玉蟾宫 悬线法

题目大意:给出一张地图,求出这张地图中最大的子矩阵,使得这个子矩阵不包含字母'R'. 思路:简单的悬线法求最大子矩阵,还是不带权值的,很好求.好久没写悬线了,复习一下. CODE: #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define MAX 1010 using namespace std; int m,n; bool map[MAX][M

【BZOJ-3039&amp;1057】玉蟾宫&amp;棋盘制作 悬线法

3039: 玉蟾宫 Time Limit: 2 Sec  Memory Limit: 128 MBSubmit: 753  Solved: 444[Submit][Status][Discuss] Description 有一天,小猫rainbow和freda来到了湘西张家界的天门山玉蟾宫,玉蟾宫宫主蓝兔盛情地款待了它们,并赐予它们一片土地.这片土地被分成N*M个格子,每个格子里写着'R'或者'F',R代表这块土地被赐予了rainbow,F代表这块土地被赐予了freda.现在freda要在这里卖

BZOJ_3039_玉蟾宫_(动态规划+悬线法)

描述 http://www.lydsy.com/JudgeOnline/problem.php?id=3039 n*m的矩阵由R和F组成,求全是F的子矩阵的大小的三倍. 分析 悬线法: 浅谈用极大化思想解决最大子矩形问题--王知昆 l[x][y]表示点(x,y)在它那一行最多能扩展到左边的位置. r[x][y]表示点(x,y)在它那一行最多能扩展到右边的位置. 每一行分别预处理l与r. 在做dp的时候:如果点(x,y)可以取,那么h[x][y]=h[x-1][y]+1,l[x][y]=max(l