N = int(input()) from collections import defaultdict # 使用defaultdict来存储立方和对应的三元组列表 cube_sums = defaultdict(list) # 生成所有可能的b, c, d，满足2 <= b <= c <= d <= N，并记录它们的立方和 for b in range(2, N + 1): for c in range(b, N + 1): for d in range(c, N + 1): total = b ** 3 + c ** 3 + d ** 3 cube_sums[total].append((b, c, d)) # 收集所有符合条件的四元组 results = [] for a in range(2, N + 1): a_cubed = a ** 3 if a_cubed in cube_sums: for triple in cube_sums[a_cubed]: results.append((a,) + triple) # 按a从小到大输出，每个四元组中的b,c,d已满足顺序 for quad in results: print(f"{quad[0]}^3 = {quad[1]}^3 + {quad[2]}^3 + {quad[3]}^3")