For enabling efficient, large-scale coordination of unmanned aerial vehicles (UAVs) under the labeled setting, in this work, we develop the first polynomial time algorithm for the reconfiguration of many moving bodies in three-dimensional spaces, with provable 1.x asymptotic makespan optimality guarantee under high robot density. More precisely, on an m1 × m2 × m3 grid, m1 ≥ m2 ≥ m3, our method computes solutions for routing up to m 1 m 2 m 3 3 uniquely labeled robots with uniformly randomly distributed start and goal configurations within a makespan of m1 + 2m2 + 2m3 + o(m1), with high probability. Because the makespan lower bound for such instances is m1 + m2 + m3 − o(m1), also with high probability, as m1 → ∞, m 1 +2m 2 +2m 3 m 1 +m 2 +m 3 optimality guarantee is achieved.3 ], yielding 1.x optimality. In contrast, it is well-known that multi-robot path planning is NP-hard to optimally solve. In numerical evaluations, our method readily scales to support the motion planning of over 100, 000 robots in 3D while simultaneously achieving 1.x optimality. We demonstrate the application of our method in coordinating many quadcopters in both simulation and hardware experiments.