LA 3695 Distant Galaxy (部分枚举+扫描法)

LA 3695 Distant Galaxy

题目大意:

一个平面上给出n(1<=n<=100)个点,找到一个矩形,使得矩形边界上的点最多.

题目分析:

若采用最简单的暴力,枚举四条边界位置+统计答案,时间复杂度为O(n^5).
若采用部分枚举思想,只O(n^2)枚举上下界,那么可以采用扫描法O(n),总时间复杂度为O(n^3),可行.
对于已确定的上下界ymin和ymax.需要用到三个数组lft,on,on2.lft[i]表示i号竖线以左(不包括当前i号竖线)在ymin和ymax上的点个数;on[i]表示i号竖线上点个数,但不包括ymin和ymax上的点;on2[i]与on[i]类似,也是表示i号竖线上的点个数,但是要包括ymin和ymax上的点.然后采用扫描法计算出当前情况下的最优解.

代码:

#include<cstdio>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn=100+10;

struct Point {
    int x,y;
    bool operator < (const Point& rhs) const {
        return x<rhs.x;
    }
    void input() {
        scanf("%d%d",&x,&y);
    }
}P[maxn];

int y[maxn],n;//y 所有点的纵坐标,用于离散化 
int lft[maxn],on[maxn],on2[maxn];
//在部分枚举确定上下界后 lft[i]表示i号竖线以左(不包括i)在边界上的点个数 
//on[i]表示i号竖线上点个数,不包括上下界上的点 on2[i]与on[i]类似,但是要包括上下界上的点 

int solve()
{
    sort(P,P+n); 
    sort(y,y+n); 
    int k=unique(y,y+n)-y,ans=0;//纵坐标离散 
    if(k<=2) return n;//若y只有两种,一定可以全部取完 
    for(int i=0;i<k;i++)
        for(int j=i+1;j<k;j++) {
            int ymin=y[i],ymax=y[j],m=0;
            for(int p=0;p<n;p++) {//找到所有竖线,递推更新lft,on,on2 
                if(p==0||P[p].x!=P[p-1].x) {//新的竖线 
                    ++m;
                    on[m]=on2[m]=0;
                    lft[m]=lft[m-1]+on2[m-1]-on[m-1];//lft[m]是不包括m竖线的 
                }
                if(P[p].y>ymin&&P[p].y<ymax) ++on[m];//on和on2的区别 
                if(P[p].y>=ymin&&P[p].y<=ymax) ++on2[m];
            }
            if(m<=2) return n;//若x只有两种,一定可以全部取完 
            int maxs=0;//maxs=max{on[p]-lft[p]|1<=p<=m} 即maxs表示左边的最大值 
            for(int p=1;p<=m;p++) { 
                ans=max(ans,lft[p]+on2[p]+maxs);
                maxs=max(maxs,on[p]-lft[p]);
            }
        }
    return ans;
}

int main()
{
    int kase=0;
    while(scanf("%d",&n)==1&&n) {
        for(int i=0;i<n;i++) P[i].input(),y[i]=P[i].y;
        printf("Case %d: %d\n",++kase,solve());
    }
    return 0;
}