Learning Generation Orders for Masked Discrete Diffusion Models via Variational Inference