20200508 SCOI模拟T2（分层图最短路）

T2 [APIO2015]雅加达的摩天楼

P3645 [APIO2015]雅加达的摩天楼
题目描述
印尼首都雅加达市有 N 座摩天楼，它们排列成一条直线，我们从左到右依次将它们编号为 0 到 N−1。除了这 N 座摩天楼外，雅加达市没有其他摩天楼。
有 M 只叫做 “doge” 的神秘生物在雅加达市居住，它们的编号依次是 0 到 M−1。编号为 i 的 doge 最初居住于编号为 $B_i$ 的摩天楼。每只 doge 都有一种神秘的力量，使它们能够在摩天楼之间跳跃，编号为 i 的 doge 的跳跃能力为 $P_i（P_i>0）$。
在一次跳跃中，位于摩天楼 b 而跳跃能力为 p 的 doge 可以跳跃到编号为 b−p （如果 $0\le b-p<N$）或 b+p （如果 $0\le b+p<N$）的摩天楼。
编号为 0 的 doge 是所有 doge 的首领，它有一条紧急的消息要尽快传送给编号为 1 的 doge。任何一个收到消息的 doge 有以下两个选择:
跳跃到其他摩天楼上；
将消息传递给它当前所在的摩天楼上的其他 doge。
请帮助 doge 们计算将消息从 0 号 doge 传递到 1 号 doge 所需要的最少总跳跃步数，或者告诉它们消息永远不可能传递到 1 号 doge。

输入格式
输入的第一行包含两个整数 N 和 M。
接下来 M 行，每行包含两个整数 $B_i$ 和 $P_i$。

输出格式
输出一行，表示所需要的最少步数。如果消息永远无法传递到 1 号 doge，输出 −1。

输入输出样例
输入
5 3
0 2
1 1
4 1
输出
5

说明/提示
【样例解释】
下面是一种步数为 5 的解决方案：
0 号 doge 跳跃到 2 号摩天楼，再跳跃到 4 号摩天楼（2 步）。
0 号 doge 将消息传递给 2 号 doge。
2 号 doge 跳跃到 3 号摩天楼，接着跳跃到 2 号摩天楼，再跳跃到 1 号摩天楼（3 步）。
2 号 doge 将消息传递给 1 号 doge。

【数据范围】
所有数据都保证 $0\le B_i<N,1\le N\le 30000,1\le P_i\le 30000,2\le M\le 30000$

思路：
考虑一种建图方法，对于每一个 doge，向他可以到的所有点连边，边权为他到这个点的步数，在图上跑最短路即可
发现这样建图时间要炸，考虑优化，建立分层图
跳跃能力小于 $\sqrt{n}$ 的 doge，对于每一种跳跃能力建立 n 个点，相邻的能到达的点之间连一条边权为 1 的边，对于每一个点，向原图对应点连一条边权为 0的边，对于每一个 doge，向他的跳跃能力的一层的对应点连一条边权为 0 的边。
跳跃能力大于 $\sqrt{n}$ 的 doge，暴力建图
复杂度证明：对于第一种边，边数小于 $n\sqrt{n}$，对于第二种边，边数小于 $n\sqrt{n}$。
嫖张图理解一下

代码：

#include<bits/stdc++.h>
using namespace std;
#define in Read()
#define re register

inline char ch(){
	static char buf[1<<21],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;
}

inline int in{
	int s=0,f=1;char x;
	for(x=ch();x<'0'||x>'9';x=ch())	if(x=='-')	f=-1;
	for( ;x>='0'&&x<='9';x=ch())	s=(s*10)+(x&15);
	return f==1?s:-s;	
}

const int A=1e7+5;
const int INF=1e9;
int n,m,blo,all;
inline int d(int x,int dep){
	return dep*n+x;
}
int head[A],tot_road=0;
struct Road{
	int nex,to,w;
}road[2*A];
inline void ljb(int x,int y,int w){
	road[++tot_road]={head[x],y,w};head[x]=tot_road;
}
int root,last;

inline void prepare(){
	for(re int i=1;i<=blo;++i)
		for(re int j=0;j<n;++j){
			if(i+j<n)	ljb(d(j,i),d(j+i,i),1),ljb(d(j+i,i),d(j,i),1);
			ljb(d(j,i),j,0);
		}
	return;
}

inline void build(int x,int val){
	int now=x,num=1;
	now+=val;
	while(now<n){
		ljb(x,now,num);
		now+=val,num++;
	}
	now=x,num=1;
	now-=val;
	while(now>=0){
		ljb(x,now,num);
		now-=val,num++;
	}
	return;
}

int dis[A];
bool ex[A];
priority_queue <pair<int,int> > q;

inline void djk(){
	fill(dis,dis+(blo+1)*n+1,INF);
	dis[root]=0;
	q.push(make_pair(0,root));
	while(!q.empty()){
		int x=q.top().second;q.pop();
		if(ex[x])	continue;
		ex[x]=1;
		for(re int y=head[x];y;y=road[y].nex){
			int z=road[y].to,w=road[y].w;
			if(dis[x]+w<dis[z]){
				dis[z]=dis[x]+w;
				q.push(make_pair(-dis[z],z));
			}
		}
	}
	return;
}

signed main(){
//	freopen("skyscraper.in","r",stdin);
//	freopen("skyscraper.out","w",stdout);
	n=in,m=in;
	blo=sqrt(n);
	prepare();
	for(re int i=0;i<m;++i){
		int u=in,v=in;
		if(v<=blo)	ljb(u,d(u,v),0);
		else	build(u,v);
		if(i==0)	root=u;
		if(i==1)	last=u;
	}
	djk();
	if(dis[last]==INF)	printf("-1\n");
	else	printf("%d\n",dis[last]);
	return 0;
}