import random
import math
N_tot = 0
N_hits = 0
centerX = 0.5
centerY = 0.5
# precision = input('input number of decimals for approximation: ')
def inside_circle(randomX,randomY):
if(math.sqrt((centerX-randomX)**2+(centerY-randomY)**2)<=0.5):
return True
else:
return False
def estimate_pi():
N_tot = 0
N_hits = 0
for i in range(0,100000):
if(inside_circle(random.random(),random.random())==True):
N_hits = N_hits+1
N_tot = N_tot+1
else:
N_tot = N_tot+1
result = N_hits/N_tot
return result
print(estimate_pi())
aW1wb3J0IHJhbmRvbQppbXBvcnQgbWF0aAoKTl90b3QgPSAwCk5faGl0cyA9IDAKY2VudGVyWCA9IDAuNQpjZW50ZXJZID0gMC41CgojIHByZWNpc2lvbiA9IGlucHV0KCdpbnB1dCBudW1iZXIgb2YgZGVjaW1hbHMgZm9yIGFwcHJveGltYXRpb246ICcpCgoKZGVmIGluc2lkZV9jaXJjbGUocmFuZG9tWCxyYW5kb21ZKToKICAgIGlmKG1hdGguc3FydCgoY2VudGVyWC1yYW5kb21YKSoqMisoY2VudGVyWS1yYW5kb21ZKSoqMik8PTAuNSk6CiAgICAgICAgcmV0dXJuIFRydWUKICAgIGVsc2U6CiAgICAgICAgcmV0dXJuIEZhbHNlCgpkZWYgZXN0aW1hdGVfcGkoKToKICAgIE5fdG90ID0gMAogICAgTl9oaXRzID0gMAoKICAgIGZvciBpIGluIHJhbmdlKDAsMTAwMDAwKToKICAgICAgICBpZihpbnNpZGVfY2lyY2xlKHJhbmRvbS5yYW5kb20oKSxyYW5kb20ucmFuZG9tKCkpPT1UcnVlKToKICAgICAgICAgICAgTl9oaXRzID0gTl9oaXRzKzEKICAgICAgICAgICAgTl90b3QgPSBOX3RvdCsxCiAgICAgICAgZWxzZToKICAgICAgICAgICAgTl90b3QgPSBOX3RvdCsxCgogICAgcmVzdWx0ID0gTl9oaXRzL05fdG90CiAgICByZXR1cm4gcmVzdWx0CgpwcmludChlc3RpbWF0ZV9waSgpKQ==