Variational and auxiliary field Monte Carlo for the Hubbard and Hubbard-Holstein models: an accurate finite-size scaling and a "sign problem" solution