# expdist3b
# Fig. 22.7  (version of Fig 22.6) 
# input: 10000

import math, random

n = int(input("n=? "))
sum = 0
for i in range(n): 
  while True:
    x = random.random() 
    y = random.random() 
    if x*x + y*y <= 1:
      break
  while True:
    z = random.random() 
    u = random.random() 
    if z*z + u*u <= 1:
      break
  if random.randint(0,1) == 0:
    z = -z
  if random.randint(0,1) == 0:
    u = -u
  dist = math.sqrt((x-z)*(x-z) + (y-u)*(y-u))
  sum += dist
print(f"exp. dist. v2 == {sum/n:.6f}")