BQM is unconstrained by definition. I need to submit a knapsack problem with two constraints (equality constraints + inequality constraints) to a quantum annealer (DWaveSampler). How would you set the problem as a BQM to be sampled fully quantum (no hybrid solver)? My understanding is that Pyqubo can convert equality constraints to bqm, but it is not clear to me if it can handle inequality constraints too.
More in general, what would be the best strategy for solving a (large!) constrained optimization problem with quantum annealing?