第53章 靈感
這種時刻劉浩然見過很多次，數學家面對難題時的那種執拗，那種不解決問題誓不罷休的勁頭。
他沒再勸，只是又熱了一杯牛奶放在了肖宿桌上。
凌晨一點，肖宿突然放下粉筆，走到窗前。
窗外是寂靜的校園，只有路燈還孤獨地亮著。
他盯著黑暗中的某一點，眼神空洞，像是在看很遠的地方。
劉浩然知道，這是肖宿在「放空」。
當思考進入死胡同時，他會暫時跳出數學的世界，讓潛意識去工作。
時間一分一秒過去。
辦公室里的鐘表滴答聲格外清晰。
凌晨兩點十七分，肖宿突然轉過身，眼睛裡閃著光。
「我明白了。」
「明白什麼？」
「我們一直在用錯誤的範式思考。」
肖宿走回白板前，語速快而清晰。
「現有的隨機微分方程數值解法，都是基於『離散化時間，逐步推進』的思路。但在流形上，這種思路天然受限，因為每一步推進都要保證不離開流形。」
他擦掉白板上一半的內容，開始重新書寫：
「正確的思路應該是：不是離散化時間，而是離散化路徑空間。」
「什麼？」劉浩然沒聽懂。
「考慮從初始點到終點的所有可能路徑組成的空間，這叫做路徑空間。」
肖宿邊寫邊解釋，「在流形上，路徑空間本身也是一個無窮維流形。我們要找的解，就是這個無窮維流形上的一個點。」
他在白板上畫了兩個示意圖。
第一個是傳統的思路，一條曲線，被切成許多小段，每段用直線近似。
第二個是他提出的新思路，一個空間，裡面有很多條曲線，其中一條被特別標出。
「傳統方法是在近似那條曲線，」肖宿指著第一個圖，「而我的想法是，直接在路徑空間裡搜索那條曲線。」
他繼續寫公式，動作流暢得像是已經演練過無數遍：
「路徑空間M上的泛函S[γ]=∫L(γ,γ')dt，其中L是拉格朗日量。我們要找的是使得S取極值的路徑，這就是變分原理。」
「但路徑空間是無窮維的，如何搜索？」
劉浩然問出了關鍵問題。
「用蒙特卡洛方法。」
肖宿說，「但不是傳統的蒙特卡洛，而是在路徑空間上的蒙特卡洛，也就是路徑積分蒙特卡洛。」
他在白板上寫下一串新的公式：
「構造一個馬爾可夫鏈，狀態空間是M，平穩分布是exp(-βS[γ])/Z。當β→∞時，這個分布集中在S的最小值點附近。通過模擬這個馬爾可夫鏈，我們就能以高概率找到最優路徑。」
劉浩然盯著那些公式，腦子裡嗡嗡作響。
這個想法太大膽了，把量子力學中的路徑積分方法，搬到隨機微分方程的數值求解中。
但仔細想想，卻又出奇地自然。
路徑積分本來就是處理無窮維空間問題的利器，而隨機微分方程的解恰恰就是一條路徑。
「這……這能實現嗎？」
他聲音有些發顫。
「能。」
肖宿肯定地說，「只需要設計合適的提議分布和接受概率，保證馬爾可夫鏈在流形上轉移。技術細節很多，但原則上沒有障礙。」
他開始詳細推導。
白板上的公式像泉水一樣湧出，一個接一個，邏輯嚴密，環環相扣。
劉浩然看著，感覺自己不是在見證一個想法的誕生，而是在見證一座數學建築的落成，從地基到框架，從結構到裝飾，每一個部分都精確而優美。
凌晨四點，推導完成。
肖宿放下粉筆，退後一步，審視整個白板。
從左到右，從上到下，從問題的提出到最終的算法設計，完整而自洽。
那些曾經困擾他兩周的障礙，此刻都變成了通往終點的階梯。
他長舒一口氣，閉上眼睛。
劉浩然能看到他臉上的疲憊，但更明顯的是那種解決問題後的滿足感，純粹的，不摻雜任何雜質的滿足。
「完成了？」劉浩然輕聲問。
「完成了。」肖宿睜開眼睛，眼神清澈，「理論框架已經完整。剩下的就是實現和測試。」
窗外，天色已經開始泛白。
新的一天即將開始。
他沒再勸，只是又熱了一杯牛奶放在了肖宿桌上。
凌晨一點，肖宿突然放下粉筆，走到窗前。
窗外是寂靜的校園，只有路燈還孤獨地亮著。
他盯著黑暗中的某一點，眼神空洞，像是在看很遠的地方。
劉浩然知道，這是肖宿在「放空」。
當思考進入死胡同時，他會暫時跳出數學的世界，讓潛意識去工作。
時間一分一秒過去。
辦公室里的鐘表滴答聲格外清晰。
凌晨兩點十七分，肖宿突然轉過身，眼睛裡閃著光。
「我明白了。」
「明白什麼？」
「我們一直在用錯誤的範式思考。」
肖宿走回白板前，語速快而清晰。
「現有的隨機微分方程數值解法，都是基於『離散化時間，逐步推進』的思路。但在流形上，這種思路天然受限，因為每一步推進都要保證不離開流形。」
他擦掉白板上一半的內容，開始重新書寫：
「正確的思路應該是：不是離散化時間，而是離散化路徑空間。」
「什麼？」劉浩然沒聽懂。
「考慮從初始點到終點的所有可能路徑組成的空間，這叫做路徑空間。」
肖宿邊寫邊解釋，「在流形上，路徑空間本身也是一個無窮維流形。我們要找的解，就是這個無窮維流形上的一個點。」
他在白板上畫了兩個示意圖。
第一個是傳統的思路，一條曲線，被切成許多小段，每段用直線近似。
第二個是他提出的新思路，一個空間，裡面有很多條曲線，其中一條被特別標出。
「傳統方法是在近似那條曲線，」肖宿指著第一個圖，「而我的想法是，直接在路徑空間裡搜索那條曲線。」
他繼續寫公式，動作流暢得像是已經演練過無數遍：
「路徑空間M上的泛函S[γ]=∫L(γ,γ')dt，其中L是拉格朗日量。我們要找的是使得S取極值的路徑，這就是變分原理。」
「但路徑空間是無窮維的，如何搜索？」
劉浩然問出了關鍵問題。
「用蒙特卡洛方法。」
肖宿說，「但不是傳統的蒙特卡洛，而是在路徑空間上的蒙特卡洛，也就是路徑積分蒙特卡洛。」
他在白板上寫下一串新的公式：
「構造一個馬爾可夫鏈，狀態空間是M，平穩分布是exp(-βS[γ])/Z。當β→∞時，這個分布集中在S的最小值點附近。通過模擬這個馬爾可夫鏈，我們就能以高概率找到最優路徑。」
劉浩然盯著那些公式，腦子裡嗡嗡作響。
這個想法太大膽了，把量子力學中的路徑積分方法，搬到隨機微分方程的數值求解中。
但仔細想想，卻又出奇地自然。
路徑積分本來就是處理無窮維空間問題的利器，而隨機微分方程的解恰恰就是一條路徑。
「這……這能實現嗎？」
他聲音有些發顫。
「能。」
肖宿肯定地說，「只需要設計合適的提議分布和接受概率，保證馬爾可夫鏈在流形上轉移。技術細節很多，但原則上沒有障礙。」
他開始詳細推導。
白板上的公式像泉水一樣湧出，一個接一個，邏輯嚴密，環環相扣。
劉浩然看著，感覺自己不是在見證一個想法的誕生，而是在見證一座數學建築的落成，從地基到框架，從結構到裝飾，每一個部分都精確而優美。
凌晨四點，推導完成。
肖宿放下粉筆，退後一步，審視整個白板。
從左到右，從上到下，從問題的提出到最終的算法設計，完整而自洽。
那些曾經困擾他兩周的障礙，此刻都變成了通往終點的階梯。
他長舒一口氣，閉上眼睛。
劉浩然能看到他臉上的疲憊，但更明顯的是那種解決問題後的滿足感，純粹的，不摻雜任何雜質的滿足。
「完成了？」劉浩然輕聲問。
「完成了。」肖宿睜開眼睛，眼神清澈，「理論框架已經完整。剩下的就是實現和測試。」
窗外，天色已經開始泛白。
新的一天即將開始。