Статья Автор: Лебедев Дмитрий

О представлении числа в виде суммы квадратов

Сумма двух квадратов есть число, значит некоторые числа можно представить в виде суммы двух квадратов 
Такие числа называют "двуквадратными"...
А если в виде суммы трех квадратов то "трехквадратными" ....
А если в виде суммы четырех квадратов, то....... - названия нет, поскольку .... существует теорема Лагранжа "О представлении любого числа в виде суммы четырех квадратов целых чисел"
Остановимся пока на проблеме "двуквадратных чисел".
Возможно Вы вспомните, что есть особое подмножество квадратов, которые можно представить в виде суммы двух других натуральных квадратов
52 = 42 +32; 132 = 122 + 52 и т.д. это 

 

Попробуем найти ответ  на вопрос 
Какие числа ГАРАНТИРОВАННО не являются "двуквадратными" 
Идею ответа можно найти, если внимательно посмотреть на таблицу умножения и ответить на вопрос
Натуральное число в разряде единиц имеет цифру 2. Может ли оно быть квадратом?
Рассмотрим таблицу
x % 4  x2 % 4
0 0
1 1
2 0
3 1
  Значит x2 + y2 по модулю 4 могут иметь остатки 0 (0+0); 1 (0 +1) и 2(1+1)
  Таким образом, если число имеет вид 4k + 3, то оно не может быть "двуквадратным"
   Если число имеет вид 4k, то  оно может быть предствалено только в виде чётных квадратов, то есть 
   оно "двуквадратное" <=> тогда число k двуквадратное
   Если число имеет вид 4k+2, то оно состваное и может быть представлено в виде нечетный квадратов
   Если число имеет вид 4k +1, то здесь можно выделить класс простых чисел, для которых есть утверждение, названнове "Рождественской теоремой Ферма" 
"Рождественская теорема Ферма" утверждает, что любое простое число, представимое в виде 4k+1 является "двуквадратным". Существует даже формула для разложения (вряд ли мы сможет её использовать sad). 
Не сложно убедиться, что для составных чисел представление в виде 4k + 1 не гарантирует "двуквадратности". Самым легким примером будет число 21.
А как быть с составными числами?

Пусть \(P = N \cdot M; N = a^2+b^2; M = c^2 + d^2\)
Тогда \(P = (a^2+b^2)\cdot(c^2 + d^2) = (..I..)^2 + (. .II..)^2 \) 
но это невозможно coolcoolcool

"Гораций, много в мире есть того,
Что вашей философии не снилось."
Воспользуемся тождеством Брахмагупты — Фибоначчи (или Диофанта

\((a^2+b^2)\cdot(c^2 + d^2) = (ad - bc)^2 + (ac + bd)^2 = (ac - bd)^2 + (ad + bc)^2\)
(желающие могут увидеть в этом тождестве магию векторных произведений / скалярного и косого)
Чтобы понять/запомнить как получаются такие тождества посмотрим на цепочку преобразований
\( 1=1\\ 1\cdot1 = 1 (sin^2\alpha + cos^2\alpha)\cdot(sin^2\beta + cos^2\beta) = (sin^2(\alpha-\beta) + cos^2(\alpha-\beta))\\ (sin^2\alpha + cos^2\alpha)\cdot(sin^2\beta + cos^2\beta) = (sin\alpha\cdot cos\beta-cos\alpha\cdot sin\beta)^2 + (cos\alpha\cdot cos\beta+sin\alpha\cdot sin\beta)^2 \)

 

Теперь становиться понятным утверждение теоремы Ферма-Эйлера

Натуральное число представимо в виде суммы двух квадратов (целых чисел) тогда и только тогда, когда ни одно простое число вида 4k+3
не входит в его разложение на простые множители в нечётной степени.
ОООк
Теперь можно и приступить к предствавлению числа в виде суммы квадратов двух целых чисел.

Проверяем число P на "двуквадратность"
  1. Находим P % 4 и если равно 3 выдаем False и завершаем работу
  2. Производим факторизацию числа P - раскладываем на простые множетели в соответствующих степенях.
  3. Из разложения выделяем "максимальный квадрат" и получаем, выражение вида
    \(P = Q^2\cdot p_1\cdot p_2\dots\cdot p_k, \ все\ p_i - \ простые\)
  4. Если среди \(p_1,\ p_2\dots p_k\) есть числа вида 4k + 3, выдаем False и завершаем работу
  5. Каждое из чисел pi предстваляем как "двуквадратное", используя тождество Диофанта последовательно "соединяем" предстваления в одно
  6. Умножаем числа предстваления на Q и получаем ответ.
Пример 1. P = 1580 = 22*5*79. 79 = 4*19 + 3 => число 1580 не является "двуквадратным"
Пример 2. P = 27144 = 2* 3* 13 * 29 = 6* (1+ 12) * (3+ 22) * (5+ 22) = 6* ((3 * 1 - 2 * 1)2 + (3 * 1 + 2 * 1)2) * (5+ 22) =
=6* (5+ 12) * (5+ 22) = 62 * ((5 * 2 - 5* 1 )2 + (5 * 5 + 2 * 1)2)  62 * ( 52 + 272) = 302 + 1622 = 900 + 26244 = 27144
Самым сложным остался вопрос предстваления в "двуквадратном" виде простого числа вида 4k+1
Какие будут предложения?
 

\(A = a_1, \dots, a_n,\dots\)Представление числа в виде суммы двуж квадратов.
 
Пропустить Навигационные Ссылки.
Чтобы оставить комментарий нужна авторизация
Печать