import numpy
a=input("k=")
k=int(a)
print ("圆周率的值为")
i=0
S=0
for i in range (k):
a=numpy.math.factorial(4*i)
b=numpy.math.factorial(i)
S+=a*(26390*i+1103)/(pow(b,4)*pow(396,4*i))
A=1/(2*(2**(1/2))/(99**2)*S)
print (A)
import random
N=input("投掷次数为:")
n=int(N)
i=0
k=0
for i in range(n):
x=random.uniform(0,1)
y=random.uniform(0,1)
if x**2+y**2<1:#在四分之一圆内的点
k=k+1
print ("pi=",4*k/n)