Random Batch Ewald and Sum-of-Gaussians Methods for High-Scalable MD simulations