[学习笔记]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......)其中要求所有元素有序,并且支持基本的操作(删除,添加,查找......) ...
随机推荐
- 字符串方法 split() & replace()
split() 语法:stringObject.split(separator) 功能:把一个字符串分割成字符串数组 返回值:Array 说明:separator 是必须的,分隔符. var str= ...
- GCC选项之-M
大多数的C/C++编译器都支持一个“-M”的选项,即自动找寻源文件中包含的头文件.举个例子,比如mian.c包含有如下头文件. #include <stdio.h> #include &q ...
- tms web core 里面调用pascal 过程。
procedure show(s:string);begin showmessage(s);end; procedure TForm3.WebButton1Click(Sender: TObject ...
- Oracle导入的常见语句
登录sql > sqlplus / as sysdba 创建表空间sql > create tablespace TABLESPACE datafile 'e:\tables1.dbf' ...
- openstack的Host Aggregates和Availability Zones
1.关系 Availability Zones 通常是对 computes 节点上的资源在小的区域内进行逻辑上的分组和隔离.例如在同一个数据中心,我们可以将 Availability Zones 规划 ...
- Linux基础笔记
基本命令: man man 查看手册 (1)pwd 显示路径 (2)cd 切换路径 (3)ls -a;-l;-h;-i;-d (隐藏的:长列表显示:human显示:索引号:只显示列表) (4)c ...
- UVaLive 3357 Pinary (Fib数列+递归)
题意:求第 k 个不含前导 0 和连续 1 的二进制串. 析:1,10,100,101,1000,...很容易发现长度为 i 的二进制串的个数正好就是Fib数列的第 i 个数,因为第 i 个也有子问题 ...
- 上传文件夹+php
最近公司做工程项目,实现文件夹上传 网上找了很久,发现网上很多代码大都存在很多问题,不过还是让我找到了一个符合要求的项目. 对项目的文件夹上传功能做出分析,找出文件夹上传的原理,对文件夹的传输模式深入 ...
- what is a resolver
resolver [rɪ'zɒlvə] 解析器 ViewResolver The ViewResolver provides a mapping between view names and actu ...
- winSockets编程（二）socket函数
初始化DLL之后,接着创建套接字,通过socket()和WSASocket()函数实现此功能. SOCKET socket( int af, int type, int protocol ); af: ...