Fast and accurate simulation of nuclear quantum effects in ab-initio molecular dynamics by a generalized Langevin equation