Get an approximation of π using probabilities of coprime integers

coprime_pi.py

from random import randint
import math

def gcd(x,y):
	while y!= 0:
		(x, y) = (y, x % y)
	return x
	
rolls = 1000
max_int = 10000000000

coprime_count = 0
for i in range(0,rolls):
	a = randint(0, max_int)
	b = randint(0, max_int)
	if gcd(a, b) == 1:
		coprime_count += 1

probability = coprime_count / rolls

pi = math.sqrt(6/probability)

print("π approximate: ", pi)