Abstract 1 Introduction 2 Preliminaries 3 Constraints for tighter separation bounds 4 Unlabeled MRMP with improved separation trade-offs 5 The Exodus algorithm: reducing the robots-separation bound 6 Discussion and future work References

Unlabeled Multi-Robot Motion Planning with Improved Separation Trade-Offs

Tsuri Farhana ORCID The Stein Faculty of Computer and Information Science, Ben-Gurion University of the Negev, Beer Sheva, Israel    Omrit Filtser ORCID Department of Mathematics and Computer Science, The Open University of Israel, Ra’anana, Israel    Shalev Goldshtein ORCID Department of Mathematics and Computer Science, The Open University of Israel, Ra’anana, Israel
Abstract

We study unlabeled multi-robot motion planning for unit-disk robots in a polygonal environment. Although the problem is hard in general, polynomial-time solutions exist under appropriate separation assumptions on start and target positions. Solovey et al. (RSS’15) provide a near-optimal solution assuming that start/target positions must have pairwise distance at least 4, and at least 52.236 from obstacles. This raises the question of whether polynomial-time algorithms can be obtained in even more densely packed environments.

In this paper we present a generalized algorithm that achieve different trade-offs on the robots-separation and obstacles-separation bounds, all significantly improving upon the state of the art. Specifically, we obtain polynomial-time constant-approximation algorithms to minimize the total path length when (i) the robots-separation is 223 and the obstacles-separation is 123, or (ii) the robots-separation is 3.291 and the obstacles-separation 1.354. Additionally, we introduce a different strategy yielding a polynomial-time solution when the robots-separation is only 2, and the obstacles-separation is 3. Finally, we show that without any robots-separation assumption, obstacles-separation of at least 1.5 may be necessary for a solution to exist.

Keywords and phrases:
multi-robot motion planning
Funding:
Omrit Filtser: This research was supported by the ISRAEL SCIENCE FOUNDATION (grant No. 2135/24).
Shalev Goldshtein: This research was supported by the ISRAEL SCIENCE FOUNDATION (grant No. 2135/24).
Copyright and License:
[Uncaptioned image] © Tsuri Farhana, Omrit Filtser, and Shalev Goldshtein; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Computational geometry
Related Version:
Full Version: https://arxiv.org/abs/2603.19502 [10]
Acknowledgements:
We thank Dan Halperin for introducing us to the problem and for helpful discussions. We are also grateful to the anonymous reviewers for their insightful comments, which helped improve the structure and presentation of the paper.
Editors:
Hee-Kap Ahn, Michael Hoffmann, and Amir Nayyeri

1 Introduction

In the classic multi-robot motion-planning (MRMP) problem, a set of robots operate in a common workspace, and need to reach a set of target positions. More precisely, given a set of m starting positions and a set of m target positions, the goal is to compute a motion strategy for the robots to reach the targets while avoiding collisions with the boundary of the workspace and with other robots. In this paper we consider unit-disk robots that are moving in a polygonal workspace in the plane. If each robot is assigned to a specific target, the problem is referred to as labeled MRMP, and otherwise it is unlabeled MRMP.

The MRMP problem and its variants have been widely investigated, both from theoretical and practical perspectives (see, e.g., [3, 4, 5, 6, 9, 11] for some of the more recent works on the problem). Very recently, both the labeled and unlabeled variants of MRMP for unit disks in polygons with holes were shown to be PSPACE-hard [1]. Many other variants of the problem are known to be hard as well, see e.g. [8, 12, 14, 16].

Interestingly, the separation between the robots plays a key role in the difficulty of the problem. Adler, de Berg, Halperin, and Solovey [2] show that in a simple polygonal workspace with n vertices and m unit-disk robots, when the robots-separation, denoted ρ, is 4 (i.e., all the start and target positions are at least distance 4 apart), a solution for unlabeled MRMP always exists, and can be found in O~(mn+m2) time. In [7] this result is further improved by distinguishing between monochromatic robots-separation (i.e. the distance between two start or two target positions), and bichromatic robots-separation (i.e. the distance between a start position and a target position). The algorithmic approach in both papers [7, 2] is similar and does not provide any reasonable approximation on the length of the solution, i.e., the total length of the paths in the motion plan.

Solovey, Yu, Zamir, and Halperin [15] present a different approach for unlabeled MRMP, with optimality guarantee even in polygons with holes, but at the cost of adding an additional assumption on the distance between the obstacles to any start and target position. Specifically, they consider the case when the robots-separation is ρ=4, and the obstacles-separation, denoted ω, is 52.236. Under these assumptions, they present an algorithm that runs in O~(m4+n2m2) time and returns a solution of total length at most OPT+4m. Here and throughout the paper, OPT denotes the minimum possible sum of lengths of the paths in a collision-free motion plan.

Notice that when assuming a large robots-separation (ρ4), the overall “density” of the robots inside the workspace becomes small, which may require the workspace to be much larger than the number of robots operating in it. Also, assuming a large obstacles-separation requires adding an additional buffering “layer” around the obstacles. This motivates the question of whether we can find weaker assumptions on either the robots-separation or the obstacles-separation, that are sufficient for having a polynomial time algorithm for unlabeled MRMP. Surprisingly, in this paper we show that both the robots-separation and obstacles-separation assumptions can be significantly relaxed, while still enabling a polynomial-time algorithm that achieves a constant-factor approximation of the solution length.

In previous works, namely [2, 7, 15], the algorithmic approach uses a monotone strategy, where robots move to their target position one by one according to some order, while other robots must remain at their start/target positions. The main result in our paper uses a weakly-monotone strategy. In a weakly-monotone motion plan, each start/target position p is associated with a revolving area - a disk of radius r>1 that contains the unit disk centered at p and is empty of obstacles and other unit disks centered at other start/target positions. The robots move to their target position one by one according to some order, while other robots must remain inside the revolving area associated with their current position. Such a weakly-monotone strategy was used, e.g., in [4, 13] for a variant of labeled MRMP with unit-disk robots.

Our contribution.

In this paper we present two different algorithms for unlabeled MRMP, that provide tighter trade-offs between robot-separation and obstacles-separation assumptions. In particular, for a parameter 0ε1, a polygonal environment of complexity n, and m unit-disc robots, our first algorithm computes in O~(m4+n2m2) time a collision-free motion plan of total length O(OPT), assuming that ρ=max{42ε,2+ε} and ω=max{1+ε,(33ε)2+1} (see Theorem 9). In Table 1 we highlight three interesting values of ε that minimizes either ρ,ω or the guaranteed length of the solution. This algorithm can be viewed as a generalization of the algorithm in [15], and indeed, as a consequence of our analysis we positively resolve a conjecture made in [15], by showing that ω=1.614 is sufficient when ρ=4. Moreover, in the full version of our paper [10] we show that a monotone motion plan (such as the one produced by the algorithm in [15]) may not always exist when ω<1.614, and therefore to achieve smaller separation bounds one must use a different strategy. Indeed, as mentioned above, our motion plan is weakly-monotone. Interestingly, we also show that a weakly-monotone motion plan may not always exist when ω<1.354.

Note that our first algorithm does not cover the case of ρ=2, which is the optimal (monochromatic) robots-separation. We therefore present a second algorithm using a completely different approach, which given a simple polygon with n vertices and m unit-disk robots, computes in O(m3+mn2) time a collision-free motion plan of total length OPT+O(m2), assuming that ρ=2 and ω=3 (see Theorem 10). Finally, we show that ω=1.5 may be necessary for a solution to exist.

Table 1: Summary of algorithms and their separation assumptions.
𝝎 𝝆 Solution length Running time Notes Reference
1 4 no guarantee O~(mn+m2) Simple polygons [2, 7]
5 4 OPT+4m O~(m4+n2m2) [15]
1.614 4 OPT+4m O~(m4+n2m2) Theorem 9
1.354 3.291 O(OPT) O~(m4+n2m2) Theorem 9
123 223 O(OPT) O~(m4+n2m2) Theorem 9
3 2 OPT+O(m2) O~(m3+mn2) Simple polygons Theorem 10

Technical ideas and overview of the paper

Our approach in the first algorithm generalizes the strategy of Solovey et al. [15] for unlabeled MRMP, which is briefly described as follows. Given a set S of start positions, a set T of target positions, and a polygonal workspace 𝒲, first compute an optimal set Γ of paths assigning each starting point to a target point (a perfect matching). This set of paths minimizes the sum of geodesic distances, when ignoring the other robots in the workspace. Then, one can show that when assuming ρ=4 and ω=5, there is always at least one target which is a so-called standalone goal, that is, a target t that does not block any other path in Γ. Such a target is useful because we can place a robot there and then treat it as an obstacle. However, there might be some other robot that blocks the path that leads to t, in which case we can make a switch: we push the last robot that blocks this path to t. This switching process is possible because of the assumption ω=5. Then, run the algorithm recursively, but after removing the start and target positions that where satisfied in the previous recursive step, and updating the workspace to treat the robot at t as an obstacle.

Our weakly-monotone motion plan is obtained by inserting a very important modification to the algorithm of [15]: instead of choosing a standalone target in each iteration, we relax the requirement and require only an “almost standalone” target – a target that may interrupt other paths, but only by a small amount. In our strategy, we distinguish between “blocking” and “interrupting” positions. Instead of always switching paths when another robot is standing in the way, we allow robots that interrupt a path by only a small amount to move away from the path. A robot that interrupts a given path may move slightly in its close neighborhood (a disk of radius equal to the amount of interruption plus one) in order to clear the way. In contrast, a robot that blocks a path may not have enough wiggle room, and we will switch the paths.

As suggested above, this combined strategy can be viewed as a generalization of [15]: we introduce another parameter, called the “overlap parameter”, which is used to determine whether a position is blocking a path or only interrupting it. This parameter allows us to define a set of constraints on ρ and ω, that allows (i) the existence of an interrupting target, (ii) the option to switch paths in the case when a robot is blocking a path, and (iii) the existence of a sufficiently large free neighborhood for an interrupting robot to clear the way. It turns out that there is an interesting trade-off between these constraints, leading to the main result of our paper (Theorem 9).

Note that the concept of interrupting positions raises some issues that were not present in the algorithm of [15]. For example, in each iteration of the recursive algorithm we need to take into account that the robot placed on a target may still interrupt other paths. In [15] this is solved simply by removing a unit disk around the interrupting target from the workspace. However, in our case, such a robot may still need to move in order to clear the way for other robots. Moreover, if we simply treat it as an obstacle, the assignment of paths in the next iteration may be too large. We solve this issue by observing that it is actually enough to remove a smaller disk around a settled target when computing the paths, while still allowing it to move within its small neighborhood when planning the motion in the following iterations. Another issue is that there may be two or more robots that interrupt a path and may need to move simultaneously. Here we describe a simple motion with respect to the ray connecting the robot moving on the path with the interrupting robots.

In Section 3, we systematically define and analyze the constraints that are required for each part of this strategy, each constraint depends on the overlap parameter. In Section 4 we present the complete generalized algorithm and prove its correctness. Then, we show how to choose an overlap parameter and adjust constraints to get the minimum possible obstacles-separation and the minimum possible robots-separation.

In our first algorithm, to clear a path, we only move robots that intersect that path. Such a strategy cannot work when the robots-separation is exactly 2, because a robot that intersect a path may be block by other robots that do not intersect that path, and thus will not be able to move and clear the path. Therefore, to achieve a separation bound of 2 between the robots, we present in Section 5 a completely different strategy, which we call the Exodus algorithm. In this strategy, we also iteratively pick a path between start and target positions that were not yet satisfied, and move a robot along that path, however, in contrast to the first algorithm, we clear the path completely by also moving robots that do not intersect the path. In simple polygons, this strategy allows us to assume robots-separation of only 2, with obstacles-separation of 3, and it applies also for labeled MRMP (see Theorem 10).

Note that the monochromatic separation must be at least 2, so the robots-separation is almost optimal (the bichromatic separation can be 0). In Section 6 we preset lower bounds on the obstacles-separation, and we leave open the question of what is the upper bound on ω when robots-separation is optimal.

2 Preliminaries

We consider a set of m unit-disk robots, moving in a polygonal workspace 𝒲2 (which may be simple or cluttered with polygonal obstacles), whose overall number of edges is n. Below, we mostly follow the notations in [7] and [15].

The obstacle space 𝒪 is the complement of the workspace 𝒲, that is, 𝒪=2𝒲. We refer to the points in 𝒲 as positions, and say that a robot is at position x𝒲 when its center is positioned at the point x. For a point x2 and a radius r+, we define Dr(x) to be the open disk of radius r centered at x. The unit-disk robots are defined to be open sets, and therefore two robots collide if and only if the distance between their positions is strictly less than 2. In addition, a robot collides with the obstacle space 𝒪 if and only if its center is at distance strictly less than 1 from 𝒪. The set of all positions in 𝒲 where a unit-disk robot does not collide with the obstacle space is called the free space, denoted ={x2D1(x)𝒪=}. Note that the free space is a closed set, and denote by the boundary of .

In the multi-robot motion-planning problem, we are given a workspace 𝒲 with n edges and two sets S={s1,,sm} and T={t1,,tm} of points in the free space . The goal is to plan a collision-free motion for m robots, each positioned on a different starting point in S, such that by the end of the motion every target position in T is occupied by some robot. When the robots are distinguishable (i.e. labeled), the robot that starts at position si is required to end up at the target ti. Otherwise, when the robots are indistinguishable (i.e., unlabeled), it does not matter which robot ends up at which target position, as long as all targets are occupied at the end of the motion.

Formally, we wish to find a set Γ of continuous paths γ1,,γm, where γi:[0,1] for every 1im, and such that i=1mγi(0)=S and i=1mγi(1)=T. In addition, the paths need to be collision-free, i.e., at any point in time t[0,1], there are no i,j such that γi(t) and γj(t) are at distance strictly less than 2. Denote the size of a solution Γ by |Γ|=γiΓ|γi|, where |γ| is the length of a path γ in the L2 norm. In the optimization version of the problem, we are interested in finding a solution Γ which minimizes |Γ|, i.e., minimizing the sum of lengths of the paths.

Separation assumptions.

We denote by ω the obstacles-separation assumption on an instance of MRMP, i.e., we assume that for every pST and every point x𝒪, it holds that xpω. We denote by ρ the robots-separation assumption, i.e., we assume that for every p1,p2ST it holds that p1p2ρ.

Overlapping, blocking and interrupting positions

Let p,q be two points in the free space. We define the overlap of p and q to be max{0,2pq} (see Figure 1). If 𝒟1(p)𝒟1(q)=, then their overlap is 0. Otherwise, we say that p,q are ε-overlapping for ε=2pq, and note that pq=2ε.

Figure 1: Overlapping, blocking, and interrupting positions.

We say that a position p is ε-blocking a path γ if there exists 0t1 such that the overlap of p and γ(t) is strictly larger than ε. Blocking robots (i.e. robots that are placed on blocking positions) may be impossible to clear from the path, and we will have to make a “switch”. We say that a point p is ε-interrupting a path γ if for every 0t1, the overlap of p and γ(t) is at most ε (that is, it is not ε-blocking). Interrupting robots (i.e. robots that are placed on interrupting positions) will have enough wiggle room to move away from γ and “clear” the path. Notice that if for every 0t1 the overlap of p and γ(t) is 0, then for every 0t1 it holds that γ(t)p>2, and therefore in this case we say that p is 0-interrupting γ. We say that a path γ is ε-blocked if there exists sS that ε-blocks it, and that γ is ε-interrupted if no sS is ε-blocking it.

The overlap parameter.

For the analysis in Section 3, we fix a parameter 0ε1 which we call the overlap parameter. In these sections, for simplicity of the presentation, we sometimes omit ε from the above notions and just say “blocking” or “interrupting”. In our analysis, we introduce a set of constrains described as functions of ε, which are required for the different parts of our algorithm to work. In Section 4 we show how to choose ε in order to optimize different separation bounds.

3 Constraints for tighter separation bounds

To give an intuition for the required constraints, we first provide a high-level overview of our strategy. The algorithm is recursive, and its input is the overlap parameter ε, a set S of starting positions, a set T of target positions, and a workspace 𝒲. In each recursive iteration, one robot is moving to a target position, then S,T and 𝒲 are updated accordingly.

The first step is to compute an optimal-assignment path set, as it is defined in [15]:

Definition 1.

Let Γ be a set of m geodesic paths γ1,,γm, where γi:[0,1] for every 1im, and such that i=1mγi(0)=S, i=1mγi(1)=T. We say that Γ is an optimal-assignment path set for S,T,𝒲 if it minimizes i=1m|γi| over all such path sets.

Note that an optimal assignment path set Γ is not necessarily a feasible solution, because it ignores possible collisions between robots. This also means that the sum of paths in Γ is at most the length of an optimal solution to our MRMP problem. Let Γ={γ1,,γm} be an optimal-assignment path set for S,T,𝒲.

In Section 3.1 we present a useful lemma showing that a blocking robot can enter the geodesic path that it blocks. This lemma will be useful throughout the analysis, and it requires one constraint on the obstacles-separation, and another on the robots-separation.

The next step is to find a target tT which does not block any path in Γ (other than the path γ from sS to t). Then we can move a robot from s to t along γ, and it will not block other paths in the following iterations. Note that t may still interrupt other paths. In Section 3.2 we show that such a target always exists, when assuming an additional constraint on the robots-separation.

If there is a position sS, which blocks the path γ, we need to make a switch, i.e., we construct a new path in from s to t, which is not blocked by any other robot. In Section 3.3 we show how to construct such a switch path γ, under all the constraints defined earlier. If there is no such blocking position s, we simply set γγ.

We now have a path γ from some s′′S (which may be either s or s) to the target t, which in not blocked by any robots, but may be interrupted. In Section 3.4 we show how to move a robot from s′′ to t along γ, while all interrupting robots (either on start or target positions) move slightly in order to clear the path. For this we introduce two additional constraints, one on the obstacles-separation, and one on the robots-separation.

Finally, we prepare the input for the next recursive iteration, by removing s′′ from S, t from T, and D1ε(t) from 𝒲. Removing D1ε(t) from 𝒲 ensures that the robots placed on t does not block any path that will be computed in any of the next iterations.

In Section 4 we present the complete algorithm based on the building blocks presented in this section, prove its correctness and analyze its running time and approximation factor.

3.1 Geodesics and blocking positions in the free space

We begin by introducing the constraints under which a blocking robot positioned on some pST can enter a geodesic path γ:[0,1] between some start and target positions, via a straight line segment connecting it to the path.

Figure 2: A geodesic path γ in from a start position γ(0)S to a target position γ(1)T is illustrated in red. The obstacles are in gray, and the trace of γ is highlighted in yellow. A position p is blocking/interrupting γ. The point a is the closest point to p on γ.

A geodesic path between two positions x,y is a shortest length path γ:[0,1] where γ(0)=x and γ(1)=y. Denote by 𝒟1(γ)=0t1𝒟1(γ(t)) the trace of γ (see Figure 2). For 0t1t21, denote by γ[t1,t2] the subpath of γ between γ(t1) and γ(t2).

In the full version of the paper [10], we provide some very useful geometric properties of geodesic paths in the free space, which result in the following constraints. Note that we have one constraint on the obstacles-separation, and another on the robots-separation.

Constraint 1.

ω(33ε)2+1

Constraint 2.

​​111Note that later we introduce two more constraints on ρ that together make Constraint 2 redundant, however, Lemma 2 holds already for this weaker constraint and may be of independent interest. ρ122+(332ε)2

Let γ:[0,1] be a geodesic path. We say that a point a on γ is locally closest to p if there exists δ>0 such that Dpa(p)γ[aδ,a+δ]=a.

Lemma 2.

Let γ:[0,1] be a geodesic path such that γ(0),γ(1)ST, and let pST be a position which is ε-blocking γ. Let a=γ(t) be a point on γ which is locally closest to p and such that pa<2ε. If Constraints 1 and 2 hold, then pa¯.

Intuitively, Lemma 2 states that a robot blocking γ and positioned on p can move to a along the path pa¯, as long as no other robots are standing in its way (see Figure 2).

3.2 Existence of an 𝜺-interrupting target

We now define formally a generalization of the standalone goal defined in [15] (Definition 3). Our definition depends on the overlap parameter ε, and it is equivalent to a standalone goal for ε=0. Let Γ={γ1,,γm} be an optimal-assignment path set for S,T,𝒲.

Definition 3.

A target position tT is an ε-interrupting target, if for every path γΓ such that γ(1)t, it holds that t is not ε-blocking γ.

To show that an ε-interrupting target always exists, we need an additional constraint.

Constraint 3.

ρ42ε

The following theorem is a generalization of Theorem 4 in [15]. The proof follows a similar logic, however, we need to take into account some issues that arise from the generalized definition and Constraints 1, 2, and 3. See the full version [10] for a complete proof.

Theorem 4.

Let Γ be an optimal-assignment path set for S,T,𝒲. If Constraints 1, 2, and 3 hold, then there exists an ε-interrupting target.

3.3 Dealing with a blocked geodesic path

Let Γ={γ1,,γm} be an optimal-assignment path set for S,T,𝒲, and assume w.l.o.g. that γi(0)=si for every 1im. Consider a geodesic path γi from si to tj, and assume that is it ε-blocked. Our goal in this section is to describe the switching process: we find another start point skS and a path γk from sk to tj such that γk is not ε-blocked. Moreover, we describe a path in from si to t, where t=γk(1) (see Figure 3). The sum of lengths of these new paths will not be much longer than the original.

Figure 3: sk is the last to ε-block γi. The switch paths are γi (dashed blue) and γk (red).

For a position skS that ε-blocks γi, denote by 0wi,k1 the largest value such that γi(wi,k) is the point on γi locally closest to sk. We say that skS is the last to ε-block γi if for any sxS{sk} that ε-blocks γi, we have wi,kwi,x (i.e. the last locally closest point to sk on γi is the closest to γi(1) along γi).

Let sk be the last to ε-block γi, and set w=wi,k. Consider the paths γk=skγi(w)¯γi[w,1] and γi=γi[0,w]γi(w)sk¯γk (where γk is the geodesic path from sk to t, see Figure 3). We call γi and γk the switch paths of i and k, respectively. Notice that skγi(w)<2ε, and therefore |γk|+|γi|<|γk|+|γi|+42ε. The following theorem will help us in showing that the switch paths are not blocked (see [10] for a full proof).

Theorem 5.

Assume that Constraints 1, 2, and 3 hold. Then (i) both γi and γk are in , (ii) skγi(w)¯ is not ε-blocked by any pTS{sk}, (iii) γk is not ε-blocked by any sjS{sk}, and (iv) if tj is not ε-blocking γk then it does not ε-block γi.

3.4 Dealing with robots interrupting a path

In this subsection, we discuss how to move a robot along an ε-interrupted path. We describe simple motion paths for the interrupting robots to clear the way, and we analyze the extra distance traveled by the interrupting robots. We note that shorter paths can be shown for the interrupting robots, however, as this does not effect the correctness of the algorithm, we chose to keep this part simple and describe a simple (yet rather short) motion path.

As mentioned at the beginning of this section, our motion plan consists of a weakly-monotone collision free path set, that is, robots moves to targets according some ordering one by one; when one robot travels from its starting position to some target, the other robots must stay in a small disk centered at their current position (which is either a start or a target position). To ensure that an interrupting robot is able to clear the way while remaining in a disk of radius 1+ε around a start/target position, we introduce two additional constraints.

Constraint 4.

ω1+ε

Constraint 5.

ρ2+ε

Consider a robot A traveling in constant speed along a path γA from start position sS to some target position tT. Let pST be either a start or a target position, occupied by another robot B, which is ε-interrupting γA. We describe a motion plan for B by constructing a clearance path γ~B:[0,1]2, along which the robot B will travel with constant speed, as follows.

Figure 4: Left: the robot B is initially positioned at p. In red is a maximal subpath of γA where B needs to clear the way. Middle, right: the robot B move inside 𝒟ε(p) while maintaining distance exactly 2 from the robot A.

Let γA[w1,w2] be a maximal subpath where γA(w)p2 for every w[w1,w2], that is, when A travels along γ[w1,w2], the robot B must clear the way. Note that there could be more than one such maximal subpath for B. We define the corresponding subpath γ~B:[w1,w2]2 for B using a polar coordinate system, where p is the pole. Denote by rx the radius and by θx the angle of a point x, and let γ~B(w):=(2rγA(w),θγA(w)+π).

The next lemma shows that the clearance path of B is feasible and relatively short. See [10] for the full proof.

Lemma 6.

γ~B(w)[w1,w2] is a continuous path that starts and ends at p, such that (i) γ~B[w1,w2]Dε(p), and (ii) |γ~B[w1,w2]||γA[w1,w2]| for 0ε1.

Recall that B may interrupt γA in more than one maximal subpath, so overall we define the clearance path γ~B:[0,1]2 as

γ~B(w)={p,if pγA(w)2,(2rγA(w),θγA(w)+π),if pγA(w)<2.

In addition, we construct such a clearance path for any other robot that interrupts γA. In the full version [10] we prove the following theorem, which shows that clearance paths are collision-free, assuming Constraints 1, 2, 3, 4, and 5.

Theorem 7.

Let A be a robot traveling in constant speed along a path γA from start position sS to some target position tT, and let {γ~BBA} be the set of clearance paths defined for the other robots. If Constraints 1, 2, 3, 4, and 5 hold for 0ε1, then the motion plan that corresponds to the paths γA{γ~BBA} is collision-free.

4 Unlabeled MRMP with improved separation trade-offs

In this section we present and analyze the algorithm based on the tools developed in the previous section. The input for the algorithm is an overlap parameter ε, a workspace 𝒲, a set S of starting positions and a set T of target positions, such that Constraints 1, 2, 3, 4, and 5 hold. In addition, we assume that each connected component of contains an equal number of start and target positions (as otherwise, there is no solution).

For the simplicity of the analysis, in this section we provide an iterative version of the algorithm. In the i’th iteration, for 1im, the algorithm constructs a collision-free motion plan for all the robots over the time interval [i1,i], in which a single robot is moving from a start position to a target position along an ε-interrupted path (Section 3.3), and the other robots may move in their small neighborhood to clear the way (Section 3.4). The output path for a robot A will then be the concatenation γA:[0,m] of the paths γA:[i1,i] constructed for it in all iterations. We then update S and T according to the path that was chosen, and remove from 𝒲 a disk of radius 1ε around the target.

At the beginning of this section we gave a high level description of our algorithm, and a more detailed pseudo code is provided in Algorithm 1. Below, we prove that this algorithm returns a collision-free motion plan, with total length O(OPT), where OPT is the optimal size of a solution to the problem.

Algorithm 1 Motion Plan for Unlabeled Disks.
Theorem 8.

Consider an input 𝒲, S, T, and 0ε1 such that Constraints 1, 2, 3, 4, and 5 hold. If each connected component of contains an equal number of start and target positions, then Algorithm 1 is guaranteed to return a collision-free motion plan Γ.

Proof.

First, observe that the solution returned by the algorithm is a path set that matches every start position in S to a unique target in T, because in each iteration we find a path between some sS and tT and then remove them from S and T for the next iteration.

Consider the first iteration of the algorithm. By Theorem 4, an ε-interrupting target tj exists, and let γi be the path from si to tj. If γi is ε-blocked, let skS be the last to block γi, and by Theorem 5(iii), the switch path γk is in , and it is not blocked by any other robot. Otherwise, γi is not blocked. In any case we choose an ε-interrupted path that leads to tj, and by Theorem 7 the motion plan computed in this iteration is collision-free.

In the following iteration, after removing D1ε(tj) from 𝒲, Constraints 1, 2, 3, 4, and 5 still hold for the updated sets S and T: the robots-separation does not change, so Constraints 2, 3, and 5 still hold, and specifically by Constraint 3 for every p1,p2ST we have p1p242ε. Therefore, for every pST{tj} and xD1ε(tj) we have px3ε, which satisfies both Constraint 1 (pxρε) and Constraint 4 (px1+ε), for any 0ε1.

However, we also need to take into account the robot that is now positioned on tj. The main observation that we will use is that tj does not block any path γ in the free space of 𝒲D1ε(tj); the reason is that the distance between any point on aγ and a point pD1ε(tj) is at least 1, and therefore the tja2ε which means that the overlap between tj and a is at most ε.

Consider an iteration a>1 of the algorithm, and denote by 𝒲a,Sa,Ta the sets 𝒲,S,T at the beginning of iteration a. The first step is to compute an optimal-assignment path set Γa. Let Γa1 be the optimal-assignment path set for 𝒲a1,Sa1,Ta1. At the beginning of iteration a, we have 𝒲a=𝒲a1𝒟1ε(tj) and Ta=Ta1{tj}, where tj is an ε-interrupting target in Γa1. Since tj does not block any path in Γa1, the disk 𝒟1ε(tj) does not overlap with the trace of any path in Γa1{γi}. Therefore, any path γΓa{γi} is in the free space of 𝒲a. If γi was not blocked in iteration a1, then Sa=Sa1{si}, and thus Γa{γi} corresponds to a matching of Sa and Ta, which implies that Γa exists. Otherwise, Sa=Sa1{sk}, where skSa1 is the last to ε-block γi. Let t be the target that was matched to it in Γa1, and let γkΓa1 be the path from sk to t. Observe that Γa1{γi,γk} corresponds to a matching between Sa{si,sk} and Ta1{tk,t}, and thus to have a matching between Sa and Ta we need to find a path in the free space of 𝒲a between si and t. We show that the switch path γi from si to t as defined in Section 3.3, has this property. Recall that γi is a concatenation of γ[1,w], γi(w)sk¯, and γk. The paths γi,γk are in Γa1 and therefore lie in the free space of 𝒲a1. By Theorem 5(i), the entire path γi is in the free space of 𝒲a1. We thus need to show that the disk 𝒟1ε(tj) does not overlap with the trace of γi. Indeed, since tj is an ε-interrupting target, it does not block γk, and therefore by Theorem 5(iv) it is also not blocking γi.

The next step is to find an ε-interrupting target tj. For this, notice that Theorem 4 still holds, because it considers only start and target positions in Sa and Ta. Let 𝒜 be the set of robots that are already lying on target positions. Notice that we can also apply Theorem 5 on 𝒲a,Sa,Ta and Γa, and we just need to make sure that no robot in 𝒜 can block the switch path γk. As mentioned above, this is true in general, that is, for any path γΓa from sSa to tTa, no robot in 𝒜 can block γ, because the overlap between any such robot and any point on γ is at most ε. In other words, robots in 𝒜 can only interrupt paths in Γa.

The final step is using Theorem 7 to construct the clearance paths. Notice that this construction applies for any pST, and therefore it can be applied for the ε-interrupted path, even if there were robots in all the start and target positions, except for the endpoints of the path.

Lastly, notice that the concatenation of the paths that correspond to a robot A from all the iterations of the algorithm results in a continuous path. This is because each clearance path starts and ends at the same start/target position.

In the full version of the paper [10] we prove the following theorem.

Theorem 9.

Given a polygonal workspace 𝒲 with n edges, a set S of m starting positions and a set T of m target positions, there is an algorithm that computes in O~(m4+n2m2) time a collision-free motion plan of length O(OPT), where OPT is the length of an optimal solution, or reports that such a motion plan does not exist, assuming that ρ=max{42ε,2+ε} and ω=max{1+ε,(33ε)2+1}.

The length of our solution.

Note that in fact we bound the size of the solution with the term (1+c)(OPT+(42ε)m) where c is the maximum number of start/target positions that may collide a unit disk under our robots-separation constraints. For example, when the robots-separation is 2 then c=6. We discuss this further in the full version of the paper [10].

Optimizing the separation bounds.

We now optimize the overlap parameter over Constraints 1, 2, 3, 4, and 5, to achieve either the smallest obstacles-separation or the smallest robots-separation possible. Recall that Constraints 1 and 4 specify the required obstacles-separation bounds, while Constraints 2, 3, and 5 specify the required robots-separation bounds. We highlight three interesting values of ε:

  • For ε=0 the motion plan is monotone, and the length of the solution is at most OPT+4m. In this case, we obtain ω=13631.614 and ρ=4.

  • For ε>0, the motion plan is weakly-monotone.

  • For ε=1563130.354, the obstacles-separation is minimized, and we obtain ω=1+1563131.354 and ρ=421563133.291.

  • For ε=23, the robots-separation is minimized, and we obtain ω=123 and ρ=223.

Lower bounds.

In the full version of the paper [10] we present lower bounds on the obstacles-separation when the motion plan is required to be (weakly) monotone. We provide an instances of the problem in which (i) ω1.614 and no monotone motion plan exists, and (ii) ω1.354 and no weakly-monotone motion plan exists. This shows that to obtain lower obstacles-separation assumptions, a completely different approach is required.

5 The Exodus algorithm: reducing the robots-separation bound

In this section we introduce a different algorithmic strategy for MRMP which we call the Exodus Algorithm. This strategy is more global in the sense that it also moves robots that do not directly block a path, and therefore it allows for robots-separation of 2, while requiring obstacles-separation of only 3. Whereas the previous algorithm relies on local path-clearing and larger robots-separation, Exodus replaces these local adjustments with a synchronized outward movement of all robots in order to clear the trace of a geodesic path, so that the selected robot can reach its target via this path. Note that in this section, we assume that the workspace 𝒲 is a simple polygon. Due to space considerations, we give an overview of the algorithm and the main result here, see [10] for a formal description and proofs.

Algorithm overview.

Given the sets S,T as before, the first step of the algorithm is to compute an optimal-assignment path set Γ={γ1,,γm} for S,T,𝒲. Assume w.l.o.g. that γi(0)=si and γi(1)=ti for every 1im. The algorithm then proceeds iteratively. In each iteration, we select a target position tiT which is not yet occupied by a robot. We first extend the path γi so that its start and end points are on the boundary of . The extended path partitions into maximally connected sub-regions of , which we call pockets. We then split each pocket using the angle bisectors corresponding to the reflex subchain of γi that defines the pocket. We call the cells in this refined partition Exodus cells. Each Exodus cell contains a single straight-line edge of γi. All the other robots are then moved outwards from the extended path by two units (which is possible because ω=3), and the movement direction for each robot is determined by the Exodus cell in which it lies. This process clears the trace of γi, allowing a robot positioned on si to move along γi until reaching ti, while all the other robots remain fixed at their offset locations. After the robot reaches its target ti, all the other robots return to their original positions by applying the inverse movement (which is possible because ρ=2, i.e., no start position is too close to ti). This iterative “corridor-opening” process motivates the name Exodus Algorithm. The procedure repeats until all robots reach their assigned targets. In each iterative step, m1 robots are moving 4m units each, while the selected robot moves along a shortest path to its assigned target (as before, |Γ|OPT). Since there are exactly m iterations, we obtain the following theorem.

Theorem 10.

Given a simple polygonal workspace 𝒲 with n edges, a set S of m starting positions and a set T of m target positions, such that ρ=2 and ω=3, the Exodus algorithm computes in O~(m3+mn2) time a collision-free motion plan of size at most OPT+4m24m, where OPT is the length of an optimal solution.

Figure 5: The extended geodesic γ^ is in red, and its trace is in blue. In pink is . Each colored region represents a distinct pocket induced by γ^, and the orange rays are the angle bisectors generated at the corresponding reflex vertices of . These bisectors subdivide each pocket into Exodus cells, within which a single direction vector is illustrated by the black arrows.
 Remark 11 (Labeled MRMP).

In the labeled setting, where each start position is already associated with a specific target, we do not compute an optimal-assignment path set. Instead, we simply compute the geodesic paths between each start and its corresponding target using Wang’s Euclidean shortest-path algorithm [17], and then apply the same iterative procedure described above. In this case, the running time becomes O(m2logn+mn2), as no matching step is required. The motion of the robots during the corridor-opening iterations remains unchanged, and the size of the solution remains at most OPT+4m24m.

6 Discussion and future work

In this paper we provide multi-robot motion plans that require tighter assumptions on the density of the input setting. Note that for the Exodus algorithm we require ρ=2. This separation bound is tight for monochromatic robots-separation, however, the bichromatic robots-separation can potentially be 0. This raises the question of what is the minimum obstacles-separation for which we can remove the robots-separation assumption completely.

Figure 6: Unlabeled MRMP with obstacles-separation of 1.5ε, for which no solution exists.

In the full version of the paper [10] we also present lower bounds for the amount of obstacles-separation needed so that a solution to the MRMP problem always exists, without any assumptions on the separation between robots. In the unlabeled case, we prove that obstacles-separation of at least 1.5 is sometimes needed (see Figure 6), and in the labeled case, we observe that obstacles-separation of at least 2 is sometimes necessary (see Figure 7).

Figure 7: Labeled MRMP with obstacles-separation of 2ε, for which no solution exists.

References

  • [1] Mikkel Abrahamsen, Kevin Buchin, Maike Buchin, Linda Kleist, Maarten Löffler, Lena Schlipf, André Schulz, and Jack Stade. Reconfiguration of unit squares and disks: Pspace-hardness in simple settings. In Oswin Aichholzer and Haitao Wang, editors, 41st International Symposium on Computational Geometry, SoCG 2025, Kanazawa, Japan, June 23-27, 2025, volume 332 of LIPIcs, pages 1:1–1:18. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2025. doi:10.4230/LIPIcs.SOCG.2025.1.
  • [2] Aviv Adler, Mark de Berg, Dan Halperin, and Kiril Solovey. Efficient multi-robot motion planning for unlabeled discs in simple polygons. IEEE Transactions on Automation Science and Engineering, 12(4):1309–1317, 2015. doi:10.1109/TASE.2015.2470096.
  • [3] Pankaj K. Agarwal, Mark de Berg, Benjamin Holmgren, Alex Steiger, and Martijn Struijs. Optimal motion planning for two square robots in a rectilinear environment. In Oswin Aichholzer and Haitao Wang, editors, 41st International Symposium on Computational Geometry, SoCG 2025, Kanazawa, Japan, June 23-27, 2025, volume 332 of LIPIcs, pages 5:1–5:17. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2025. doi:10.4230/LIPIcs.SOCG.2025.5.
  • [4] Pankaj K. Agarwal, Tzvika Geft, Dan Halperin, and Erin Taylor. Multi-robot motion planning for unit discs with revolving areas. Computational Geometry, 114:102019, 2023. doi:10.1016/j.comgeo.2023.102019.
  • [5] Pankaj K. Agarwal, Dan Halperin, Micha Sharir, and Alex Steiger. Near-optimal min-sum motion planning for two square robots in a polygonal environment. In David P. Woodruff, editor, Proceedings of the 2024 ACM-SIAM Symposium on Discrete Algorithms, SODA 2024, Alexandria, VA, USA, January 7-10, 2024, pages 4942–4962. SIAM, 2024. doi:10.1137/1.9781611977912.176.
  • [6] Pankaj K. Agarwal, Benjamin Holmgren, and Alex Steiger. Near-optimal min-sum motion planning in a planar polygonal environment. CoRR, abs/2510.21639, 2025. doi:10.48550/arXiv.2510.21639.
  • [7] Bahareh Banyassady, Mark de Berg, Karl Bringmann, Kevin Buchin, Henning Fernau, Dan Halperin, Irina Kostitsyna, Yoshio Okamoto, and Stijn Slot. Unlabeled Multi-Robot Motion Planning with Tighter Separation Bounds. In Xavier Goaoc and Michael Kerber, editors, 38th International Symposium on Computational Geometry (SoCG 2022), volume 224 of Leibniz International Proceedings in Informatics (LIPIcs), pages 12:1–12:16, Dagstuhl, Germany, 2022. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.SoCG.2022.12.
  • [8] Thomas Brocken, G. Wessel van der Heijden, Irina Kostitsyna, Lloyd E. Lo-Wong, and Remco J. A. Surtel. Multi-robot motion planning of k-colored discs is pspace-hard. In Martin Farach-Colton, Giuseppe Prencipe, and Ryuhei Uehara, editors, 10th International Conference on Fun with Algorithms, FUN 2021, Favignana Island, Sicily, Italy, May 30 - June 1, 2021, volume 157 of LIPIcs, pages 15:1–15:16. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2021. doi:10.4230/LIPIcs.FUN.2021.15.
  • [9] Guillermo Esteban, Dan Halperin, Víctor Ruíz, Vera Sacristán, and Rodrigo I. Silveira. Shortest coordinated motion for square robots. In Pat Morin and Subhash Suri, editors, Algorithms and Data Structures - 18th International Symposium, WADS 2023, Montreal, QC, Canada, July 31 - August 2, 2023, Proceedings, volume 14079 of Lecture Notes in Computer Science, pages 430–443. Springer, 2023. doi:10.1007/978-3-031-38906-1_28.
  • [10] Tsuri Farhana, Omrit Filtser, and Shalev Goldshtein. Unlabeled multi-robot motion planning with improved separation trade-offs. CoRR, abs/2603.19502, 2026. arXiv:2603.19502.
  • [11] Sándor P. Fekete, Ramin Kosfeld, Peter Kramer, Jonas Neutzner, Christian Rieck, and Christian Scheffer. Coordinated motion planning: Multi-agent path finding in a densely packed, bounded domain. In Julián Mestre and Anthony Wirth, editors, 35th International Symposium on Algorithms and Computation, ISAAC 2024, Sydney, Australia, December 8-11, 2024, volume 322 of LIPIcs, pages 29:1–29:15. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2024. doi:10.4230/LIPIcs.ISAAC.2024.29.
  • [12] John E Hopcroft, Jacob Theodore Schwartz, and Micha Sharir. On the complexity of motion planning for multiple independent objects; pspace-hardness of the" warehouseman’s problem". The international journal of robotics research, 3(4):76–88, 1984.
  • [13] Israela Solomon and Dan Halperin. Motion planning for multiple unit-ball robots in d. In Marco Morales, Lydia Tapia, Gildardo Sánchez-Ante, and Seth Hutchinson, editors, Algorithmic Foundations of Robotics XIII, Proceedings of the 13th Workshop on the Algorithmic Foundations of Robotics, WAFR 2018, Mérida, Mexico, December 9-11, 2018, volume 14 of Springer Proceedings in Advanced Robotics, pages 799–816. Springer, 2018. doi:10.1007/978-3-030-44051-0_46.
  • [14] Kiril Solovey and Dan Halperin. On the hardness of unlabeled multi-robot motion planning. In Lydia E. Kavraki, David Hsu, and Jonas Buchli, editors, Robotics: Science and Systems XI, Sapienza University of Rome, Rome, Italy, July 13-17, 2015, 2015. doi:10.15607/RSS.2015.XI.046.
  • [15] Kiril Solovey, Jingjin Yu, Or Zamir, and Dan Halperin. Motion planning for unlabeled discs with optimality guarantees. In Lydia E. Kavraki, David Hsu, and Jonas Buchli, editors, Robotics: Science and Systems XI, Sapienza University of Rome, Rome, Italy, July 13-17, 2015, 2015. doi:10.15607/RSS.2015.XI.011.
  • [16] Paul G. Spirakis and Chee-Keng Yap. Strong np-hardness of moving many discs. Inf. Process. Lett., 19(1):55–59, 1984. doi:10.1016/0020-0190(84)90130-3.
  • [17] Haitao Wang. A new algorithm for euclidean shortest paths in the plane. In Samir Khuller and Virginia Vassilevska Williams, editors, STOC ’21: 53rd Annual ACM SIGACT Symposium on Theory of Computing, Virtual Event, Italy, June 21-25, 2021, pages 975–988. ACM, 2021. doi:10.1145/3406325.3451037.