史上最贱的数学题
我正在有選擇性地拿過去Andrew和Richard Guy研究的一些立體問題做消遣。數值結果真是令人嘆為觀止。（來自MathOverflow的評論）
這就是Allan MacLeod，一位退休的數學家，在幾年前偶然發現的方程，它真的令人嘆為觀止。老實說，我也算是大風大浪見的多了，但這么精妙的丟番圖方程（注：有一個或者幾個變量的整系數方程，它們的求解僅僅在整數范圍內進行。最后這個限制使得丟番圖方程求解與實數范圍方程求解有根本的不同。也叫不定方程。）還是第一次見。
在我碰到這道題之前，它已經被某人心懷惡意地發布在網絡上，成為流行的朋友圈圖片，肆意捉弄那些老實人（Scridhar，這個人是不是你？）。我根本沒意識到我偶然看到的這道題到底是個什么樣的怪物。它長這個樣：
你可能已經在朋友圈看到過很多這樣的圖了，它們一般都是標題黨的垃圾：什么“95%的麻省理工畢業生無法解決的問題”，這個“問題”要么很空洞，要么偷換概念，要么就是不重要的腦筋急轉彎。
但這個問題不是。這張圖片就是一個精明的，或者說陰險的圈套。大概99.999995%的人根本沒有任何機會解決它，甚至包括一大批頂級大學非數論方向的數學家。它的確是可解的，但那真的真的不得了的難。
（順便說一句。發布的人實際上不是Scridhar，或者說不能怪他。）
你可能會這樣想，如果所有的嘗試都失敗了，我們還可以直接用電腦計算大力出奇跡。這年頭，寫個電腦程序解決這種形式簡單的方程真是太容易了，只要它真的有答案，那電腦最終一定會找出來。但很抱歉，大錯特錯。用電腦暴力計算在這里毫無用處。
如果不把Quora的讀者都當作橢圓曲線的入門者的話，我不知道怎么才能寫出適合的答案。我在這能做的只是一個簡要的概覽。主要參考文獻是最近Bremmer和MacLeod2014年在《數學和信息學年鑒（Annales Mathematicae etInformaticae）》上發表的一篇名為《一個不一般的立體代表性問題（An unusual cubic representationproblem）》的精彩論文。? ? ? ?
讓我們開始吧。
我們求解的是這個方程的整數解
?(為了與論文的變量名相適應，我把蘋果、香蕉和菠蘿修改過來了)
面對任何方程，你需要做的第一步是嘗試并確定問題背景。這到底被劃歸到哪一類問題？嗯，我們被要求找到整數解，所以這是一個數論問題。就題而言，方程涉及有理函數（多項式除多項式的函數形式），但很顯然我們可以用通分移項的方法化成一個多項式函數，所以我們實際上解得是一個丟番圖方程（?Diophantine equation）。正數解的要求有一點不同尋常，接下來我們會看到這個要求會讓問題變得多么難。
現在，我們有了多少變量？這個問題看起來很蠢：很明顯，我們有三個變量，分別是a、b、c。讓我們慢一點來。一個科班出身的數論學家第一眼就能察覺到，這個方程是齊次的。這意味著如果（a，b，c）是方程的一個特解的話，那（7a，7b，7c）都是它的解。你能看出為什么嗎？給每一個變量乘一個常數沒有改變方程的結構（7只是一個例子），因為分子分母全部都約掉了。
這意味著這個方程看上去像是三維的，但它實際上只有兩維。在幾何學中，它對應著一個面（一個三元方程一般定義一個兩維的面。一般來說，k個n元方程定義一個d維的流形，d=n-k）。這個面是由一條過原點的線旋轉形成的，可以通過截取的單平面來理解。這是一條投影曲線。
在大多數初等的情形，這種降維可以這么解釋：無論解是什么，我們都可以分為兩類，c=0的情形和c≠0的情形。第一類僅僅涉及兩個變量（所以自然是二維的），而第二類情形我們可以對所有解同時除以c并得到一個c=1的解（注：在上上一段，我們已經說明了這樣一組解也一定是方程的解）。因此我們可以在c=1的情況下尋找a和b的有理數解，只要乘以一個公分母，就得到了a，b，c的正數解。一般來說，齊次方程的整數解對應一個低一個維度的非齊次方程的有理數解。
接下來的問題是：這個方程的次數是什么？次數指的是各項中最高的冪次，對于涉及多個變量相乘的項，冪次就是各變量冪次之和。舉個例子，如果某項為a2 bc4 ，那此項的次數就是7=2+1+4 。
丟番圖方程在不同次數難度完全不一樣，寬泛地說：
一次的非常簡單。
二次的也被理解得非常透徹，一般能用相對初等的方法解決。
三次的就是滿山滿海的深奧理論和數不勝數的開放問題。
四次的，嗯，真的真的很難。
我們這個方程是三次的。為什么？嗯，去分母之后就很顯然了：
即使沒有合并同類項，你也可以明白地看到次數為3：沒有超過三個變量的乘積，最后我們得到的是類似a^3 、b^2 c、abc這樣的項，而沒有冪次超過3的。合并同類項后，方程整理如下：
a^3+b^3+c^3?3(a^2b+ab^2+a^2c+ac^2+b^2c+bc^2)?5abc=0
你可能會反對這樣的變形：因為這樣獲得的解可能恰好使某個分母等于0，使得原方程沒有意義。這是對的，我們的新方程的確有些解不與原方程對應。但這是好事（@F91）。這個多項式形式給原方程打上了一些補丁使得它便于處理；對于我們找到的任何特解，只需要代入原方程檢驗一下分母等不等于0就可以了。
事實上，多項式方程很容易處理。比如說，?a=?1?，b=1,?c=0。這是好事：我們有了有理數解，或者說有理點。這意味著我們的立體方程（3維）實際上是個橢圓曲線。
當你發現這個方程是橢圓曲線時，你會喜出望外，然后悲從中來（注：這里不是大家熟悉的圓錐曲線中的橢圓，而是域上虧格為1的光滑射影曲線。對于特征不等于2的域，它的仿射方程可以寫成：y^2=x^3+ax^2+bx+c。復數域上的橢圓曲線為虧格為1的黎曼面。Mordell證明了整體域上的橢圓曲線是有限生成交換群，這是著名的BSD猜想的前提條件。阿貝爾簇是橢圓曲線的高維推廣。By 百度百科。），因為你發現橢圓曲線問題是個龐然大物（學渣哇的一聲哭出來）。這個方程是一個展現橢圓曲線理論強大的經典案例，證明它可以被用來尋找一些爆難問題的解。
我們需要做的第一件事把橢圓曲線化成魏爾斯特拉斯（注：Weierstrass，提起他最著名的成就就是嚴密化微積分的ε-δ語言）形式。這是一個長得像這樣的等式：
? ? ? y2=x3+ax+b
或者有時候也會化成
? ? ? y2=x3+ax2+bx+c
(這被稱為長魏爾斯特拉斯形式。它并不是嚴格必需的，但有時候會帶來一些便利)
眾所周知，任何橢圓曲線都可以化成這種形式（在特征為2或者3的域特別基礎，如果你研究特征特別小的域，那結果就不一樣了，我們此處不作討論）。如果想講清楚怎么把橢圓曲線化成這種形式，那可就是長篇大論了（學渣的碎碎念：我信我信）。你只需要知道，這種變形是完全機械的操作（關鍵在于方程至少存在一個有理數點，而我們已經確定了一個有理數點）。現在有若干計算機函數包可以輕而易舉地幫你搞定這件事。
但即使你不知道如何完成變換，驗證它也是很容易的，或者說至少是機械的。對于我們而言，需要的變換由令人生畏的公式導出。
我知道這看上去就像隨意的巫毒把戲（注：巫毒，是目前最為人熟悉的非洲信仰，在西方文化中就是神秘力量的象征符號，可以類比國人心中的毒盅、趕尸和降頭），但請相信我它不是。一旦你完成了這些變形，沉悶但異常直白的代數計算可以證明它是對的。
y2=x3+109x2+224x
這個方程盡管看起來很原方程長得不怎么像，但確是如假包換的忠實模型。在圖像上它長成這樣，一條有著兩個實部的經典橢圓曲線：
右邊的“魚尾”連續延伸至正負無窮。左邊的封閉橢圓曲線也將給我們帶來解決問題的驚喜。給定這個方程的任意解（x，y），你都可以通過下面的等式還原所求的a，b，c：
你需要記住，三元組（a：b：c）是用投影曲線理解的——無論你從這些方程中獲得什么數值，你都可以隨意乘上一個你想要的常數。
對于我們展示的兩個圖像，無論是從a，b，c到x，y還是反過來，都可以證明這兩個方程從數論的角度是等價的：一個方程的有理數解可以導出另一個方程的有理數解。專業術語叫做雙向有理等價（birational equivalence），而這個概念在代數幾何里面是一個非常基本的。如我們之前注意到的那樣，可能存在一些不相互對應的特殊點，而情形是a+b，a+c 或者b+c恰好等于0 。這是構造雙有理等價的必要代價，而不需要對此有任何擔心。
讓我們來看看手里的這個例子。它的橢圓曲線存在一個很好的有理數點：x=?100,?y=260。可能找到這個點不太容易，但檢驗它在曲線上就很簡單了：直接代入原方程檢驗等式兩邊是否相等（我不是隨機摸的點，但各位不用關心這個問題）。我們可以簡單地驗證a，b，c代入的結果。
我們得到了a=2/7,b=?1/14,c=11/14，既然我們可以隨意乘以一個公分母，那我們就可以變形為a=4,b=?1,c=11.
代入原方程，的確4/（?1+11）?1/（4+11）+11/（4?1）=4
你可以很容易地驗證。這就是我們原方程的一個簡單整數解——但很遺憾，不是正整數解。找到這個解用手算不太容易，但用一點耐心即使不用計算機也不算太難。它將成為我們找到正數解的緣起之地。
現在，一旦你在橢圓曲線上找到了有理數點，如P（-100,260），你就可以利用弦切技巧進行加法，生成其它的有理數點（有理數的加法是封閉的，有理數加有理數還是有理數）。
圖解：橢圓曲線上點的加法
在任何情形下，在一個域（實數域R或者有理數域Q）中給定一個方程，你可以把解視為位于R2或者Q2的點（來自R2或者Q2的投影），而相加律就是弦—切結構的變形：想要對兩個點A和B做加法，構造一條過二點的直線（弦），當A，B重合時，這就是曲線的切線。找到直線與曲線的第三個交點P，對O和P重復上述操作，再次得到的交點就是A+B。當O點被選為無窮遠處的點（一般都這么處理），圖像就如上所示（注：至于O點是什么，這就涉及群論和更深奧的橢圓曲線知識，懂的自然懂，不懂的我也講不懂，因為我也不懂）。更詳細的見原作者的Quora回答previous，再詳細的請去翻代數幾何。
一開始，我們可以通過作P點的切線，找到它和曲線再次相交的點，以此增加P點的值。結果看上變得有點嚇人P+P=2P=(8836/25,?950716/125)。
同樣的，這個新的點也對應一組a，b，c的值,(a,b,c)=(9499,?8784,5165)。
這個解用手算就很困難，但用電腦就是小意思了。然而，它還不是正的。
當然，困難嚇不倒我們，我們繼續計算3P=2P+P，操作方法就是連接P和2P找到與曲線的第三個交點再與O點相連找到第四個交點。同樣的，我們計算a，b，c，然而還是同樣的，結果不是正數。以此類推，計算4P，5P等等等等。直到我們計算到9P。
9P=(-66202368404229585264842409883878874707453676645038225/13514400292716288512070907945002943352692578000406921,
58800835157308083307376751727347181330085672850296730351871748713307988700611210/1571068668597978434556364707291896268838086945430031322196754390420280407346469)
很明顯這不是人算的了，但交給機器，這也就是9次簡單的幾何程序迭代。對應的a，b，c值也很恐怖：
a=154476802108746166441951315019919837485664325669565431700026634898253202035277999,
b=36875131794129999827197811565225474825492979968971970996283137471637224634055579,
c=4373612677928697257861252602371390152816537558161613618621437993378423467772036
這些是80位數！你不可能通過暴力計算找到一個80位數（注：簡單的算術題，三個1080的數，總共的組合數就是10240，神威太湖之光的峰值計算能力為12.5億億次每秒，折算不過1018 次/s，至少需要10222秒，大約10214年，更震撼的寫法就是1億億億億億億億億億億億億億億億億億億億億億億億億億億年）！但無論它看上去怎么不可思議，但這些數值代回原方程，的確等于4：
讓我們回到理論本身再探討一下。定義在有理數上的橢圓曲線存在一個階（rank），它表示我們最開始至少需要知道多少個有理點才能通過弦切方法找到曲線上所有的有理數點。我們這條橢圓曲線的階等于1，說明雖然它上面有無窮多個有理點但都是由一個有理點生成的，而這個點不是別的恰好就是我們最開始的那個P點（-100,260）。
計算階數并找到這樣的一個生成子的算法非同尋常，但SageMath（現在叫CoCalc）只需要幾行代碼1秒鐘以內就搞定了。你可以查看我的代碼（here）,它從頭開始再現了整個解法，當然其中有Sage內置的橢圓曲線處理方法。
在我們看來，P點位于曲線的橢圓部分，而其它的mP（m為正整數）點也一樣。它們會逐漸跑遍整個橢圓并最終均勻地分布在整個曲線上。而我們是很幸運的，因為只有很少一部分橢圓能產生a，b，c的正數解：它們是下面這張圖加粗的部分（引自Bremmer和MacLeod的論文）。
P，2P等點并不在黑色加粗的部分，但9P恰好在，使我們得到一個80位的正整數解。
Bremmer和MacLeod還研究了如果我們把等式右邊的4換成其它的東西會怎么樣。如果你覺得我們的解太大了，那是因為你還沒見識到把4換成178的結果。那就不僅僅是80位了，你需要398,605,460位數。對，你沒看錯，那個解就是這么大。如果你試試896，位數就飆升到數萬億位了。沒錯，數萬億位的解，屬于這個看上去人畜無害的方程。
上述的丟番圖方程就是一個系數很小但整數解位數巨大的駭人案例。它不僅僅是令人生畏的符號，而是一項意義深遠的研究。希爾伯特第十大問題的否證陳述意味著，隨著系數逐漸增大，解的增長將變為一個不可計算的方程——因為如果它是可計算的，那我們就能得到一個解開丟番圖方程的簡單算法，而事實上并沒有，無論是簡單的還是復雜的。這項研究展現了與那個問題的某種聯系：4->80位，178->數億位，896->數萬億位，讓我們瞥見那個怪異的、不可計算的函數的一貌。稍稍把我們的方程改動一下，解就會迅速增長到蓋過我們這個“可憐的”、“渺小的”宇宙的任何事物。
何其美妙、何其揶揄的小小方程！
————
編輯?∑Pluto
來源：數學建模
更多精彩：
?泰勒定理的奇聞軼事
?丘成桐：漫談微分幾何
?Leibniz 如何想出微積分？（一）
?線性相關和秩的物理意義
?數學史上你認為最丑陋的公式是什么？
?陶哲軒談什么是好的數學
?田淵棟：數學的用處（下篇）
?你絕對沒想過原來數學家這么流氓，一言不合就進行暴力證明
?世界上最牛的五篇博士論文
?數學中有哪些巧合讓人眼前一亮？
?算法立功！清華畢業教授美國被搶車，警察無能為力自己用“貪心算法”找回
?學術史上的奇文：怎樣用數學抓獅子
?臺大教授的反思：最難的一課 我們卻沒教給學生
?麻省理工學院(MIT)研究生學習指導—— 怎樣做研究生
?分享 數學，常識和運氣 ——投資大師詹姆斯·西蒙斯2010年在MIT的講座
算法數學之美微信公眾號歡迎賜稿
稿件涉及數學、物理、算法、計算機、編程等相關領域，經采用我們將奉上稿酬。
投稿郵箱：math_alg@163.com
總結
- 上一篇： 霍金承认有鬼神，为什么霍金和牛顿这些人最
- 下一篇： 德媒看2019年的世界：西方在争吵，中国