使用精度为 10^-3 和置信度为 99% 的标准蒙特卡洛方法计算 π 的计算成本是多少？

Question

我想知道如何解决这个问题，因为我使用了蒙特卡洛方法，但我没有得到正确的结果。

import numpy as np 

N = 10**-3
X = np.random.rand(N)
Y = np.random.rand(N)
pi = 4.0 / N * (X**2 + Y**2 <= 1.0).sum()

print(pi)