While you could just use a tool like quadprog with continuous varibles in the interval [0,1], then round the result at the end, but there is no assurance the result is then optimal. And if you do round the solution at the end, it will probably no longer even satisfy the eqaulity constraints. So you would just be wasting your time.
Therefore GA is probably your main choice, with binary variables and a quadratic objective.
