The semiclassical Schrodinger equation with multiscale and random potentials often appears when studying electron dynamics in heterogeneous quantum systems. As time evolves, the wavefunction develops high-frequency oscillations in both the physical space and the random space, which poses severe challenges for numerical methods. We propose a multiscale reduced basis method, where we construct multiscale reduced basis functions using an optimization method and the proper orthogonal decomposition method in the physical space and employ the quasi-Monte Carlo method in the random space. Our method is verified to be efficient: the spatial grid size is only proportional to the semiclassical parameter and (under suitable conditions) almost first order convergence rate is achieved in the random space with respect to the sample number.