This article explores the comparative strengths of Higher-Order Unconstrained Binary Optimization (HUBO) and Quadratic Unconstrained Binary Optimization (QUBO) models in the context of probability optimization using tensor sampling techniques. HUBO can represent interactions beyond pairwise terms without introducing auxiliary variables, while QUBO benefits from a simpler algebraic structure and a wide range of established solvers. We combine a theoretical analysis of formulation expressiveness and reduction overhead with empirical experiments on synthetic polynomial objectives, RSA factorization, and Max-Cut instances. The results show that native HUBO formulations are often preferable when the original problem contains high-order interactions, because reducing such problems to QUBO increases the effective dimension seen by tensor samplers. At the same time, QUBO remains competitive for naturally quadratic problems. These findings provide practical guidance for selecting an optimization framework that balances model fidelity, dimension growth, and solver performance.
This article explores the comparative strengths of Higher-Order Unconstrained Binary Optimization (HUBO) and Quadratic Unconstrained Binary Optimization (QUBO) models in the context of probability optimization using tensor sampling techniques. HUBO can represent interactions beyond pairwise terms without introducing auxiliary variables, while QUBO benefits from a simpler algebraic structure and a wide range of established solvers. We combine a theoretical analysis of formulation expressiveness and reduction overhead with empirical experiments on synthetic polynomial objectives, RSA factorization, and Max-Cut instances. The results show that native HUBO formulations are often preferable when the original problem contains high-order interactions, because reducing such problems to QUBO increases the effective dimension seen by tensor samplers. At the same time, QUBO remains competitive for naturally quadratic problems. These findings provide practical guidance for selecting an optimization framework that balances model fidelity, dimension growth, and solver performance.
The increasing complexity of industrial scheduling and transport routing problems motivates the study of alternative optimization formulations and computational paradigms. In this work, we study how higher-order unconstrained binary optimization (HUBO) formulations of such problems map onto quantum optimization workflows in both noisy and fault-tolerant regimes. We consider three representative logistics and manufacturing use cases and formulate each as a HUBO problem. This captures process intricacies, such as highly correlated assembly-line scheduling rules, which are difficult to express faithfully with the standard quadratic (QUBO) form, while at the same time reducing the number of binary variables required in the quantum mapping, thus lowering qubit demand. We compare the HUBO formulations with corresponding QUBO encodings, highlighting a key trade-off: while HUBO reduces qubit requirements through compact binary encoding, it introduces higher-order interaction terms that increase circuit depth, limiting feasibility on current quantum hardware. The proposed formulations are validated using classical solvers across several problem instances and benchmark small routing problem instances using bias-field digitized counterdiabatic quantum optimization in classical simulation. We complement these results with a resource and scalability analysis, focusing on the capacitated vehicle routing problem as a representative large-scale industrial use case. Our analysis indicates that while HUBO formulations offer advantages in qubit scaling compared to QUBO encodings, their practical implementation is constrained by gate fidelity, coherence, and circuit depth, making hybrid quantum-classical workflows and early fault-tolerant quantum hardware the most plausible settings for their practical use.
Quantum approaches to combinatorial optimization problems (COPs) are often limited by the resource demands of Quadratic Unconstrained Binary Optimization (QUBO) encodings, which enlarge circuits through penalty terms and increase qubit and gate counts. We show that Higher-Order Unconstrained Binary Optimization (HUBO) enables a more resource-efficient formulation. Our method systematically constructs HUBO Hamiltonians and, compared to a QUBO formulation in benchmarks on Gate Assignment (GAP), Maximum k-Colorable Subgraph (MkCS), and Integer Programming (IP) problems, significantly reduces qubit requirements and decreases total CNOT gate counts by at least 89.6% for all tested instances. These results highlight HUBO as a practical alternative for quantum optimization on near-term devices. To promote adoption, we release an open-source Python library that automates HUBO model construction, extends beyond the examples presented in this work, and broadens access to resource-efficient quantum optimization.
We present a quantum feature-selection framework based on a higher-order unconstrained binary optimization (HUBO) formulation that explicitly incorporates multivariate dependencies beyond standard quadratic encodings. In contrast to QUBO-based approaches, the proposed model includes one-, two-, and three-body interaction terms derived from mutual-information measures, enabling the objective function to capture feature relevance, pairwise redundancy, and higher-order statistical structure within a unified energy model. To suppress trivial all-selected solutions, we further include structured linear penalties that promote sparsity while preserving informative variables. The resulting HUBO instances are optimized with digitized counterdiabatic quantum optimization on IonQ Forte and compared against noiseless quantum simulation as well as two classical dimensionality-reduction baselines: SelectKBest based on mutual information and principal component analysis (PCA). We evaluate the proposed workflow on two benchmark classification datasets, namely the Gallstone dataset and the Spambase dataset, and analyze both predictive performance and selected-subset structure. The results show good qualitative agreement between hardware executions and noiseless simulations, supporting the feasibility of implementing higher-order feature-selection Hamiltonians on current trapped-ion processors. In addition, the quantum approach yields competitive classification performance while producing compact and informative feature subsets, highlighting the potential of higher-order quantum optimization for machine-learning preprocessing tasks.
The Probabilistic Optimizer with TEnsor Sampling (PROTES) is a gradient-free optimization method for multidimensional arrays and discretized multivariate functions based on the tensor train (TT) format. This novel tensor-network methodology demonstrates significant promise for complex Quadratic Unconstrained Binary Optimization (QUBO) problems. We establish the first application of PROTES to machine learning model training via QUBO formulations and provide the inaugural comparative analysis against D-Wave quantum annealing. Our results reveal a critical performance dichotomy: For smallscale, structurally complex QUBO instances, PROTES achieves superior solution accuracy with enhanced computational stability compared to quantum annealing, while maintaining competitive efficiency. Conversely, Quantum Annealing retains scalability advantages for larger problem dimensions. These findings delineate precise operational domains for tensor-network versus quantum optimization paradigms, establishing PROTES as the preferred solver for compact, high-complexity QUBO problems. This research provides fundamental insights for strategic solver selection and advances the engineering of hybrid quantumclassical optimization frameworks.
This document details the methodology and steps taken to convert Higher Order Unconstrained Binary Optimization (HUBO) models into Quadratic Unconstrained Binary Optimization (QUBO) models. The focus is primarily on prime factorization problems; a critical and computationally intensive task relevant in various domains including cryptography, optimization, and number theory. The conversion from Higher-Order Binary Optimization (HUBO) to Quadratic Unconstrained Binary Optimization (QUBO) models is crucial for harnessing the capabilities of advanced computing methodologies, particularly quantum computing and DYNEX neuromorphic computing. Quantum computing offers potential exponential speedups for specific problems through its intrinsic parallelism capabilities. Conversely, DYNEX neuromorphic computing enhances efficiency and accelerates the resolution of intricate, pattern-oriented tasks by simulating memristors in GPUs, employing a highly decentralized approach, via Blockchain technology. This transformation enables the exploitation of these cutting-edge computing paradigms to address complex optimization challenges effectively. Through detailed explanations, mathematical formulations, and algorithmic strategies, this document aims to provide a comprehensive guide to understanding and implementing the conversion process from HUBO to QUBO. It underscores the importance of such transformations in making prime factorization computationally feasible on both existing classical computers and emerging computing technologies.
Ising machines, including quantum annealing machines, are promising next-generation computers for combinatorial optimization problems. Due to hardware limitations, most Ising-type hardware can only solve functions expressed in linear or quadratic terms of binary variables. Therefore, problems with higher-order terms require order reduction, which increases the number of variables and constraints and may degrade solution quality. This study evaluates the effectiveness of directly solving such problems without order reduction by using the polynomial unconstrained binary optimization (PUBO) solver implemented in Toshiba's Simulated Quantum-inspired Bifurcation Machine+ (SQBM+). We compare its performance against a conventional quadratic unconstrained binary optimization (QUBO) solver on the same hardware platform. As benchmarks, we use a vehicle routing problem with distance balancing and the low autocorrelation binary sequence (LABS) problem, both of which naturally include higher-order interactions. Results show that the PUBO solver consistently achieves superior solution quality and stability compared to its QUBO counterpart, while maintaining comparable computational time and requiring no order-reduction compilation, indicating the potential advantages of directly handling higher-order terms in practical optimization problems.