This paper develops an algorithmic solution using Ising machines to solve large-scale higher-order binary optimization (HOBO) problems with inequality constraints for resource optimization in wireless communications systems. Quadratic unconstrained binary optimization (QUBO) aims to solve a special category of these problems widely encountered in engineering and science. To solve QUBO instances, specialized Ising machines have been designed, while sophisticated quantum annealing algorithm and quantum-inspired classical heuristics have been developed. However, the application of QUBO in wireless communications has limited practical interest mainly due to the complexity of resource optimization problems which are often characterized by high-order polynomial terms and strict inequality constraints. To overcome these bottlenecks and take advantage of recent advancements in Ising machines, in this paper, we propose an iterative algorithmic solution to solve HOBO problems, which is based on the augmented Lagrangian method to handle constraints. Specifically, Taylor expansion is employed to approximate higher-order polynomials to quadratic ones in the augmented Lagrangian function, which enables the solution of a single QUBO problem at each iteration without auxiliary variables. As an illustrative case study, we consider the problem of phase optimization in a simultaneous wireless information and power transfer system, where a reconfigurable intelligent surface with 1-bit phase resolution is used to facilitate information/energy transfer. Simulation results verify that the proposed algorithm achieves satisfactory performance and outperforms heuristic benchmark schemes.
This paper develops an algorithmic solution using Ising machines to solve large-scale higher-order binary optimization (HOBO) problems with inequality constraints for resource optimization in wireless communications systems. Quadratic unconstrained binary optimization (QUBO) aims to solve a special category of these problems widely encountered in engineering and science. To solve QUBO instances, specialized Ising machines have been designed, while sophisticated quantum annealing algorithm and quantum-inspired classical heuristics have been developed. However, the application of QUBO in wireless communications has limited practical interest mainly due to the complexity of resource optimization problems which are often characterized by high-order polynomial terms and strict inequality constraints. To overcome these bottlenecks and take advantage of recent advancements in Ising machines, in this paper, we propose an iterative algorithmic solution to solve HOBO problems, which is based on the augmented Lagrangian method to handle constraints. Specifically, Taylor expansion is employed to approximate higher-order polynomials to quadratic ones in the augmented Lagrangian function, which enables the solution of a single QUBO problem at each iteration without auxiliary variables. As an illustrative case study, we consider the problem of phase optimization in a simultaneous wireless information and power transfer system, where a reconfigurable intelligent surface with 1-bit phase resolution is used to facilitate information/energy transfer. Simulation results verify that the proposed algorithm achieves satisfactory performance and outperforms heuristic benchmark schemes.
Quantum annealing (QA) offers a physics-inspired approach for such subproblems by mapping them onto an Ising Hamiltonian, which is algebraically equivalent to a Quadratic Unconstrained Binary Optimization (QUBO) model executable on quantum annealers or hybrid solvers. However, their practical employment in wireless systems critically depends on the specific problem formulation, on the particular variable encoding, constraint handling, solver selection, hardware embedding, and fair benchmarking. This paper develops wireless-oriented QUBO modeling guidelines for QA-assisted resource allocation. We first present a beginner-friendly toy example in which a terrestrial-aerial resource-allocation problem with a fractional signal-to-interference-plus-noise ratio (SINR) maximization objective is reformulated into a two-bit QUBO via Boolean interpolation, converted into an Ising representation, and decoded back into a wireless scheduling decision. The same scheduling is applicable to UAV, aircraft, satellite, and terrestrial links in SAGIN systems. We then review the QA foundations and solver ecosystems, including binary quadratic, discrete quadratic, constrained quadratic, and hybrid models from the perspective of NG wireless optimization. Furthermore, we discuss both binary and discrete encoding, penalty-based constraint conversion, embedding overhead, solution decoding, and feasibility verification. A QUBO Compatibility Index is introduced as a lightweight diagnostic metric for assessing the associated QUBO formulation readiness, rather than a predictor of runtime or solution quality. Finally, a case study on max-min uplink power control in cell-free massive MIMO illustrates how a fairness-oriented wireless problem can be discretized, constrained, converted into the concept of a QUBO, and evaluated against classical baselines. A host of key challenges and future directions are also discussed.
This study proposes a novel method for simplifying inequality constraints in Higher-Order Binary Optimization (HOBO) formulations. The proposed method addresses challenges associated with Quadratic Unconstrained Binary Optimization (QUBO) formulations, specifically the increased computational complexity and reduced solution accuracy caused by the introduction of slack variables and the resulting growth in auxiliary qubits. By efficiently integrating constraints, the method enhances the computational efficiency and accuracy of both quantum and classical solvers. The effectiveness of the proposed approach is demonstrated through numerical experiments applied to combinatorial optimization problems. The results indicate that this method expands the applicability of quantum algorithms to high-dimensional problems and improves the practicality of classical optimization solvers for optimization problems involving inequality constraints.
Ising machines, including quantum annealing machines, are promising next-generation computers for combinatorial optimization problems. However, due to hardware limitations, most Ising-type hardware can only solve objective functions expressed in linear or quadratic terms of binary variables. Therefore, problems with higher-order terms require an order-reduction process, which increases the number of variables and constraints and may degrade solution quality. In this study, we evaluate the effectiveness of directly solving such problems without order reduction by using a high-performance simulated annealing-based optimization solver capable of handling polynomial unconstrained binary optimization (PUBO) formulations. We compare its performance against a conventional quadratic unconstrained binary optimization (QUBO) solver on the same hardware platform. As benchmarks, we use the low autocorrelation binary sequence (LABS) problem and the vehicle routing problem with distance balancing, 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 potential advantages of directly handling higher-order terms in practical optimization problems.
Recent hardware advances in quantum and quantum-inspired annealers promise substantial speedup for solving NP-hard combinatorial optimization problems compared to general-purpose computers. These special-purpose hardware are built for solving hard instances of Quadratic Unconstrained Binary Optimization (QUBO) problems. In terms of number of variables and precision of these hardware are usually resource-constrained and they work either in Ising space or in Boolean space. Many naturally occurring problem instances are higher-order in nature. The known method to reduce the degree of a higher-order optimization problem uses Rosenberg's polynomial. The method works in Boolean space by reducing the degree of one term by introducing one extra variable. In this work, we prove that in Ising space the degree reduction of one term requires the introduction of two variables. Our proposed method of degree reduction works directly in Ising space, as opposed to converting an Ising polynomial to Boolean space and applying previously known Rosenberg's polynomial. For sparse higher-order Ising problems, this results in a more compact representation of the resultant QUBO problem, which is crucial for utilizing resource-constrained QUBO solvers.
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.
An Ising machine (IM), as a type of analog computer tailored for tackling intractable combinatorial optimization problems, has attracted remarkable attention in recent years. In contrast to the blossoming field of bespoke IM hardware, developing metaheuristics from IMs remains largely uninvestigated. Here, we propose a physics-inspired evolutionary computation paradigm, termed the Ising optimization framework (IOF); it comprises a unique Ising algorithm and a hybrid annealing scheme, which together are well-suited for solving quadratic unconstrained binary optimization (QUBO) problems embedded in Ising system energy. The Ising algorithm leverages a set of iterated self-mapping functions to evolve an Ising-spin swarm, enabling efficient energy minimization in artificial Ising systems while mitigating detrimental chaos. Complementing the algorithm, a hybrid annealing scheme integrating singular value dropout, bifurcation control, and a nudging strategy, is devised to augment the overall optimization capacity. The effectiveness of the IOF is validated on various Ising and Max-cut problems with decision variables ranging from 625 to 5000 in number. In comparison to four major types of methods for solving QUBOs, including IM simulations, nature-inspired algorithms, a state-of-the-art heuristic, and the commercial solver Gurobi, the IOF consistently demonstrates notable optimization quality and computational efficiency. This paper provides a theoretical foundation and practical guidelines for bridging Ising-inspired approaches with evolutionary computation, offering an evolutionary perspective on Ising optimizations and suggesting a fertile avenue for future research and application.