看看就想到了二分答案,然後用前綴和驗證一下。
再一看數據:他喵的居然是10000,n*nlogn絕對爆炸。
經過一會兒的思索,我們發現草只分布在最多500個格子裏。
那一定有很多很多行和列是沒有草的,而我們需要做的就是去掉這些行列,這就是離散化。
(感覺很nb的樣子)
這樣的話,我們一開始就要存下每個有草的位置,然後排個序,找到每個草的相對位置,建立一個二維矩陣,並對這個矩陣求矩陣前綴和。預處理結束。接下來就是二分答案。對於每個長度,我們枚舉正方形左上角的點,然後找到他可以夠到的範圍,看是否滿足條件就行了。
PS:所有的原數據我們都存在了B數組裏面,而離散化後的值就是每個數值的相對排序。
#include <bits/stdc++.h>
using namespace std;
const int maxn=10005;
int c,n,x[maxn],y[maxn],b[maxn];
struct node
{
int x,y;
}a[maxn];int tot=0;
int s[maxn][maxn];
void lisan()
{
sort(b+1,b+1+tot);
tot=unique(b+1,b+1+tot)-b-1;
for(int i=1;i<=n;i++)
{
int tx=lower_bound(b+1,b+tot+1,a[i].x)-b;//相同的第一個
int ty=lower_bound(b+1,b+tot+1,a[i].y)-b;
s[tx][ty]++;
}
b[++tot]=10001;//
for(int i=1;i<=tot;i++)
{
for(int j=1;j<=tot;j++) s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+s[i][j];
}
}
int check(int x)
{
int p=upper_bound(b+1,b+tot+1,b[tot]-x+1)-b-1;
for(int i=1;i<=p;i++)
{
for(int j=1;j<=p;j++)
{
int tx=upper_bound(b+1,b+tot+1,b[i]+x-1)-b-1;//相同的最後一個
int ty=upper_bound(b+1,b+tot+1,b[j]+x-1)-b-1;
if(s[tx][ty]-s[i-1][ty]-s[tx][j-1]+s[i-1][j-1]>=c) return true;
}
}
return false;
}
int main()
{
scanf("%d%d",&c,&n);
for(int i=1;i<=n;i++)
{
scanf("%d%d",&a[i].x,&a[i].y);
b[++tot]=a[i].x;b[++tot]=a[i].y;
}
lisan();
int l=1,r=10000;int ans=0;
while(l<=r)
{
int mid=(l+r)/2;
if(check(mid))
{
ans=mid;r=mid-1;
}
else l=mid+1;
}
printf("%d",ans);
return 0;
}