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

Question

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.2, 0.7],
        [0.1, 0.2, 0.7]],
       [[0.3, 0.2, 0.5],
        [0.4, 0.1, 0.5]]])

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