def maybeRightX(maybeX, polys):
    sum = 0
for i in range(len(polys)):
    sum += polys[i] * maybeX ** i
return sum

def solve(y, polys):
    lo = 0
hi = y
while lo <= hi:
    mid = (lo + hi) // 2
SymPy cannot solve polynomial equation of 4th order