统计模拟综合

def solve():
    import numpy as np
    from scipy import stats
    np.random.seed(____)
    n_sim = ____
    n_sample = ____
    alpha = ____
    rejections = ____
    for _ in range(____):
        sample = np.random.normal(____, ____, n_sample)
        _, p = stats.ttest_1samp(____)
        if p < alpha:
            rejections += ____
    rate = rejections / n_sim
    print(____)