Multivariate Pólya distribution