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

[arcondello]

Right now, to remove multiple variables from a QM, one uses a loop like

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.