AtCoder Beginner Contest 075 D - Axis-Parallel Rectangle

https://beta.atcoder.jp/contests/abc075/tasks/abc075_d

题意：

给出坐标平面上n个点的坐标，要求找到一个面积最小的矩形使得这个矩形的边界加上内部的点的数量大于等于k。

思路：

由于坐标过大，所以离散化。

离散化之后用前缀和，但是Orz求前缀和的时候写错了。

枚举左下角和右上角的坐标，分别为(i,j)和(s,e)。

那么点的数量为sum[s][e] - sum[s][j-1] - sum[i-1][q] + sum[i-1][j-1]，切记，切记，画图画图。

代码：

 #include <stdio.h>
 #include <string.h>
 #include <algorithm>
 #include <map>
 using namespace std;
 
 struct node
 {
     int x,y;
 } a[];
 
 map<int,int> mpx,mpy;
 map<int,int> invx,invy;
 
 int sum[][];
 int b[][];
 
 bool cmp1(node aa,node bb)
 {
     return aa.x < bb.x;
 }
 
 bool cmp2(node aa,node bb)
 {
     return aa.y < bb.y;
 }
 
 int main()
 {
     int n,k;
 
     scanf("%d%d",&n,&k);
 
     for (int i = ;i <= n;i++)
     {
         scanf("%d%d",&a[i].x,&a[i].y);
     }
 
     long long ans;
 
     sort(a+,a+n+,cmp1);
 
     for (int i = ;i <= n;i++)
     {
         invx[i] = a[i].x;
         mpx[a[i].x] = i;
     }
 
     long long dx = a[n].x - a[].x;
 
     sort(a+,a+n+,cmp2);
 
     for (int i = ;i <= n;i++)
     {
         invy[i] = a[i].y;
         mpy[a[i].y] = i;
     }
 
     long long dy = a[n].y - a[].y;
 
     ans = dx * dy;
 
     for (int i = ;i <= n;i++)
     {
         int x = a[i].x,y = a[i].y;
 
         b[mpx[x]][mpy[y]] = ;
     }
 
     for (int i = ;i <= n;i++)
         for (int j = ;j <= n;j++)
         {
             sum[i][j] = sum[i-][j] + sum[i][j-] - sum[i-][j-] + b[i][j];
         }
 
     for (int i = ;i <= n;i++)
         for (int j = ;j <= n;j++)
         {
             for (int s = i;s <= n;s++)
                 for (int q = j;q <= n;q++)
                 {
                     if (sum[s][q] - sum[s][j-] - sum[i-][q] + sum[i-][j-] >= k)
                     {
                         long long dx = invx[s] - invx[i];
                         long long dy = invy[q] - invy[j];
 
                         long long tmp = abs(dx) * abs(dy);
 
                         ans = min(ans,tmp);
                     }
                 }
         }
 
     printf("%lld\n",ans);
 
     return ;
 }