---
title: Cardinality-Constrained Portfolio Selection
url: https://www.emergentmind.com/topics/cardinality-constrained-portfolio-selection
type: topic
---

# Cardinality-Constrained Portfolio Selection

Searching arXiv for recent and foundational papers on cardinality-constrained portfolio selection to ground the article in current research.
Cardinality-constrained portfolio selection is the class of portfolio optimization problems in which the number of active positions is explicitly limited, either globally or within groups such as sectors, while the portfolio still satisfies risk, return, budget, and often trading constraints. In the classical long-only mean–variance setting, a representative formulation is
\[
\min_{x \in \mathbb{R}^n}\quad x^\top A x - \tau \mu^\top x
\quad \text{subject to}\quad
e^\top x = 1,\;
x \ge 0,\;
\|x\|_0 \le k,
\]
where \(A \succeq 0\) is the covariance matrix, \(\mu\) the expected return vector, and \(k\) the maximum number of held assets [2309.16004]. More general formulations incorporate sector-level sparsity, transaction costs, benchmark-relative risk, tracking objectives, dynamic factor structure, and even higher-order interaction terms [1810.10563], [1404.3286], [2503.18609], [2502.17915], [2603.15947]. The central technical feature is that the \(\ell_0\)-type restriction turns an otherwise convex portfolio problem into a nonconvex, combinatorial, and typically NP-hard optimization problem [1810.10563], [2309.16004].

## 1. Problem class and mathematical structure

At its core, the topic concerns sparse portfolio design under explicit support constraints. In the simplest Markowitz-style case, sparsity is imposed by \(\|x\|_0 \le k\), with full investment and nonnegativity ensuring a long-only portfolio on the simplex [2309.16004]. In more general settings, the weights are partitioned into blocks \(w=[w_1,\dots,w_m]\), and the feasible set may require both sector-weight bounds and groupwise cardinality limits:
\[
\Omega = \Big\{p_i \le \mathbf{1}^\top w_i \le q_i,\; \|w_i\|_0 \le k_i,\; i=1,\dots,m\Big\},
\]
so that sparsity is controlled simultaneously at group and total-portfolio levels [1810.10563].

A practically important extension is single-period rebalancing with transaction costs and benchmark-relative risk. In that setting, the investor starts from a current portfolio \(P\), trades through buy and sell vectors \(x^b,x^s\), and obtains a new portfolio \(x\) satisfying
\[
P + x^b - x^s = x,
\]
while minimizing benchmark-relative variance
\[
(x-x^*)^\top Q(x-x^*)
\]
subject to a net expected return constraint after transaction costs and a cardinality requirement \(\sum_j z_j=\text{card}\), with linking bounds \(a_j z_j \le x_j \le b_j z_j\) [1404.3286]. In index tracking, the same sparsity principle appears in a different objective: the portfolio is chosen to minimize tracking error relative to an index, optionally with an additive enhancement target, while holding only \(n\) stocks from a much larger index universe [2503.18609].

This suggests that “cardinality-constrained portfolio selection” is not a single model but a family of sparse allocation problems whose common element is exact control over the support of the portfolio. The risk functional may be variance, tracking-error variance, CVaR, or a higher-order objective, but the nonconvexity arises from the same discrete support selection mechanism [1810.10563], [2603.15947].

## 2. Mixed-integer formulations and combinatorial hardness

The standard exact modeling device is to introduce binary selection variables. In the mean–variance case, the support restriction can be written through \(z_i \in \{0,1\}\) with
\[
x_i \le M z_i,\qquad \sum_i z_i \le k,
\]
or, more tightly, through lower and upper bound linking constraints \(a_j z_j \le x_j \le b_j z_j\) [2309.16004], [1404.3286]. In benchmark-relative and rebalancing models, the resulting problem is a mixed-integer quadratic program with a quadratic objective, linear constraints, and binary variables [1404.3286]. In the minimum-variance setting, the classical Big-\(M\) formulation is
\[
\begin{aligned}
\min_{\beta,s}\quad & \beta^\top \Sigma \beta\\
\text{s.t.}\quad & \mathbf{1}^\top \beta = 1,\;
\mathbf{1}^\top s \le k,\;
-Ms_j \le \beta_j \le Ms_j,\;
s_j \in \{0,1\},
\end{aligned}
\]
which explicitly couples support selection and continuous weights [2505.10099].

The source of hardness is the support search. Even the simplest constraint \(\|w\|_0 \le k\) converts a convex problem on the simplex into subset selection over a nonconvex union of coordinate subspaces [1810.10563]. The combinatorial scale is immediate: choosing 10 out of 30 assets already yields \(\binom{30}{10} \approx 3\times 10^8\) subsets [1810.10563]. In a CAPM-based industry universe with \(n \approx 94\) and \(K=10\), the support count \(\binom{94}{10}\) is described as astronomically large, rendering exhaustive search infeasible [2603.15652].

This computational barrier is amplified by realistic constraints. Transaction costs, lower and upper bounds, groupwise sparsity, tracking-error restrictions, and dynamic state dependence enlarge the continuous part of the problem while leaving the discrete search intact [1404.3286], [2502.17915]. The exact MIQP route therefore remains conceptually clean but can become impractical as \(n\), \(k\), or the number of auxiliary variables grows [2309.16004], [2505.10099].

A common misconception is that sparsity can be enforced adequately by an \(\ell_1\) constraint. Under simplex constraints, however, \(\|w\|_1=1\) already holds, so an \(\ell_1\) restriction is ineffective as a cardinality surrogate in long-only portfolio models [1810.10563].

## 3. Continuous reformulations and exact-penalty approaches

A major line of work replaces the discrete support decision by a continuous reformulation while preserving exact cardinality semantics. One route is difference-of-convex programming. For the rebalancing model with transaction costs and binary variables \(z\), the integrality condition is encoded by the concave penalty
\[
\alpha(z)=\sum_{j=1}^n z_j(1-z_j),
\]
which is nonnegative on the relaxed feasible set and vanishes if and only if \(z_j \in \{0,1\}\) for all \(j\) [1404.3286]. With a sufficiently large exact-penalty parameter \(\theta\), the original MIQP is equivalent to the continuous DC program
\[
\min_{(x,x^b,x^s,z)\in A}
\Big[(x-x^*)^\top Q(x-x^*) + \theta \sum_{j=1}^n z_j(1-z_j)\Big],
\]
where \(A\) is the convex relaxation of the original feasible set [1404.3286]. The associated DC Algorithm linearizes the convex part \(h(z)=\theta\sum_j z_j(z_j-1)\) and solves a sequence of convex quadratic programs.

Another route introduces an auxiliary sparse variable. In the relaxed framework for mean–variance and CVaR models, one keeps a continuous portfolio variable \(w\in\Delta_1\), introduces \(v\in\Omega\) carrying the exact \(\ell_0\) restrictions, and couples the two through a quadratic penalty:
\[
\min_{w,v}\quad f(w)+\frac{\nu}{2}\|w-v\|^2
\quad\text{s.t.}\quad
w\in\Delta_1,\; v\in\Omega.
\]
Here the nonconvexity is confined to \(v\), while \(w\) remains in a convex simplex [1810.10563]. Projection onto \(\Omega\) is tractable because within each sector the Euclidean projection reduces to keeping the \(k_i\) largest components and projecting the retained block onto a simplex or interval-simplex set [1810.10563].

A third reformulation uses selector variables in \([0,1]\) and a bilinear complementarity-type constraint. For a generic convex cardinality-constrained problem with nonnegative variables, the exact equivalence
\[
\|x\|_0 \le K
\iff
\exists y\in[0,1]^n:
\sum_i y_i \le K,\;
\sum_i x_i - \sum_i y_i x_i \le 0
\]
yields a continuous reformulation in \((x,y)\) [1906.00504]. Linearizing the bilinear term \(\sum_i y_i x_i\) produces a successive convex approximation scheme in which each step solves a convex subproblem; under the paper’s assumptions, KKT points of the reformulation with \(\|x^*\|_0=K\) are local optimizers of the original problem [1906.00504].

These constructions share a common principle: preserve the exact support logic, but move the hard combinatorics into penalties, projections, or complementarity structures that admit continuous optimization substeps.

## 4. Decomposition, projections, and scalable first-order methods

A different research direction seeks scalability by separating support choice from weight optimization or by working with continuous relaxations whose geometry encourages binary solutions.

Penalty decomposition is one example. For the cardinality-constrained mean–variance problem, the split formulation
\[
\min_{x,y}\quad x^\top A x - \tau \mu^\top x
\quad\text{s.t.}\quad
e^\top x=1,\;
y\ge 0,\;
\|y\|_0\le k,\;
x-y=0
\]
is penalized with an \(\ell_2\) coupling term
\[
q_\rho(x,y)=x^\top A x-\tau\mu^\top x+\rho\|x-y\|_2^2.
\]
For fixed \(\rho\), block coordinate descent has closed-form updates: the \(x\)-step is an equality-constrained convex quadratic with analytic solution, and the \(y\)-step is the hard-thresholding projection \(y_*=H_k(x^+)\) [2309.16004]. The outer penalty sequence drives \(x\) and \(y\) together, and accumulation points are proved to be local minimizers of the original mean–variance problem under the paper’s conditions [2309.16004].

Column-generation-type methods address cardinality indirectly by solving repeated restricted master problems over a candidate set of assets. In the benchmark-relative Markowitz extension with active-weight, sector, market-cap quintile, beta, active-share, and tracking-error considerations, the algorithm repeatedly solves a convex quadratic problem on a restricted set \(S\) of candidate assets, then uses dual information to compute a marginal effect
\[
\delta_i = m_i - k_i
\]
for assets outside \(S\), where \(m_i\) is the direct contribution to the objective and \(k_i\) is the dual-cost term from the constraints [1812.00093]. Assets with the most negative \(\delta_i\) are added, low-weight or zero-weight names are removed, and the candidate-set size enforces the portfolio cardinality heuristically [1812.00093].

Recent scalable first-order work reformulates sparse minimum-variance selection through a Boolean relaxation on the simplex. Defining
\[
\Sigma_t = T_t \Sigma T_t + \delta(I-T_t^2),
\qquad
f_\delta(t) = -t^\top \Sigma_t^{-1} t,
\]
with \(t \in [0,1]^p\) and \(T_t=\operatorname{diag}(t)\), one obtains a continuous objective that agrees exactly with the binary subset objective on \(\{0,1\}^p\) [2505.10099]. The parameter \(\delta\) controls a convex-to-concave transition: for small \(\delta\), \(f_\delta\) is strictly convex on a truncated domain; for sufficiently large \(\delta\), it becomes strictly concave [2505.10099]. A Frank–Wolfe continuation scheme then starts from the simplex center, follows a geometric grid in \(\delta\), and is driven toward binary extreme points corresponding to sparse portfolios.

This suggests a methodological continuum. At one end are exact MIQP formulations; at the other are purely continuous paths that use geometry, penalties, or projections to recover discrete supports without branching. The relative appeal depends on the scale of \(p\), the desired exactness, and the structure of the risk model.

## 5. Heuristics, metaheuristics, and repair-based methods

Metaheuristics remain prominent because they can directly explore the combinatorial support space without exact branching. In benchmark studies of mean–variance portfolios with exact cardinality \(K\) and buy-in bounds \(\ell_i z_i \le w_i \le d_i z_i\), Asexual Reproduction Optimization encodes a portfolio by a chromosome containing \(K\) asset indices and \(K\) corresponding weights [2101.03312]. Cardinality is built into the representation, while budget and bound feasibility are restored by a repair procedure derived from Chang et al. [2101.03312]. On five standard benchmark datasets with 31, 85, 89, 98, and 225 assets, the reported mean percentage error relative to the unconstrained efficient frontier is lower for ARO than for GA, SA, TS, and PSO on four of the five datasets, and the average error is reduced by approximately 20 percent relative to the minimum average error of those alternatives [2101.03312].

In index tracking, another heuristic strategy avoids the NP-hard joint problem by separating asset pre-selection from weight estimation. Eight pre-selection procedures are studied, combining forward selection or backward elimination with OLS or LAD regression, with or without an intercept [2503.18609]. The preferred variant is BE-OLS(n), and for S&P 500 tracking the paper reports that out-of-sample tracking errors are roughly proportional to \(1/\sqrt{\text{cardinality}}\) [2503.18609]. For enhanced index tracking, by contrast, cardinalities of the order 10 to 20 are reported as most effective [2503.18609].

Metaheuristics also require repair operators. A covariance-aware repair framework for long-only portfolios with \(\|w\|_0 \le K\), simplex, and box constraints introduces the distance
\[
d_\Omega(a,b)=\sqrt{(a-b)^\top \Omega (a-b)},
\]
whose square equals the tracking-error variance between portfolios \(a\) and \(b\) [2512.19986]. The proposed Covariance-Aware Simplex Projection first selects the active set using volatility-normalized scores \(|z_i|/\sqrt{\Omega_{ii}}\), then projects onto the selected simplex by minimizing \((w-z)^\top\Omega(w-z)\) rather than Euclidean distance [2512.19986]. On S&P 500 data from 2020–2024, CASP-Basic is reported to deliver materially lower portfolio variance than standard Euclidean repair, with most of the variance reduction driven by volatility-normalized selection and an additional consistent improvement from the covariance-aware projection [2512.19986].

A separate controversy concerns higher-order portfolio objectives. In native cubic cardinality-constrained portfolio optimization, where the objective contains quadratic Markowitz terms plus cubic three-way sector co-movement penalties, SA and tabu search are applied only after Rosenberg quadratization expands the variable count from \(n\) to \(5n\) and adds penalty terms for auxiliary consistency and cardinality [2603.15947]. Under matched 60-second CPU budgets, Hyper-Adaptive Momentum Dynamics operates directly on the native cubic objective with exact cardinality-preserving projection and iterated local search, and the paper reports substantially lower decoded native cubic objective values than SA and tabu search at \(n=200,300,500,1000\) [2603.15947]. The paper’s decoded-feasibility analysis shows that SA can satisfy all exact cardinality and Rosenberg auxiliary constraints yet still decode to a native objective 80–88% worse than HAMD at \(n=200\), which the authors interpret as a surrogate-distortion effect rather than simple infeasibility [2603.15947].

## 6. Extensions: CVaR, tracking, tangent portfolios, dynamic factors, and global optimization

Cardinality constraints are now used well beyond static variance minimization. In CVaR optimization, sparsity can be handled by introducing both a sparse auxiliary portfolio \(v\) and an auxiliary hinge variable \(u\), leading to a relaxed problem in \((w,u,v)\) solved by PALM with FISTA-style acceleration [1810.10563]. For small instances, the paper reports that the method finds a very good local minimum for the cardinality-constrained CVaR model, while for larger dimensions it yields feasible portfolios nearly as efficient as their unconstrained counterparts [1810.10563].

In tangent portfolio optimization, sparsity is imposed directly on the Sharpe-ratio maximizer:
\[
\max_w \frac{\mu^\top w}{\sqrt{w^\top \Sigma w}}
\quad\text{s.t.}\quad
\operatorname{Card}(w)\le k.
\]
A Cholesky-based heuristic exploits the characterization that Sharpe ratios are ordered by angles in the transformed space \(L_\Sigma^\top w\) [2502.11701]. The surrogate problem is solved by selecting the \(k\) largest absolute components of \(L_\Sigma^\top \hat w\), where \(\hat w\) is the unconstrained tangency portfolio, and then re-optimizing on the reduced universe [2502.11701]. The paper further introduces a covariance-based diagonal-dominance score and reports a strong positive correlation, about 0.8265, between this score and heuristic performance relative to CPLEX [2502.11701].

Dynamic settings generalize the support constraint across time. In a multi-period mean–variance problem under dynamic factor models, cardinality appears as a cone-type constraint
\[
\mathcal{A}_t=
\left\{
\boldsymbol\pi\in\mathbb{R}^N:
\sum_{i=1}^N |\operatorname{sign}(\pi_i)| \le N_{\text{active}}
\right\},
\]
which limits the number of active risky positions each period [2502.17915]. Despite the non-convexity, the paper derives a semi-analytical optimal policy as a piecewise linear feedback rule in wealth, with factor dependence embedded in slope vectors \(\mathbf{k}_t^\pm(\mathbf{s}_t)\) obtained from backward stochastic optimization recursions [2502.17915]. In a Markov regime-switching illustration with no-short-selling and cardinality at most 2 out of 4 assets, the active asset set changes with regime and sign regime, showing that sparse dynamic portfolios can exhibit discrete support switching over time [2502.17915].

At the other extreme, globally optimal sparse quadratic optimization can be pursued by interval branch-and-bound. For the generic convex quadratic problem
\[
\min_{x\in B}\; \frac12 x^\top Q x + q^\top x + c
\quad\text{s.t.}\quad
\|x\|_0 \le k,
\]
with \(Q\succeq 0\) and box constraints \(B\), an interval branch-and-bound method branches only at 0 and encodes each coordinate’s status by flags indicating fixed zero, undecided, or definitely nonzero [2504.04043]. The paper states that the proposed algorithm can also accommodate linear inequality constraints while maintaining global convergence, so a budget constraint can be incorporated by linearization into inequalities [2504.04043]. In regression benchmarks, the method is competitive with GUROBI’s quadratic mixed-integer solver and scales to dimensions up to \(p=2000\) for small \(k\), which suggests applicability to sparse portfolios when exact support control is required [2504.04043].

A complementary theoretical line studies relaxations of quadratically constrained cardinality minimization. For the problem
\[
\min_x C(x)
\quad\text{s.t.}\quad
(x-c)^\top Q(x-c)\le \gamma,
\]
where \(C(x)\) counts nonzero components, continuous relaxations derived from mixed-integer formulations are shown to have arbitrarily poor approximation ratios, whereas optimized diagonal relaxations often give much stronger lower bounds and greatly reduce branch-and-bound complexity [1210.5225]. In portfolio terms, this means that risk ellipsoids close to diagonal, diagonally dominant, or nearly coordinate-aligned structures are especially favorable for diagonal relaxation quality [1210.5225].

## 7. Computational trade-offs, misconceptions, and open directions

Several broad conclusions emerge. First, exact cardinality control matters. Regularization and surrogate penalties may induce sparsity, but they generally do not let the investor prescribe the support size directly [2309.16004]. This is why many recent methods favor exact \(\ell_0\) constraints handled through penalties, projections, split variables, or discrete neighborhoods rather than through tuning-only surrogates [1810.10563], [2505.10099].

Second, the main algorithmic trade-off is between global optimality and scalability. MIQP and global branch-and-bound methods can produce exact or provably optimal solutions, but their cost rises rapidly with problem size and support complexity [2504.04043], [2309.16004]. Continuous methods such as DCA, PALM, SCA, Boolean relaxation, and penalty decomposition typically converge only to stationary points or local minimizers, yet they can solve medium- and large-scale instances in fractions of a second to a few seconds and often match exact solvers closely on benchmark sets [1404.3286], [1810.10563], [1906.00504], [2505.10099].

Third, sparsity interacts strongly with model specification. For pure index tracking, larger cardinality improves out-of-sample tracking error, transaction volume, and return–risk ratios, approaching the index as \(n\) rises [2503.18609]. For enhancement, small cardinalities preserve alpha better [2503.18609]. For tangent portfolios, low correlation and diagonal dominance improve the quality of Cholesky-based support selection [2502.11701]. For cubic models, surrogate quadratization can distort the search landscape severely [2603.15947]. For dynamic factor models, cardinality constraints produce regime-dependent support switching, but the optimal policy retains a tractable piecewise linear feedback form [2502.17915].

A plausible implication is that future progress will depend less on a single dominant algorithmic paradigm than on matching the support-handling mechanism to the surrounding portfolio model. Static mean–variance, CVaR, tracking, dynamic factor allocation, and higher-order co-movement models all impose cardinality in different ways, and the most effective methods exploit that structure directly rather than treating sparsity as a generic afterthought.

Source: https://www.emergentmind.com/topics/cardinality-constrained-portfolio-selection