Quantum effects on virial coefficients: numerical approach using centroids