【题意】
给一个N*M的数矩阵
现在求一个子矩阵 要求子矩阵中最大值与最小值的差<=C。而且子矩阵的宽度(横)不超过100(长(竖)没有限制)。 求子矩阵的最大面积。
【输入格式】
第一行两个整数 M(左右方向),N(上下方向)和 C (N,M<=500 0<=C<= 10 )
接下来 N行 每行M个数 每个数(-30000~30000)
【输出格式】
子矩阵的最大面积
【样例输入】
10 15 4
41 40 41 38 39 39 40 42 40 40
39 40 43 40 36 37 35 39 42 42
44 41 39 40 38 40 41 38 35 37
38 38 33 39 36 37 32 36 38 40
39 40 39 39 39 40 40 41 43 41
39 40 41 38 39 38 39 39 39 42
36 39 39 39 39 40 39 41 40 41
31 37 36 41 41 40 39 41 40 40
40 40 40 42 41 40 39 39 39 39
42 40 44 40 38 40 39 39 37 41
41 41 40 39 39 40 41 40 39 40
47 45 49 43 43 41 41 40 39 42
42 41 41 39 40 39 42 40 42 42
41 44 49 43 46 41 42 41 42 42
45 40 42 42 46 42 44 40 42 41
【样例输出】
35
题解:
我的方法复杂度似乎是对的,但被卡了1个点?我太菜了……我的方法是先预处理出
mx[i][j][k]
和
mn[i][j][k]
,分别表示第
i
行第然而我还是被卡。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=505;
int m,n,c,a[maxn][maxn],start[maxn];
int mx[maxn][maxn][105];
int mn[maxn][maxn][105];
bool mark[maxn];
struct Queue{int mx,mn;}q[maxn];
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
int main()
{
// freopen("airport9.in","r",stdin);
m=read();n=read();c=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]=read(),mx[i][j][1]=mn[i][j][1]=a[i][j];
for(int k=2;k<=min(m,35);k++)
for(int i=1;i<=n;i++)
for(int j=1;j+k-1<=m;j++)
mx[i][j][k]=max(mx[i][j][k-1],a[i][j+k-1]),
mn[i][j][k]=min(mn[i][j][k-1],a[i][j+k-1]);
int ans=1;
q[0].mx=-999999999;
q[0].mn=999999999;
for(int i=1;i<=m;i++)start[i]=1;
for(int k=1;k<=min(m,35);k++)
{
for(int j=1;j+k-1<=m;j++)
{
if(mark[j])continue;
int head=1,tail=1,i;
for(i=start[j];i<=n;i++)
if(mx[i][j][k]-mn[i][j][k]<=c)break;
if(i==n+1){mark[j]=true;continue;}
start[j]=i;
if(mx[i][j][k]-mn[i][j][k]<=c)
{
q[1].mx=mx[i][j][k];q[1].mn=mn[i][j][k];
ans=max(ans,k);
for(int l=i+1;l<=n;l++)
{
if(mx[l][j][k]-mn[l][j][k]>c)
{
if((n-l)*k<=ans)break;
head=1;tail=0;continue;
}
while(head<=tail&&max(mx[l][j][k],q[head].mx)-min(mn[l][j][k],q[head].mn)>c)
{
head++;
q[head].mx=max(q[head].mx,q[head-1].mx);
q[head].mn=min(q[head].mn,q[head-1].mn);
}
if(head<=tail)ans=max(ans,k*(tail-head+2));
tail++;
q[tail].mx=mx[l][j][k];
q[tail].mn=mn[l][j][k];
q[head].mx=max(q[head].mx,q[tail].mx);
q[head].mn=min(q[head].mn,q[tail].mn);
}
}
}
}
printf("%d\n",ans);
}