From Sets to Multisets: Provable Variational Inference for Probabilistic Integer Submodular Models