Задача: для целых значений I, J, K, L найдите N первых четверок чисел, чтобы выполнялось равенство: i^3-j^3=(k^2+l^2)^2 - C (СИ)

Формулировка задачи:

Никак не могу решить, потому что не подступиться к задаче. Программу вполне по силам написать, да только не представляю, с чего начинать решение. Искать сами числа вполне допустимо с нуля, потому что дан весь диапазон int. По идее вывод программы должен выглядеть как-то так: N=30 1. I=1 J=2 K=3 L=4 2. I=5 J=6 K=7 L=8 ...

Решение задачи: «Задача: для целых значений I, J, K, L найдите N первых четверок чисел, чтобы выполнялось равенство: i^3-j^3=(k^2+l^2)^2»

Листинг программы

n = 0;
for(I=0; ; I++) {
  for(J=0; J<=I; J++) {
    rk = I*I*I - J*J*J;
    for(K=0; K*K <=rk; K++) {
       for(L=K; L*L <= rk -K*K; L++) {
          // Проверка
          // В случае успеха n++
          // Если L!=K еще раз n++
          // Если n >=N выход из всех циклов

Объяснение кода листинга программы

В данном коде:

Переменная n инициализируется нулем и используется для подсчета количества найденных чететок чисел.
В цикле for перебираются все возможные значения I.
В каждой итерации цикла for перебираются все возможные значения J от 0 до I.
Вычисляется значение выражения i^3-j^3=(k^2+l^2)^2.
В цикле for перебираются все возможные значения K от 0 до sqrt(rk).
Для каждого значения K в цикле for перебираются все возможные значения L от K до sqrt(rk-K*K).
Проверяется условие, если оно выполняется, то увеличивается значение n на единицу.
Если L не равно K, то значение n увеличивается еще на единицу.
Если n больше или равно N, то все циклы выходят.