USACO 6.1.2 A Rectangular Barn

题目大意：

给出r*c的矩阵中p个障碍点的位置，求能放置矩形的最大面积。

简单题解：
令h[i][j]表示(i,j)到竖直向上的最近障碍点（或顶端）的线段的长度，left[i][j]、right[i][j]分别表示该线段向左、右延伸的长度。则max{h[i][j]*(left[i][j]+right[i][j]-1)}即为答案。
使用动态规划求h[i][j]、left[i][j]、right[i][j]。

我的代码：

/*
ID:t-x.h1
LANG:C++
TASK:rectbarn
*/
#include<cstdio>
#define up(x,y) if((y)>(x)) (x)=(y)
FILE *fi=fopen("rectbarn.in","r"),*fo=fopen("rectbarn.out","w");
const int MAXn=3000+9,INF=99999999;
bool map[MAXn][MAXn];
int left[MAXn],right[MAXn],h[MAXn];
inline int min(int a,int b)
{
return a?(a<b?a:b):b;
}
int main()
{
int r,c,p,i,j,k,ans=0;
fscanf(fi,"%d%d%d",&r,&c,&p);
for(i=1;i<=p;++i)
{
fscanf(fi,"%d%d",&j,&k);
map[j][k]=1;
}
for(i=1;i<=r;++i)
{
for(k=0,j=1;j<=c;++j)
if(map[i][j])
h[j]=0,left[j]=k=0;
else
++h[j],left[j]=min(left[j],++k);
for(k=0,j=c;j;--j)
{
right[j]=map[i][j]?k=0:min(right[j],++k);
up(ans,h[j]*(left[j]+right[j]-1));
}
}
fprintf(fo,"%d\n",ans);
fclose(fi);
fclose(fo);
return 0;
}