UVALive 2689 极大子正方形

      本题思路详见论文《浅谈用极大化思想解决最大子矩形问题》的算法一，我们的子正方形一定在子矩形的内部，且边长为矩形长宽的一个较小值。

#include<cstdio>
#include<algorithm>
#define maxn 109
using namespace std;
int ansl,ansx,ansy;
struct node
{
    int x,y;
    bool operator<(const node &rhs)const
    {
        return x<rhs.x;
    }
}p[maxn];
bool cmp(node a,node b)
{
    return a.y<b.y;
}
void update(int x1,int x2,int y1,int y2)
{
    if(ansl<min(x2-x1,y2-y1))
    {
        ansl=min(x2-x1,y2-y1);
        ansx=x1;
        ansy=y1;
    }
}
int main()
{
    int tt;
    scanf("%d",&tt);
    while(tt--)
    {
        int n,W,H;
        ansl=1,ansx=ansy=0;
        scanf("%d%d%d",&n,&W,&H);
        bool f0=0,fn=0;
        for(int i=0;i<n;i++)
        {
            scanf("%d%d",&p[i].x,&p[i].y);
            if(p[i].x==0&&p[i].y==0)f0=1;
            if(p[i].x==W&&p[i].y==H)fn=1;
        }
        if(!f0)p[n++]=(node){0,0};  //把边界作为特殊的点插入点集
        if(!fn)p[n++]=(node){W,H};
        sort(p,p+n);
        for(int i=0;i<n;i++)
        {
            int X=p[i].x,Y=p[i].y,L=0,R=H;
            for(int j=i+1;j<n;j++)
            {
                if(p[j].x==X)continue;
                update(X,p[j].x,L,R);
                if(p[j].y<Y)
                    L=max(L,p[j].y);
                else if(p[j].y>Y)
                    R=min(R,p[j].y);
                else
                    break;
            }
        }
        for(int i=n-1;i>=0;i--)
        {
            int L=0,R=H,X=p[i].x,Y=p[i].y;
            for(int j=i-1;j>=0;j--)
            {
                if(X==p[j].x)continue;
                update(p[j].x,X,L,R);
                if(p[j].y>Y)
                    R=min(R,p[j].y);
                else if(p[j].y<Y)
                    L=max(L,p[j].y);
                else
                    break;
            }
        }
        sort(p,p+n,cmp);
        for(int i=1;i<n;i++)
        {
            if(p[i].y!=p[i-1].y)
                update(0,W,p[i-1].y,p[i].y);
        }
        printf("%d %d %d\n",ansx,ansy,ansl);
        if(tt)printf("\n");
    }
}