uva 1382 - Distant Galaxy

点击打开链接uva 1382

题意：给出平面上的n个点，找出一个矩形，使得边界上含有尽量多的点
思路：
1 很清楚，如果输入的n个点在同一行或者同一列的话那么ans = n。还有一种情况就是n个点的横坐标和纵坐标只有2种，那么这种情况ans = n。
2 对于这一题我们考虑的是枚举矩形的上下边界（纵坐标），然后利用其它的方法求左右边界，见下图

3 对于竖线i，我们用left[i]表示竖线左边位于上下边界的点数（不包括位于竖线i）, on[i]表示竖线上位于上下边界之间的点数（和on2[i]的区别就是on[i]不统计位于上下边界的点数），这样给定左右边界i和j的话，矩形边界上的点数为left[j]+on2[j]+on[i]-left[i],当有边界j确定的时候，on[i]-left[i]要最大
4 那么我们枚举完上下边界后，我们利用O(n)的时间去求left[] , on[] , on2[]，然后枚举有边界j，维护最大的on[i]-left[i]

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int MAXN = 110;
struct Node{
    int x;
    int y;
    bool operator<(const Node& s)const{
        return x < s.x; 
    }
};
Node node[MAXN];
int n , numy[MAXN];

int solve(){
    sort(node , node+n);
    sort(numy , numy+n);
    int num = unique(numy , numy+n)-numy;
    if(num <= 2)//如果纵坐标最多只有2个那么ans = n
       return n;
    //枚举上下界
    int miny , maxy , ans;
    int left[MAXN] , on[MAXN] , on2[MAXN];
    ans = 0;
    for(int i = 0 ; i < num ; i++){
       for(int j = i+1 ; j < num ; j++){
           miny = numy[i];
           maxy = numy[j];
           //求left , on , on2数组;
           int k = -1;
           memset(left , 0 , sizeof(left));
           memset(on , 0 , sizeof(on));
           memset(on2 , 0 , sizeof(on2));
           for(int t = 0 ; t < n ; t++){
              if(!t || node[t].x != node[t-1].x){
                 k++;
                 left[k] = k == 0 ? 0 : left[k-1]+on2[k-1]-on[k-1]; 
              }
              if(node[t].y > miny && node[t].y < maxy)
                 on[k]++;
              if(node[t].y >= miny && node[t].y <= maxy)
                 on2[k]++;
           }
           if(k <= 1)//如果横坐标最多只有2个那么ans = n
              return n;
           int Max = 0;
           for(int t = 0 ; t <= k ; t++){
              ans = max(ans , left[t]+on2[t]+Max); 
              Max = max(Max , on[t]-left[t]); 
           }
       }
    }
    return ans;
}

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