【poj3141】 Distant Galaxy

http://poj.org/problem?id=3141 (题目链接)

题意：给出平面上n个点，找出一个矩形，使边界上包含尽量多的点。

solution
　　不难发现，除非所有输入点都在同一行或同一列上，最优矩形的4条边上都至少有一个点。这样的话，我们可以枚举四条边穿过的点，然后统计点数。
　　考虑部分枚举，只枚举矩形上下界，用其它方法确定左右界。
　　设一条竖线i，用left[i]表示竖线左边位于上下界上单点数（不统计位于该竖线上的点），on[i]和on2[i]表示竖线上位于上下边界之间的点数（on[i]不统计位于上下界上的点，而on2[i]要统计）。这样，当左右边界分别为i，j时，矩形边界上的点数为left[j]-left[i]+on[i]+on2[j]。当右边界j确定时，on[i]-left[i]应最大。
　　枚举完上下边界后，我们花O(n)的时间按照从左到右的顺序for一遍所有的的点，计算left,on[i],on2[i]，然后枚举右边界j，同时维护on[i]-left[i]的最大值。

代码：

// poj3141
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 2147483640
#define Pi 3.1415926535898
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;

struct point {int x,y;}a[10010];

int y[10010],on[10010],on2[10010],L[10010],n;

bool cmp(point a,point b) {return a.x<b.x;}
int main() {
    int T=0;
    while (scanf("%d",&n)!=EOF && n) {
        T++;
        printf("Case %d: ",T);
        for (int i=0;i<n;i++) {scanf("%d%d",&a[i].x,&a[i].y);y[i]=a[i].y;}
        sort(a,a+n,cmp);
        sort(y,y+n);
        int m=unique(y,y+n)-y;
        if (m<=2) {printf("%d\n",n);continue;}
        int ans=0;
        for (int l=0;l<m;l++)
            for (int r=l+1;r<m;r++) {
                int ymin=y[l],ymax=y[r],k=0;
                memset(L,0,sizeof(L));
                for (int i=0;i<n;i++) {
                    if (i==0 || a[i].x!=a[i-1].x) {
                        k++;
                        on[k]=on2[k]=0;
                        L[k]=L[k-1]+on2[k-1]-on[k-1];
                    }
                    if (a[i].y>ymin && a[i].y<ymax) on[k]++;
                    if (a[i].y>=ymin && a[i].y<=ymax) on2[k]++;
                }
                int M=0;
                for (int j=1;j<=k;j++) {
                    ans=max(ans,L[j]+on2[j]+M);
                    M=max(M,on[j]-L[j]);
                }
            }
        printf("%d\n",ans);
    }
    return 0;
}