XSY3382

二分$c$，問題變成能否用一個(gè)長(zhǎng)$c$寬$c$的矩形框住至少$k$個(gè)點(diǎn)。

二維數(shù)點(diǎn)問題考慮用掃描線解決。將所有點(diǎn)按照$x$從小到大排序。

枚舉一段$x$坐標(biāo)相差不超過$c$的點(diǎn)（雙指針推進(jìn)），初始想法是根據(jù)這些點(diǎn)的$y$值建一棵權(quán)值線段樹，維護(hù)$y$值在某段區(qū)間內(nèi)的點(diǎn)的個(gè)數(shù)，然后判斷能否找到一段長(zhǎng)度為$c$的區(qū)間，包含了至少$k$個(gè)點(diǎn)。

但發(fā)現(xiàn)很難進(jìn)行下去，所以換一個(gè)想法，定義$a_i$為$y$值在區(qū)間$[i,i+c]$內(nèi)的點(diǎn)的個(gè)數(shù)，根據(jù)$y$值建權(quán)值線段樹，維護(hù)$a$的最大值，判斷$max\{a\}$是否大于等于$k$。

題外話：我是大傻逼，想到正解結(jié)果算錯(cuò)復(fù)雜度又否掉了。。。

#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=100010; struct Point{ll x,y;friend bool operator < (Point a,Point b){return a.x<b.x;} }p[N]; ll y[N];int tot; int n,k; ll l,r,ans,mid; int tag[N<<2],mx[N<<2]; int lft[N]; void add(int u,int l,int r,int ql,int qr,int v){if(ql<=l&&r<=qr){mx[u]+=v;tag[u]+=v;return;}int mid=(l+r)>>1;if(ql<=mid) add(u<<1,l,mid,ql,qr,v);if(qr>mid) add(u<<1|1,mid+1,r,ql,qr,v);mx[u]=max(mx[u<<1],mx[u<<1|1])+tag[u]; } bool check(ll c){memset(tag,0,sizeof(tag));memset(mx,0,sizeof(mx));for(int i=1,l=1;i<=tot;i++){while(y[i]-y[l]>c) l++;lft[i]=l;}for(int i=1,l=1,r=0;i<=n;i++){while(r<n&&p[r+1].x-p[i].x<=c){r++;add(1,1,tot,lft[p[r].y],p[r].y,1);} while(l<r&&p[l].x<p[i].x){add(1,1,tot,lft[p[l].y],p[l].y,-1);l++;} if(mx[1]>=k) return 1;if(r==n) return 0;}return 0; } int main(){scanf("%d%d",&n,&k);for(int i=1;i<=n;i++){scanf("%lld%lld",&p[i].x,&p[i].y);y[++tot]=p[i].y;}sort(p+1,p+n+1);sort(y+1,y+n+1);tot=unique(y+1,y+n+1)-y-1;for(int i=1;i<=n;i++){p[i].y=lower_bound(y+1,y+tot+1,p[i].y)-y;}l=0;r=4e9;while(l<=r){mid=(l+r)>>1;if(check(mid)){ans=mid;r=mid-1;}else l=mid+1;}printf("%lld\n",ans);return 0; }

總結(jié)

以上是生活随笔為你收集整理的[XSY3382] 专家系统（二分+线段树）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。