MathLabs

Số học và Lý thuyết số

Phương trình Pell

Phương trình x2−dy2=1x^2-dy^2=1, có nghiệm nguyên liên hệ với phân số liên tục.

Trực giácTrực giác: điểm nguyên trên hyperbol xấp xỉ d\sqrt{d}

Tìm số nguyên x,yx,y với x2−2y2=1x^2-2y^2=1. Thử các giá trị nhỏ: (x,y)=(3,2)(x,y)=(3,2) thỏa, vì 9−8=19-8=1. Biến đổi lại, (xy)2=2+1y2\left(\frac{x}{y}\right)^2=2+\frac{1}{y^2}, nên x/y=3/2=1.5x/y=3/2=1.5 đã là một xấp xỉ đáng kinh ngạc của 2≈1.41421\sqrt{2}\approx 1.41421. Các phương trình dạng x2−dy2=1x^2 - d y^2 = 1 — gọi là phương trình Pell — hóa ra không thể tách rời khỏi lý thuyết phân số liên tục và các xấp xỉ hữu tỉ tốt nhất.

Đồ thị đường đi của các giản phân với một đỉnh được tô nổi bật là nghiệm cơ bản của phương trình Pell.
Các giản phân của d\sqrt{d} vẽ dưới dạng đồ thị đường đi: mỗi đỉnh là một giản phân pk/qkp_k/q_k, đỉnh tô nổi bật là nghiệm cơ bản.

Nâng caoĐịnh nghĩa, nghiệm cơ bản và cấu trúc

Định nghĩa: Phương trình Pell và nghiệm cơ bản

Với số nguyên dương dd không phải số chính phương, x2−dy2=1x^2 - d y^2 = 1 được gọi là phương trình Pell. Trong các nghiệm nguyên dương (x,y)(x,y) của nó, nghiệm có xx nhỏ nhất (tương đương yy nhỏ nhất) là nghiệm cơ bản (x1,y1)(x_1,y_1).

x2−dy2=1,d∈Z>0 not a perfect squarex^2 - d y^2 = 1,\qquad d\in\mathbb{Z}_{>0}\ \text{not a perfect square}

Phân tích vế trái thành (x−yd)(x+yd)=1(x-y\sqrt d)(x+y\sqrt d)=1 biến phương trình Pell thành một phát biểu về vành Z[d]\mathbb{Z}[\sqrt d]: nghiệm tương ứng với phần tử có chuẩn 11, và nhân hai phần tử như vậy cho một phần tử khác. Đặc biệt, một khi biết (x1,y1)(x_1,y_1), mọi nghiệm tiếp theo có được bằng cách lấy lũy thừa: xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n.

(x1+y1d)(x1−y1d)=1 ⟹ xn+ynd=(x1+y1d)n(x_1+y_1\sqrt d)(x_1-y_1\sqrt d)=1 \ \Longrightarrow\ x_n+y_n\sqrt d = (x_1+y_1\sqrt d)^n
Nghiệm cơ bản với dd nhỏ
ddPhân số liên tục của d\sqrt dNghiệm cơ bản (x1,y1)(x_1,y_1)
22[1;2‾][1;\overline{2}](3,2)(3,2)
33[1;1,2‾][1;\overline{1,2}](2,1)(2,1)
55[2;4‾][2;\overline{4}](9,4)(9,4)
77[2;1,1,1,4‾][2;\overline{1,1,1,4}](8,3)(8,3)

Nâng caoĐịnh lý và chứng minh

Với mọi số nguyên dương dd không phải số chính phương, x2−dy2=1x^2 - d y^2 = 1 có nghiệm với x,yx,y là số nguyên dương.

Vì sao đúng?

Hoàn toàn không hiển nhiên rằng hyperbol x2−dy2=1x^2-dy^2=1, chắc chắn có điểm thực, phải đi qua một điểm nguyên — định lý này đảm bảo điều đó luôn đúng, với mọi dd không chính phương, nhờ một lập luận chuồng bồ câu khéo léo trên các xấp xỉ hữu tỉ.

Chứng minh

Theo định lý xấp xỉ Dirichlet, với mọi số nguyên Q>0Q>0 tồn tại số nguyên p,qp,q với 1≤q≤Q1\le q\le Q và ∣d−pq∣<1qQ\left|\sqrt d - \frac{p}{q}\right|<\frac{1}{qQ}. Cho Q→∞Q\to\infty sinh ra vô hạn cặp (p,q)(p,q) với ∣d−pq∣<1q2\left|\sqrt d-\frac{p}{q}\right|<\frac{1}{q^2}.

Với mỗi cặp như vậy, ∣p−qd∣<1q|p-q\sqrt d|<\frac{1}{q}, nên ∣p2−dq2∣=∣p−qd∣⋅∣p+qd∣<1q(2qd+1q)<2d+1|p^2-dq^2|=|p-q\sqrt d|\cdot|p+q\sqrt d|<\frac{1}{q}\left(2q\sqrt d+\frac{1}{q}\right)<2\sqrt d+1. Vậy p2−dq2p^2-dq^2 chỉ nhận một trong hữu hạn giá trị nguyên trong (−2d−1, 2d+1)(-2\sqrt d-1,\,2\sqrt d+1), trong khi có vô hạn cặp (p,q)(p,q).

Theo nguyên lý chuồng bồ câu, một số nguyên khác không cố định kk trong khoảng hữu hạn đó thỏa p2−dq2=kp^2-dq^2=k với vô hạn cặp (p,q)(p,q). Trong vô hạn cặp này, lại theo chuồng bồ câu, vô hạn cặp có cùng lớp dư p≡p0, q≡q0(mod∣k∣)p\equiv p_0,\ q\equiv q_0\pmod{|k|}.

Lấy hai cặp phân biệt (p1,q1)≠(p2,q2)(p_1,q_1)\neq(p_2,q_2) như vậy với p12−dq12=p22−dq22=kp_1^2-dq_1^2=p_2^2-dq_2^2=k và cùng lớp dư theo môđun ∣k∣|k|. Đặt x+yd=(p1+q1d)(p2−q2d)kx+y\sqrt d = \frac{(p_1+q_1\sqrt d)(p_2-q_2\sqrt d)}{k}; khai triển cho thấy x=p1p2−dq1q2kx=\frac{p_1p_2-dq_1q_2}{k} và y=p1q2−p2q1ky=\frac{p_1q_2-p_2q_1}{k} thực sự là số nguyên chính vì cùng lớp dư theo môđun ∣k∣|k|, và tính nhân của chuẩn N(a+bd)=a2−db2N(a+b\sqrt d)=a^2-db^2 cho x2−dy2=k⋅kk2=1x^2-dy^2=\frac{k\cdot k}{k^2}=1. Vì (p1,q1)≠(p2,q2)(p_1,q_1)\neq(p_2,q_2) nhưng chúng cho cùng tỉ số ở giới hạn, kiểm tra được y≠0y\neq 0, và thay (x,y)(x,y) bằng (∣x∣,∣y∣)(|x|,|y|) (vẫn là nghiệm, vì chỉ xuất hiện bình phương) cho một nghiệm nguyên dương.

Nếu (x1,y1)(x_1,y_1) là nghiệm cơ bản của x2−dy2=1x^2 - d y^2 = 1, thì mọi nghiệm nguyên dương (x,y)(x,y) bằng (xn,yn)(x_n,y_n) với n≥1n\ge 1 nào đó, trong đó xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n.

Vì sao đúng?

Nó nói rằng vô hạn nghiệm không phải một tập rải rác bí ẩn mà là một dãy kiểu cấp số nhân hoàn toàn dự đoán được, sinh ra bằng cách "nhân" lặp đi lặp lại nghiệm nhỏ nhất — đưa một tìm kiếm vô hạn về việc chỉ cần tìm một số.

Chứng minh

Trước hết kiểm tra (xn,yn)(x_n,y_n) định nghĩa bởi xn+ynd=(x1+y1d)nx_n+y_n\sqrt d=(x_1+y_1\sqrt d)^n đúng là nghiệm với mọi nn: lấy liên hợp, xn−ynd=(x1−y1d)nx_n-y_n\sqrt d=(x_1-y_1\sqrt d)^n, nên xn2−dyn2=(xn+ynd)(xn−ynd)=[(x1+y1d)(x1−y1d)]n=(x12−dy12)n=1n=1x_n^2-dy_n^2=(x_n+y_n\sqrt d)(x_n-y_n\sqrt d)=\left[(x_1+y_1\sqrt d)(x_1-y_1\sqrt d)\right]^n=(x_1^2-dy_1^2)^n=1^n=1.

Giờ giả sử (x,y)(x,y) là một nghiệm nguyên dương bất kỳ không có dạng này; vì xn→∞x_n\to\infty khi n→∞n\to\infty, tồn tại duy nhất nn với xn+ynd≤x+yd<xn+1+yn+1d=(xn+ynd)(x1+y1d)x_n+y_n\sqrt d \le x+y\sqrt d < x_{n+1}+y_{n+1}\sqrt d = (x_n+y_n\sqrt d)(x_1+y_1\sqrt d).

Chia cả hai vế cho (xn+ynd)(x_n+y_n\sqrt d), tức nhân với nghịch đảo (xn−ynd)(x_n-y_n\sqrt d) (hợp lệ vì xn2−dyn2=1x_n^2-dy_n^2=1): đặt x′+y′d=(x+yd)(xn−ynd)x'+y'\sqrt d = (x+y\sqrt d)(x_n-y_n\sqrt d). Khi đó 1≤x′+y′d<x1+y1d1\le x'+y'\sqrt d < x_1+y_1\sqrt d, và x′2−dy′2=(x2−dy2)(xn2−dyn2)=1⋅1=1x'^2-dy'^2=(x^2-dy^2)(x_n^2-dy_n^2)=1\cdot 1=1, nên (x′,y′)(x',y') cũng là nghiệm của phương trình Pell.

Một tính toán ngắn dùng x′+y′d≥1x'+y'\sqrt d\ge 1 và x′2−dy′2=1x'^2-dy'^2=1 cho thấy x′≥1x'\ge 1 và y′≥0y'\ge 0 (một nghiệm với x′+y′d≥1x'+y'\sqrt d\ge 1 nhưng y′<0y'<0 sẽ buộc x′>x1x'>x_1, mâu thuẫn với x′+y′d<x1+y1dx'+y'\sqrt d<x_1+y_1\sqrt d kết hợp với phương trình chuẩn). Nếu y′>0y'>0, thì (x′,y′)(x',y') là nghiệm dương với x′+y′d<x1+y1dx'+y'\sqrt d<x_1+y_1\sqrt d, mâu thuẫn với tính nhỏ nhất của nghiệm cơ bản. Vậy y′=0y'=0, buộc x′=1x'=1, tức x+yd=xn+yndx+y\sqrt d=x_n+y_n\sqrt d, nên (x,y)=(xn,yn)(x,y)=(x_n,y_n) suy cho cùng — mâu thuẫn với giả thiết.

Vậy mọi nghiệm dương chính xác là một (xn,yn)(x_n,y_n) nào đó, chứng minh nghiệm cơ bản sinh ra toàn bộ tập nghiệm.

Nâng caoỨng dụng thực tiễn và Ví dụ minh họa

Phương trình Pell không chỉ là một câu đố: nó chi phối các xấp xỉ hữu tỉ tốt nhất dùng trong thiết kế bánh răng cơ khí, số học của nó nằm sau các thuật toán phân tích thừa số nguyên cổ điển, và việc tính nghiệm cơ bản của nó với dd lớn chính là bài toán tính toán mà thuật toán lượng tử của Hallgren giải được trong thời gian đa thức — mà chưa có phương pháp cổ điển hiệu quả nào được biết.

Ví dụ: Tìm nghiệm cơ bản với d=2d=2

Dùng phân số liên tục của 2\sqrt{2} để tìm nghiệm cơ bản của x2−2y2=1x^2-2y^2=1, rồi sinh ra nghiệm tiếp theo.

Lời giải

Phân số liên tục của 2\sqrt2 là [1;2‾]=1+12+12+⋯[1;\overline{2}]=1+\cfrac{1}{2+\cfrac{1}{2+\cdots}}, với các giản phân 1, 32, 75, 1712, 4129,…1,\ \frac{3}{2},\ \frac{7}{5},\ \frac{17}{12},\ \frac{41}{29},\dots

Kiểm tra giản phân 32\frac{3}{2}: 32−2⋅22=9−8=13^2-2\cdot 2^2=9-8=1. Đây là giản phân nhỏ nhất thỏa mãn, nên nghiệm cơ bản là (x1,y1)=(3,2)(x_1,y_1)=(3,2).

Dùng công thức truy hồi xn+ynd=(x1+y1d)nx_n+y_n\sqrt{d} = (x_1+y_1\sqrt{d})^n với n=2n=2: x2+y22=(3+22)2=9+122+8=17+122x_2+y_2\sqrt2=(3+2\sqrt2)^2=9+12\sqrt2+8=17+12\sqrt2, nên (x2,y2)=(17,12)(x_2,y_2)=(17,12).

Kiểm chứng: 172−2⋅122=289−288=117^2-2\cdot 12^2=289-288=1, xác nhận nghiệm tiếp theo, khớp với giản phân 1712\frac{17}{12} tìm được ở trên — đúng như lý thuyết dự đoán.

Ví dụ: Xấp xỉ kỹ thuật: thiết kế tỉ số bánh răng gần 2\sqrt{2}

Một cơ cấu cần hai bánh răng ăn khớp có tỉ số răng xấp xỉ 2\sqrt2 càng gần càng tốt bằng số răng nguyên nhỏ, để sai số vẫn nhỏ sau nhiều vòng quay. Dùng các giản phân của phương trình Pell để chọn số răng và chặn sai số.

Lời giải

Từ nghiệm cơ bản (3,2)(3,2) của x2−2y2=1x^2-2y^2=1, tỉ số 3/2=1.53/2=1.5 cho ∣3/2−2∣≈0.0858|3/2-\sqrt2|\approx 0.0858 — dùng được, nhưng thô với máy móc cần độ chính xác.

Dùng giản phân tiếp theo từ (17,12)(17,12) (tìm được bằng cách bình phương 3+223+2\sqrt2): cặp bánh răng với 1717 và 1212 răng cho tỉ số 17/12≈1.4166717/12\approx 1.41667, và ∣17/12−2∣≈0.00245|17/12-\sqrt2|\approx 0.00245, chính xác hơn gần 3535 lần so với thiết kế 3/23/2 mà chỉ tăng số răng vừa phải.

Chặn tổng quát cho mọi giản phân pk/qkp_k/q_k của một phân số liên tục là ∣pkqk−2∣<1qk2\left|\frac{p_k}{q_k}-\sqrt2\right|<\frac{1}{q_k^2}, nên khi kỹ sư chuyển sang nghiệm Pell tiếp theo (x3,y3)(x_3,y_3) (có được từ (3+22)3=99+702(3+2\sqrt2)^3=99+70\sqrt2, tức 99/7099/70) số răng tăng lên 7070 nhưng sai số giảm xuống dưới 1/702≈0.00021/70^2\approx 0.0002.

Điều này minh họa sự đánh đổi kỹ thuật cơ bản hiện rõ trong phương trình Pell: mỗi nghiệm liên tiếp (xn,yn)(x_n,y_n) đánh đổi việc tăng chi phí sản xuất (nhiều răng hơn, tức yny_n lớn hơn) lấy một sai số xấp xỉ giảm nhanh có thể định lượng, cho phép người thiết kế chọn đúng điểm trên đường cong này phù hợp ngân sách và yêu cầu độ chính xác.

Nghiệm cơ bản của x2−2y2=1x^2-2y^2=1 là gì?

Nếu (x1,y1)(x_1,y_1) là nghiệm cơ bản, nghiệm (x2,y2)(x_2,y_2) được suy ra thế nào?

x2−3y2=−1x^2-3y^2=-1 có nghiệm nguyên không?

Thuật toán Hallgren giải phương trình Pell trong thời gian đa thức bằng loại tính toán nào?

Tài liệu tham khảo

  1. Sean Hallgren (2007). Polynomial-time quantum algorithms for Pell's equation and the principal ideal problem · DOI:10.1145/1206035.1206039
  2. Hendrik W. Lenstra Jr. (2002). Solving the Pell Equation