Цепные дроби для sqrtd и квадратичных иррациональностей

Рассмотрим число [math]\alpha=[\sqrt{D}]+\sqrt{D}[/math]. Заметим, что оно приведённое [math]\alpha\gt 1, [\sqrt{D}]-\sqrt{D}\in(-1;0)[/math]. Тогда сразу следуют следующие утверждения

• число [math][\sqrt{D}]+\sqrt{D}[/math] представимо в виде чисто периодической цепной дроби.
• [math]\sqrt{D}[/math] представимо в виде цепной дроби из [math]a_0[/math] и периода.
• [math]\sqrt{D}=[\sqrt{D}]+\sqrt{D}-a_0[/math] значит [math]\sqrt{D}=\langle a_0, \overline{a_1,\cdots, a_n, 2a_0} \rangle[/math].

Теорема (Лагранж):

Число[math]\alpha[/math] представимо в виде периодической цепной дроби тогда и только тогда, когда [math]\alpha[/math] квадратичная иррациональность.

Доказательство:

[math]\triangleright[/math]

[math]\Rightarrow[/math]. [math]\alpha=\langle a_0,a_1,\cdots,\overline{a_k,\cdots a_n}\rangle[/math], тогда введём [math]\alpha_k=\langle \overline{a_k,\cdots, a_n}\rangle[/math]. Тогда [math]\alpha_k=\langle a_k,\cdots, a_n, \overline{\alpha_k} \rangle[/math]. [math]\alpha_k=\frac{P_n'\alpha_k+P_{n-1}'}{Q_n'\alpha_k+Q_{n-1}'}\Rightarrow Q_n'\alpha_k^2+(P_n'+Q_{n-1}')\alpha_k+P_{n-1}'=0[/math] Поэтому [math]\alpha_k[/math] квадратичная иррациональность, так как иррационально и удовлетворяет уравнению с целыми коэффициентами. Аналогично получим, что [math]\alpha = \frac{P_k\alpha_k+P_{k-1}}{Q_k\alpha_k+Q_{k-1}}[/math]. Поэтому и [math]\alpha[/math] квадратичная иррациональность. [math]\Leftarrow[/math]. Пусть [math]a\alpha^2+b\alpha+c=0[/math]. Разложим [math]\alpha[/math] в цепную дробь и для [math]\forall k:\alpha=\frac{P_k\alpha_k+P_{k-1}}{Q_k\alpha_k+Q_{k-1}}[/math]. Подставим в уравнение и заменим коэффициенты [math]A_k\alpha_k^2+B_k\alpha_k+C_k=0[/math]. Где [math]A_k=aP_k^2+bP_kQ_k+cQ_k^2[/math], [math]B_k=2aP_kP_{k-1}+bP_{k-1}Q_k+bP_kQ_{k-1}+2cQ_kQ_{k-1}[/math] и [math]C_k=aP_{k-1}^2+bP_{k-1}Q_{k-1}+cQ_{k-1}^2[/math]. Вычислим и упростим дискриминант и получим: [math]B_k^2-4A_kC_k=b^2-4ac[/math]. Ограничим [math]A_k,B_k,C_k[/math]. По тому, что [math]~|\alpha-\frac{P_k}{Q_k}|\lt \frac{1}{Q_k^2}[/math] имеем </tex>\frac{P_k}{Q_k}=\alpha-\frac{\epsilon}{Q_k^2},\epsilon\in(-1;1)</tex>.

[math]\triangleleft[/math]