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");
}
}