luogu 5471 [NOI2019]弹跳 KDtree + Dijkstra

题目链接  

第一眼就是 $KDtree$ 优化建图
然而，空间只有 $128mb$，开不下  
时间不吃紧，考虑直接跑 $Dijkstra$
$Dijkstra$ 中存储的是起点到每个输入时给出的矩阵的最短距离
当取出堆顶时就将这个矩阵中所有点 "裂开"，并更新每一个小点的答案
如果该点在之前已经被一个最短路更短的矩阵更新过了，就不扩展它
否则，扩展每一个分裂出的点，将新的矩阵加入优先队列中即可
一个重要的优化就是每次 "分裂" 一个矩阵时要在 $KDtree$ 中将对应节点都删掉，加快下一次查询 

#include<bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin) 
#define maxn 2000000  
#define inf 100000000
#define lson t[now].ch[0]
#define rson t[now].ch[1] 
#define ll long long            
using namespace std;  

struct Edge {  
    int t,x1,y1,x2,y2;  
    Edge(int t=0,int x1=0,int y1=0,int x2=0,int y2=0):t(t),x1(x1),y1(y1),x2(x2),y2(y2){} 
}; 
struct Data {     
    int ch[2],minv[2],maxv[2],siz,id,p[2];                      
}t[maxn],arr[maxn];   
struct P { 
    int dis,id;  
    P(int dis=0,int id=0):dis(dis),id(id) {}    
    bool operator<(P b) const {
        return b.dis<dis;    
    }
};  
vector<int>G[maxn];  
priority_queue<P>Q;    
vector<Edge>edges;     

int d,cnt,len;  
int dis[maxn],answer[maxn],nd[maxn],vis[maxn];         
void Min(int &a,int b) { 
    if(b<a) a=b;  
}
void Max(int &a,int b) { 
    if(b>a) a=b;  
}
bool cmp(Data i,Data j) {  
    return i.p[d]==j.p[d]?i.p[d^1]<j.p[d^1]:i.p[d]<j.p[d];   
}  
void pushup(int now) {
    t[now].minv[0]=t[now].minv[1]=inf,t[now].maxv[0]=t[now].maxv[1]=-inf; 
    for(int i=0;i<2;++i) {
        if(lson) {
            Min(t[now].minv[0],t[lson].minv[0]),Min(t[now].minv[1],t[lson].minv[1]),
            Max(t[now].maxv[0],t[lson].maxv[0]),Max(t[now].maxv[1],t[lson].maxv[1]);   
        }
        if(rson) {
            Min(t[now].minv[0],t[rson].minv[0]),Min(t[now].minv[1],t[rson].minv[1]),
            Max(t[now].maxv[0],t[rson].maxv[0]),Max(t[now].maxv[1],t[rson].maxv[1]);   
        }
    }         
}
void build(int l,int r,int &now,int o) {
    now=++cnt;    
    if(l==r) {
        t[now].siz=1, t[now].id=arr[l].id;        
        t[now].p[0]=t[now].minv[0]=t[now].maxv[0]=arr[l].p[0];  
        t[now].p[1]=t[now].minv[1]=t[now].maxv[1]=arr[l].p[1];       
        return;  
    }
    int mid=(l+r)>>1;    
    d=o,nth_element(arr+l,arr+mid,arr+1+r,cmp);      
    build(l,mid,lson,o^1);         
    if(r>mid) build(mid+1,r,rson,o^1);       
    pushup(now),t[now].siz=t[lson].siz+t[rson].siz;               
}
void dfs(int now,int x1,int y1,int x2,int y2) {
    if(t[now].minv[0]>y1||t[now].maxv[0]<x1||t[now].minv[1]>y2||t[now].maxv[1]<x2 || !t[now].siz) return;    
    if(t[now].id) {  
        nd[++len]=t[now].id;       
        t[now].siz=0;            
        return; 
    }  
    if(lson) dfs(lson,x1,y1,x2,y2); 
    if(rson) dfs(rson,x1,y1,x2,y2);  
    t[now].siz=t[lson].siz+t[rson].siz;           
}
void debug(int now) {
    if(!now) return; 
    debug(lson), printf("%d %d %d %d\n",t[now].minv[0],t[now].maxv[0],t[now].minv[1],t[now].maxv[1]),debug(rson);   
}
int main(){
    //   setIO("input");  
    int n,m,w,h,rt; 
    scanf("%d%d%d%d",&n,&m,&w,&h); 
    for(int i=1;i<=n;++i)  scanf("%d%d",&arr[i].p[0],&arr[i].p[1]),arr[i].id=i;           
    build(1,n,rt,0);                               
    for(int i=1;i<=m;++i) {
        int a,b,c,d,e,f;         
        scanf("%d%d%d%d%d%d",&a,&b,&c,&d,&e,&f), edges.push_back(Edge(b,c,d,e,f)), G[a].push_back(i-1);            
    }                                                      
    answer[1]=0,vis[1]=1;                                
    for(int i=0;i<G[1].size();++i) {
        int to=G[1][i]; 
        dis[to]=edges[to].t, Q.push(P(edges[to].t, to));         
    }
    while(!Q.empty()) { 
        P e=Q.top(); Q.pop();
        int cn=e.dis;   
        len=1,dfs(rt, edges[e.id].x1,edges[e.id].y1,edges[e.id].x2,edges[e.id].y2);          
        for(int i=1;i<=len;++i) {   
            int cur=nd[i]; 
            if(vis[cur]) continue;  
            vis[cur]=1, answer[cur]=cn;       
            for(int j=0;j<G[cur].size();++j) {   
                int to=G[cur][j];                    
                dis[to] = cn + edges[to].t;                                           
                Q.push(P(dis[to], to));            
            }
        }
    } 
    for(int i=2;i<=n;++i) printf("%d\n",answer[i]);     
    return 0; 
}