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 n thành tổng của k phần tử lấy từ một tập A các số nguyên (chẳng hạn các số nguyên tố, hoặc các luỹ thừa bậc k đúng). Mã hoá A 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 α trong [0,1), lập tổng mũ FA(α)=∑a∈Ae(aα), trong đó e(x):=e2πix. Khi α quét từ 0 đến 1, FA(α)k dao động; gần một số ít điểm "cộng hưởng" — các số hữu tỉ a/q với mẫu số q 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 θ sẽ di chuyển điểm e2πiθ/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(α)k khi θ/360=α chạy khắp [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≥0, đặt FA(α)=∑a∈Ae(aα) với e(x):=e2πix. Với n≥0 và k≥1, gọi rk(n)=#{(a1,…,ak)∈Ak:a1+⋯+ak=n} là số cách biểu diễn có thứ tự của n thành tổng k phần tử của A.
FA(α)=a∈A∑e(aα),e(x):=e2πix
Nâng FA lên luỹ thừa k rồi khai triển, hệ số của e(nα) "lẽ ra" phải đúng bằng rk(n); định danh trích xuất dưới đây (Định lý 1) làm rõ điều này: rk(n)=∫01FA(α)ke(−nα)dα. Ở đây FA(α)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ố, n, 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(α)ke(−nα)dα
Cung lớn so với cung nhỏ
Đặc điểm
Cung lớn (gần a/q, q 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/q với q≤Q
Phần còn lại của [0,1) sau khi bỏ hết các cung lớn
Độ lớn của FA(α)
Gần giá trị lớn nhất tầm thường ∣A∣: các số hạng cộng dồn
Được kỳ vọng nhỏ hơn nhiều so với ∣A∣: các số hạng triệt tiêu
Với mọi tập hữu hạn A⊂Z≥0, số nguyên k≥1 và n≥0, rk(n)=∫01FA(α)ke(−nα)dα, trong đó rk(n) là số bộ k phần tử có thứ tự lấy từ A có tổng bằng n.
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α=1 nếu m=0, và =0 nếu m=0: viết e(mα)=cos(2πmα)+isin(2πmα), khi m=0 đây là trọn vẹn một số nguyên chu kỳ của sóng sin/cos trên [0,1) nên tích phân bằng 0; khi m=0 hàm dưới dấu tích phân là hằng số 1.
Bây giờ khai triển luỹ thừa k của hàm sinh: FA(α)k=(∑a∈Ae(aα))k=∑(a1,…,ak)∈Ake((a1+⋯+ak)α), một tổng hữu hạn trên mọi bộ k phần tử có thứ tự lấy từ A, nhóm theo số mũ m=a1+⋯+ak.
Nhân cả hai vế với e(−nα) rồi lấy tích phân từng số hạng trên [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α.
Theo hệ thức trực giao, mỗi số hạng bên phải bằng 1 đúng khi a1+⋯+ak=n và bằng 0 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=n, tức là rk(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 α và mọi số nguyên Q≥1, tồn tại số hữu tỉ a/q với 1≤q≤Q, gcd(a,q)=1, sao cho α−qa≤q(Q+1)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) thành các cung lớn — những khoảng ngắn quanh các xấp xỉ tốt a/q với q nhỏ, nơi FA lớn — và các cung nhỏ, phần còn lại.
Chứng minh
Xét Q+2 số 0,{α},{2α},…,{Qα},1, trong đó {x} là phần thập phân của x; tất cả đều nằm trong [0,1]. Chia [0,1] thành Q+1 khoảng con bằng nhau, mỗi khoảng dài 1/(Q+1): [0,Q+11),[Q+11,Q+12),….
Ta có Q+2 số nhưng chỉ Q+1 khoảng con, nên theo nguyên lý chuồng bồ câu, hai trong số các số {jα} và {iα} (với 0≤i<j≤Q, cho phép i=0 nên {iα}=0) rơi vào cùng một khoảng con, do đó sai khác nhau ít hơn 1/(Q+1): ∣{jα}−{iα}∣<Q+11.
Đặt q=j−i, vậy 1≤q≤Q. Vì {jα}−{iα}=(jα−iα)−(⌊jα⌋−⌊iα⌋)=qα−a với số nguyên a=⌊jα⌋−⌊iα⌋, ta được ∣qα−a∣<Q+11, tức là α−qa<q(Q+1)1.
Cuối cùng, nếu gcd(a,q)=d>1, chia cả a và q cho d: phân số thu được có mẫu số còn nhỏ hơn và khoảng cách tới α bằng hoặc nhỏ hơn, nên ta có thể giả sử 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)∼4n31exp(π32n) 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)) đượ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 n.
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}. Dùng định danh trích xuất, tính r2(9), số cặp có thứ tự (a,b)∈A2 với a+b=9.
Lời giải
Theo Định lý 1, r2(9)=∫01FA(α)2e(−9α)dα với FA(α)=∑a=18e(aα); 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) với a,b∈{1,…,8} và a+b=9: (1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1). Mỗi toạ độ đầu từ 1 đến 8 xác định duy nhất b=9−a, và b luôn rơi lại vào {1,…,8} (vì 1≤a≤8⇒1≤9−a≤8), nên cả 8 giá trị của a đều hợp lệ.
Vậy r2(9) bằng 8. 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 A 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 và Q=5 để tìm một phân số a/q, 1≤q≤5, thoả cận của Định lý 2, rồi kiểm tra bằng số.
Lời giải
Ta muốn a/q với 1≤q≤5 và ∣2−a/q∣<1/(q⋅6) (lấy Q+1=6 trong cận của định lý). Thử q=5: số nguyên gần 52≈7.0711 nhất là a=7, cho 7/5=1.4.
Kiểm tra cận: ∣2−7/5∣=∣1.41421…−1.4∣≈0.01421, trong khi định lý hứa hẹn 1/(5⋅6)=1/30≈0.0333; quả thật 0.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/5 này thực chất là giản phân của liên phân số 2=[1;2,2,2,…] 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 57 là một chứng nhân hợp lệ.
Trong định danh trích xuất rk(n)=∫01FA(α)ke(−nα)dα, vì sao tích phân này bằng đúng rk(n)?
Với A={1,2,…,8}, r2(9) — số cặp có thứ tự (a,b)∈A2 với a+b=9 — bằng bao nhiêu?
Theo định lý xấp xỉ Dirichlet, với mọi số thực α và số nguyên Q≥1, đ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)∼4n31exp(π32n) cho hàm phân hoạch giúp ước lượng điều gì?