It is well-known that solving the exact optimization problem in (REF ) is NP-hard [1]} in general. For instance, a popular exact method called leaps and bounds [2]} is practically useful only for values of \(p\)  smaller than 30, as for example is the case when using the R package leaps [3]} or mplot [4]}.
