思路详见 王知昆《浅谈用极大化思想解决最大子矩形问题》
写得很详细（感谢~....）
因为不太会用递推，所以用了第一种方法，时间复杂度是O（n^2)，n为枚举的点数，对付这题绰绰有余
思路也很简单（建议一定看原文）
先根据x排序
之后两重循环，枚举i后的每一个点j到i可以形成的矩形面积
怎么求这个矩形面积呢？
非常简单，miny,maxy,分别表示纵坐标的上下界
如果枚举的点j比i的y大，那么就修改上界，反之，修改下界（具体的，可以看论文中的图，更直观些）
这里需要注意两个遗漏的地方（论文中也有特别提到）
就是子矩形左右边界是否与大矩形有重合的问题
由于我写的代码，是在算完i到j的矩形面积之后，再更新上下界的
所以，在第二重循环结束之后，在算一次：s=(lx-x)*(maxy-miny)即可
这是从左到右枚举的情况（可解决右边界重合）
同理，再从右到左枚举一遍，即可解决左边界重合的情况
（mark:我并没有解决，左右边界同时重合的情况，但是还是A了...数据太弱？）
ps:两边同时重合的判断也很简单，分别扩展一次即可（这里就不完善了QAQ...因为懒）
附上代码：
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
int lx,ly,n,s,maxs=0;
struct node{
int x,y;
}f[5001];
bool cmp(const node &a,const node &b){//按照x排序，其次按y排序（其实与y如何排序没有太大关系
if(a.x<b.x) return 1;
else if(a.x==b.x && a.y>b.y) return 1;
else return 0;
}
int main(){
//freopen("data.txt","r",stdin);
scanf("%d%d",&lx,&ly);
scanf("%d",&n);
if(n==0){//特判
cout<<lx*ly;
return 0;
}
for(int i=1;i<=n;i++){
scanf("%d%d",&f[i].x,&f[i].y);
}
sort(f+1,f+n+1,cmp);
int dx,dy;
for(int i=1;i<=n;i++){//从左到右枚举
s=0;
int miny=0,maxy=ly;//上下界初始化
for(int j=i+1;j<=n;j++){
dx=f[j].x-f[i].x;
dy=maxy-miny;
s=dx*dy;
if(s>maxs) maxs=s;//根据论文中的图解模拟即可，需要在算完S之后更新上下界
if(f[j].y>=f[i].y && f[j].y<maxy) maxy=f[j].y;
if(f[j].y<=f[i].y && f[j].y>miny) miny=f[j].y;
}
dx=lx-f[i].x;
dy=maxy-miny;//由于已更新了最后扫到的点，最后不需要再更新
s=dx*dy;
if(s>maxs) maxs=s;
}
for(int i=n;i>=1;i--){//从右到左枚举同上思路（其实可以写一个函数减少代码量QAQ）
s=0;
int miny=0,maxy=ly;
for(int j=i-1;j>=1;j--){
dx=f[i].x-f[j].x;
dy=maxy-miny;
s=dx*dy;
if(s>maxs) maxs=s;
if(f[j].y>=f[i].y && f[j].y<maxy) maxy=f[j].y;
if(f[j].y<=f[i].y && f[j].y>miny) miny=f[j].y;
}
dx=f[i].x-0;
dy=maxy-miny;
s=dx*dy;
if(s>maxs) maxs=s;
}
cout<<maxs;
return 0;
}