Метод полного перебора - Pascal ABC

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

Помогите!!! Напишите программу, которая определяет, сколько точек с целочисленными координатами содержится в сфере радиуса R с центром в начале системы координат. Подразумевается, что R – натуральное число, 1≤R≤30. Расстояние d между точкой с координатами (x, y, z) и началом системы координат определяется по формуле d=√(x^2+y^2+z^2)

Решение задачи: «Метод полного перебора»

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

var
  r, n, k, x, y, z: longint;
 
begin
  write('R = ');
  readln(r);
  for x := 1 to r do
    for y := 1 to r do
      for z := 1 to r do
        if x * x + y * y + z * z <= r * r
          then inc(n);
  for x := 1 to r do
    for y := 1 to r do
      if x * x + y * y <= r * r
        then inc(k);
  n := 8 * (n + k + r) + 1;
  write('n = ', n)
end.

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

В этом коде используется метод полного перебора для поиска всех возможных комбинаций трех чисел (x, y, z), удовлетворяющих условию, что их сумма не превышает заданное число r. В первой части кода:

Задается значение переменной r.
Для каждой комбинации трех чисел (x, y, z) от 1 до r:
- Проверяется, удовлетворяет ли сумма квадратов этих чисел условию xx + yy + zz <= rr.
- Если условию удовлетворяет, то к переменной n прибавляется единица. Во второй части кода:
Для каждой комбинации двух чисел (x, y) от 1 до r:
- Проверяется, удовлетворяет ли сумма квадратов этих чисел условию xx + yy <= r*r.
- Если условию удовлетворяет, то к переменной k прибавляется единица. В конце кода:
Значение переменной n увеличивается на 8 разность между n и k, и на единицу.
Выводится значение переменной n.