Unlabeled Multi-Robot Motion Planning with Improved Separation Trade-Offs
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 , and at least 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 and the obstacles-separation is , or (ii) the robots-separation is and the obstacles-separation . Additionally, we introduce a different strategy yielding a polynomial-time solution when the robots-separation is only , and the obstacles-separation is . Finally, we show that without any robots-separation assumption, obstacles-separation of at least may be necessary for a solution to exist.
Keywords and phrases:
multi-robot motion planningFunding:
Omrit Filtser: This research was supported by the ISRAEL SCIENCE FOUNDATION (grant No. 2135/24).Copyright and License:
2012 ACM Subject Classification:
Theory of computation Computational geometryAcknowledgements:
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 NayyeriSeries and Publisher:
Leibniz International Proceedings in Informatics, Schloss Dagstuhl – Leibniz-Zentrum für Informatik
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 starting positions and a set of 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 vertices and unit-disk robots, when the robots-separation, denoted , is (i.e., all the start and target positions are at least distance apart), a solution for unlabeled MRMP always exists, and can be found in 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 , and the obstacles-separation, denoted , is . Under these assumptions, they present an algorithm that runs in time and returns a solution of total length at most . 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 (), 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 is associated with a revolving area - a disk of radius that contains the unit disk centered at 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 , a polygonal environment of complexity , and unit-disc robots, our first algorithm computes in time a collision-free motion plan of total length , assuming that and (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 is sufficient when . 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 , 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 .
Note that our first algorithm does not cover the case of , which is the optimal (monochromatic) robots-separation. We therefore present a second algorithm using a completely different approach, which given a simple polygon with vertices and unit-disk robots, computes in time a collision-free motion plan of total length , assuming that and (see Theorem 10). Finally, we show that may be necessary for a solution to exist.
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 of start positions, a set 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 and , there is always at least one target which is a so-called standalone goal, that is, a target 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 , in which case we can make a switch: we push the last robot that blocks this path to . This switching process is possible because of the assumption . 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 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 , 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 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 , with obstacles-separation of , and it applies also for labeled MRMP (see Theorem 10).
Note that the monochromatic separation must be at least , so the robots-separation is almost optimal (the bichromatic separation can be ). 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 unit-disk robots, moving in a polygonal workspace (which may be simple or cluttered with polygonal obstacles), whose overall number of edges is . Below, we mostly follow the notations in [7] and [15].
The obstacle space is the complement of the workspace , that is, . We refer to the points in as positions, and say that a robot is at position when its center is positioned at the point . For a point and a radius , we define to be the open disk of radius centered at . 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 . In addition, a robot collides with the obstacle space if and only if its center is at distance strictly less than 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 . 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 edges and two sets and of points in the free space . The goal is to plan a collision-free motion for robots, each positioned on a different starting point in , such that by the end of the motion every target position in is occupied by some robot. When the robots are distinguishable (i.e. labeled), the robot that starts at position is required to end up at the target . 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 , where for every , and such that and . In addition, the paths need to be collision-free, i.e., at any point in time , there are no such that and are at distance strictly less than . Denote the size of a solution by , where is the length of a path in the 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 and every point , it holds that . We denote by the robots-separation assumption, i.e., we assume that for every it holds that .
Overlapping, blocking and interrupting positions
Let be two points in the free space. We define the overlap of and to be (see Figure 1). If , then their overlap is . Otherwise, we say that are -overlapping for , and note that .
We say that a position is -blocking a path if there exists such that the overlap of and 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 is -interrupting a path if for every , the overlap of and 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 the overlap of and is , then for every it holds that , and therefore in this case we say that is -interrupting . We say that a path is -blocked if there exists that -blocks it, and that is -interrupted if no is -blocking it.
The overlap parameter.
For the analysis in Section 3, we fix a parameter 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 of starting positions, a set of target positions, and a workspace . In each recursive iteration, one robot is moving to a target position, then , 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 geodesic paths , where for every , and such that , . We say that is an optimal-assignment path set for if it minimizes 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 be an optimal-assignment path set for .
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 which does not block any path in (other than the path from to ). Then we can move a robot from to along , and it will not block other paths in the following iterations. Note that 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 , which blocks the path , we need to make a switch, i.e., we construct a new path in from to , 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 , we simply set .
We now have a path from some (which may be either or ) to the target , which in not blocked by any robots, but may be interrupted. In Section 3.4 we show how to move a robot from to 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 from , from , and from . Removing from ensures that the robots placed on 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 can enter a geodesic path between some start and target positions, via a straight line segment connecting it to the path.
A geodesic path between two positions is a shortest length path where and . Denote by the trace of (see Figure 2). For , denote by the subpath of between and .
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.
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.
Let be a geodesic path. We say that a point on is locally closest to if there exists such that .
Lemma 2.
Let be a geodesic path such that , and let be a position which is -blocking . Let be a point on which is locally closest to and such that . If Constraints 1 and 2 hold, then .
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 . Let be an optimal-assignment path set for .
Definition 3.
A target position is an -interrupting target, if for every path such that , it holds that is not -blocking .
To show that an -interrupting target always exists, we need an additional constraint.
Constraint 3.
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 . If Constraints 1, 2, and 3 hold, then there exists an -interrupting target.
3.3 Dealing with a blocked geodesic path
Let be an optimal-assignment path set for , and assume w.l.o.g. that for every . Consider a geodesic path from to , and assume that is it -blocked. Our goal in this section is to describe the switching process: we find another start point and a path from to such that is not -blocked. Moreover, we describe a path in from to , where (see Figure 3). The sum of lengths of these new paths will not be much longer than the original.
For a position that -blocks , denote by the largest value such that is the point on locally closest to . We say that is the last to -block if for any that -blocks , we have (i.e. the last locally closest point to on is the closest to along ).
Let be the last to -block , and set . Consider the paths and (where is the geodesic path from to , see Figure 3). We call and the switch paths of and , respectively. Notice that , and therefore . 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 and are in , (ii) is not -blocked by any , (iii) is not -blocked by any , and (iv) if is not -blocking then it does not -block .
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 around a start/target position, we introduce two additional constraints.
Constraint 4.
Constraint 5.
Consider a robot traveling in constant speed along a path from start position to some target position . Let be either a start or a target position, occupied by another robot , which is -interrupting . We describe a motion plan for by constructing a clearance path , along which the robot will travel with constant speed, as follows.
Let be a maximal subpath where for every , that is, when travels along , the robot must clear the way. Note that there could be more than one such maximal subpath for . We define the corresponding subpath for using a polar coordinate system, where is the pole. Denote by the radius and by the angle of a point , and let .
The next lemma shows that the clearance path of is feasible and relatively short. See [10] for the full proof.
Lemma 6.
is a continuous path that starts and ends at , such that (i) , and (ii) for .
Recall that may interrupt in more than one maximal subpath, so overall we define the clearance path as
In addition, we construct such a clearance path for any other robot that interrupts . 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 be a robot traveling in constant speed along a path from start position to some target position , and let be the set of clearance paths defined for the other robots. If Constraints 1, 2, 3, 4, and 5 hold for , then the motion plan that corresponds to the paths 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 of starting positions and a set 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 ’th iteration, for , the algorithm constructs a collision-free motion plan for all the robots over the time interval , 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 will then be the concatenation of the paths constructed for it in all iterations. We then update and according to the path that was chosen, and remove from a disk of radius 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 , where OPT is the optimal size of a solution to the problem.
Theorem 8.
Consider an input , , , and 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 to a unique target in , because in each iteration we find a path between some and and then remove them from and for the next iteration.
Consider the first iteration of the algorithm. By Theorem 4, an -interrupting target exists, and let be the path from to . If is -blocked, let be the last to block , and by Theorem 5(iii), the switch path is in , and it is not blocked by any other robot. Otherwise, is not blocked. In any case we choose an -interrupted path that leads to , and by Theorem 7 the motion plan computed in this iteration is collision-free.
In the following iteration, after removing from , Constraints 1, 2, 3, 4, and 5 still hold for the updated sets and : the robots-separation does not change, so Constraints 2, 3, and 5 still hold, and specifically by Constraint 3 for every we have . Therefore, for every and we have , which satisfies both Constraint 1 () and Constraint 4 (), for any .
However, we also need to take into account the robot that is now positioned on . The main observation that we will use is that does not block any path in the free space of ; the reason is that the distance between any point on and a point is at least , and therefore the which means that the overlap between and is at most .
Consider an iteration of the algorithm, and denote by the sets at the beginning of iteration . The first step is to compute an optimal-assignment path set . Let be the optimal-assignment path set for . At the beginning of iteration , we have and , where is an -interrupting target in . Since does not block any path in , the disk does not overlap with the trace of any path in . Therefore, any path is in the free space of . If was not blocked in iteration , then , and thus corresponds to a matching of and , which implies that exists. Otherwise, , where is the last to -block . Let be the target that was matched to it in , and let be the path from to . Observe that corresponds to a matching between and , and thus to have a matching between and we need to find a path in the free space of between and . We show that the switch path from to as defined in Section 3.3, has this property. Recall that is a concatenation of , , and . The paths are in and therefore lie in the free space of . By Theorem 5(i), the entire path is in the free space of . We thus need to show that the disk does not overlap with the trace of . Indeed, since is an -interrupting target, it does not block , and therefore by Theorem 5(iv) it is also not blocking .
The next step is to find an -interrupting target . For this, notice that Theorem 4 still holds, because it considers only start and target positions in and . Let be the set of robots that are already lying on target positions. Notice that we can also apply Theorem 5 on and , and we just need to make sure that no robot in can block the switch path . As mentioned above, this is true in general, that is, for any path from to , 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 .
The final step is using Theorem 7 to construct the clearance paths. Notice that this construction applies for any , 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 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 edges, a set of starting positions and a set of target positions, there is an algorithm that computes in time a collision-free motion plan of length , where OPT is the length of an optimal solution, or reports that such a motion plan does not exist, assuming that and .
The length of our solution.
Note that in fact we bound the size of the solution with the term where 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 then . 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 the motion plan is monotone, and the length of the solution is at most . In this case, we obtain and .
-
For , the motion plan is weakly-monotone.
-
For , the obstacles-separation is minimized, and we obtain and .
-
For , the robots-separation is minimized, and we obtain and .
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) and no monotone motion plan exists, and (ii) 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 , while requiring obstacles-separation of only . 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 as before, the first step of the algorithm is to compute an optimal-assignment path set for . Assume w.l.o.g. that and for every . The algorithm then proceeds iteratively. In each iteration, we select a target position which is not yet occupied by a robot. We first extend the path 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 that defines the pocket. We call the cells in this refined partition Exodus cells. Each Exodus cell contains a single straight-line edge of . All the other robots are then moved outwards from the extended path by two units (which is possible because ), and the movement direction for each robot is determined by the Exodus cell in which it lies. This process clears the trace of , allowing a robot positioned on to move along until reaching , while all the other robots remain fixed at their offset locations. After the robot reaches its target , all the other robots return to their original positions by applying the inverse movement (which is possible because , i.e., no start position is too close to ). This iterative “corridor-opening” process motivates the name Exodus Algorithm. The procedure repeats until all robots reach their assigned targets. In each iterative step, robots are moving units each, while the selected robot moves along a shortest path to its assigned target (as before, ). Since there are exactly iterations, we obtain the following theorem.
Theorem 10.
Given a simple polygonal workspace with edges, a set of starting positions and a set of target positions, such that and , the Exodus algorithm computes in time a collision-free motion plan of size at most , where OPT is the length of an optimal solution.
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 , 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 .
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 . This separation bound is tight for monochromatic robots-separation, however, the bichromatic robots-separation can potentially be . This raises the question of what is the minimum obstacles-separation for which we can remove the robots-separation assumption completely.
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 is sometimes needed (see Figure 6), and in the labeled case, we observe that obstacles-separation of at least is sometimes necessary (see Figure 7).
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 . 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.
