uva 1382 Distant Galaxy

https://vjudge.net/problem/UVA-1382

题意：

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

思路：

参考训练指南。

首先如果枚举起起点，终点，再进行统计的话，复杂度显然太大了，因为坐标的范围是10的9次方级别的。

所以，又想到扫描线的方法。

我们首先对数组排个序，枚举上下边界，之后，用三个数组，left[i]表示i这条线左边（不包括这条线）的在上下边界的点的数量，on[i]表示在i这条线上的位于上下边界之间的点的数量（不包括上下边界），on2[i]表示，

在i这条线上的位于上下边界之间的点的数量（包括上下边界），我们从左往右扫一遍，维护这三个数组的信息。

left[i] = left[i-1] + on2[i] - on[i-1],在i-1这条边的左边的在上下边界上的点的数量等于它上一条边的左边的在上下边界上的点的数量加上上一条边的在上下边界上的点的数量，on2[i-1] - on[i-1]就是位于上下边界的点的数量。

ans = max(left[i] + on2[i]+on[k] - left[k])(k < i)，我们要维护的就是on[k] - left[k]的最大值，在从左往右扫的时候维护。

一些细节的处理需要注意，比如点都在一条边或者两条边上。

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;

struct node
{
    int x,y;
} a[105];

int n,m;

int y[105];
int left[105],on[105],on2[105];

bool cmp(node aa,node bb)
{
    return aa.x < bb.x;
}

int solve(void)
{
    sort(a,a+n,cmp);

    int ans = 0;

    sort(y,y+n);

    m = unique(y,y+n) - y;

    if (m <= 2) return n;

    for (int i = 0;i < m;i++)
        for (int j = i + 1;j < m;j++)
    {
        int ymin = y[i],ymax = y[j];

        int num = 0;

        for (int k = 0;k < n;k++)
        {
            if (k == 0 || a[k].x != a[k-1].x)
            {
                num++;

                on[num] = on2[num] = 0;

                left[num] = num == 0 ? 0 : left[num-1] + on2[num-1] - on[num-1];
            }


            if (a[k].y > ymin && a[k].y < ymax) on[num]++;
            if (a[k].y >= ymin && a[k].y <= ymax) on2[num]++;
        }

        if (num <= 2) return n;

        int M = 0;

        for (int k = 1;k <= num;k++)
        {
            ans = max(ans,M + left[k] + on2[k]);
            M = max(M,on[k] - left[k]);
        }
    }

    return ans;
}

int main()
{
    int kase = 0;

    while (scanf("%d",&n) != EOF)
    {
        if (!n) break;

        for (int i = 0;i < n;i++)
        {
            scanf("%d%d",&a[i].x,&a[i].y);

            y[i] = a[i].y;
        }

        int ans = solve();

        printf("Case %d: %d\n",++kase,ans);
    }

    return 0;
}