Спасибо за ответ.
Я нашла монографию Венкова "Элементарная теория чисел", там описана теория тернарных форм, в частности, формула для количества представлений натурального числа в виде суммы квадратов трех целых чисел. Я использовала ее, когда выводила формулу для числа решений уравнения

,

,

при заданном натуральном

. Знаю еще, что тернарными формами в 30-40-х годах прошлого века занимался Линник, но пока ничего подходящего в его работах не обнаружила.
У меня ситуация, имхо, не проще. Да, есть хорошая формула для представления натурального числа в виде суммы четырех квадратов целых чисел. Но в моем-то случае все четыре числа обязательно натуральные и попарно различные, нули исключены.