[NOI2019] 弹跳

题意：

给你平面上的$n$个点，共有$m$个弹跳装置。

每个弹跳装置可以从点$p_{i}$以$t_{i}$的代价跳到矩形$(L_{i},D_{i}),(R_{i},U_{i})$中的任何一个点。

现在请你对于每座城市求出从1号点跳到它的最小代价。

$n\leq 70000,m\leq 150000$。

题解：不把边建出来的$KD-tree$优化建图。

看一眼就知道$KD-tree$优化建图，但如果把所有边都建出来就爆炸了。

设原图上的点是实点，$KD-tree$上的点（代表一个实点和一个矩形）是虚点。

那么在$Dijkstra$到每个点的时候：

• 若是实点，在$KD-tree$上查找能连的虚点/实点并向其连带权边。
• 若是虚点，向它的两个儿子和它对应的实点连权为0的边。

时间复杂度$O(m\sqrt{n})$，空间复杂度$O(m\sqrt{n})$（实际上只有优先队列可能达到这个空间，其他都是$O(n)$）。

没了。不知道为什么有人写线段树。

代码：

#include<bits/stdc++.h>
#define maxn 200005
#define maxm 500005
#define inf 0x7fffffff
#define ll long long
#define rint register int
#define debug(x) cerr<<#x<<": "<<x<<endl
#define fgx cerr<<"--------------"<<endl
#define dgx cerr<<"=============="<<endl

using namespace std;
struct Point{int x[],iid;}p[maxn];
int mx[maxn][],mn[maxn][],id[maxn],ls[maxn],rs[maxn];
int tot,np,dis[maxn<<],vis[maxn<<],tim[maxn];
int L[maxn],R[maxn],D[maxn],U[maxn];
struct node{int u,d;bool operator<(const node &b)const{return d>b.d;}};
priority_queue<node> q; vector<int> vc[maxn];

inline int read(){
    int x=,f=; char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
    for(;isdigit(c);c=getchar()) x=x*+c-'';
    return x*f;
}

inline node md(int u,int d){node res;res.u=u,res.d=d;return res;}
inline void upd(int u,int v){if(dis[u]>v){dis[u]=v,q.push(md(u,v));}}
inline bool ok(int xa,int ya,int xb,int yb){return xb<=xa&&ya<=yb;}
inline bool cmp(Point a,Point b){return (!np)?(a.x[]<b.x[]):(a.x[]<b.x[]);}
inline void pushup(int k){
    for(rint i=;i<;i++){
        mx[k][i]=mn[k][i]=p[id[k]].x[i];
        if(ls[k]) mx[k][i]=max(mx[k][i],mx[ls[k]][i]),mn[k][i]=min(mn[k][i],mn[ls[k]][i]);
        if(rs[k]) mx[k][i]=max(mx[k][i],mx[rs[k]][i]),mn[k][i]=min(mn[k][i],mn[rs[k]][i]);
    }
}
inline int build(int l,int r,int type){
    if(l>r) return ;
    int mid=l+r>>,k=++tot;
    id[k]=mid,np=type;
    nth_element(p+l,p+mid,p+r+,cmp);
    ls[k]=build(l,mid-,type^);
    rs[k]=build(mid+,r,type^);
    pushup(k); return k;
}
inline void add(int l,int r,int d,int u,int c,int k){
    if(mx[k][]<l || mn[k][]>r || mx[k][]<d || mn[k][]>u) return;
    if(ok(mn[k][],mx[k][],l,r)&&ok(mn[k][],mx[k][],d,u)){upd(k,c);return;}
    if(ok(p[id[k]].x[],p[id[k]].x[],l,r)&&ok(p[id[k]].x[],p[id[k]].x[],d,u)) upd(p[id[k]].iid,c);
    if(ls[k]) add(l,r,d,u,c,ls[k]); if(rs[k]) add(l,r,d,u,c,rs[k]);
} 

inline void Dijkstra(int s,int n){
    memset(vis,,sizeof(vis));
    memset(dis,,sizeof(dis));
    dis[s]=,q.push(md(s,));
    while(!q.empty()){
        rint u=q.top().u; q.pop();
        if(vis[u]) continue; vis[u]=;
        if(u<=n) for(rint i=;i<vc[u].size();i++)
            add(L[vc[u][i]],R[vc[u][i]],D[vc[u][i]],U[vc[u][i]],dis[u]+tim[vc[u][i]],n+);
        else upd(p[id[u]].iid,dis[u]),upd(ls[u],dis[u]),upd(rs[u],dis[u]);
    }
}

int main(){
    freopen("jump.in","r",stdin);
    freopen("jump.out","w",stdout);
    int n=read(),m=read(),w=read(),h=read();
    for(rint i=;i<=n;i++)
        p[i].x[]=read(),p[i].x[]=read(),p[i].iid=i;
    tot=n,build(,n,);
    for(rint i=;i<=m;i++){
        int p=read();tim[i]=read();
        L[i]=read(),R[i]=read();
        D[i]=read(),U[i]=read();
        vc[p].push_back(i);
    }
    Dijkstra(,n);
    for(rint i=;i<=n;i++) printf("%d\n",dis[i]);
    return ;
}

D2T1