Generalization of Lieb's variational principle to Bogoliubov-Hartree-Fock theory