Кое-что добавлю.
Если четверка

- решение системы, то четверки


- тоже решения. Поэтому каждое решение заключено в арифметическую прогрессию между двумя "нулевыми" решениями. Пример для

:

Для простых

количество искомых решений в такой прогрессии равно

(если условиться

), а количество возможных прогрессий -

. Почему это происходит - не очень понятно, тут связь с количеством "нулевых" решений, которые имеют вид

. Можно отсюда вывести выражение для общего числа решений системы при простых

:

Доказать не берусь, но если формула верна, верны также следующие утверждения:
- Для попарно вз. простых

в натуральных числах разрешимо ровно одно из двух уравнений:

- В случае простых

для

однозначно определено

противоположной четности такое, что

. Тогда если

, то для чисел

разрешимо уравнение (1), в противном случае - уравнение (2). Сумма всех

, не превышающих

, для которых разрешимо уравнение (1)

. Сумма всех чисел, не превышающих

, для которых разрешимо ур. (2)

.
Или я сильно ошибаюсь.