Bài toán đàn bò của Archimedes và phương trình Pell
Cạm Bẫy của Những Con Số Ngây Thơ và Thách Thức Archimedes
Vào thế kỷ 3 TCN, nhà bác học vĩ đại Archimedes xứ Syracuse đã gửi một bài thơ gồm 44 câu thơ bằng tiếng Hy Lạp cổ đến nhà toán học Eratosthenes xứ Cyrene. Bài thơ chứa đựng một bài toán đố được gọi là Bài toán Đàn bò của Archimedes (Problema Bovinum), thách thức giới trí thức đương thời tính toán số lượng bò trong đàn bò của Thần Mặt Trời Apollo.
Đề bài đưa ra một hệ các điều kiện tỉ lệ giữa các con bò đực và bò cái thuộc 4 màu sắc: trắng, đen, vàng và đốm. Nhìn thoáng qua, đây có vẻ như một bài toán đố dân gian vui nhộn. Tuy nhiên, khi biến đổi các điều kiện ràng buộc hình học (bò trắng và bò đen xếp thành một hình vuông, bò vàng và bò đốm xếp thành một hình tam giác), bài toán dẫn trực tiếp đến một phương trình Diophantine bậc hai có dạng:
$$x^2 – d y^2 = 1$$
Trực giác ban đầu khiến người ta nghĩ rằng đáp án sẽ là một vài nghìn hay vài triệu con bò. Thế nhưng, câu trả lời thực sự đã làm rúng động toàn bộ lịch sử toán học: Số lượng bò nhỏ nhất thỏa mãn bài toán Archimedes là một con số khổng lồ chứa tới 206.545 chữ số!
Nhà toán học Peter Koymans (Đại học Leiden) từng bình luận trên Quanta Magazine: “Để in ra nghiệm nhỏ nhất của bài toán Archimedes, bạn cần tới 50 trang giấy. Theo một nghĩa nào đó, đây là một trò đùa khổng lồ mà Archimedes đã để lại cho nhân loại.”
Phát Biểu Đầy Đủ Bài Toán Đàn Bò của Archimedes
Để thấy được đỉnh cao tư duy và sự tinh quái của Archimedes, hãy cùng xem xét phát biểu toán học đầy đủ của bài toán này.
Gọi $W, B, Y, D$ lần lượt là số lượng bò đực các màu Trắng (White), Đen (Black), Vàng/Nâu (Yellow), và Đốm (Dappled). Gọi $w, b, y, d$ lần lượt là số lượng bò cái cùng các màu Trắng, Đen, Vàng, và Đốm.
Phần 1: Hệ Ràng Buộc Tỉ Lệ Số Lượng (7 Phương Trình Tuyến Tính)
Archimedes đưa ra 7 điều kiện tỉ lệ giữa các loại bò đực và bò cái như sau:
Số bò đực trắng bằng số bò đực vàng cộng thêm $(\frac{1}{2} + \frac{1}{3})$ số bò đực đen: $$W = Y + \left(\frac{1}{2} + \frac{1}{3}\right) B = Y + \frac{5}{6} B$$
Số bò đực đen bằng số bò đực vàng cộng thêm $(\frac{1}{4} + \frac{1}{5})$ số bò đực đốm: $$B = Y + \left(\frac{1}{4} + \frac{1}{5}\right) D = Y + \frac{9}{20} D$$
Số bò đực đốm bằng số bò đực vàng cộng thêm $(\frac{1}{6} + \frac{1}{7})$ số bò đực trắng: $$D = Y + \left(\frac{1}{6} + \frac{1}{7}\right) W = Y + \frac{13}{42} W$$
Số bò cái trắng bằng $(\frac{1}{3} + \frac{1}{4})$ tổng đàn bò đen (cả đực và cái): $$w = \left(\frac{1}{3} + \frac{1}{4}\right) (B + b) = \frac{7}{12} (B + b)$$
Số bò cái đen bằng $(\frac{1}{4} + \frac{1}{5})$ tổng đàn bò đốm (cả đực và cái): $$b = \left(\frac{1}{4} + \frac{1}{5}\right) (D + d) = \frac{9}{20} (D + d)$$
Số bò cái đốm bằng $(\frac{1}{5} + \frac{1}{6})$ tổng đàn bò vàng (cả đực và cái): $$d = \left(\frac{1}{5} + \frac{1}{6}\right) (Y + y) = \frac{11}{30} (Y + y)$$
Số bò cái vàng bằng $(\frac{1}{6} + \frac{1}{7})$ tổng đàn bò trắng (cả đực và cái): $$y = \left(\frac{1}{6} + \frac{1}{7}\right) (W + w) = \frac{13}{42} (W + w)$$
Nếu chỉ giải Phần 1 (hệ 7 phương trình 8 ẩn trên $\mathbb{Z}^+$), nghiệm nguyên dương nhỏ nhất cho tổng số đàn bò đã lên tới 50.389.082 con (trong đó $W = 10.366.482$, $B = 7.460.514$, $D = 7.358.060$, $Y = 4.149.387$).
Phần 2: Điều Kiện Ràng Buộc Hình Học (2 Điều Kiện Bẻ Khóa)
Thách thức thực sự nằm ở 2 điều kiện bổ sung được Archimedes cài cắm ở cuối bài thơ:
Điều kiện Hình Vuông: Tổng số bò đực trắng và bò đực đen khi đứng chung tạo thành một hình vuông hoàn chỉnh (tức là một số chính phương $m^2$): $$W + B = m^2$$
Điều kiện Hình Tam Giác: Tổng số bò đực vàng và bò đực đốm khi đứng chung tạo thành một hình tam giác (tức là một số tam giác $\frac{n(n+1)}{2}$): $$D + Y = \frac{n(n+1)}{2}$$
Khi kết hợp Phần 1 và Phần 2, sau các bước thế đại số phức tạp để triệt tiêu ẩn, bài toán Archimedes dẫn trực tiếp về phương trình Pell đại diệt vong:
$$x^2 – 410,286,423,278,424 , y^2 = 1$$
hoặc sau khi rút gọn hệ số chính phương:
$$u^2 – 4,729,494 , v^2 = 1$$
Tên Gọi Nhầm Lẫn Lịch Sử và Đồng Nhất Thức Brahmagupta
Phương trình $x^2 – d y^2 = 1$ (với $d$ là số nguyên dương không chính phương) ngày nay được cả thế giới biết đến với tên gọi Phương trình Pell. Tuy nhiên, đây là một trong những sự nhầm lẫn lịch sử nổi tiếng nhất của toán học. Nhà toán học vĩ đại Leonhard Euler vào thế kỷ 18 đã gán nhầm tên của nhà toán học người Anh John Pell (1611–1685) cho phương trình này, dù Pell hầu như không đóng góp công trình nào liên quan. Người mà Euler muốn nhắc đến thực chất là William Brouncker (1620–1684).
Trên thực tế, phương trình này đã được nghiên cứu sâu sắc từ hơn một thiên niên kỷ trước đó tại Ấn Độ. Vào năm 628 CN, nhà toán học Brahmagupta trong tác phẩm kinh điển Brahmasphutasiddhanta đã phát minh ra một hằng đẳng thức đại số kỳ diệu được gọi là Phép hợp thành Samasa:
Đồng nhất thức Brahmagupta (Samasa Identity): Nếu $(x_1, y_1)$ là nghiệm của $x^2 – d y^2 = k_1$ và $(x_2, y_2)$ là nghiệm của $x^2 – d y^2 = k_2$, thì: $$(x_1^2 – d y_1^2)(x_2^2 – d y_2^2) = (x_1 x_2 + d y_1 y_2)^2 – d(x_1 y_2 + x_2 y_1)^2 = k_1 k_2$$
Đồng nhất thức này chứng minh rằng nếu ta có hai nghiệm $(x_1, y_1)$ và $(x_2, y_2)$ của phương trình $x^2 – d y^2 = 1$ ($k_1 = k_2 = 1$), ta có thể “nhân” chúng lại với nhau để sinh ra một nghiệm mới $(x_3, y_3)$ hoàn toàn nguyên:
Phát kiến của Brahmagupta chính là tiền thân của Cấu trúc Nhóm Đơn vị (Unit Group) trong đại số hiện đại. Nghiệm nguyên của phương trình Pell chính là các đơn vị trong vành số nguyên $\mathbb{Z}[\sqrt{d}]$ của trường số vô tỷ $\mathbb{Q}(\sqrt{d})$.
Phân Số Liên Tục: Chìa Khóa Khai Mở Bức Tường Lập Phương
Dù Brahmagupta đã tìm ra cách sinh ra vô số nghiệm từ một nghiệm ban đầu, câu hỏi hóc húa nhất vẫn tồn tại: Làm sao để tìm thấy nghiệm nguyên nhỏ nhất (Nghiệm cơ bản – Fundamental Solution)?
Đến thế kỷ 12, nhà toán học Ấn Độ Bhāskara II đã sáng tạo ra thuật toán vòng lặp Chakravala. Đến năm 1768, nhà toán học Pháp Joseph-Louis Lagrange mới đưa ra bản chứng minh hoàn chỉnh đầu tiên khẳng định phương trình $x^2 – d y^2 = 1$ luôn luôn có vô số nghiệm nguyên với mọi số nguyên không chính phương $d > 0$, bằng cách sử dụng Phân số liên tục (Continued Fractions).
Căn bậc hai của một số không chính phương $d$ luôn biểu diễn được dưới dạng phân số liên tục chu kỳ:
Các phân số tiệm cận $\frac{p_n}{q_n}$ của phân số liên tục này chính là những mảnh ghép cung cấp chính xác nghiệm nguyên cho phương trình Pell!
Hình 1: Cấu trúc bề mặt hyperbol $x^2 – 2y^2 = 1$. Các điểm màu đỏ đại diện cho các nghiệm nguyên phân bố rời rạc trên nhánh hyperbol.
Minh Họa Tính Toán Từng Bước với $d = 2$ và $d = 3$
Để thấy rõ cơ chế hoạt động kỳ diệu của Phân số liên tục và Đồng nhất thức Brahmagupta, hãy cùng thực hiện các bước tính toán số học cụ thể với $d = 2$:
Bài toán: Tìm các nghiệm nguyên của $x^2 – 2 y^2 = 1$
Bước 1: Khai triển Phân số liên tục của $\sqrt{2}$ Ta có $\sqrt{2} \approx 1.41421356\dots$ Khai triển phân số liên tục chu kỳ của $\sqrt{2}$ là: $$\sqrt{2} = [1; \overline{2}] = 1 + \frac{1}{2 + \frac{1}{2 + \frac{1}{2 + \dots}}}$$
Bước 2: Tính các Phân số Tiệm cận $\frac{p_n}{q_n}$
Tiệm cận $n=0$: $$\frac{p_0}{q_0} = 1 = \frac{1}{1} \implies 1^2 – 2(1)^2 = -1 \quad \text{(Nghiệm của phương trình Pell âm)}$$
Tiệm cận $n=1$: $$\frac{p_1}{q_1} = 1 + \frac{1}{2} = \frac{3}{2} \implies 3^2 – 2(2)^2 = 9 – 8 = 1 \quad \text{(Thành công!)}$$ Ta thu được Nghiệm Cơ Bản Nhỏ Nhất: $(x_1, y_1) = (3, 2)$.
Bước 3: Sử dụng Đồng nhất thức Brahmagupta để nhân nghiệm Áp dụng Phép hợp thành Samasa cho hai nghiệm $(3, 2)$ và $(3, 2)$: $$x_2 = x_1 x_1 + d y_1 y_1 = 3 \times 3 + 2 \times (2 \times 2) = 9 + 8 = 17$$ $$y_2 = x_1 y_1 + x_2 y_1 = 3 \times 2 + 3 \times 2 = 6 + 6 = 12$$ Thử lại: $17^2 – 2(12)^2 = 289 – 2(144) = 289 – 288 = 1$! Ta có nghiệm thứ hai: $(17, 12)$.
Bước 4: Sinh ra Nghiệm thứ ba $(x_3, y_3)$ Nhân nghiệm $(17, 12)$ với nghiệm cơ bản $(3, 2)$: $$x_3 = 17 \times 3 + 2 \times (12 \times 2) = 51 + 48 = 99$$ $$y_3 = 17 \times 2 + 12 \times 3 = 34 + 36 = 70$$ Thử lại: $99^2 – 2(70)^2 = 9801 – 2(4900) = 9801 – 9800 = 1$! Ta có nghiệm thứ ba: $(99, 70)$.
Nhờ cấu trúc nhóm này, từ duy nhất một nghiệm cơ bản nhỏ $(3,2)$, ta có thể nhân lên vô tận để tạo ra chuỗi nghiệm $(1,0), (3,2), (17,12), (99,70), (577,408), \dots$ chạy ra vô cực.
Cơn Bùng Nổ Chữ Số: Từ $d = 61$ Đến 206.545 Chữ Số
Tuy nhiên, sự êm đềm của $d = 2$ chỉ là một ngoại lệ hiếm hoi. Khi hệ số $d$ tăng lên, độ dài của nghiệm cơ bản không tăng theo cấp số cộng hay cấp số nhân thông thường, mà bùng nổ theo cấp số nhân tầng.
Vào năm 1657, Pierre de Fermat đã gửi một thách đố tới các nhà toán học châu Âu, yêu cầu giải phương trình Pell với $d = 61$:
$$x^2 – 61 y^2 = 1$$
Trực giác ngây thơ khiến nhiều người tin rằng nghiệm sẽ chỉ nằm ở hàng chục hay hàng trăm. Nhưng kết quả thực tế đã khiến tất cả bàng hoàng. Nghiệm nguyên nhỏ nhất cho $d = 61$ đòi hỏi những con số dài tới 10 và 11 chữ số:
Sự bùng nổ này đạt đến đỉnh cao kinh hoàng trong Bài toán Đàn bò của Archimedes. Khi quy đổi các điều kiện của Archimedes về phương trình Pell, ta thu được phương trình:
$$x^2 – 410,286,423,278,424 , y^2 = 1$$
Vào năm 1880, nhà toán học Đức A. Amthor mới lần đầu tiên tính toán thành công dung lượng nghiệm cho bài toán Archimedes. Nghiệm nhỏ nhất đòi hỏi số lượng bò là:
$$\text{Tổng số bò} \approx 7.760271 \times 10^{206,544}$$
Số lượng bò này chứa tới 206.545 chữ số — lớn hơn tổng số nguyên tử trong toàn bộ vũ trụ quan sát được (vốn chỉ khoảng $10^{80}$)!
Hình 2: Sự bùng nổ độ dài nghiệm nguyên. Nghiệm của bài toán Archimedes vượt xa tổng số nguyên tử trong vũ trụ
Các vấn đề mở
Phương trình Pell không dừng lại ở một câu đố lịch sử cổ đại, mà là cánh cửa dẫn trực tiếp vào trung tâm của Lý thuyết Số học Đại số hiện đại (Algebraic Number Theory).
Đặc biệt, vào tháng 8 năm 2022, hai nhà toán học Peter Koymans và Carlo Pagano đã tạo nên một đột phá chấn động khi chứng minh thành công Giả thuyết Stevenhagen (1993) về Phương trình Pell Âm ($x^2 – d y^2 = -1$), giải quyết trọn vẹn câu hỏi: Với tần suất chính xác bao nhiêu % thì phương trình Pell âm có nghiệm nguyên?
Dưới đây là các kết quả tổng hợp liên quan đến Phương trình Pell:
Chủ đề
Hiện trạng
Trạng thái Phương trình Dương ($x^2 – dy^2 = 1$)
Đã giải quyết triệt để 100%. Lagrange (1768) chứng minh luôn luôn có vô số nghiệm nguyên cho mọi $d$ không chính phương.
Đột phá Phương trình Âm ($x^2 – dy^2 = -1$)
Đã chứng minh Giả thuyết Stevenhagen (1993) bởi Peter Koymans & Carlo Pagano (Quanta Magazine / arXiv:2201.13424). Tần suất có nghiệm của các số $d$ thỏa điều kiện là chính xác $\approx 58.0578%$.
Công cụ Toán học cốt lõi
Nhóm lớp (Class Groups), Định lý Đơn vị Dirichlet, Regulator $R_d = \ln(\epsilon_d)$, Luật Tương hỗ Rédei, Bất biến Alexander Smith (2017/2022).
Ứng dụng Mật mã học Hậu Lượng tử
Nhóm Lớp Lý Tưởng (Ideal Class Groups) áp dụng trong thuật toán mật mã CSIDH (Mật mã đồng xạ hậu lượng tử) và Hàm Tạo Độ Trễ Kiểm Chứng (VDFs) trong công nghệ Blockchain (Chia, Ethereum).
Kiểm chứng Tự động bằng AI
Chuẩn hóa định lý Lagrange, Phân số liên tục và Nhóm Đơn vị trong thư viện Mathlib của ngôn ngữ kiểm chứng Lean 4.
Từ thách thức 44 câu thơ của Archimedes đến những thuật toán mật mã hậu lượng tử của thế kỷ 21, Phương trình Pell là minh chứng rực rỡ cho thấy một cấu trúc toán học cổ đại có thể tiếp tục tỏa sáng và định hình tương lai công nghệ của nhân loại.
Comments
So empty here ... leave a comment!