Scalable estimation and inference for massive linear mixed models with crossed random effects