MathLabs

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

Phương pháp vòng tròn

Kỹ thuật giải tích dùng tích phân Fourier trên đường tròn đơn vị để đếm nghiệm của các phương trình cộng tính.

Trực giácTrực giác: đếm bằng cách đi vòng quanh một đường tròn

Giả sử bạn muốn đếm số cách viết một số nguyên nn thành tổng của kk phần tử lấy từ một tập AA các số nguyên (chẳng hạn các số nguyên tố, hoặc các luỹ thừa bậc kk đúng). Mã hoá AA bằng một "sóng" chạy đúng một vòng quanh đường tròn đơn vị: với mỗi số thực α\alpha trong [0,1)[0,1), lập tổng mũ FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha), trong đó e(x):=e2πixe(x) := e^{2\pi i x}. Khi α\alpha quét từ 00 đến 11, FA(α)kF_A(\alpha)^k dao động; gần một số ít điểm "cộng hưởng" — các số hữu tỉ a/qa/q với mẫu số qq nhỏ — các số hạng của tổng thẳng hàng và cộng dồn (gọi là cung lớn), trong khi ở hầu hết những chỗ còn lại các pha chỉ theo hướng gần như ngẫu nhiên và triệt tiêu lẫn nhau (cung nhỏ). Phương pháp vòng tròn biến bức tranh hình học này thành một công thức chính xác rồi ước lượng công thức đó trên từng cung.

Đường tròn đơn vị tương tác cho thấy một điểm ở góc theta, dùng để dựng tổng mũ của phương pháp vòng tròn.
Kéo θ\theta sẽ di chuyển điểm e2πiθ/360e^{2\pi i\theta/360} quanh đường tròn đơn vị; phương pháp vòng tròn lấy tích phân các điểm như vậy nhân với FA(α)kF_A(\alpha)^k khi θ/360=α\theta/360=\alpha chạy khắp [0,1)[0,1).

Đại họcĐịnh nghĩa: hàm sinh trên đường tròn

Định nghĩa: Tổng mũ và số cách biểu diễn

Với tập hữu hạn A⊂Z≥0A\subset\mathbb Z_{\ge0}, đặt FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha) với e(x):=e2πixe(x) := e^{2\pi i x}. Với n≥0n\ge0 và k≥1k\ge1, gọi rk(n)=#{(a1,…,ak)∈Ak:a1+⋯+ak=n}r_k(n)=\#\{(a_1,\dots,a_k)\in A^k : a_1+\cdots+a_k=n\} là số cách biểu diễn có thứ tự của nn thành tổng kk phần tử của AA.

FA(α)=∑a∈Ae(aα),e(x):=e2πixF_A(\alpha) = \sum_{a \in A} e(a\alpha), \qquad e(x) := e^{2\pi i x}

Nâng FAF_A lên luỹ thừa kk rồi khai triển, hệ số của e(nα)e(n\alpha) "lẽ ra" phải đúng bằng rk(n)r_k(n); định danh trích xuất dưới đây (Định lý 1) làm rõ điều này: rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha. Ở đây FA(α)kF_A(\alpha)^k chính là sóng trong hình trên, và tích phân sẽ lọc ra đúng một tần số, nn, mà ta quan tâm — đó là toàn bộ nội dung của phương pháp: thay một bài toán đếm bằng một bài toán ước lượng tích phân.

rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha
Cung lớn so với cung nhỏ
Đặc điểmCung lớn (gần a/qa/q, qq nhỏ)Cung nhỏ (phần còn lại)
Vị tríMột khoảng ngắn quanh mỗi số hữu tỉ a/qa/q với q≤Qq\le QPhần còn lại của [0,1)[0,1) sau khi bỏ hết các cung lớn
Độ lớn của FA(α)F_A(\alpha)Gần giá trị lớn nhất tầm thường ∣A∣|A|: các số hạng cộng dồnĐược kỳ vọng nhỏ hơn nhiều so với ∣A∣|A|: các số hạng triệt tiêu
Vai trò trong ước lượngCho ra số hạng chính (một "chuỗi kỳ dị")Phải được chặn trên và gộp vào số hạng sai số

Đại họcHai định lý nền tảng

Với mọi tập hữu hạn A⊂Z≥0A\subset\mathbb Z_{\ge0}, số nguyên k≥1k\ge1 và n≥0n\ge0, rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha, trong đó rk(n)r_k(n) là số bộ kk phần tử có thứ tự lấy từ AA có tổng bằng nn.

Vì sao đúng?

Điều này cho phép ta thay một bài toán đếm tổ hợp thuần tuý bằng câu hỏi về độ lớn của một tích phân giải tích: nếu chứng minh được tích phân dương, thì chắc chắn tồn tại một cách biểu diễn.

Chứng minh

Trước hết, dùng hệ thức trực giao: ∫01e(mα) dα\int_0^1 e(m\alpha)\, d\alpha=1=1 nếu m=0m=0, và =0=0 nếu m≠0m\neq 0: viết e(mα)=cos⁡(2πmα)+isin⁡(2πmα)e(m\alpha)=\cos(2\pi m\alpha)+i\sin(2\pi m\alpha), khi m≠0m\neq0 đây là trọn vẹn một số nguyên chu kỳ của sóng sin/cos trên [0,1)[0,1) nên tích phân bằng 00; khi m=0m=0 hàm dưới dấu tích phân là hằng số 11.

Bây giờ khai triển luỹ thừa kk của hàm sinh: FA(α)k=(∑a∈Ae(aα))k=∑(a1,…,ak)∈Ake((a1+⋯+ak)α)F_A(\alpha)^k=\Big(\sum_{a\in A}e(a\alpha)\Big)^k=\sum_{(a_1,\dots,a_k)\in A^k} e\big((a_1+\cdots+a_k)\alpha\big), một tổng hữu hạn trên mọi bộ kk phần tử có thứ tự lấy từ AA, nhóm theo số mũ m=a1+⋯+akm=a_1+\cdots+a_k.

Nhân cả hai vế với e(−nα)e(-n\alpha) rồi lấy tích phân từng số hạng trên [0,1)[0,1) (hợp lệ vì tổng hữu hạn nên tích phân và tổng giao hoán được): ∫01FA(α)ke(−nα) dα=∑(a1,…,ak)∈Ak∫01e((a1+⋯+ak−n)α) dα\int_0^1 F_A(\alpha)^k e(-n\alpha)\,d\alpha=\sum_{(a_1,\dots,a_k)\in A^k}\int_0^1 e\big((a_1+\cdots+a_k-n)\alpha\big)\,d\alpha.

Theo hệ thức trực giao, mỗi số hạng bên phải bằng 11 đúng khi a1+⋯+ak=na_1+\cdots+a_k=n và bằng 00 trong mọi trường hợp khác. Vậy toàn bộ tổng suy biến về đúng số bộ có a1+⋯+ak=na_1+\cdots+a_k=n, tức là rk(n)r_k(n). Điều này chứng minh định danh một cách chính xác tuyệt đối, không hề xấp xỉ: đây là một đẳng thức, không phải một công thức tiệm cận.

Với mọi số thực α\alpha và mọi số nguyên Q≥1Q\ge1, tồn tại số hữu tỉ a/qa/q với 1≤q≤Q1\le q\le Q, gcd⁡(a,q)=1\gcd(a,q)=1, sao cho ∣α−aq∣≤1q(Q+1)\left|\alpha - \dfrac{a}{q}\right| \le \dfrac{1}{q(Q+1)}.

Vì sao đúng?

Điều này bảo đảm mọi điểm trên đường tròn đều nằm gần một số hữu tỉ có mẫu số nhỏ, và đó chính là điều cho phép ta chia [0,1)[0,1) thành các cung lớn — những khoảng ngắn quanh các xấp xỉ tốt a/qa/q với qq nhỏ, nơi FAF_A lớn — và các cung nhỏ, phần còn lại.

Chứng minh

Xét Q+2Q+2 số 0,{α},{2α},…,{Qα},10,\{\alpha\},\{2\alpha\},\dots,\{Q\alpha\},1, trong đó {x}\{x\} là phần thập phân của xx; tất cả đều nằm trong [0,1][0,1]. Chia [0,1][0,1] thành Q+1Q+1 khoảng con bằng nhau, mỗi khoảng dài 1/(Q+1)1/(Q+1): [0,1Q+1),[1Q+1,2Q+1),…[0,\tfrac1{Q+1}), [\tfrac1{Q+1},\tfrac2{Q+1}),\dots.

Ta có Q+2Q+2 số nhưng chỉ Q+1Q+1 khoảng con, nên theo nguyên lý chuồng bồ câu, hai trong số các số {jα}\{j\alpha\} và {iα}\{i\alpha\} (với 0≤i<j≤Q0\le i<j\le Q, cho phép i=0i=0 nên {iα}=0\{i\alpha\}=0) rơi vào cùng một khoảng con, do đó sai khác nhau ít hơn 1/(Q+1)1/(Q+1): ∣{jα}−{iα}∣<1Q+1|\{j\alpha\}-\{i\alpha\}|<\tfrac{1}{Q+1}.

Đặt q=j−iq=j-i, vậy 1≤q≤Q1\le q\le Q. Vì {jα}−{iα}=(jα−iα)−(⌊jα⌋−⌊iα⌋)=qα−a\{j\alpha\}-\{i\alpha\} = (j\alpha - i\alpha) - (\lfloor j\alpha\rfloor - \lfloor i\alpha\rfloor) = q\alpha - a với số nguyên a=⌊jα⌋−⌊iα⌋a=\lfloor j\alpha\rfloor-\lfloor i\alpha\rfloor, ta được ∣qα−a∣<1Q+1|q\alpha - a| < \tfrac{1}{Q+1}, tức là ∣α−aq∣<1q(Q+1)\left|\alpha-\dfrac aq\right| < \dfrac{1}{q(Q+1)}.

Cuối cùng, nếu gcd⁡(a,q)=d>1\gcd(a,q)=d>1, chia cả aa và qq cho dd: phân số thu được có mẫu số còn nhỏ hơn và khoảng cách tới α\alpha bằng hoặc nhỏ hơn, nên ta có thể giả sử gcd⁡(a,q)=1\gcd(a,q)=1 mà không mất tính tổng quát. Chứng minh hoàn tất; đây là lập luận chuồng bồ câu hữu hạn, mang tính dựng, không cần dựa vào bất kỳ ước lượng chưa được chứng minh nào.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Ngoài lý thuyết số, chính ý tưởng "cộng dồn các sóng và tìm cộng hưởng" cũng là nguyên lý hoạt động của biến đổi Fourier rời rạc dùng khắp trong xử lý tín hiệu và kỹ thuật điện, và công thức tiệm cận p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) mà Hardy và Ramanujan rút ra bằng một phiên bản sơ khai của phương pháp vòng tròn (áp dụng cho hàm phân hoạch p(n)p(n)) được dùng trong vật lý thống kê để ước lượng số vi trạng thái — do đó là entropy — của một hệ các kích thích boson không phân biệt được ở tổng năng lượng cố định nn.

Ví dụ: Một tính toán thu nhỏ của phương pháp vòng tròn

Cho A={1,2,…,8}A=\{1,2,\dots,8\}. Dùng định danh trích xuất, tính r2(9)r_2(9), số cặp có thứ tự (a,b)∈A2(a,b)\in A^2 với a+b=9a+b=9.

Lời giải

Theo Định lý 1, r2(9)=∫01FA(α)2e(−9α) dαr_2(9)=\int_0^1 F_A(\alpha)^2 e(-9\alpha)\,d\alpha với FA(α)=∑a=18e(aα)F_A(\alpha)=\sum_{a=1}^{8}e(a\alpha); không cần thực sự tính tích phân bằng giải tích — định danh đã bảo đảm nó bằng đúng số đếm trực tiếp, nên ta có thể tính số đếm bằng tổ hợp rồi tin vào định danh.

Liệt kê mọi cặp có thứ tự (a,b)(a,b) với a,b∈{1,…,8}a,b\in\{1,\dots,8\} và a+b=9a+b=9: (1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1)(1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1). Mỗi toạ độ đầu từ 11 đến 88 xác định duy nhất b=9−ab=9-a, và bb luôn rơi lại vào {1,…,8}\{1,\dots,8\} (vì 1≤a≤8⇒1≤9−a≤81\le a\le8 \Rightarrow 1\le 9-a\le8), nên cả 88 giá trị của aa đều hợp lệ.

Vậy r2(9)r_2(9) bằng 88. Ví dụ nhỏ này chính là định danh của Định lý 1 đang hoạt động: tích phân giải tích và số đếm tổ hợp, theo cách dựng, là một và cùng một con số — nội dung thực sự của phương pháp chỉ xuất hiện khi AA trở thành một tập vô hạn hoặc đang lớn dần (như tập các số nguyên tố) và ta phải ước lượng thay vì liệt kê.

Ví dụ: Xấp xỉ hữu tỉ tốt nhất qua định lý Dirichlet

Áp dụng định lý xấp xỉ Dirichlet với α=2\alpha=\sqrt2 và Q=5Q=5 để tìm một phân số a/qa/q, 1≤q≤51\le q\le5, thoả cận của Định lý 2, rồi kiểm tra bằng số.

Lời giải

Ta muốn a/qa/q với 1≤q≤51\le q\le5 và ∣2−a/q∣<1/(q⋅6)|\sqrt2-a/q|<1/(q\cdot6) (lấy Q+1=6Q+1=6 trong cận của định lý). Thử q=5q=5: số nguyên gần 52≈7.07115\sqrt2\approx7.0711 nhất là a=7a=7, cho 7/5=1.47/5=1.4.

Kiểm tra cận: ∣2−7/5∣=∣1.41421…−1.4∣≈0.01421|\sqrt2-7/5|=|1.41421\ldots-1.4|\approx0.01421, trong khi định lý hứa hẹn 1/(5⋅6)=1/30≈0.03331/(5\cdot6)=1/30\approx0.0333; quả thật 0.01421<0.03330.01421<0.0333, nên cận được thoả một cách thoải mái, đúng như bảo đảm.

Phân số 7/57/5 này thực chất là giản phân của liên phân số 2=[1;2,2,2,… ]\sqrt2=[1;2,2,2,\dots] sau hai bước, đó là lý do nó xấp xỉ tốt đến vậy — chứng minh chuồng bồ câu của Dirichlet không cần liên phân số để hoạt động, nhưng luôn cho ra các xấp xỉ tốt tương đương. Vậy 75\dfrac{7}{5} là một chứng nhân hợp lệ.

Trong định danh trích xuất rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha, vì sao tích phân này bằng đúng rk(n)r_k(n)?

Với A={1,2,…,8}A=\{1,2,\dots,8\}, r2(9)r_2(9) — số cặp có thứ tự (a,b)∈A2(a,b)\in A^2 với a+b=9a+b=9 — bằng bao nhiêu?

Theo định lý xấp xỉ Dirichlet, với mọi số thực α\alpha và số nguyên Q≥1Q\ge1, điều gì được bảo đảm tồn tại?

Trong vật lý thống kê, công thức tiệm cận Hardy–Ramanujan p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) cho hàm phân hoạch giúp ước lượng điều gì?

Tài liệu tham khảo

  1. G. H. Hardy, S. Ramanujan (1918). Asymptotic formulae in combinatory analysis · DOI:10.1112/plms/s2-17.1.75
  2. J. Bourgain, C. Demeter, L. Guth (2016). Proof of the main conjecture in Vinogradov's Mean Value Theorem for degrees higher than three · DOI:10.4007/annals.2016.184.2.7 · arXiv:1512.01565
  3. B. Green, T. Tao (2008). The primes contain arbitrarily long arithmetic progressions · DOI:10.4007/annals.2008.167.481 · arXiv:math/0404188