Improve the performance of removing many variables from a C++ quadratic model

Right now, to remove multiple variables from a QM, one uses a loop like
```C++
for (auto& v : to_remove) {
expression.remove_variable(v);
}
```
it would nice to be able to do so in bulk, and we can save some performance in the reindexing.