# randCoprime.py 
# Fig. 6.8 

import random
from egcd import gcd

M = 1000000
N = 3000000
count = 0
for i in range(M): 
  a = random.randint(1,N)
  b = random.randint(1,N)
  if gcd(a, b) == 1:
    count += 1

print(f"{count/M:.8f}")