#!/usr/bin/env pypy3 import math M = 10 ** 5 X = 10 ** 18 INV = -1 def calc_d(n): p = n * (n + 2) q, r = divmod(p, 4) if r != 0: return INV s = math.sqrt(q) if s * s != q: return INV else: return q def main(): g = 0 for n in range(X): d = calc_d(n) if d == INV: continue print(d) g += 1 if g >= M: break if __name__ == '__main__': main()