雅加达的摩天楼
题意描述:
- 有$N$座摩天楼，从左到右依次编号为$0$到$N-1$。
- 有$M$个信息传递员，编号依次为$0$到$M-1$。编号为i的传递员最初在编号为$B_i$的摩天楼，邮递员可以在摩天楼之间跳跃(向前或者向后)，编号为$i$的传递员跳跃能力为$P_i$。
- 当传递员到达一个摩天楼他可以执行两种操作
- 1: 跳跃到其他摩天楼上。
- 2: 将消息传递给当前摩天楼的其他传递员。
- 最终要将编号为$0$的摩天楼信息传递到编号为$1$的摩天楼地方。
输入格式:
- 第一行输入$N$和$M$。$(1\leq n \leq 310^4， 1\leq m \leq 310^4)$
- 接下来第$2$到$M+1$，每行包含两个整数$B_i$和$P_i$。
输出格式:
- 输出一行，表示最小步数，如果无法到达则输出$-1$。
解题思路：
- 分析完题目后可以想到是最短路，对于每一个传递员都可以跳到别处，可以以此建立起点和终点，边权为1。
- 考虑暴力加边的做法，但是如果$p$很小的时候，最坏情况下有$n^2$条边，任何最短路算法都无法通过。
- 优化建图，采用分块的思想: 将$n$点合并为$\sqrt{n}$个块
- 对于$P_i\leq\sqrt{n}$的传递员：
- 对于某一栋摩天大楼，就把他从一个点转化为一栋真的摩天大楼。
- 对于每栋楼建$\sqrt{n}$层，每层有$n$个点，第一层代表$p=1$的走法，第二层代表$p=2$的走法,....。其中第$0$代表原先的那个点。
- 然后暴力加边。第$i$层代表$P=i$的情况，间距为$i$的点相互连长度为$1$的边。再把这些辅助点全部连向底边。
- 对于每个$P_i>\sqrt{n}$的点向其能达到的所有点的第$0$层加边。
- 总边数不会超过$N\sqrt{N}$。
- 建图完毕，跑最短路。
代码：
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = (3e4 + 10)*500;
const int maxm = (3e4 + 10)*500;
int n, m;
int s, t; //起点和终点
int B[maxn], P[maxn]; //起始于B, 跳跃能力为P
int block;
int head[maxn], nex[maxn<<1], ver[maxn<<1], edge[maxn<<1], tot;
inline void add_edge(int x, int y, int z)
{
ver[++tot] = y; edge[tot] = z;
nex[tot] = head[x]; head[x] = tot;
}
inline int get_block(int x, int y){
return x*n + y;
}
int dist[maxn];
bool v[maxn];
void SPFA()
{
memset(dist, 0x3f, sizeof(dist));
dist[s] = 0; queue<int> q;
q.push(s); v[s] = 1;
while(q.size())
{
int x = q.front(); q.pop();
v[x] = 0;
for(int i = head[x]; i; i = nex[i])
{
int y = ver[i], z = edge[i];
if(dist[y] > dist[x] + z)
{
dist[y] = dist[x] + z;
if(!v[y])
{
q.push(y);
v[y] = 1;
}
}
}
}
printf("%d\n", dist[t] == 0x3f3f3f3f ? -1 : dist[t]);
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = 1, x, y; i <= m; i++)
{
scanf("%d%d", &x, &y);
B[i] = x + 1, P[i] = y;
}
//起始位置 终止位置 块的大小
s = B[1], t = B[2]; block = min((int)sqrt(n), 100);
//每一个点转化为一栋楼, 有sqrt(n)层
for(int i = 1; i <= block; i++)
for(int j = 1; j <= n; j++)
{
add_edge(get_block(i, j), j, 0); //楼底到每一层连边
if(j <= n - i) //j条一次跳出范围了
{
//每一层间隔为i的点 跳一次边权为1
add_edge(get_block(i,j), get_block(i,j)+i, 1);
add_edge(get_block(i,j)+i, get_block(i,j), 1);
}
}
for(int i = 1; i <= m; i++) //考虑每一个doge的跳跃能力
{
// 对于每一个传递员, 从楼底到Pi层对应连边
if(P[i] <= block)
add_edge(B[i], get_block(P[i], B[i]), 0);
else
{ //如果大于一个块了 那就暴力加边 因为有sqrt个块 所以加sqrt条边
for(int j = 1; B[i] + j * P[i] <= n; j++)
add_edge(B[i], B[i] + j*P[i], j); //向前跳
for(int j = 1; B[i] - j * P[i] >= 1; j++)
add_edge(B[i], B[i] - j*P[i], j); //向后跳
} //这里最大可能加N*sqrt(N)条边
} SPFA(); //优化后边数有 N*sqrt(N)的规模
return 0;
}