# 10**7ということは1周しかできない、2重ループは不可 # 因数分解できる形に変形 # 相加相乗平均、二乗式への式変形 # 3項なら2項にして考察 # 0 <= x <= y <= z <= Nとして考える # 真ん中全探索y決め打ち # 1つを固定してあとは因数分解で決めるのでどうだ def divisors(n): lower_divisors , upper_divisors = [], [] i = 1 while i*i <= n: if n % i == 0: lower_divisors.append(i) if i != n // i: upper_divisors.append(n//i) i += 1 return lower_divisors + upper_divisors[::-1] def main(): N = int(input()) ans = 0 ans_list = [] for z in range(0, N+1): # 因数分解形 # (X+z)(Y+z) - z**2 = N # (X+z)(Y+z) = N+z**2 divs = divisors(N+z**2) #print('z', z, 'divs', divs) for d1 in divs: d2 = (N+z**2)//d1 if d1-z >= 0 and d2-z >= 0: ans += 1 ans_list.append((d1-z, d2-z, z)) print(len(ans_list)) for X, Y, Z in ans_list: print(X, Y, Z) main()