[NOI2019]弹跳，K-D Tree+Dijkstra

正题

      Portal

      好卡常啊，我的常数太不优秀了。

      刷KDT的时候看到有一道有意思的NOI题，就顺便来做了。

      前面不卡空间的分数还是挺好想的，用KDT来维护这个区间连边即可，树上的连边也很显然。

      发现只有56分，因为只有128mb。

      但是发现在Dijkstra的时候，我们每一个点都只会松弛一遍，对于一个在KDT上的点，它的出边最多只有三条。

      对于一个原矩阵的点，他可能通过跳板跳到另一个节点，那么松弛的时候相当于查询这些子矩阵对应的节点。

      那么就不用连边，时间复杂度还是一样的。

      加了读优和ZZY神仙的Ofast，终于在LGOJ不开O2的情况下跑过去了。

#pragma GCC optimize("Ofast")
#include
using namespace std;

const int N=70010,M=150010;
int cmp=0;
struct node{
	int x[2],id;
	bool operator<(const node q)const{
		return x[cmp] E[N];
struct Heap{
	int x,d;
	bool operator<(const Heap q)const{
		return d>q.d;
	}
};
priority_queue Q;
int lim[2][2];
int dis[N<<1];
bool tf[N<<1];
int n,m,w,h,num=0,rt,a,b;

inline char nc(){
    static char buf[200000000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,200000000,stdin),p1==p2)?EOF:*p1++;
}
inline void read(int&sum){
    char ch=nc();sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
}

inline void ins(int x,int y,int c){
	if(dis[y]>dis[x]){
		dis[y]=dis[x];
		Q.push((Heap){y,dis[y]});
	}
};

#define ls tr[now].ls
#define rs tr[now].rs

inline void update(int now){
	if(ls){
		tr[now].mmin[0]=min(tr[now].mmin[0],tr[ls].mmin[0]);
		tr[now].mmin[1]=min(tr[now].mmin[1],tr[ls].mmin[1]);
		tr[now].mmax[0]=max(tr[now].mmax[0],tr[ls].mmax[0]);
		tr[now].mmax[1]=max(tr[now].mmax[1],tr[ls].mmax[1]);
	}
	if(rs){
		tr[now].mmin[0]=min(tr[now].mmin[0],tr[rs].mmin[0]);
		tr[now].mmin[1]=min(tr[now].mmin[1],tr[rs].mmin[1]);
		tr[now].mmax[0]=max(tr[now].mmax[0],tr[rs].mmax[0]);
		tr[now].mmax[1]=max(tr[now].mmax[1],tr[rs].mmax[1]);
	}
}

void build(int &now,int l,int r,int type){
	now=++num;cmp=type;
	int mid=(l+r)/2;
	nth_element(s+1+l,s+1+mid,s+1+r);
	tr[now].p=s[mid];
	tr[now].mmin[0]=tr[now].mmax[0]=s[mid].x[0];
	tr[now].mmin[1]=tr[now].mmax[1]=s[mid].x[1];
	if(l