from tools import number_theory _eps = 0.0000000000001 def n_peak(p): for n in range(1, p): right = n ** 3 + p * (n ** 2) m = int(right ** (1 / 3) + _eps) print(n, m, right - m ** 3) print(n_peak(19))