Генерация PMF этой настраиваемой случайной величины с использованием Python

Question

Пусть (X_1,...X_20)~Multinomial(500,(p_1,...,p_20)). И пусть V=(X_1-25)^2+...+(X_20-25)^2. Как найти PMF для V, используя Python? Любая помощь будет оценена.