---
title: Bi-Objective Traveling Thief Problem (BI-TTP)
url: https://www.emergentmind.com/topics/bi-objective-traveling-thief-problem-bi-ttp
type: topic
---

# Bi-Objective Traveling Thief Problem (BI-TTP)

The Bi-Objective Traveling Thief Problem (BI-TTP) is a bi-objective variant of the Traveling Thief Problem in which routing and packing must be optimized simultaneously under explicit interdependence between the traveling salesperson problem and the 0-1 knapsack problem. In the cited literature, a solution consists of a tour together with a packing plan, and the Pareto front is the set of non-dominated trade-offs between the chosen objectives [1802.02434][2002.04303][2011.05081][2603.18038]. A central structural feature is that knapsack weight accumulated along the route reduces travel speed, so item-selection decisions alter the effective cost of subsequent edges. The literature also shows that BI-TTP is not represented by a single universally fixed objective pair: one line of work uses a reward-versus-weight formulation, whereas others use profit-versus-travel-time formulations. This suggests that BI-TTP is best understood as a family of closely related bi-objective formulations built on the same route-packing coupling.

## 1. Formal definitions and alternative objective formulations

A standard BI-TTP instance specifies a set of cities, a set of items distributed over those cities, a knapsack capacity, distances, and the maximal and minimal velocities. In one formulation, a feasible solution is a pair \((\pi,z)\) where \(\pi=\langle \pi_1,\dots,\pi_n\rangle\) is a permutation of the cities with \(\pi_1=1\), and \(z\in\{0,1\}^m\) is a binary picking vector satisfying the capacity constraint \(\sum_{j=1}^m w_j z_j \le W\). The cumulative knapsack weight after visiting city \(\pi_i\) is
\[
\omega(i,\pi,z)=\sum_{k=1}^i\sum_{j\in I:\,\mathrm{loc}(j)=\pi_k} w_j z_j,
\]
and the total travel time is
\[
f_2(\pi,z)=\sum_{i=1}^{n-1}\frac{d(\pi_i,\pi_{i+1})}{v_{\max}-\bigl(v_{\max}-v_{\min}\bigr)\,\omega(i,\pi,z)/W}
+\frac{d(\pi_n,\pi_1)}{v_{\max}-\bigl(v_{\max}-v_{\min}\bigr)\,\omega(n,\pi,z)/W}.
\]
The corresponding profit objective is
\[
f_1(\pi,z)=\sum_{j=1}^m p_j z_j.
\]
This yields the bi-objective optimization problem \(\max (f_1(\pi,z),-f_2(\pi,z))\), or equivalently \(\min (f_1,f_2)\) after sign conversion for profit [2011.05081].

A closely related formulation uses the same route-packing structure but writes the objectives as total travel time and item profit directly. For a tour \(\mathcal A=(\alpha_1,\dots,\alpha_N)\) and packing \(\mathcal Z\), the cumulative weight is
\[
W_i(\mathcal A,\mathcal Z)=\sum_{j=1}^i \sum_{k=1}^{B_{\alpha_j}} w_{\alpha_j,k}\,z_{\alpha_j,k},
\]
the travel speed on leg \(i\to i+1\) is
\[
v_{i,i+1}=v_{\max}-\frac{W_i}{W}(v_{\max}-v_{\min}),
\]
the travel-time objective is
\[
f_1(\mathcal A,\mathcal Z)=\sum_{i=1}^N \frac{d_{\alpha_i,\alpha_{i+1}}}{v_{\max}-(W_i/W)(v_{\max}-v_{\min})},
\]
and the profit objective is \(P(\mathcal A,\mathcal Z)\), or in minimization form \(f_2(\mathcal A,\mathcal Z)=-P(\mathcal A,\mathcal Z)\) [2603.18038].

Wu et al. introduced a different bi-criteria formulation in which the first objective is a net reward
\[
f(\pi,\rho)=\text{profit}-R\cdot \text{travel time},
\]
and the second objective is the total picked weight
\[
\varphi(\rho)=\sum w\,\rho,
\]
subject to \(\varphi(\rho)\le C\). In that setting, BI-TTP seeks the Pareto set of all non-dominated \((\pi,\rho)\) such that no other feasible solution can improve \(f\) without worsening \(\varphi\), or vice versa [1802.02434].

## 2. Coupling of routing and packing

The essential difficulty of BI-TTP is the interweaving of its two classical subproblems. If the knapsack is ignored, the problem reduces to a traveling salesperson problem over the cities. If the tour is fixed and travel-time effects are ignored, the problem reduces to a 0-1 knapsack problem over the items. In BI-TTP, however, each stolen item increases the current knapsack weight, which decreases the current speed and therefore increases the travel time of later legs [2002.04303].

This dependence is formalized through a velocity function of the current load. One formulation writes
\[
v(W)=
\begin{cases}
v_{\max}-\dfrac{W}{Q}(v_{\max}-v_{\min}) & W\le Q,\\[4pt]
v_{\min} & W>Q \quad (\text{infeasible}),
\end{cases}
\]
so a packing choice does not merely contribute profit; it also modifies the cost structure of the remainder of the route [2002.04303]. In the profit-versus-time formulations, this coupling is visible in the denominator of each edge-traversal term, because the effective speed depends on the cumulative weight after previously visited cities [2011.05081][2603.18038].

A recurring source of confusion is the assumption that one may optimize the tour and the packing nearly independently and then combine the two results. The cited formulations directly contradict that view. The best tour under constant speed need not remain best once speed becomes load-dependent, and a packing plan that is attractive in a static knapsack model may become inferior once the carrying cost over remaining distance is accounted for. This is why the literature consistently treats BI-TTP as a multi-component optimization problem rather than as a simple juxtaposition of TSP and KP [1802.02434][2002.04303].

## 3. Exact dynamic programming and indicator-based evolutionary search

For a fixed tour \(\pi\), Wu et al. consider the Packing-While-Travelling (PWT) subproblem: choose the packing \(\rho\) to maximize \(f(\pi,\rho)\) subject to the capacity constraint. Building on Neumann et al. (2017), they use an exact dynamic-programming scheme that not only returns one optimal packing but naturally constructs the entire Pareto front of the fixed-tour subproblem [1802.02434].

With the items ordered along the tour, the dynamic-programming table is
\[
B[k,w] = \max\{\text{total reward for choosing among items }1,\dots,k\text{ with total weight exactly }w\}.
\]
The recurrence compares excluding and including the \(k\)-th item:
\[
B[k,w]=\max\Bigl\{B[k-1,w],\; B[k-1,w-w_k]+p_k-\Delta_k(w)\Bigr\},
\]
where \(\Delta_k(w)\) accounts for the additional slowdown induced by picking item \(k\). After each row is computed, dominated cells are discarded by pruning. The naïve version runs in \(O(mC)\) time and space; with pruning, the practical performance is often much better [1802.02434].

The final non-dominated pairs \((w,B[m,w])\) form the DP front \(\tau_\pi\) for tour \(\pi\). Wu et al. use these fronts inside an indicator-based evolutionary algorithm in which the individual is the tour \(\pi\), not the full pair \((\pi,\rho)\). For every tour in the population, the DP subroutine produces a full Pareto front. The union of all fronts is
\[
\Phi=\bigcup_{\pi\in\Pi}\tau_\pi,
\]
and its Pareto-optimal surface is \(\omega=\mathrm{Pareto}(\Phi)\). Tour quality is then defined by its contribution to that global surface through two indicators: Loss-of-Surface-Contribution (LSC) and Loss-of-Hypervolume (LHV). The general indicator template is
\[
\mathcal I(\tau)=1-\frac{\mathcal M(\omega\setminus \tau)}{\mathcal M(\omega)}.
\]
For surface contribution,
\[
SC(\tau)=\frac{|\omega\cap\tau|}{|\omega|},
\qquad
\mathrm{LSC}(\tau)=1-SC(\omega\setminus \tau),
\]
and for hypervolume,
\[
\mathrm{LHV}(\tau)=1-\frac{HV(\omega\setminus\tau)}{HV(\omega)}.
\]
Survivor selection repeatedly removes the tour with smallest indicator value, while parent selection biases mating toward tours whose DP fronts make larger contributions to the global surface [1802.02434].

This architecture is distinctive because the evolutionary layer perturbs only the tour component, while the packing trade-offs for each route are handled exactly by dynamic programming. A plausible implication is that the method avoids one common source of noise in multi-objective search: the need to assess a route using only one sampled packing configuration rather than its full fixed-route Pareto structure.

## 4. Scalarization methods: weighted sums and \(\varepsilon\)-constraints

Another major line of work treats BI-TTP through scalarization. In the weighted-sum method, a parameter \(\alpha\in[0,1]\) controls the balance between profit and travel time:
\[
F_\alpha(\pi,z)=\alpha\,f_1(\pi,z)-(1-\alpha)\,R\,f_2(\pi,z).
\]
Varying \(\alpha\) over \([0,1]\) yields different trade-offs, with \(\alpha=0\) emphasizing a pure TSP tour and \(\alpha=1\) recovering the KP-only perspective. The method uses randomized versions of existing heuristics: the Chained Lin-Kernighan heuristic for tours, a randomized packing heuristic with item scores
\[
s_j=\frac{p_j^a}{(w_j)^b(d_j)^c},
\]
and exploitation through restricted 2-opt moves and probabilistic bit flips [2011.05081].

The same paper states the classical limitations of weighted-sum decomposition: an inability to find Pareto points in non-convex regions of the front and a tendency to produce unevenly spaced solutions when the front is curved. Because the BI-TTP subproblems are solved heuristically rather than exactly, the algorithm may still stumble upon some non-convex front points by virtue of heuristic imperfections [2011.05081]. This is an explicit methodological caveat rather than a contradiction.

A more recent scalarization strategy reformulates BI-TTP by the \(\varepsilon\)-constraint method. First, the extreme profits \(g_{\min}\) and \(g_{\max}\) are determined, and the interval \([g_{\min},g_{\max}]\) is divided into \(S\) segments
\[
\varepsilon_0=g_{\min}<\varepsilon_1<\dots<\varepsilon_S=g_{\max}.
\]
For each band \(s\), the solver minimizes travel time subject to a profit interval:
\[
\min_{\mathcal A,\mathcal Z} f_1(\mathcal A,\mathcal Z)
\quad\text{s.t.}\quad
\varepsilon_s\le f_2(\mathcal A,\mathcal Z)\le \varepsilon_{s+1},\;
W_N(\mathcal A,\mathcal Z)\le W.
\]
The resulting fractional objective is rewritten with auxiliary variables \(b_i>0\), transformed into a Quadratic Unconstrained Binary Optimization model, and solved with quantum annealing, after which a tailored heuristic refinement called LEA (Later and Enough Accumulation) is applied [2603.18038].

The main algorithmic families in the cited literature can be organized as follows.

| Approach | Core mechanism | Reported emphasis |
|---|---|---|
| Indicator-based EA + DP | Exact PWT DP front inside a bi-objective EA | Surface contribution and hypervolume guidance |
| Weighted-sum method | Randomized scalarization over \(\alpha\) | Broad heuristic Pareto approximation |
| NDS-BRKGA | Random-key GA with non-dominated sorting | Elite survival and customized repair |
| \(\varepsilon\)-constraint + QA | Profit bands, QUBO reformulation, QA, LEA | Broad Pareto front and time efficiency |

## 5. Population-based genetic algorithms and customized representations

The non-dominated sorting based customized random-key genetic algorithm (NDS-BRKGA) represents BI-TTP by a real-valued genotype \(x\in[0,1]^{(n-1+m)}\), partitioned into tour genes and packing genes. Decoding sorts the tour genes to obtain a permutation of cities \(2,\dots,n\), prepends city \(1\), and converts each packing gene to a binary decision through a threshold at \(0.5\) [2002.04303].

Initialization incorporates domain knowledge from both subproblems. The algorithm solves the TSP subproblem with LKH for \(t\) seconds to obtain a near-optimal tour \(\pi^\*\) and its symmetric reverse, and solves the KP subproblem by a two-stage Greedy+DP heuristic, controlled by a parameter \(\delta\), to obtain a packing plan \(z^\*\). Candidate \((\pi,z)\) solutions are built by combining these structures, filtered by non-domination, encoded back into random keys, and mixed with uniformly random individuals [2002.04303]. The paper explicitly argues that domain-knowledge seeding accelerates convergence.

Capacity feasibility is enforced by a custom repair operator. If a decoded packing is overweight, the operator scans the tour from end to start and removes items whose collection contributes to the overweight, preferentially dropping items collected latest and thereby maximizing early-tour speed. After repair, the corresponding keys are reset to reflect the feasible packing [2002.04303].

The evolutionary core follows BRKGA-style biased crossover:
\[
c_i=
\begin{cases}
a_i & \text{w.p. }\rho_e,\\
b_i & \text{w.p. }1-\rho_e,
\end{cases}
\]
where \(a\) is an elite parent and \(b\) is drawn from the full population. Diversity is maintained through mutant individuals sampled uniformly in the genotype space. Bi-objective survival uses the NSGA-II mechanism: fast non-dominated sorting, crowding distance in normalized 2-D objective space, and truncation by least-crowded individuals when a front must be cut [2002.04303].

The parameter study spans 3072 parameter combinations. The reported best hypervolumes arise with \(N=500\ldots1000\), \(N_e\approx0.4\ldots0.6N\), small \(N_m\), \(\rho_e\approx0.6\ldots0.8\), \(\alpha\ge 0.1\), and regular local search with \(\omega\approx10\ldots50\) [2002.04303]. In competition results, the preliminary “jomar” entry won 1st place on seven of nine instances at EMO-2019, and the refined NDS-BRKGA placed 2nd overall at GECCO-2019 [2002.04303].

## 6. Benchmarks, empirical results, and open directions

The empirical literature uses several benchmark sets. Wu et al. evaluate on the standard TTP benchmark families eil51, eil76, and eil101, each with three item-correlation types: uncorrelated, strongly bounded, and similar weights. Their parameter setting is population size \(\mu=50\), maximum generations \(\alpha=20\,000\), five independent runs per setting, and 16 configurations arising from two indicators \(\{\mathrm{LHV},\mathrm{LSC}\}\) crossed with eight parent-selection rules. Performance is measured by the global hypervolume of the final union-surface and the total reward \(f(\pi,\rho)\) of the best single solution, compared with the single-objective state-of-the-art MA2B. Across all 16 settings, harmonic rank-based selection on LHV and best-half on LSC performed best, fitness-proportionate selection was a safe default, the BI-TTP approach consistently found larger hypervolumes than uniform-random selection with \(p\ll .001\) under Welch’s \(t\)-test, and the best BI-TTP variants outperformed MA2B on 8 of 9 test instances [1802.02434].

The weighted-sum study evaluates nine medium and large BI-TTP instances from the EMO-2019 and GECCO-2019 competitions together with 960 single-objective TTP instances for parameter tuning and validation. Its metrics are hypervolume of the Pareto set, number of non-dominated points, single-objective TTP score, runtime, and the curve of hypervolume over time. Against NDSBRKGA on 960 instances, WSM achieved a significantly higher average hypervolume in 82.2% of cases, was worse in 15.0%, and tied in 2.8% under a Wilcoxon test with \(p<0.05\). On the six \(fnl4461_*\) and \(pla33810_*\) instances, WSM attained the top hypervolume, improving the best known by up to 1.7%, and it matched or exceeded those of 21 leading TTP algorithms on 379 of 960 benchmark instances in the single-objective setting [2011.05081].

The NDS-BRKGA study uses nine benchmark instances from Polyakovskiy et al. (2014), with sizes ranging from \(a280\) to \(pla33810\), and evaluates performance by normalized two-dimensional hypervolume. It reports that LKH needs only about 300 seconds to reach at most a 0.15% gap on the TSP subproblem, and that the GH+DP knapsack solver with \(\delta=5\cdot 10^4\) finds optimal KP solutions in acceptable time. Even though BI-TTP does not explicitly maximize the single-objective TTP score, many of the reported bi-objective solutions strictly dominate the best known single-objective solutions on the two smallest instances [2002.04303].

The \(\varepsilon\)-constraint plus quantum-annealing study uses benchmark instances from Polyakovskiy et al. (2014): \(ch150\) variants with \(N=150\) and \(a280\) variants with \(N=280\). It compares against NSGA-II, U-NSGA-III, and MOEA/D, using normalized hypervolume and total runtime, including classical CPU time, QPU access, and LEA. The reported results are hypervolume values around \(0.7\) to \(0.9\) across all instances in a single hybrid run, approximately \(1.3\) to \(1.4\) seconds of QPU time for 10 subproblems with less than \(0.5\%\) LEA overhead, overall runtimes such as about \(724\) seconds for \(ch150\_n149\_bsc\) and about \(4\,144\) seconds for \(a280\_n279\_bsc\), and up to \(9\times\) speed-up relative to the cited baselines while obtaining strictly better or comparable hypervolume [2603.18038].

Across these studies, several limitations and open directions recur. The weighted-sum method inherits the classical difficulty of covering strictly non-convex front regions under exact solution, though heuristic search can mitigate that effect [2011.05081]. The DP-based hybrid may require more aggressive pruning or approximation when \(C\) is very large [1802.02434]. The random-key GA notes the computational expense of repeated LKH runs and local search on very large instances, and possible bias from initialization by single-objective subproblem solvers [2002.04303]. The quantum-annealing formulation is motivated by severe scalability issues in conventional methods and by the complexity of reformulating fractional objectives into QUBO-compatible form [2603.18038]. Future directions explicitly proposed in the literature include richer tour-encoding operators such as edge-assembly crossover, multi-level DP/EAs, per-instance algorithm configuration, alternative scalarization techniques including \(\varepsilon\)-constraint and Tchebycheff, and extensions to dynamic, multi-depot, and other interwoven multi-component optimization problems [1802.02434][2002.04303][2011.05081].

In aggregate, the BI-TTP literature portrays the problem as a stringent testbed for hybrid optimization: exact dynamic programming embedded in evolutionary selection, randomized heuristic scalarization, customized non-dominated sorting genetic algorithms, and quantum-classical \(\varepsilon\)-constraint pipelines all exploit the same route-packing interdependence but expose different trade-offs in exactness, diversity, scalability, and runtime.

Source: https://www.emergentmind.com/topics/bi-objective-traveling-thief-problem-bi-ttp