[学习笔记]K-D Tree
以前其实学过的但是不会拍扁重构……所以这几天学了一下 \(K-D\ Tree\) 的正确打开姿势。
\(K\) 维 \(K-D\ Tree\) 的单次操作最坏时间复杂度为 \(O(k\times n^{1-\frac 1k})\)
1、[Violet]天使玩偶/SJY摆棋子
二维 \(K-D\ Tree\)。
不过要暴力重构。。。我终于会了，但不开 \(O_2\) 过不去。。。
\(Code\ Below:\)
// luogu-judger-enable-o2
#include <bits/stdc++.h>
using namespace std;
const int maxn=500000+10;
const int inf=0x3f3f3f3f;
const double alpha=0.75;
int n,m,D,rt,tot,ans,cnt,rub[maxn<<1],top;
struct node{
int d[2];
}a[maxn<<1];
struct KD_Tree{
int d[2],Max[2],Min[2],ch[2],siz;
inline void init(){
d[0]=d[1]=Max[0]=Max[1]=Min[0]=Min[1]=ch[0]=ch[1]=siz=0;
}
inline void get(node a){
Max[0]=Min[0]=d[0]=a.d[0];
Max[1]=Min[1]=d[1]=a.d[1];
}
}t[maxn<<1];
inline bool operator < (const node &a,const node &b){
return a.d[D]<b.d[D];
}
inline int read(){
register int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return (f==1)?x:-x;
}
inline int newnode(){
int x=top?rub[top--]:++cnt;
t[x].init();return x;
}
inline void update(int x,int y){
t[x].Max[0]=max(t[x].Max[0],t[y].Max[0]);
t[x].Max[1]=max(t[x].Max[1],t[y].Max[1]);
t[x].Min[0]=min(t[x].Min[0],t[y].Min[0]);
t[x].Min[1]=min(t[x].Min[1],t[y].Min[1]);
}
inline void pushup(int x){
t[x].siz=t[t[x].ch[0]].siz+t[t[x].ch[1]].siz+1;
if(t[x].ch[0]) update(x,t[x].ch[0]);
if(t[x].ch[1]) update(x,t[x].ch[1]);
}
int build(int l,int r,int now){
int mid=(l+r)>>1,x=newnode();D=now;
nth_element(a+l,a+mid,a+r+1);
t[x].get(a[mid]);
if(l<mid) t[x].ch[0]=build(l,mid-1,now^1);
if(mid<r) t[x].ch[1]=build(mid+1,r,now^1);
pushup(x);return x;
}
inline void del(int x){
if(t[x].ch[0]) del(t[x].ch[0]);
a[++tot]=(node){t[x].d[0],t[x].d[1]};
rub[++top]=x;
if(t[x].ch[1]) del(t[x].ch[1]);
}
inline void check(int &x,int now){
if(1.0*max(t[t[x].ch[0]].siz,t[t[x].ch[1]].siz)>1.0*alpha*t[x].siz)
tot=0,del(x),x=build(1,tot,now);
}
void insert(int &x,node c,int now){
if(!x){
x=newnode();
t[x].get(c);
return ;
}
if(c.d[now]<=t[x].d[now]) insert(t[x].ch[0],c,now^1);
else insert(t[x].ch[1],c,now^1);
pushup(x);check(x,now);
}
inline int dis(int x,int X,int Y){
return max(t[x].Min[0]-X,0)+max(X-t[x].Max[0],0)+max(t[x].Min[1]-Y,0)+max(Y-t[x].Max[1],0);
}
void query(int x,int X,int Y){
int L=abs(t[x].d[0]-X)+abs(t[x].d[1]-Y),d[2],k;
d[0]=t[x].ch[0]?dis(t[x].ch[0],X,Y):inf;
d[1]=t[x].ch[1]?dis(t[x].ch[1],X,Y):inf;
ans=min(ans,L);k=d[0]>d[1];
if(d[k]<ans) query(t[x].ch[k],X,Y);
if(d[k^1]<ans) query(t[x].ch[k^1],X,Y);
}
int main()
{
n=read(),m=read();tot=n;
for(int i=1;i<=n;i++) a[i].d[0]=read(),a[i].d[1]=read();
rt=build(1,n,0);
int op,x,y;
while(m--){
op=read(),x=read(),y=read();
if(op==1) insert(rt,(node){x,y},0);
else ans=inf,query(rt,x,y),printf("%d\n",ans);
}
return 0;
}
\(CDQ\) 分治解法：
我们做 \(4\) 次 \(CDQ\) 分治，每次只计算左下角的点，将问题转换成在树状数组上查最小。
\(Code\ Below:\)
// luogu-judger-enable-o2
#include <bits/stdc++.h>
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int maxn=500000+10;
const int maxm=1000000+10;
const int inf=0x3f3f3f3f;
int n,m,X,Y,c[maxm],ans[maxn],cnt,tot;
struct Query{
int op,x,y,id;
}p[maxn<<1],q[maxn<<1],t[maxn<<1];
namespace IO{
#define gc() (iS==iT?(iT=(iS=ibuff)+fread(ibuff,1,SIZ,stdin),(iS==iT?EOF:*iS++)):*iS++)
const int SIZ=1<<21|1;
char *iS,*iT,ibuff[SIZ],obuff[SIZ],*oS=obuff,*oT=oS+SIZ-1,fu[110],c;int fr;
inline void out(){
fwrite(obuff,1,oS-obuff,stdout);
oS=obuff;
}
template <class T>
inline void read(T &x){
x=0;T y=1;
for(c=gc();(c>'9'||c<'0')&&c!='-';c=gc());
c=='-'?y=-1:x=(c&15);
for(c=gc();c>='0'&&c<='9';c=gc()) x=x*10+(c&15);
x*=y;
}
template <class T>
inline void print(T x,char text='\n'){
if(x<0) *oS++='-',x*=-1;
if(x==0) *oS++='0';
while(x) fu[++fr]=x%10+'0',x/=10;
while(fr) *oS++=fu[fr--];
*oS++=text;out();
}
}
inline void update(int x,int y){
for(;x<=Y;x+=lowbit(x)) c[x]=max(c[x],y);
}
inline int query(int x){
int ans=0;
for(;x;x-=lowbit(x)) ans=max(ans,c[x]);
return ans;
}
inline void clear(int x){
for(;x<=Y;x+=lowbit(x)){
if(c[x]) c[x]=0;
else break;
}
}
void CDQ(int l,int r){
if(l==r) return ;
int mid=(l+r)>>1;
CDQ(l,mid);CDQ(mid+1,r);
int u=l,v=mid+1,j=l,tmp;
for(int i=mid+1;i<=r;i++)
if(q[i].op==2){
for(;j<=mid&&q[j].x<=q[i].x;j++)
if(q[j].op==1) update(q[j].y,q[j].x+q[j].y);
tmp=query(q[i].y);
if(tmp) ans[q[i].id]=min(ans[q[i].id],q[i].x+q[i].y-tmp);
}
for(int i=l;i<j;i++)
if(q[i].op==1) clear(q[i].y);
j=l;
while(u<=mid&&v<=r){
if(q[u].x<=q[v].x) t[j++]=q[u++];
else t[j++]=q[v++];
}
while(u<=mid) t[j++]=q[u++];
while(v<=r) t[j++]=q[v++];
for(int i=l;i<=r;i++) q[i]=t[i];
}
void del(){
int X=0,Y=0,cnt=0;
for(int i=1;i<=tot;i++)
if(q[i].op==2) X=max(X,q[i].x),Y=max(Y,q[i].y);
for(int i=1;i<=tot;i++)
if(q[i].x<=X&&q[i].y<=Y) q[++cnt]=q[i];
CDQ(1,cnt);
}
int main()
{
IO::read(n),IO::read(m);
int op,x,y;
for(int i=1;i<=n;i++){
IO::read(x),IO::read(y);
x++;y++;
p[++tot]=(Query){1,x,y,0};
X=max(X,x);Y=max(Y,y);
}
for(int i=1;i<=m;i++){
IO::read(op),IO::read(x),IO::read(y);
x++;y++;
if(op==1) p[++tot]=(Query){1,x,y,0};
else p[++tot]=(Query){2,x,y,++cnt};
X=max(X,x);Y=max(Y,y);
}
X++;Y++;
for(int i=1;i<=cnt;i++) ans[i]=inf;
for(int i=1;i<=tot;i++) q[i]=p[i];
del();
for(int i=1;i<=tot;i++) q[i]=p[i],q[i].x=X-p[i].x;
del();
for(int i=1;i<=tot;i++) q[i]=p[i],q[i].y=Y-p[i].y;
del();
for(int i=1;i<=tot;i++) q[i]=p[i],q[i].x=X-p[i].x,q[i].y=Y-p[i].y;
del();
for(int i=1;i<=cnt;i++) IO::print(ans[i]);
return 0;
}
2、[CQOI2016]K远点对
我们用一个小根堆来计算第 \(k\) 小。
\(Code\ Below:\)
#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=100000+10;
const int inf=0x3f3f3f3f;
int n,k,D,rt,cnt;
priority_queue<ll,vector<ll>,greater<ll> > pq;
struct node{
int d[2];
}a[maxn];
struct KD_Tree{
int d[2],Max[2],Min[2],ch[2],siz;
inline void init(){
d[0]=d[1]=Max[0]=Max[1]=Min[0]=Min[1]=ch[0]=ch[1]=siz=0;
}
inline void get(node a){
Max[0]=Min[0]=d[0]=a.d[0];
Max[1]=Min[1]=d[1]=a.d[1];
}
}t[maxn];
inline bool operator < (const node &a,const node &b){
return a.d[D]<b.d[D];
}
inline int read(){
register int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return (f==1)?x:-x;
}
inline void update(int x,int y){
t[x].Max[0]=max(t[x].Max[0],t[y].Max[0]);
t[x].Max[1]=max(t[x].Max[1],t[y].Max[1]);
t[x].Min[0]=min(t[x].Min[0],t[y].Min[0]);
t[x].Min[1]=min(t[x].Min[1],t[y].Min[1]);
}
inline void pushup(int x){
t[x].siz=t[t[x].ch[0]].siz+t[t[x].ch[1]].siz+1;
if(t[x].ch[0]) update(x,t[x].ch[0]);
if(t[x].ch[1]) update(x,t[x].ch[1]);
}
int build(int l,int r,int now){
int mid=(l+r)>>1,x=++cnt;D=now;
nth_element(a+l,a+mid,a+r+1);
t[x].get(a[mid]);
if(l<mid) t[x].ch[0]=build(l,mid-1,now^1);
if(mid<r) t[x].ch[1]=build(mid+1,r,now^1);
pushup(x);return x;
}
inline ll dis(int x,int X,int Y){
return max(1ll*(t[x].Min[0]-X)*(t[x].Min[0]-X),1ll*(t[x].Max[0]-X)*(t[x].Max[0]-X))
+max(1ll*(t[x].Min[1]-Y)*(t[x].Min[1]-Y),1ll*(t[x].Max[1]-Y)*(t[x].Max[1]-Y));
}
inline ll dist(int x,int X,int Y){
return 1ll*(t[x].d[0]-X)*(t[x].d[0]-X)+1ll*(t[x].d[1]-Y)*(t[x].d[1]-Y);
}
void query(int x,int X,int Y){
ll L=dist(x,X,Y),d[2];int k;
d[0]=t[x].ch[0]?dis(t[x].ch[0],X,Y):-inf;
d[1]=t[x].ch[1]?dis(t[x].ch[1],X,Y):-inf;
if(L>pq.top()) pq.pop(),pq.push(L);
k=d[1]>d[0];
if(d[k]>pq.top()) query(t[x].ch[k],X,Y);
if(d[k^1]>pq.top()) query(t[x].ch[k^1],X,Y);
}
int main()
{
n=read(),k=read();
for(int i=1;i<=n;i++) a[i].d[0]=read(),a[i].d[1]=read();
rt=build(1,n,0);
for(int i=1;i<=k*2;i++) pq.push(0ll);
for(int i=1;i<=n;i++) query(rt,a[i].d[0],a[i].d[1]);
printf("%lld\n",pq.top());
return 0;
}
3、[国家集训队]JZPFAR
我们在 \(K\) 远点对的基础上再记一个 \(minid\)，一样暴力查询。
\(Code\ Below:\)
#include <bits/stdc++.h>
#define ll long long
#define pli pair<ll,int>
#define mp make_pair
#define F first
#define S second
using namespace std;
const int maxn=100000+10;
const int inf=0x3f3f3f3f;
int n,m,D,rt,cnt;
priority_queue<pli,vector<pli>,greater<pli> > pq;
struct node{
int d[2],id;
}a[maxn];
struct KD_Tree{
int d[2],id,ID,Max[2],Min[2],ch[2],siz;
inline void init(){
d[0]=d[1]=id=ID=Max[0]=Max[1]=Min[0]=Min[1]=ch[0]=ch[1]=siz=0;
}
inline void get(node a){
Max[0]=Min[0]=d[0]=a.d[0];
Max[1]=Min[1]=d[1]=a.d[1];
id=ID=a.id;
}
}t[maxn];
inline bool operator < (const node &a,const node &b){
return (a.d[D]!=b.d[D])?a.d[D]<b.d[D]:a.id>b.id;
}
inline int read(){
register int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return (f==1)?x:-x;
}
inline void update(int x,int y){
t[x].Max[0]=max(t[x].Max[0],t[y].Max[0]);
t[x].Max[1]=max(t[x].Max[1],t[y].Max[1]);
t[x].Min[0]=min(t[x].Min[0],t[y].Min[0]);
t[x].Min[1]=min(t[x].Min[1],t[y].Min[1]);
}
inline void pushup(int x){
t[x].siz=t[t[x].ch[0]].siz+t[t[x].ch[1]].siz+1;
t[x].ID=min(t[x].ID,min(t[t[x].ch[0]].ID,t[t[x].ch[1]].ID));
if(t[x].ch[0]) update(x,t[x].ch[0]);
if(t[x].ch[1]) update(x,t[x].ch[1]);
}
int build(int l,int r,int now){
int mid=(l+r)>>1,x=++cnt;D=now;
nth_element(a+l,a+mid,a+r+1);
t[x].get(a[mid]);
if(l<mid) t[x].ch[0]=build(l,mid-1,now^1);
if(mid<r) t[x].ch[1]=build(mid+1,r,now^1);
pushup(x);return x;
}
inline ll dis(int x,int X,int Y){
return max(1ll*(t[x].Min[0]-X)*(t[x].Min[0]-X),1ll*(t[x].Max[0]-X)*(t[x].Max[0]-X))
+max(1ll*(t[x].Min[1]-Y)*(t[x].Min[1]-Y),1ll*(t[x].Max[1]-Y)*(t[x].Max[1]-Y));
}
inline ll dist(int x,int X,int Y){
return 1ll*(t[x].d[0]-X)*(t[x].d[0]-X)+1ll*(t[x].d[1]-Y)*(t[x].d[1]-Y);
}
void query(int x,int X,int Y){
ll L=dist(x,X,Y),d[2];int k;
d[0]=t[x].ch[0]?dis(t[x].ch[0],X,Y):-inf;
d[1]=t[x].ch[1]?dis(t[x].ch[1],X,Y):-inf;
if((L>pq.top().F)||(L==pq.top().F&&t[x].id<-pq.top().S)) pq.pop(),pq.push(mp(L,-t[x].id));
k=d[1]>d[0];
if((d[k]>pq.top().F)||(d[k]==pq.top().F&&t[t[x].ch[k]].ID<-pq.top().S)) query(t[x].ch[k],X,Y);
if((d[k^1]>pq.top().F)||(d[k^1]==pq.top().F&&t[t[x].ch[k^1]].ID<-pq.top().S)) query(t[x].ch[k^1],X,Y);
}
int main()
{
n=read();
for(int i=1;i<=n;i++) a[i].d[0]=read(),a[i].d[1]=read(),a[i].id=i;
rt=build(1,n,0);
m=read();
int x,y,k;
while(m--){
x=read(),y=read(),k=read();
for(int i=1;i<=k;i++) pq.push(mp(0ll,-inf));
query(rt,x,y);
printf("%d\n",-pq.top().S);
while(!pq.empty()) pq.pop();
}
return 0;
}
4、[SDOI2010]捉迷藏
最远点对最近点对问题，经典 \(K-D\ Tree\) 问题。不过最近点对要除去自己。
\(Code\ Below:\)
// luogu-judger-enable-o2
#include <bits/stdc++.h>
using namespace std;
const int maxn=100000+10;
const int inf=0x3f3f3f3f;
int n,D,rt,ans=inf,ans1,ans2;
struct node{
int d[2];
}a[maxn];
inline bool operator < (const node &a,const node &b){
return a.d[D]<b.d[D];
}
struct KD_Tree{
int d[2],Max[2],Min[2],ch[2];
inline void init(){
d[0]=d[1]=Max[0]=Max[1]=Min[0]=Min[1]=ch[0]=ch[1]=0;
}
inline void get(node a){
Max[0]=Min[0]=d[0]=a.d[0];
Max[1]=Min[1]=d[1]=a.d[1];
}
}t[maxn];
inline int read(){
register int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return (f==1)?x:-x;
}
inline void update(int x,int y){
t[x].Max[0]=max(t[x].Max[0],t[y].Max[0]);
t[x].Max[1]=max(t[x].Max[1],t[y].Max[1]);
t[x].Min[0]=min(t[x].Min[0],t[y].Min[0]);
t[x].Min[1]=min(t[x].Min[1],t[y].Min[1]);
}
inline void pushup(int x){
if(t[x].ch[0]) update(x,t[x].ch[0]);
if(t[x].ch[1]) update(x,t[x].ch[1]);
}
int build(int l,int r,int now){
int mid=(l+r)>>1,x=mid;D=now;
nth_element(a+l,a+mid,a+r+1);
t[x].init();t[x].get(a[mid]);
if(l<mid) t[x].ch[0]=build(l,mid-1,now^1);
if(mid<r) t[x].ch[1]=build(mid+1,r,now^1);
pushup(x);return x;
}
inline int maxdis(int x,int X,int Y){
return max(abs(t[x].Min[0]-X),abs(t[x].Max[0]-X))+max(abs(t[x].Min[1]-Y),abs(t[x].Max[1]-Y));
}
inline int mindis(int x,int X,int Y){
return max(t[x].Min[0]-X,0)+max(X-t[x].Max[0],0)+max(t[x].Min[1]-Y,0)+max(Y-t[x].Max[1],0);
}
inline int dis(int x,int X,int Y){
return abs(t[x].d[0]-X)+abs(t[x].d[1]-Y);
}
void query_max(int x,int X,int Y){
int L=dis(x,X,Y),d[2],k;
d[0]=t[x].ch[0]?maxdis(t[x].ch[0],X,Y):-inf;
d[1]=t[x].ch[1]?maxdis(t[x].ch[1],X,Y):-inf;
ans1=max(ans1,L);k=d[1]>d[0];
if(d[k]>ans1) query_max(t[x].ch[k],X,Y);
if(d[k^1]>ans1) query_max(t[x].ch[k^1],X,Y);
}
void query_min(int x,int X,int Y){
int L=dis(x,X,Y),d[2],k;
d[0]=t[x].ch[0]?mindis(t[x].ch[0],X,Y):inf;
d[1]=t[x].ch[1]?mindis(t[x].ch[1],X,Y):inf;
if(L) ans2=min(ans2,L);k=d[1]<d[0];
if(d[k]<ans2) query_min(t[x].ch[k],X,Y);
if(d[k^1]<ans2) query_min(t[x].ch[k^1],X,Y);
}
int main()
{
n=read();
for(int i=1;i<=n;i++) a[i].d[0]=read(),a[i].d[1]=read();
rt=build(1,n,0);
for(int i=1;i<=n;i++){
ans1=0;ans2=inf;
query_max(rt,a[i].d[0],a[i].d[1]);
query_min(rt,a[i].d[0],a[i].d[1]);
ans=min(ans,ans1-ans2);
}
printf("%d\n",ans);
return 0;
}
5、Generating Synergy
带下传标记的 \(K-D\ Tree\)。修改的时候暴力修改，查询的时候要将其祖先的标记下传。
\(Code\ Below:\)
// luogu-judger-enable-o2
#include <bits/stdc++.h>
using namespace std;
const int maxn=100000+10;
const int inf=0x3f3f3f3f;
int n,D,rt,ans=inf,ans1,ans2;
struct node{
int d[2];
}a[maxn];
inline bool operator < (const node &a,const node &b){
return a.d[D]<b.d[D];
}
struct KD_Tree{
int d[2],Max[2],Min[2],ch[2];
inline void init(){
d[0]=d[1]=Max[0]=Max[1]=Min[0]=Min[1]=ch[0]=ch[1]=0;
}
inline void get(node a){
Max[0]=Min[0]=d[0]=a.d[0];
Max[1]=Min[1]=d[1]=a.d[1];
}
}t[maxn];
inline int read(){
register int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return (f==1)?x:-x;
}
inline void update(int x,int y){
t[x].Max[0]=max(t[x].Max[0],t[y].Max[0]);
t[x].Max[1]=max(t[x].Max[1],t[y].Max[1]);
t[x].Min[0]=min(t[x].Min[0],t[y].Min[0]);
t[x].Min[1]=min(t[x].Min[1],t[y].Min[1]);
}
inline void pushup(int x){
if(t[x].ch[0]) update(x,t[x].ch[0]);
if(t[x].ch[1]) update(x,t[x].ch[1]);
}
int build(int l,int r,int now){
int mid=(l+r)>>1,x=mid;D=now;
nth_element(a+l,a+mid,a+r+1);
t[x].init();t[x].get(a[mid]);
if(l<mid) t[x].ch[0]=build(l,mid-1,now^1);
if(mid<r) t[x].ch[1]=build(mid+1,r,now^1);
pushup(x);return x;
}
inline int maxdis(int x,int X,int Y){
return max(abs(t[x].Min[0]-X),abs(t[x].Max[0]-X))+max(abs(t[x].Min[1]-Y),abs(t[x].Max[1]-Y));
}
inline int mindis(int x,int X,int Y){
return max(t[x].Min[0]-X,0)+max(X-t[x].Max[0],0)+max(t[x].Min[1]-Y,0)+max(Y-t[x].Max[1],0);
}
inline int dis(int x,int X,int Y){
return abs(t[x].d[0]-X)+abs(t[x].d[1]-Y);
}
void query_max(int x,int X,int Y){
int L=dis(x,X,Y),d[2],k;
d[0]=t[x].ch[0]?maxdis(t[x].ch[0],X,Y):-inf;
d[1]=t[x].ch[1]?maxdis(t[x].ch[1],X,Y):-inf;
ans1=max(ans1,L);k=d[1]>d[0];
if(d[k]>ans1) query_max(t[x].ch[k],X,Y);
if(d[k^1]>ans1) query_max(t[x].ch[k^1],X,Y);
}
void query_min(int x,int X,int Y){
int L=dis(x,X,Y),d[2],k;
d[0]=t[x].ch[0]?mindis(t[x].ch[0],X,Y):inf;
d[1]=t[x].ch[1]?mindis(t[x].ch[1],X,Y):inf;
if(L) ans2=min(ans2,L);k=d[1]<d[0];
if(d[k]<ans2) query_min(t[x].ch[k],X,Y);
if(d[k^1]<ans2) query_min(t[x].ch[k^1],X,Y);
}
int main()
{
n=read();
for(int i=1;i<=n;i++) a[i].d[0]=read(),a[i].d[1]=read();
rt=build(1,n,0);
for(int i=1;i<=n;i++){
ans1=0;ans2=inf;
query_max(rt,a[i].d[0],a[i].d[1]);
query_min(rt,a[i].d[0],a[i].d[1]);
ans=min(ans,ans1-ans2);
}
printf("%d\n",ans);
return 0;
}
[学习笔记]K-D Tree的更多相关文章
- [学习笔记]Dsu On Tree
[dsu on tree][学习笔记] - Candy? - 博客园 题单: 也称:树上启发式合并 可以解决绝大部分不带修改的离线询问的子树查询问题 流程: 1.重链剖分找重儿子 2.sol:全局用桶 ...
- 决策树学习笔记（Decision Tree）
什么是决策树? 决策树是一种基本的分类与回归方法.其主要有点事模型具有可得性,分类速度快.学习时,利用训练数据,根据损失函数最小化原则建立决策树模型:预测时,对新数据,利用决策树模型进行分类. 决策树 ...
- [学习笔记] Uplift Decision Tree With KL Divergence
Uplift Decision Tree With KL Divergence Intro Uplift model 我没找到一个合适的翻译,这方法主要应用是,探究用户在给予一定激励之后的表现,也就是 ...
- 学习笔记——k近邻法
对新的输入实例,在训练数据集中找到与该实例最邻近的\(k\)个实例,这\(k\)个实例的多数属于某个类,就把该输入实例分给这个类. \(k\) 近邻法(\(k\)-nearest neighbor, ...
- 【学习笔记】K-D tree 区域查询时间复杂度简易证明
查询算法的流程 如果查询与当前结点的区域无交集,直接跳出. 如果查询将当前结点的区域包含,直接跳出并上传答案. 有交集但不包含,继续递归求解. K-D Tree 如何划分区域 可以借助下文图片理解. ...
- R语言学习笔记—K近邻算法
K近邻算法(KNN)是指一个样本如果在特征空间中的K个最相邻的样本中的大多数属于某一个类别,则该样本也属于这个类别,并具有这个类别上样本的特性.即每个样本都可以用它最接近的k个邻居来代表.KNN算法适 ...
- [学习笔记]k短路
A*:我已经忘了怎么写了,反正n=30,m=1000都能卡掉... 正解:可持久化左偏树+堆维护可能集合 原论文:http://www.docin.com/p-1387370338.html 概括: ...
- [ML学习笔记] 决策树与随机森林（Decision Tree&Random Forest）
[ML学习笔记] 决策树与随机森林(Decision Tree&Random Forest) 决策树 决策树算法以树状结构表示数据分类的结果.每个决策点实现一个具有离散输出的测试函数,记为分支 ...
- 树上启发式合并（dsu on tree）学习笔记
有丶难,学到自闭 参考的文章: zcysky:[学习笔记]dsu on tree Arpa:[Tutorial] Sack (dsu on tree) 先康一康模板题吧:CF 600E($Lomsat ...
- 珂朵莉树(Chtholly Tree)学习笔记
珂朵莉树(Chtholly Tree)学习笔记 珂朵莉树原理 其原理在于运用一颗树(set,treap,splay......)其中要求所有元素有序,并且支持基本的操作(删除,添加,查找......) ...
随机推荐
- 初识kbmmw 中的ORM
在kbmmw 5.02.1 中,加入了ORM 的功能(这里可能和其他语言的定义不完全一样),我们就简单的认为 它就是一个类与数据库的转换吧.今天就先介绍一下如何通过kbmmw 的ORM 功能,实现类与 ...
- Java语法基础课 原码 反码 补码
原码就是符号位加上真值的绝对值, 即用第一位表示符号, 其余位表示值. 反码的表示方法是:正数的反码是其本身:负数的反码是在其原码的基础上, 符号位不变,其余各个位取反. 补码的表示方法是在反码的基础 ...
- boost--signal
1.signals2库 signals2库实现了线程安全的观察者模式,在signals2中观察者模式被称为信号/插槽(signals/slots),它是一种函数回调机制.一个信号可以关联一个或多个插槽 ...
- C语言常用数据类型说明
1.取值范围: short一般占两个字节,取值范围:-32768 - 32767 int一般占两个或四个字节,取值范围:-2147483648 - 2147483647 unsigned int ...
- 使用Docker镜像
1 使用Docker镜像 1.1 获取镜像 命令格式:docker pull NAME[:TAG] NAME为镜像仓库的名称 TAG为镜像的标签(表示版本号) 描述一个镜像需要包括:名称+ ...
- jQuery length 和 size()区别
jQuery length和size()区别总结如下: 1.length是属性,size()是方法. 2.如果你只是想获取元素的个数,两者效果一样既 $("img").length ...
- 20155205 信息安全技术第二次实验 Windows口令破解
20155205 信息安全技术第二次实验 Windows口令破解 实验原理 一.口令破解方法 口令破解主要有两种方法:字典破解和暴力破解. 字典破解是指通过破解者对管理员的了解,猜测其可能使用某些信息 ...
- DS8800后端的光纤通道交换式互连方式
DS8800 使用SAS 硬盘.使用了FC 到SAS 转换,光纤通道交换技术被用于DS8800 后端. FC 技术是普遍用于在一个光纤通道仲裁环路(Fibre Channel Arbitrated L ...
- JAVA所属公司与非盈利组织
版权 现在是Oracle公司的 Apache 负责Java发展的,重要的非盈利组织,主要产品包括Struts.Tomcat
- 【repost】Chrome 控制台console的用法
下面我们来看看console里面具体提供了哪些方法可以供我们平时调试时使用. 目前控制台方法和属性有: ["$$", "$x", "dir" ...