Zenn
Open16

ペル方程式

Toru3Toru3

x2ny2=1x^2-ny^2=1をペル方程式という。
(x,y)=(1,0)(x, y) = (1, 0)は自明な解である。
nnが平方数でない場合は非自明な解を持つ。
以下nnは平方数でない正の整数とする。

右辺の111-144, 4-4などとしたものを考えることもある。

Toru3Toru3

ペル方程式は基本解をx1+y1nx_1+y_1\sqrt{n}とすると全ての解はxk+ykn=(x1+y1n)kx_k + y_k \sqrt{n} = (x_1+y_1\sqrt{n})^kと表せる。

n\sqrt{n}の連分数展開を使うとペル方程式の基本解が求まる。

Toru3Toru3

n\sqrt{n}連分数展開を[a0;a1,a2,,ar][a_0; \overline{a_1, a_2, \dots, a_r}][1]とする。

pkqk:=[a0;a1,a2,,ak1]\frac{p_k}{q_k}:=[a_0; a_1, a_2, \dots, a_{k-1}]

とする。

  • rrが偶数のとき
    • x2ny2=1x^2-ny^2=-1は解を持たない
    • (pr,qr)(p_r, q_r)x2ny2=1x^2-ny^2=1の基本解
  • rrが奇数のとき
    • (pr,qr)(p_r, q_r)x2ny2=1x^2-ny^2=-1の基本解
    • x1+y1n=(pr+qrn)2x_1+y_1\sqrt{n}=(p_r+q_r\sqrt{n})^2x2ny2=1x^2-ny^2=1の基本解
脚注
  1. n\sqrt{n}の連分数展開は必ず周期的になるのでこの形だけ考えればよい。また、a0=na_0 = \lfloor \sqrt{n} \rfloor, ar=2na_r = 2\lfloor \sqrt{n} \rfloorとなる。また、a1,a2,...,ar2,ar1a_1, a_2,..., a_{r-2}, a_{r-1}は回文となる。 ↩︎

Toru3Toru3

pkp_k, qkq_kp0=1p_0=1, q0=0q_0=0, p1=a0p_1=a_0, q1=1q_1=1, pk=ak1pk1+pk2p_k=a_{k-1}p_{k-1}+p_{k-2}, qk=ak1qk1+qk2q_k=a_{k-1}q_{k-1}+q_{k-2}という三項間漸化式で表せる。
よって、漸化式は2x2行列を使って以下のように書ける。

[pkqkpk1qk1]=[ak1110][pk1qk1pk2qk2] \begin{bmatrix} p_k & q_k \\ p_{k-1} & q_{k-1} \end{bmatrix}= \begin{bmatrix} a_{k-1} & 1\\ 1 & 0 \end{bmatrix} \begin{bmatrix} p_{k-1} & q_{k-1}\\ p_{k-2} & q_{k-2} \end{bmatrix}

繰り返し展開すると以下のようになる。

[pkqkpk1qk1]=[ak1110][ak2110][a1110][a0110] \begin{bmatrix} p_k & q_k\\ p_{k-1} & q_{k-1} \end{bmatrix}= \begin{bmatrix} a_{k-1} & 1\\ 1 & 0 \end{bmatrix} \begin{bmatrix} a_{k-2} & 1\\ 1 & 0 \end{bmatrix}\cdots \begin{bmatrix} a_1 & 1\\ 1 & 0 \end{bmatrix} \begin{bmatrix} a_0 & 1\\ 1 & 0 \end{bmatrix}
Toru3Toru3

途中の pk,qkp_k, q_kpk2nqk=(1)k1γk+1p_k^2-nq_k=(-1)^{k-1}\gamma_{k+1} となる。
ここでαk=n+βkγk\alpha_k=\dfrac{\sqrt{n}+\beta_k}{\gamma_k}, α0=n\alpha_0 = \sqrt{n}, αk+1=1αkαk\alpha_{k+1}=\dfrac{1}{\alpha_k-\lfloor\alpha_k\rfloor}

Toru3Toru3
Toru3Toru3

Amthor's solutionの項を読むとある整数mmとsquare freeな整数nn'を用いてn=m2nn=m^2n'と表せるnnの解を効率よく求めることが出来そう。

  1. nn'での基本解(x,y)(x', y')を求める
  2. x+myn=(x+yn)kx+my\sqrt{n'}=(x'+y'\sqrt{n'})^kと表せるkkを見つける
  3. 上の(x,y)(x, y)が解
    • x+myn=x+ym2n=x+ynx+my\sqrt{n'}=x+y\sqrt{m^2n'}=x+y\sqrt{n}
Toru3Toru3

Smooth numbersの項を見ると2次ふるい法と似たような手法で高速化出来そう。(あとでちゃんと読む)

Toru3Toru3

x^2-ny^2=1の解は求まるが、x^2-ny^2=-1の解は求まらなさそう

Toru3Toru3

x24ny2=±4x^2-4ny^2=\pm4の解とx2ny2=±1x^2-ny^2=\pm1の解は(x,y)(x/2,y)(x, y)\mapsto(x/2, y)で1対1対応する。[1]

脚注
  1. 整数論1 第1版のp.94 定理4.2.3 ↩︎

Toru3Toru3

n≢0(mod4)n \not\equiv 0 \pmod 4かつx2ny2=±4x^2-ny^2=\pm 4の解が存在する\Rightarrow u2nv2=±1u^2-nv^2=\pm 1の解が存在してu+vn=(x+yn2)3u+v\sqrt{n}=\left(\dfrac{x+y\sqrt{n}}{2}\right)^3

u2nv2=±1u^2-nv^2=\pm 1の解が存在する\Rightarrowx2ny2=±4x^2-ny^2=\pm 4の解が存在して(x,y)=(2u,2v)(x, y)=(2u, 2v)

Toru3Toru3

証明

x+yn2xyn2=x2y2n4=±1\dfrac{x+y\sqrt{n}}{2}\cdot\dfrac{x-y\sqrt{n}}{2}=\dfrac{x^2-y^2n}{4}=\pm 1 より(有理数)解である.

ここでyyが偶数であるとすると±4+y2\pm 4 + y^2は4の倍数となる。よってxxも偶数となる。
また、xxが偶数であるとするとx24x^2\mp 4は4の倍数となる。ここでn≢0(mod4)n \not\equiv 0 \pmod{4}よりy2y^2は少なくとも偶数である必要がある。よってyyは偶数となる。
よって、(x,y)(x, y)の組としては(偶数,偶数)(偶数, 偶数)あるいは(奇数,奇数)(奇数, 奇数)しかない。
(偶数,偶数)(偶数, 偶数)の場合、x+yn2Z[n]\dfrac{x+y\sqrt{n}}{2}\in \mathbb{Z}[\sqrt{n}]. よって(x+yn2)3Z[n]\left(\dfrac{x+y\sqrt{n}}{2}\right)^3\in \mathbb{Z}[\sqrt{n}]. 従って整数解となる。
よって、以下では(奇数,奇数)(奇数, 奇数)の場合のみ考える。

(x+yn2)3=x(x2+3y2n)8+y(y2n+3x2)8n\left(\dfrac{x+y\sqrt{n}}{2}\right)^3=\dfrac{x(x^2+3y^2n)}{8}+\dfrac{y(y^2n+3x^2)}{8}\sqrt{n}.

x2+3y2n=x2+3(x24)=4x212=4(x23)x^2+3y^2n=x^2+3(x^2\mp 4)=4x^2\mp 12=4(x^2\mp 3).
よってx(x2+3y2n)=4x(x23)x(x^2+3y^2n)=4x(x^2\mp 3).
ここでxxは奇数なのでx23x^2\mp 3は偶数. よって84x(x23)8 \mid 4x(x^2\mp 3).

y2n+3x2=(x24)+3x2=4x24=4(x21)y^2n+3x^2=(x^2\mp 4)+3x^2=4x^2\mp 4=4(x^2\mp 1).
よってy(y2n+3x2)=4y(x21)y(y^2n+3x^2)=4y(x^2\mp 1).
ここでxxは奇数なのでx21x^2\mp 1は偶数. よって84y(x21)8 \mid 4y(x^2\mp 1).

以上より(x+yn2)3=x(x2+3y2n)8+y(y2n+3x2)8nZ[n]\left(\dfrac{x+y\sqrt{n}}{2}\right)^3=\dfrac{x(x^2+3y^2n)}{8}+\dfrac{y(y^2n+3x^2)}{8}\sqrt{n}\in \mathbb{Z}[\sqrt{n}].

Toru3Toru3

(x,y)=(奇数,奇数)(x,y)=(奇数, 奇数)の場合, x21(mod4)x^2 \equiv 1 \pmod{4}, y21(mod4)y^2 \equiv 1 \pmod{4}よりn1(mod4)n\equiv 1\pmod{4}である必要がある.

ログインするとコメントできます