Добавлю маленько. Если знать, как решается уравнение с двумя переменными, то увеличение числа переменных принципиальных трудностей не создаёт. Подробности можно поискать в известных учебниках по теории чисел. Вычленить из Бухштаба вычислительный алгоритм будет для начинающего несколько трудноват. Также не искал есть ли в инете нормальные сайты с примерами. Изложу одну из возможных схем, ограничиваясь этим примером.
Итак, требуется найти все целые решения уравнения
Временно считаем

ио параметра, то есть рассматриваем уравнение

считая параметр

известным.
Так как левая часть делится на 3, то это уравнение явно неразрешимо, если

не делится на 3, иначе же

, то есть

.
Подставив в уравнение и сократив на 3, получим

Это уже стандартное уравнение общим решением:

Как итог переменные

выражены через произвольные целочисленные параметры

и

.
В общем случае, линейное уравнение с

неизвестными будет разрешимо только если правая часть делится на НОД коэффициентов левой части и общее решение в таком случае будет содержать

произвольных целочисленных параметров. Схема та же, как показано - одно из переменных ио параметра (чтобы заслужить такую честь, надо иметь соответствующий вид - что-то на что-то должно делиться), получаем уравнение с параметром от меньшего числа неизвестных, опять уменьшаем число неизвестных и т.д.