---
title: Multiobjective Temporal Shortest Path Problem
url: https://www.emergentmind.com/topics/multiobjective-temporal-shortest-path-problem
type: topic
---

# Multiobjective Temporal Shortest Path Problem

Searching arXiv for the cited papers and closely related formulations.
The multiobjective temporal shortest path problem denotes a family of path problems on temporal or time-augmented graphs in which feasible paths must respect time and are evaluated under several objectives simultaneously. In the single-source formulation, given a directed discrete-time temporal graph \(G=(V,R)\), a start node \(s\in V\), and \(p\geq 1\) objectives, the task is, for each \(v\in V\), to compute the set of nondominated images of temporal \(s\)-\(v\) paths together with one corresponding efficient path for each image; alternatively, one may have to detect the existence of an improving zero-duration cycle [2605.05954]. Closely related formulations model dynamic obstacles by safe intervals [2108.00745], encode timetable uncertainty through interval arc costs and scenario sets [2412.14654], or treat concurrent routing queries as mutually interacting through edge loads in temporal load-aware road networks [2110.09937].

## 1. Core definitions and solution concepts

A directed discrete-time temporal graph is defined as \(G=(V,R)\), where \(V\) is a finite set of nodes and \(R\subseteq V\times V\times \mathbb{Q}_{\ge 0}\times \mathbb{Q}_{\ge 0}\) is a finite set of arcs. Each arc \(r\) has a start node \(\alpha(r)\), an end node \(\omega(r)\), a start time \(\tau(r)\), and a traversal time \(\lambda(r)\). A temporal path is a sequence
\[
P=(v_0,r_1,v_1,\dots,r_k,v_k)
\]
such that \(\alpha(r_i)=v_{i-1}\), \(\omega(r_i)=v_i\), and \(\tau(r_i)+\lambda(r_i)\le \tau(r_{i+1})\) for each \(i\); a zero-arcs path \((v_0)\) is also allowed [2605.05954].

The objective model in the general single-source setting is also explicit. Each objective is a tuple
\[
(M,\le,f,\oplus,\overline 0,dir),
\]
where \(M\) is a totally ordered set, \(f:R\to M\) assigns arc-values, \(\oplus\) is associative with left-neutral \(\overline 0\), and \(dir\in\{\min,\max\}\) specifies preference direction. For \(p\) objectives, the image of a path \(P\) is
\[
\bigl(f_1(P),\dots,f_p(P)\bigr)\in \prod_{j=1}^p M_j.
\]
Weak domination \(y'\preceq y\) holds componentwise according to each objective direction, and strict domination \(y'\prec y\) means weak domination with inequality in at least one component. An \(s\)-\(v\) path is efficient if no other \(s\)-\(v\) path strictly dominates its image; the image of an efficient path is nondominated [2605.05954].

In the safe-interval formulation for dynamic obstacles, the same Pareto logic is expressed on a time-augmented graph \(G^t=G\times\{0,1,2,\dots,T\}\). For a single agent with start \(v_0\) and goal \(v_g\), the aim is to find all collision-free paths \(\pi\) that minimize \(M\) nonnegative criteria simultaneously, with cost vector
\[
c(\pi)=(c_1(\pi),c_2(\pi),\dots,c_M(\pi))\in (\mathbb{R}^+)^M.
\]
A path \(\pi^*\) is Pareto-optimal if there is no other collision-free \(\pi\) with \(c(\pi)\succeq c(\pi^*)\), and the target set is
\[
T^*=\{\pi\mid \pi \text{ is collision-free, Pareto-optimal, and cost-unique}\}.
\]
This formulation makes explicit that temporal feasibility and Pareto efficiency can be imposed simultaneously even when occupancy constraints vary over time [2108.00745].

## 2. Temporal modeling variants

The literature represented here uses several temporal abstractions, each tailored to a different source of temporal heterogeneity.

| Formulation | Temporal representation | Additional state |
|---|---|---|
| SSMTSPP | Directed discrete-time temporal graph \(G=(V,R)\) | Arc start times and traversal times |
| MO-SIPP / MO-CBS | Time-augmented graph \(G^t\) with safe intervals | Vertex occupancies and interval states |
| Collective routing | Static directed graph \(G=(V,E)\) with continuous time and discrete partition \(T\) | Edge–Load–Matrix \(L(e,\tau)\) |
| Tropical Dijkstra | Directed event-activity graph | Interval arc costs \([\ell_e,u_e]\) and scenarios |

In the load-aware road-network model, the topology remains static while costs become time- and load-dependent. The graph is \(G=(V,E)\), time is continuous \(a\in \mathbb{R}_{\ge 0}\) with a discrete partition \(T=\{\tau_1,\tau_2,\dots,\tau_T\}\), and each edge \(e\) carries a capacity \(F_e\), free-flow travel-time \(\Upsilon_e\), speed-limit \(z_e\), and length \(\delta_e\). The edge-load function \(\ell_e:E\times T\to \mathbb{N}\) records the number of vehicles scheduled to traverse edge \(e\) in interval \(\tau\), and these values are collected in the Edge–Load–Matrix (ELM), \(L(e,\tau)=\ell_e(\tau)\) [2110.09937].

In the safe-interval model, time is discrete and occupancy is externalized as dynamic obstacles. For each vertex \(u\), the occupied times are \(Occ(u)\subseteq\{0,\dots,T\}\), and safe intervals are the maximal contiguous intervals in \([0,\infty)\setminus Occ(u)\). A SIPP state is therefore not simply a vertex-time pair but a pair \((u,[t_a,t_b])\), where \([t_a,t_b]\) is one of the non-overlapping safe intervals at \(u\) [2108.00745].

In the interval-cost passenger-routing model, the graph consists of events and activities. Nodes represent departure or arrival events of a vehicle at a stop, while arcs represent a driving leg, a dwelling time, a transfer between lines, or an access/egress between a street location and a stop. Each arc \(e\) has an interval cost \([\ell_e,u_e]\subset \mathbb{R}_{\ge 0}\), and a scenario is a choice \(c_e\in[\ell_e,u_e]\) for every arc. This shifts the temporal dimension from explicit time expansion toward a scenario-parametric shortest-path model [2412.14654].

A plausible implication is that the term “multiobjective temporal shortest path problem” encompasses not a single canonical encoding but a class of models in which temporal feasibility, uncertain travel times, dynamic occupancy, or endogenous congestion determine how efficient paths are defined and computed.

## 3. Dominance structure, objective interactions, and complexity

A central distinction in the recent theory is between monotonicity and isotonicity. An objective is monotone when extending a path cannot improve its value in the preferred direction merely by appending an arc; it is isotonic when the order of two path values is preserved if both paths are extended by the same arc. Under both properties, Dijkstra-style label-setting is possible. Without them, zero-duration temporal cycles may generate infinitely many strictly better images or may have to be traversed an arbitrarily large finite number of times to obtain certain nondominated images [2605.05954].

This observation rules out a common static-graph intuition. When monotonicity or isotonicity is missing, one cannot fix in advance a small upper bound, such as \(n-1\), on the length of all efficient paths. The restricted problem variant SSMTSPP-MPL therefore imposes a maximum admissible path length \(K\) and considers only paths containing at most \(K\) arcs. The same paper establishes several sufficient conditions under which such a bound is not required: if no zero-duration cycle reachable from \(s\) exists, then \(K=|R|\) suffices; if each node enforces a positive minimum waiting time \(\Delta(v)>0\), then again \(K=|R|\) suffices; if at each node there are at most \(\kappa\) distinct images of \(s\)-\(v\) paths, both algorithms stabilize within \(m\kappa\) iterations [2605.05954].

The load-aware collective formulation introduces a different objective interaction. For \(Q\) queries \(q_k\), one seeks paths \(p_k\) that collectively optimize three desiderata: minimize total or average travel-time,
\[
\min_x \sum_{k=1}^Q a_d^k,
\]
ensure fairness through the congestion penalty
\[
\pi_k = |p_k|-|\phi_k|,
\]
and balance edge loads through
\[
\sigma(e,\tau)=\min\{\ell_e(\tau)/F_e,1\},\quad
FFCU=\frac{1}{E\,T}\sum_{e,\tau}\sigma(e,\tau),\quad
LD=\frac{1}{E\,T}\sum_{e,\tau}\mathbf{1}\{\ell_e(\tau)>0\}.
\]
A weighted-sum or lexicographic multiobjective can be formed, although the reported implementation focuses on minimizing \(\sum_k a_d^k\) while reporting fairness and load-utilization as secondary metrics [2110.09937].

The same collective routing problem is stated to be NP-hard. The argument given is that even in the static case, routing many origin-destination demands so as to minimize total travel-time under capacity constraints is known to be NP-hard, and the addition of time and dynamic congestion only enriches the problem [2110.09937]. In the interval-cost setting, complexity is instead governed by the size of the label set \(\lambda=\max_{v\in V}|\mathcal L(v)|\), yielding a worst-case running time
\[
O\bigl(|E|\cdot \lambda^2\cdot T_{\rm eval}\bigr)
\]
for tropical Dijkstra [2412.14654].

## 4. Algorithmic paradigms

When all objectives are monotone and isotonic and the first objective is earliest-arrival, a label-setting algorithm can be used. Labels are prioritized by the first objective, permanent labels are stored per node, and temporally feasible outgoing arcs are relaxed only if the resulting label is not dominated. The complexity bound reported for this setting is
\[
O\bigl(L_{\max}(1+\overline d)\log L_{\max}\bigr),
\]
where \(L_{\max}\) is the total number of permanent labels and \(\overline d=m/n\) is the average out-degree [2605.05954].

If isotonicity holds but monotonicity does not, the problem is handled by a multiobjective temporal label-correcting algorithm for SSMTSPP-MPL. It maintains, for each node \(v\) and each length \(k\), a set \(L(v,k)\) of candidate labels with path-length exactly \(k\), and dominated labels may still be discarded immediately because isotonicity guarantees that they cannot lead to new nondominated images. If neither monotonicity nor isotonicity holds, a fully general label-correcting algorithm keeps all generated labels of length at most \(K\) and removes dominated labels only once no label-sets grow any further; its worst-case running time is
\[
O\bigl(K\,(L_{\max}+|R|)\bigr),
\]
with space \(O(K\,L_{\max})\) [2605.05954].

A different low-level strategy appears in the dynamic-obstacle setting. MO-SIPP associates each state \(s=(u,[t_a,t_b])\) with a set of labels \(l=(s,\vec g,t_r)\), where \(\vec g\) is the cost-to-come vector and \(t_r\) is the arrival time. The algorithm uses vector dominance for goal pruning and a specific label-dominance relation at the same state,
\[
l \succeq_\ell l' \iff t_r(l)\le t_r(l')\ \text{and}\ 
\vec g(l)+(t_r(l')-t_r(l))\cdot \vec c_w \le \vec g(l')
\]
component-wise. The stated theorem is that, upon termination, MO-SIPP generates exactly the set \(T^*\) of all cost-unique Pareto-optimal paths from \(v_0\) to \(v_g\) in the presence of dynamic obstacles [2108.00745].

The interval-cost formulation replaces vector labels by tropical polynomials. For each node \(v\), tropical Dijkstra maintains a label set \(\mathcal L(v)\subset \mathbb T[x_1,\dots,x_m]\) of non-dominated polynomials. The path monomial of a path \(p\) is
\[
\tau(p)=\bigotimes_{e\in p}x_e=\sum_{e\in p}x_e,
\]
and the path polynomial of a path set \(\mathcal P\) is
\[
L(\mathcal P)(c)=\bigoplus_{p\in\mathcal P}\tau(p)(c)=\min_{p\in\mathcal P}c(p).
\]
Dominance is defined pointwise over all scenarios, and the algorithm returns complete or essential sets by changing the dominance relation used for pruning [2412.14654].

The load-aware road-network variant extends Dijkstra and A* by replacing static edge weights with a FIFO-compliant arrival-time function
\[
f_{ij}(a_i)=\tau_i+(a_i-\tau_i)^{\epsilon_e(\tau_i)}+\Upsilon_{ij}=a_j,
\]
where \(\tau_i=\lfloor a_i\rfloor\) and
\[
\epsilon_e(\tau)=
\begin{cases}
1, & \ell_e(\tau)\le F_e\\[4pt]
\dfrac{1}{\ell_e(\tau)-F_e}, & \ell_e(\tau)>F_e.
\end{cases}
\]
The paper states that this gives a FIFO network, implies that subpaths of shortest paths are shortest, and ensures that no waiting at a node can improve arrival time; the worst-case time for the Temporal–Load–Aware A* variant is \(O(V\log V + ET)\) [2110.09937].

## 5. Specialized formulations: safe intervals, interval scenarios, and collective reassignment

The safe-interval line of work is motivated by collision avoidance in dynamic environments. In that setting, waiting is explicitly represented as a graph action, each edge has a fixed \(M\)-vector cost, waiting for one time-step incurs \(\vec c_w\), and the objective is to enumerate all cost-unique Pareto-optimal collision-free paths for a single agent. MO-SIPP is then embedded in MO-CBS as the low-level search procedure for multi-objective conflict-based search, yielding a formulation referred to as multi-objective MAPF (MOMAPF) [2108.00745].

The interval-cost line of work starts from public transportation. The Complete Interval Shortest Path Problem asks, for each target \(t\), for the union of all \(s\to t\) shortest paths over all scenarios,
\[
\mathcal P_t^{\rm compl}
=\bigcup_{c\in[\ell,u]}
\{\text{all }s\to t\text{ shortest paths w.r.t. }c\},
\]
whereas the Essential Interval Shortest Path Problem seeks a scenario-covering minimal subset \(\mathcal P_t^{\rm ess}\subseteq \mathcal P_t^{\rm compl}\) such that for every scenario at least one path in \(\mathcal P_t^{\rm ess}\) is shortest. The paper further states that C-ISPP and E-ISPP can be cast as a multiobjective shortest-path problem with \(d=2^{|E|}\)-dimensional costs, in which C-ISPP corresponds to all weakly-efficient paths and E-ISPP to a smallest efficient cover of the Pareto front [2412.14654].

The collective routing line of work changes the problem semantics more substantially. A \(k\)-th query \(q_k=(s_k,d_k,t_k)\) departs \(s_k\) at time \(t_k\) and must reach \(d_k\) via a path \(p_k=\{(e,\tau)\}\). Because assigned paths update the ELM and therefore future travel times, shortest-path queries are no longer isolated tasks. The CS-MAT heuristic repeatedly selects, from a batch of unassigned queries, the query whose earliest-feasible load-aware shortest path gives the lowest arrival time at its destination; once that path is assigned, the ELM is updated and only queries whose previously computed best path intersected the newly assigned path in space-time are re-evaluated. The method is described as “almost embarrassingly parallel” and is reported to scale to hundreds of thousands of queries on modest clusters [2110.09937].

This suggests a useful conceptual distinction. Some formulations seek the full nondominated image set for a fixed source or source-goal pair, while others optimize a system-wide objective over many temporally interacting paths. Both remain temporal and multiobjective, but the latter couples path computations through shared network state.

## 6. Empirical findings, misconceptions, and open questions

The empirical record in these papers is heterogeneous because the tasks differ. In the collective load-aware setting, experiments on Porto and New York datasets during the AM peak report that Average Journey Time was reduced by up to \(63\%\) over naïve free-flow Dijkstra and by up to \(\sim 20\%\) over static load-aware (SLAD); absolute savings are reported as approximately \(4.8\) minutes per trip in Porto and approximately \(3.4\) minutes in New York. The same study reports lower mean and lower standard deviation of the congestion penalty \(\pi_k\) than TLAA* under high congestion, increases of up to \(30\)–\(40\%\) in Free-Flow Capacity Utilization and Load Distribution relative to baselines, and routine use of more than \(97\%\) of network capacity under heavy load. For \(50{,}000\) queries in NYC with \(6\) cores, end-to-end runtime including journey is reported as \(13.95\) minutes for CS-MAT, \(14.66\) minutes for TLAA*, and \(23.16\) minutes for Dijkstra [2110.09937].

In the safe-interval setting, MO-SIPP was compared, for single-agent performance, to a standard multi-objective A* (NAMOA*) over \(G^t\). On benchmark grids with random risk fields and \(M=2,3\) objectives, it achieved an order-of-magnitude reduction in planning time per call and similar or smaller frontier sizes because of time-dimension compression via safe intervals. The paper states, for example, a reduction from \(10\) seconds to \(0.5\) seconds on a room map with \(M=2\) [2108.00745].

In the tropical Dijkstra study, computational experiments were carried out on \(17\) condensed event-activity instances derived from the Wuppertal public-transport system and \(26\) TimPassLib instances. Small networks with \(|E|<1000\) finished in milliseconds even without transfer caps; the largest Wuppertal instance, with approximately \(13{,}000\) arcs, took approximately \(15\) minutes per source for the complete set and approximately \(10\) minutes for the essential set. Restricting to at most \(3\) transfers sped up all instances by one to two orders of magnitude, making full city-scale networks solvable in less than \(1\) minute per source. The number of complete paths per OD pair ranged from fewer than \(2\) to approximately \(1000\), essential sets were \(80\)–\(95\%\) of complete sets on Wuppertal, and empirically \(\lambda\) stayed in the low tens [2412.14654].

Two misconceptions are directly contradicted by the available results. First, a static-graph bound on efficient path length does not carry over automatically to temporal multiobjective routing; Example 1 in the label-correcting paper shows that an efficient path can require exponential length in \(|R|\) when objectives are non-monotone, and Example 2 shows that non-simple combinations of zero-duration cycles may be needed to generate a new nondominated image [2605.05954]. Second, dominated labels cannot always be discarded immediately: this is valid for isotonic objectives, but when isotonicity is absent the fully general label-correcting algorithm keeps all generated labels until stabilization before deleting dominated ones [2605.05954].

Open questions are also explicit. For tropical Dijkstra, the stated problems include tight worst-case bounds on \(\lambda\) for structured networks, integration into combined timetabling-routing optimizers such as preprocessing in modulo network-simplex timetabling, and handling truly dynamic costs or time-dependent intervals [2412.14654]. In the broader temporal multiobjective setting, the role of zero-duration cycles, the need for admissible bounds \(K\), and the distinction between monotone, isotone, and fully general objectives remain the principal structural fault lines [2605.05954].

Source: https://www.emergentmind.com/topics/multiobjective-temporal-shortest-path-problem