题意:链接
方法:最大子矩阵
解析:
做过了那么多最大子矩阵,这道题随便YY就出来了,但我有个想法啊,如果某个n*m的矩阵,有部分坏点,非坏点有权值,可正可负,则选出的无坏点子矩阵的最大和为多少?
貌似dp?不知道啊。。
先说这道题,因为坏点权值为0,可以预处理出sum[i,j]代表i,j为右下角,(1,1)为左上角的矩阵的权值和。
当然我们可以暴力枚举矩阵,不过这是4次方的过不了。
再观察发现n,m<=1000,所以我们可以选用算♂法②来搞这道题。
先预处理sum,以及le,ri,h。
然后枚举所有好点,计算出其最大子矩阵的四个顶点坐标,一个简单的容斥就随便A了。
代码:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 1010
using namespace std;
int n,m;
int map[N][N];
bool M[N][N];
int sum[N][N];
int le[N][N];
int ri[N][N];
int h[N][N];
void getsum()
{
for(int i=1;i<=n;i++)
{
int s=0;
for(int j=1;j<=m;j++)
{
s+=map[i][j];
sum[i][j]=s+sum[i-1][j];
}
}
}
void init()
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(!M[i][j])
le[i][j]=0;
else le[i][j]=le[i][j-1]+1;
}
for(int j=m;j>=1;j--)
{
if(!M[i][j])
ri[i][j]=0;
else ri[i][j]=ri[i][j+1]+1;
}
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(M[i][j]&&M[i-1][j])
{
h[i][j]=h[i-1][j]+1;
le[i][j]=min(le[i][j],le[i-1][j]);
ri[i][j]=min(ri[i][j],ri[i-1][j]);
}
}
}
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
scanf("%d",&map[i][j]);
M[i][j]=!map[i][j]?0:1;
}
}
getsum();
init();
int ans=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(M[i][j])
{
int ptx1,pty1,ptx2,pty2,ptx3,pty3,ptx4,pty4;
ptx1=i,pty1=j-le[i][j];
ptx2=i-h[i][j]-1,pty2=j-le[i][j];
ptx3=i-h[i][j]-1,pty3=j+ri[i][j]-1;
ptx4=i,pty4=j+ri[i][j]-1;
ans=max(ans,sum[ptx4][pty4]-sum[ptx1][pty1]-sum[ptx3][pty3]+sum[ptx2][pty2]);
}
}
}
printf("%d\n",ans);
}
版权声明:本文为博主原创文章,未经博主允许不得转载。
时间: 2024-10-29 08:06:40