Add `-` and `-=` operator overload to `adiar::bdd`

[SSoelvsten]

The - and -= operator adiar::bdd should be added as an overload for the adiar::bdd_diff function. Remember to add a test or two to ensure it works as intended.