负二项分布

def solve():
    import matplotlib.pyplot as plt
    import numpy as np
    from scipy import stats
    r, p = ____, ____
    x = np.arange(____)
    pmf = stats.nbinom.pmf(____)
    plt.figure(figsize=(____, ____))
    plt.bar(____)
    plt.xlabel(____)
    plt.ylabel(____)
    plt.title(____)
    plt.grid(____)
    plt.show()