PseudoBooleanOptimization.jl

julia> import Pkg; Pkg.add("PseudoBooleanOptimization")

using PseudoBooleanOptimization
const PBO = PseudoBooleanOptimization

f = PBO.PBF{Symbol,Float64}([
    [:x, :y] => 10.0,
    [:x, :z] => 15.0,
    [:y, :z] => -5.0,
])

g = PBO.PBF{Symbol,Float64}([
    :x => 1.0,
    :y => 2.0,
    :z => 3.0,
])

h = f + g

1.0x + 2.0y + 3.0z + 10.0x*y + 15.0x*z - 5.0y*z