主題
Search

一百美元,百位數字挑戰難題


一百美元,百位數字挑戰難題是由十個數值分析問題組成,發表於 2002 年 1/2 月號的 SIAM News (http://www.siam.org/siamnews/01-02/challenge.pdf)。 這些問題的提出者 Nick Trefethen 懸賞 100 美元給在 2002 年 5 月 20 日之前,對這些問題獲得最多正確位數(最多 10 位)的個人或團隊。 Trefethen 低估瞭解題者的聰明才智,20 個獨立的團隊獲得了所有 10 個問題的 10 位正確數字。 在一位匿名捐助者的介入以幫助支付超出預期的獎金後,獎金支票於 2002 年 12 月寄給了所有獲獎者。

提出的問題如下。

HundredDollar1Plot

1. lim_(epsilon->0)int_epsilon^1x^(-1)cos(x^(-1)lnx)dx 是什麼?這個問題在其原始形式中很難進行數值積分,因為原點附近存在強烈的振盪。 然而,透過代換 u=-lnx/x,它可以轉化為積分

 h_1=int_0^infty(cosu)/(u{1+[W(u)]^(-1)})du
(1)

使用振盪數值積分技術,這個積分可以相當快地收斂。

Boersma 和 Jansen 使用圍道積分將問題轉化為

 h_1=int_0^(pi/2)sin(tsint)e^(-tcost)dt+int_1^infty(cos(1/2piy))/(y^(y+1))dy,
(2)

它們每一個都收斂得相當快。

Laurie 指出,該積分可以寫成

 h_1=R[int_Cz^(i/z-1)dz],
(3)

其中 C 是從 0 到 1 的任何圍道,使得在 C 和從 0 到 1 的線段定義的區域內沒有奇點(Wagon 2004)。 例如,C=(0,1/2+i,1) 就是這樣一條圍道。

這個問題也可以透過將被積函式寫成漸近級數來解決

 xcos[xln(x^(-1))]=sum_(k=0)^infty((-1)^k)/((2k)!)(ln^(2k)x)x^(2k+1)
(4)

並逐項積分以獲得強振盪和

 h_1=sum_(k=1)^infty(-1)^(k+1)(2k)^(2k-1).
(5)

令人驚訝的是,可以使用 Wynn 的 epsilon 方法,使用適當的項數和外推度來正則化這個和,從而獲得大約 7 位正確數字。

2. 一個光子在 x-x-y 平面以速度 1 運動,在 t=0 時從 (x,y)=(0.5,0.1) 出發,方向正東。 在平面上的每個整數格點 (i,j) 周圍,都豎立了一個半徑為 1/3 的圓形鏡子。 光子在 t=10 時離原點有多遠?

3. 具有條目 a_(11)=1a_(12)=1/2a_(21)=1/3a_(13)=1/4a_(22)=1/5a_(31)=1/6 等的無限矩陣 Al^2 上的有界運算元。 ||A|| 是什麼? 這個問題等價於找到具有條目的無限矩陣的最大奇異值

 a_(ij)=2/(2-i+i^2-3j+2ij+j^2),
(6)

即,矩陣

 A=[1 1/2 1/4 ...; 1/3 1/5 ... ...; 1/6 ... ... ...; ... ... ... ...; | | | ...].
(7)
HundredDollar4Plot

4. 函式的全域性最小值是什麼

 exp(sin(50x))+sin(60e^y)+sin(70sinx)+sin(sin(80y))-sin(10(x+y))+1/4(x^2+y^2)?
(8)

(參見 Bailey 等人,2007 年,第 12 頁和 219 頁;Kampas 和 Pintér 2006 年)。 這個問題可以用 Wolfram 語言中的一行程式碼解決。

  NMinimize[f[x, y], {x, y}, Method -> {"RandomSearch",
    "SearchPoints" -> 700}, WorkingPrecision -> 20]

5. 令 f(z)=1/Gamma(z),其中 Gamma(z)伽馬函式,令 p(z) 是在單位圓盤上以上確界範數 ||·||_infty 最佳逼近 f(z)三次多項式||f-p||_infty 是什麼?

6. 一隻跳蚤從無限二維整數格點上的 (0,0) 開始,進行有偏隨機遊走:每一步,它以 1/4 的機率向北或向南跳,以 1/4+epsilon 的機率向東跳,以 1/4-epsilon 的機率向西跳。 跳蚤在其遊走過程中返回 (0, 0) 的機率為 1/2。 epsilon 是什麼?

解由求解下式給出

 (pisqrt(2+16epsilon^2-2sqrt(1-16epsilon^2)))/(K(sqrt((2sqrt(1-16epsilon^2))/(sqrt(1-16epsilon^2)-1-8epsilon^2))))=2,
(9)

其中 K(k)第一類完全橢圓積分。 等價地,它由以下方程的根給出

 sqrt(2)AGM(sqrt(1+8epsilon^2-sqrt(1-16epsilon^2)),sqrt(1+8epsilon^2+sqrt(1-16epsilon^2))=1,
(10)

其中 AGM(a,b)算術-幾何平均值

它也可以透過求解 u(epsilon)=2 給出,其中

u(epsilon)=2/piint_0^pi(dphi)/(sqrt(3-4cosphi+cos^2phi+16epsilon^2))
(11)
=2/piint_(-1)^1(dt)/(sqrt(1-t^2)sqrt(3-4t+t^2+14epsilon^2))
(12)
=(2sqrt(2))/(pisqrt(1+8epsilon^2sqrt(1-16epsilon^2)))×K(sqrt(1-(1+8epsilon^2-sqrt(1-16epsilon^2))/(1+8epsilon^2+sqrt(1-16epsilon^2)))),
(13)

(Bornemann 2002)。

7. 令 A20000×20000 矩陣,其條目在所有位置都為零,除了主對角線上的素數 2, 3, 5, 7, ..., 224737 以及所有位置 a_(ij) 中為 1 的數字,其中 |i-j|=1, 2, 4, 8, ..., 16384。 A^(-1) 的 (1, 1) 項是什麼?

這個問題可以精確求解,得到一個有理數,其分子和分母各有 97389 位數字

 h_3=(31016407491...417983612357075)/(42776629106...013006012935182)
(14)

(Wagon 2004)。

8. 一個正方形板 [-1,1]×[-1,1]溫度u=0。 在時間 t=0 時,溫度沿四個邊中的一個邊增加到 u=5,同時沿其他三個邊保持在 u=0,然後熱量根據 u_t=Deltau 流入板中。 板中心處的溫度何時達到 u=1

給出 10 位正確數字的表示式由下式給出

 h_8=-1/(pi^2)ln(-81pi^4+518400x-691200x^3+345600x^5 
 -76800x^7+6400x^9)_1,
(15)

其中 (P(x))_n 是一個多項式根。 獲得更高的精度需要使用級數的更多項。

HundredDollar9Plot

9. 積分 I(alpha)=int_0^2[2+sin(10alpha)]x^alphasin(alpha/(2-x))dx 取決於引數 alpha。 在 alpha in [0,5] 中的哪個 I(alpha) 值處,I(alpha) 達到其最大值?

透過進行變數替換 u=1/(2-x),積分可以轉化為

 I(alpha)=4sqrt(pi)Gamma(alpha)G_(2,4)^(3,0)((alpha^2)/(16)|(alpha+2)/2,(alpha+3)/2; 1/2,1/2,1,0)[sin(10alpha)+2].
(16)

透過繪圖,可以看到最大值出現在 alpha=0.78 附近,並且可以使用標準的求根技術來高精度地確定它。

10. 一個位於 10×1 矩形中心的粒子進行布朗運動(即,步長無窮小的二維隨機遊走),直到它到達邊界。 它在端點而不是側面擊中的機率是多少? 令人驚訝的是,這個問題有一個閉式解,由下式給出

h_(10)=4/pisum_(k=0)^(infty)((-1)^k)/(2k+1)sech[5pi(2k+1)]
(17)
=8/pisum_(k=0)^(infty)(-1)^ktan^(-1)[e^(-5pi(2k-1))]
(18)
=2/picos^(-1)sqrt(lambda(1/(10)i))
(19)
=2/pisin^(-1)[(3-2sqrt(2))^2(2+sqrt(5))^2×(sqrt(10)-3)^2(5^(1/4)-sqrt(2))^4],
(20)

其中 lambda(tau)橢圓 lambda 函式(參見 Bailey 等人,2007 年,第 48 頁)。

解決方案總結在下表中。

#OEISh_n
1.A1172310.3233674316
2.A1172320.9952629194
3.A1172331.274224152
4.A117234-3.306868647
5.A1172350.2143352345
6.A1172360.06191395447
7.A1172370.7250783462
8.A1172380.4240113870
9.A1172390.7859336743
10.A1172403.837587979×10^(-7)

使用 探索

參考文獻

Bailey, D. H. 和 Borwein, J. M. “實驗數學的示例問題。” 2003 年 9 月 22 日。 http://crd.lbl.gov/~dhbailey/expmath/expmath-probs.pdf.Bailey, D. H.; Borwein, J. M.; Calkin, N. J.; Girgensohn, R.; Luke, D. R.; 和 Moll, V. H. Experimental Mathematics in Action. Wellesley, MA: A K Peters, 2007.Beard, B. B.; Medley, B.; 和 van Gans, M. "The 2002 SIAM Challenge." http://www.maxwellian.demon.co.uk/~marijke/SIAM2002/.Boersma, J.; Jansen, J.; Simons, S.; 和 Steutel, F. "The SIAM 100-Dollar 100-Digit Challenge." http://www.win.tue.nl/scg/siamcontest/.Bornemann, F. "Short Remarks on the Solution of Trefethen's Hundred-Digit Challenge." 2002 年 11 月 5 日。 http://www-m3.ma.tum.de/m3old/ftp/Bornemann/pdf/short.pdf.Bornemann, F.; Lauire, D.; Wagon, S.; 和 Waldvogel, J. The SIAM 100-Digit Challenge: A Study in High-Accuracy Numerical Computing. Philadelphia, PA: SIAM, 2004. 附加材料可在 http://www-m8.ma.tum.de/m3/bornemann/challengebook/ 獲取。Borwein, J. M. "The 100 Digit Challenge: An Extended Review." Math. Intelligencer 27, 40-48, 2005.Borwein, J. 和 Bailey, D. Mathematics by Experiment: Plausible Reasoning in the 21st Century. Wellesley, MA: A K Peters, pp. 22-24, 2003.Briggs, K. "Hundred-Dollar, Hundred-Digit Challenge." http://keithbriggs.info/solutions.html.Kampas, F. J. 和 Pintér, J. D. "Configuration Analysis and Design Using Optimization Tools in Mathematica." Mathematica J. 10, 128-154, 2006.Kern, M. "Solution to the SIAM 'Hundred-Dollar, Hundred-Digit Challenge'." 報告。 2002 年 5 月。 http://www-rocq.inria.fr/~kern/Challenge/RR-challenge.pdf.Laurie, D. "Trefethen Challenge Problems." http://dip.sun.ac.za/~laurie/trefethen-challenge/.Leslie, M. (Ed.). "NetWatch: Decimal Decathlon." Science 295, 1431, 2002.Sloane, N. J. A. 序列 A117231, A117232, A117233, A117234, A117235, A117236, A117237, A117238, A117239, 和 A117240 在 “整數序列線上百科全書” 中。Trefethen, N. "A Hundred-Dollar, Hundred-Digit Challenge." SIAM News 35, No. 1, 2002 年 1/2 月。 http://www.siam.org/siamnews/01-02/challenge.pdf.Trefethen, N. "The SIAM 100-Dollar, 100-Digit Challenge." http://web.comlab.ox.ac.uk/oucl/work/nick.trefethen/hundred.html.Trefethen, N. L. "Chastened Challenge Sponsor: "I Misjudged." SIAM News 35, No. 6, 1-3, 2002 年 7/8 月。Trott, M. The Mathematica GuideBook for Programming. New York: Springer-Verlag, p. 109, 2004. http://www.mathematicaguidebooks.org/.Weisstein, E. W. "A Hundred-Dollar Challenge." Headline News, 2002 年 2 月 4 日。 https://mathworld.tw/news/2002-02-04/challenge/.Weisstein, E. W. "Hundred-Dollar Challenge Winners Announced." Headline News, 2002 年 5 月 25 日。 https://mathworld.tw/news/2002-05-25/challenge/.Wagon, S. "Solutions." http://stanwagon.com/wagon/Misc/Links/SIAMchallenge_lnk_2.html.Wagon, S. "The SIAM 100-Digit Challenge." Wolfram Technology Conference, Champaign IL, 2004. http://library.wolfram.com/infocenter/Conferences/5353/.

在 中被引用

一百美元,百位數字挑戰難題

請引用為

Weisstein, Eric W. “一百美元,百位數字挑戰難題。” 來自 Web 資源。 https://mathworld.tw/Hundred-DollarHundred-DigitChallengeProblems.html

主題分類