有没有一种有效的方法可以并行生成多项式随机变量？

numpy.random具有以下用于生成多项式随机样本的功能。

multinomial(n,p,size)

但是我想知道是否存在一种有效的方法来生成针对不同参数n和p的多项式样本。例如，

n = np.array([[10],[20]])
p = np.array([[0.1,0.2,0.7],[0.4,0.4,0.2]])

甚至对于更高的n和p像这样：

n = np.array([[[10],[20]],[[10],[20]]])
p = np.array([[[0.1,[0.1,0.7]],[[0.3,0.5],0.1,0.5]]])

我知道对于单变量随机变量，我们可以做这种事情，但是不知道如何在python中对多项式做事。