赶牛入圈.

赶牛入圈

农夫约翰希望为他的奶牛们建立一个畜栏。

这些挑剔的畜生要求畜栏必须是正方形的，而且至少要包含 C 单位的三叶草，来当做它们的下午茶。

畜栏的边缘必须与 X，Y 轴平行。

约翰的土地里一共包含 N 单位的三叶草，每单位三叶草位于一个 1×1 的土地区域内，区域位置由其左下角坐标表示，并且区域左下角的 X,Y 坐标都为整数，范围在 1 到 10000 以内。

多个单位的三叶草可能会位于同一个 1×1 的区域内，因为这个原因，在接下来的输入中，同一个区域坐标可能出现多次。

只有一个区域完全位于修好的畜栏之中，才认为这个区域内的三叶草在畜栏之中。

请你帮约翰计算一下，能包含至少 C 单位面积三叶草的情况下，畜栏的最小边长是多少。

输入格式

第一行输入两个整数 C 和 N。

接下来 N 行，每行输入两个整数 X 和 Y，代表三叶草所在的区域的 X,Y 坐标。

同一行数据用空格隔开。

输出格式

输出一个整数，代表畜栏的最小边长。

数据范围

1≤C≤500,
C≤N≤500

输入样例：

3 4
1 2
2 1
4 1
5 2

输出样例：

4

算法：二分，前缀和，离散化

题目分析：

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>

using namespace std;

typedef pair<int, int> PII;
const int N = 510, M = 1010;
PII points[N];  //总坐标数量不超过500个
vector<int> nums;    //x坐标与y坐标个数之和最多1000个
int sum[M][M];
int n, c;

int find(int x)
{
    int l = 1, r = nums.size() - 1;
    while(l < r)
    {
        int mid = l + r >> 1;
        if(nums[mid] >= x) r = mid;
        else l = mid + 1;
    }
    return r;
}

bool check(int k)
{
    for(int x1 = 1, x2 = 1; x2 < nums.size(); x2++)
    {
        while(nums[x2] + 1 - nums[x1] > k) x1++;
        for(int y1 = 1, y2 = 1; y2 < nums.size(); y2++)
        {
            while(nums[y2] + 1 - nums[y1] > k) y1++;
            if(sum[x2][y2] - sum[x1 - 1][y2] - sum[x2][y1 - 1] + sum[x1 - 1][y1 - 1] >= c)
            {
                return true;
            }
        }
    }
    return false;
}

int main()
{
    cin >> c >> n;
    nums.push_back(0);  //使之后二维离散化的矩阵下标从一开始
    for(int i = 0; i < n; i++)
    {
        int x, y;
        cin >> x >> y;
        points[i] = {x, y};
        nums.push_back(x);
        nums.push_back(y);
    }
    //离散化: 将1~10000范围内的x, y坐标离散化为1~nums.size()范围的坐标
    {
        sort(nums.begin(), nums.end());
        nums.erase(unique(nums.begin(), nums.end()), nums.end());
    }
    //求前缀和
    {
        for(int i = 0; i < n; i++)
        {
            int x = find(points[i].first), y = find(points[i].second);
            sum[x][y]++;
        }
        for(int i = 1; i < nums.size(); i++)
        {
            for(int j = 1; j < nums.size(); j++)
            {
                sum[i][j] += sum[i - 1][j] + sum[i][j - 1] - sum[i - 1][j - 1];
            }
        }
    }
    //二分判定边长是否合法
    {
        int l = 1, r = 10000;
        while(l < r)
        {
            int mid = l + r >> 1;
            if(check(mid)) r = mid;
            else l = mid + 1;
        }
        cout << r << endl;
    }
    return 0;
}