from sympy import isprime def unary(n): return "1"*(n-1)+"0" def elias_gamma(n): if n ==1: return "1" k=int(math.log(n,2)) fp=unary(1+k) sp=n-2**(k) nb=k sp=bin(sp)[2:] if len(sp)