Don't include terms with 0 coefficient

This makes things much faster, because build_BDD is faster.
Show both building the constraint and solving in the debug output