from sympy.ntheory.modular import crt N = [86812553978993, 81744303091421, 83695120256591] c = [8875674977048, 70744354709710, 29146719498409] x = crt(N, c)[0] print("Hex String:") print(round(x ** (1. /3)))