gammaTrialMethod.py

def gammaTrialMethod(n): 
    gamma = 1
    count = 0
    while n % 2 == 0: 
        count += 1
        n = n // 2
    gamma *= pow(2, ceil(count / 2))
    
    for i in range(3, ceil(math.sqrt(n)) + 1, 2): 
        count = 0
        while n % i== 0:
            count += 1
            n = n // i 
        gamma *= pow(i, ceil(count / 2))
              
    if n > 2:
        gamma *= n
    return gamma