1. 項目概述一次完整的數(shù)學建模競賽復盤去年帶隊參加美賽MCM/ICM的經(jīng)歷現(xiàn)在回想起來依然覺得收獲滿滿。當時我們選的C題題目是關(guān)于“確認關(guān)于黃蜂的傳言”本質(zhì)上是一個數(shù)據(jù)驅(qū)動的生態(tài)學與傳播動力學交叉問題。很多同學對美賽既向往又畏懼覺得它高深莫測尤其是看到“內(nèi)含完整代碼”這樣的標題可能會以為是一堆難以理解的算法堆砌。其實不然美賽的核心在于用數(shù)學語言清晰地定義現(xiàn)實問題并給出一個有說服力的解決方案。代碼只是實現(xiàn)工具思路才是靈魂。這篇記錄我會以一個過來人的身份徹底復盤我們當時解決2021年美賽C題的全過程。我不會只扔給你一堆代碼而是會詳細拆解我們是如何一步步從拿到賽題時的茫然到建立模型、求解、驗證最終完成一篇20多頁論文的。無論你是正在備賽的新手還是對數(shù)學建模感興趣的朋友都能從中看到一套可復現(xiàn)的方法論。我們會重點講清楚幾個核心問題如何從一段描述中提煉出數(shù)學模型面對不完整的數(shù)據(jù)該怎么辦如何讓論文邏輯自洽、讓評委信服以及那些在官方指導里不會寫的、我們踩過的坑和總結(jié)出的實戰(zhàn)技巧。2. 賽題核心解析與破題思路2.1 題目背景與問題重述2021年美賽C題的標題是“Confirming the Buzz about Hornets”直譯是“確認關(guān)于大黃蜂的傳言”。題目給了一個場景有傳言說在華盛頓州發(fā)現(xiàn)了一種亞洲大黃蜂Vespa mandarinia俗稱“殺人蜂”這種蜂對本地蜜蜂種群和農(nóng)業(yè)有巨大威脅。我們的任務就是分析這個傳言的真實性并評估如果傳言屬實其潛在的生態(tài)影響。題目提供了幾類數(shù)據(jù)一是關(guān)于這種亞洲大黃蜂的生物學特性如飛行速度、活動范圍、繁殖率等二是華盛頓州的地理、氣候和植被數(shù)據(jù)三是可能不可靠的公眾目擊報告。這里第一個難點就出現(xiàn)了數(shù)據(jù)不完整、不一致且存在噪聲。公眾報告可能誤報生態(tài)參數(shù)存在不確定性。美賽的經(jīng)典套路就是給你一個“臟”的現(xiàn)實問題考驗你處理不確定性和缺失信息的能力。我們的破題思路始于對問題本質(zhì)的界定。這不僅僅是一個“是否存在”的判斷題而是一個基于不確定證據(jù)的推理與預測問題。我們需要構(gòu)建一個模型能夠評估傳言可信度根據(jù)零散、可能錯誤的目擊報告結(jié)合地理和生物知識計算該蜂種在華盛頓州出現(xiàn)的概率。預測傳播動態(tài)如果存在它們會如何擴散速度和范圍如何評估生態(tài)風險擴散會對本地蜜蜂、農(nóng)作物授粉造成多大影響我們將問題分解為三個子模型信度評估模型、空間傳播模型和生態(tài)影響模型。這三個模型環(huán)環(huán)相扣前者的輸出是后者的輸入。2.2 模型選型背后的邏輯為什么選這些模型這是論文能否拿獎的關(guān)鍵必須在論文中清晰闡述你的“Why”。對于信度評估我們放棄了復雜的機器學習分類因為數(shù)據(jù)量小且標簽模糊選擇了貝葉斯推理網(wǎng)絡(luò)。這是核心亮點。貝葉斯方法的優(yōu)勢在于它能天然地處理不確定性和先驗知識。我們將“蜂群存在”作為一個待估計的概率事件將每一條目擊報告包含時間、地點、描述可信度作為證據(jù)。通過設(shè)定先驗概率基于該物種的原始分布和入侵歷史并利用證據(jù)逐步更新后驗概率。例如一個來自經(jīng)驗豐富的昆蟲學家的報告其證據(jù)強度似然比遠高于一個普通游客模糊的描述。這個模型不僅能給出一個整體的存在概率還能量化每個證據(jù)的貢獻甚至能進行反事實分析“如果某條最關(guān)鍵的報告是誤報結(jié)論會如何變化”這極大地增強了分析的深度。對于空間傳播模型常見的選擇有元胞自動機CA和偏微分方程PDE反應擴散模型。我們選擇了后者——Fisher-KPP方程。這是一個經(jīng)典的生物入侵擴散模型。原因在于CA雖然直觀但參數(shù)多物理意義不如PDE清晰。Fisher-KPP方程?u/?t D?2u ru(1-u/K)形式優(yōu)美其中u是種群密度D是擴散系數(shù)r是內(nèi)稟增長率K是環(huán)境承載量。D可以由大黃蜂的飛行能力和景觀阻力如森林、水域估算r和K來自其生物學特性及華盛頓州的氣候-植被匹配度。這個模型的優(yōu)勢在于它有豐富的數(shù)學理論支撐我們可以討論波前傳播的漸進速度v ≈ 2√(Dr)這能給出一目了然的擴散速度預測比單純的模擬更顯功力。對于生態(tài)影響模型我們采用了耦合的Lotka-Volterra競爭模型將亞洲大黃蜂和本地蜜蜂視為競爭關(guān)系。同時我們建立了一個簡化的“授粉服務-作物產(chǎn)量”關(guān)系鏈將蜜蜂種群數(shù)量的變化映射到幾種主要依賴蟲媒授粉的農(nóng)作物如蘋果、藍莓的潛在產(chǎn)量損失上。這部分的關(guān)鍵在于參數(shù)估計的合理性我們大量引用了已有的生態(tài)學研究文獻來支撐參數(shù)取值。注意模型的美觀性與實用性的權(quán)衡。美賽評委欣賞有數(shù)學深度的模型但更看重模型是否恰當?shù)亟鉀Q了問題。Fisher-KPP方程比元胞自動機“更數(shù)學”但我們必須解釋清楚參數(shù)Dr是如何從題目給出的有限數(shù)據(jù)中合理推斷出來的否則就是空中樓閣。我們的策略是用主要模型PDE展示數(shù)學能力同時用更簡單的輔助模型如基于GIS的成本加權(quán)距離分析進行交叉驗證體現(xiàn)思維的嚴謹。3. 核心模塊實現(xiàn)與代碼詳解3.1 貝葉斯信度評估模型的實現(xiàn)這部分是分析的起點我們用Python的pymc3現(xiàn)為pymc庫來實現(xiàn)概率編程。核心是構(gòu)建一個貝葉斯分層模型。首先定義關(guān)鍵隨機變量presence_prob: 全局存在概率先驗設(shè)為0.05基于“傳言”的保守估計。report_reliability[i]: 第i條報告的可信度Beta分布根據(jù)報告來源類型設(shè)定先驗如專家報告的先驗均值更高。observation[i]: 觀察到的第i條報告伯努利分布以presence_prob * report_reliability[i]為參數(shù)。import pymc3 as pm import numpy as np # 假設(shè)有n條報告reports_prior是每條報告可信度的先驗均值來自人工評估 n len(reports) with pm.Model() as bayesian_model: # 先驗蜂群存在的概率我們假設(shè)較低 presence_prob pm.Beta(presence_prob, alpha1, beta20) # 期望約0.05 # 每條報告的可信度先驗 # 例如如果報告來源是“官方機構(gòu)”alpha和beta參數(shù)使其均值接近0.9 # 如果是“社交媒體”均值可能只有0.3 reliability_priors define_priors_based_on_source(reports_source) report_reliability pm.Beta(report_reliability, alphareliability_priors[alpha], betareliability_priors[beta], shapen) # 似然觀測到的報告1為真0為假或無關(guān)。這里簡化實際應根據(jù)報告內(nèi)容建模。 # 核心邏輯一條真實的報告既要求蜂群確實存在也要求該報告本身可靠。 p_obs presence_prob * report_reliability observations pm.Bernoulli(observations, pp_obs, observedreports_observed) # 采樣推斷 trace pm.sample(2000, tune1000, return_inferencedataFalse)實操要點先驗的選擇不是隨意的。我們?yōu)閜resence_prob選擇了Beta(1,20)意味著在沒有任何證據(jù)前我們認為存在的可能性只有5%1/(120)。這是一個保守的、對傳言懷疑的先驗。在論文中我們專門用一小節(jié)論證了這個先驗的合理性并測試了不同先驗如Beta(1,1)均勻先驗的敏感性結(jié)果顯示后驗概率雖然數(shù)值有變化但定性結(jié)論概率是否超過某個閾值是穩(wěn)健的。這體現(xiàn)了分析的嚴謹性。reports_observed是我們輸入的數(shù)據(jù)。這里有一個技巧我們并沒有簡單地將報告記為“1”存在。我們對報告內(nèi)容進行了文本分析提取了關(guān)鍵詞如“體型巨大”、“橙色頭部”、“筑巢行為”將其與亞洲大黃蜂的特征進行匹配匹配度作為一個0到1的數(shù)值輸入。這比二值化處理包含了更多信息。采樣后通過pm.plot_posterior(trace)可以直觀看到presence_prob的后驗分布。我們得到的后驗概率均值達到了0.78這為“傳言可能屬實”提供了定量支持。3.2 基于Fisher-KPP方程的擴散模擬我們用有限差分法在二維網(wǎng)格上數(shù)值求解Fisher-KPP方程。網(wǎng)格對應華盛頓州地圖每個網(wǎng)格點的K環(huán)境承載量由該點的植被類型和氣候適宜度決定。import numpy as np from scipy.ndimage import distance_transform_edt def simulate_invasion(D, r, K_map, initial_map, dx, dt, total_steps): 使用顯式有限差分法求解2D Fisher-KPP方程。 參數(shù): D: 擴散系數(shù) (km^2/week) r: 內(nèi)稟增長率 (1/week) K_map: 空間承載量地圖 (2D array) initial_map: 初始種群分布 (2D array) dx: 空間步長 (km) dt: 時間步長 (week) - 必須滿足穩(wěn)定性條件 dt dx^2 / (4*D) total_steps: 總模擬步數(shù) 返回: 時間序列的種群分布 u initial_map.copy() history [u.copy()] ny, nx u.shape # 穩(wěn)定性檢查 max_dt dx**2 / (4 * D) if dt max_dt: print(f警告時間步長{dt}可能不穩(wěn)定建議小于{max_dt:.4f}) # 實踐中我們這里選擇減小dt或使用隱式方法 for step in range(total_steps): u_new u.copy() # 使用五點差分格式計算拉普拉斯項擴散 laplacian (np.roll(u, 1, axis0) np.roll(u, -1, axis0) np.roll(u, 1, axis1) np.roll(u, -1, axis1) - 4 * u) / dx**2 # Fisher-KPP方程離散化 reaction r * u * (1 - u / K_map) diffusion D * laplacian u_new u dt * (diffusion reaction) # 邊界條件假設(shè)邊界為不可逾越的障礙Neumann零通量邊界 # 通過roll操作實現(xiàn)的周期性邊界并不合適這里更佳做法是直接處理邊界網(wǎng)格。 # 為簡化示例我們使用零通量近似的簡便寫法邊界點不參與擴散計算。 u_new[0, :] u[0, :] dt * reaction[0, :] # 上邊界 u_new[-1, :] u[-1, :] dt * reaction[-1, :] # 下邊界 u_new[:, 0] u[:, 0] dt * reaction[:, 0] # 左邊界 u_new[:, -1] u[:, -1] dt * reaction[:, -1] # 右邊界 u np.clip(u_new, 0, None) # 種群密度非負 history.append(u.copy()) return np.array(history)參數(shù)估計細節(jié)擴散系數(shù) D題目給出了大黃蜂的飛行速度約40 km/年。我們將其轉(zhuǎn)換為擴散系數(shù)。對于隨機擴散D與速度v和擴散時間t的關(guān)系近似為σ2 2Dt其中σ是標準差。我們假設(shè)一個季節(jié)13周內(nèi)的擴散距離標準差約為速度乘以時間進行量綱換算后估算出D ≈ 10-20 km2/week。這是一個關(guān)鍵推導在論文附錄中我們詳細展示了計算過程。內(nèi)稟增長率 r根據(jù)其繁殖生物學蜂后產(chǎn)卵量、發(fā)育周期估算出在理想條件下每周種群增長率約為0.1-0.15。然后根據(jù)華盛頓州的氣候匹配度使用MAXENT生態(tài)位模型的思路比較原產(chǎn)地和華盛頓的氣候變量對其進行折減得到空間變化的r_map。初始地圖initial_map我們假設(shè)如果存在初始入侵點位于公眾報告最密集的聚類中心。用一個二維高斯分布模擬初始種群分布。踩坑實錄數(shù)值穩(wěn)定性。最初我們沒注意穩(wěn)定性條件dt dx^2 / (4D)導致模擬后期出現(xiàn)數(shù)值爆炸種群密度激增到天文數(shù)字。調(diào)試了很久才發(fā)現(xiàn)。心得是對于任何微分方程數(shù)值求解第一步永遠是進行量綱分析和穩(wěn)定性條件估算。后來我們改用更穩(wěn)定的隱式方法如Crank-Nicolson或自適應步長庫如solve_ivp但為了論文代碼的簡潔和可讀性最終版本保留了顯式格式并嚴格限制了步長。3.3 生態(tài)影響與風險評估耦合我們將傳播模型的輸出未來不同時間點的種群分布圖作為輸入接入生態(tài)影響模型。def assess_impact(bee_population_map, hornet_population_map, crop_map): 評估對本地蜜蜂和農(nóng)作物的影響。 參數(shù): bee_population_map: 本地蜜蜂種群密度圖 hornet_population_map: 大黃蜂種群密度圖 crop_map: 作物分布與依賴度圖值表示授粉依賴程度和經(jīng)濟效益權(quán)重 返回: competition_loss: 蜜蜂種群相對損失 crop_risk_index: 作物風險指數(shù) # 使用Lotka-Volterra競爭模型計算平衡點偏移 # 簡化計算假設(shè)蜜蜂受到直接捕殺和競爭壓力損失率與大黃蜂密度成正比 # α: 競爭系數(shù)表示單位大黃蜂密度對蜜蜂增長率的抑制 alpha 0.05 # 基于文獻的估計值 bee_growth_rate_reduction alpha * hornet_population_map # 估算蜜蜂種群損失穩(wěn)態(tài)近似 competition_loss bee_population_map * (1 - np.exp(-bee_growth_rate_reduction * 1)) # 考慮一年影響 # 計算作物風險蜜蜂損失導致授粉服務下降進而影響產(chǎn)量 # 假設(shè)產(chǎn)量損失與蜜蜂損失成正比并加權(quán)作物經(jīng)濟價值 pollination_deficit competition_loss / (bee_population_map 1e-5) # 避免除零 crop_yield_loss_ratio pollination_deficit * crop_map[pollination_dependency] crop_risk_index np.sum(crop_yield_loss_ratio * crop_map[economic_weight]) return competition_loss, crop_risk_index注意事項生態(tài)影響模型是不確定性最大的環(huán)節(jié)。競爭系數(shù)alpha、作物依賴度等參數(shù)都有很大的變化范圍。我們的策略是進行全面的敏感性分析和情景模擬。在論文中我們設(shè)置了樂觀、基準、悲觀三組參數(shù)分別運行模型給出了影響的范圍例如“預計5年內(nèi)本地蜜蜂種群可能減少10%-40%”。這種表述方式比給出一個單一的確切數(shù)字更科學、更令人信服。代碼上我們用一個循環(huán)包裹了上面的評估函數(shù)遍歷多組參數(shù)并生成一系列圖表來展示結(jié)果的范圍和分布。4. 論文寫作與結(jié)果可視化心法4.1 如何組織一篇邏輯清晰的論文美賽論文有標準的框架摘要、引言、假設(shè)、模型、求解、分析、結(jié)論等但內(nèi)在邏輯的流暢性才是高分關(guān)鍵。我們的行文邏輯如下摘要用一段話概括全部工作。我們采用了“問題-方法-關(guān)鍵結(jié)果-結(jié)論”的四段式。明確指出使用了貝葉斯模型后驗概率0.78、Fisher-KPP擴散模型預測前沿擴散速度~15 km/年和生態(tài)風險評估蜜蜂種群潛在損失10-40%并給出了核心建議加強監(jiān)測、優(yōu)先控制特定區(qū)域。引言與問題重述用自己的語言復述問題并明確列出我們要解決的幾個具體子問題即2.1中的三個。讓評委一眼就知道你理解準確且思路清晰。假設(shè)與合理性論證這是展示你科學素養(yǎng)的地方。每一條假設(shè)都要有理由。例如“假設(shè)大黃蜂的擴散主要受景觀類型影響而忽略細微地形起伏”理由是飛行高度足以克服小地形障礙且現(xiàn)有數(shù)據(jù)精度不支持更細粒度分析。為關(guān)鍵假設(shè)進行敏感性測試如改變擴散系數(shù)D看結(jié)果變化能極大增強論文的魯棒性。模型建立分小節(jié)介紹貝葉斯模型、擴散模型、影響模型。每個小節(jié)都遵循“為什么選這個模型-模型數(shù)學形式-參數(shù)如何估計-如何求解”的結(jié)構(gòu)。把公式和文字解釋穿插好避免大段純公式或純文字。模型求解與結(jié)果分析展示核心結(jié)果圖如后驗概率分布、擴散模擬動畫截圖、風險地圖。分析要深入不僅說“圖1顯示種群在擴散”還要說“擴散呈現(xiàn)明顯的各向異性東部山區(qū)速度慢于西部河谷這與我們設(shè)定的景觀阻力一致...”。靈敏度分析與模型檢驗專門一節(jié)展示改變關(guān)鍵參數(shù)的結(jié)果討論模型的局限性。我們甚至用歷史上另一種入侵物種如非洲化蜜蜂的數(shù)據(jù)做了粗略的回顧性預測檢驗以佐證模型的有效性。結(jié)論與建議總結(jié)發(fā)現(xiàn)并提出具體、可操作的建議。例如“建議在斯諾夸爾米山口設(shè)立監(jiān)測點因為模型顯示這是擴散到東部農(nóng)業(yè)區(qū)的主要通道”這比泛泛而談“加強監(jiān)測”要好得多。4.2 讓圖表“說話”的可視化技巧好的圖表能抵千言萬語。我們主要使用matplotlib和seaborn。貝葉斯結(jié)果可視化使用seaborn的kdeplot繪制先驗和后驗概率的分布對比清晰展示數(shù)據(jù)證據(jù)如何更新了我們的認知。import seaborn as sns import matplotlib.pyplot as plt fig, ax plt.subplots(1, 2, figsize(12,4)) # 繪制先驗分布 prior_samples np.random.beta(1, 20, 10000) sns.kdeplot(prior_samples, axax[0], fillTrue) ax[0].set_title(Prior Distribution of Presence Probability) ax[0].set_xlabel(Probability) ax[0].set_ylabel(Density) # 繪制后驗分布從trace中提取 posterior_samples trace[presence_prob] sns.kdeplot(posterior_samples, axax[1], fillTrue) ax[1].axvline(posterior_samples.mean(), colorred, linestyle--, labelfMean{posterior_samples.mean():.2f}) ax[1].set_title(Posterior Distribution of Presence Probability) ax[1].set_xlabel(Probability) ax[1].legend() plt.tight_layout()擴散動態(tài)可視化制作多面板圖展示不同時間點如0125年后的種群分布。使用matplotlib的imshow并配上地理底圖從GIS數(shù)據(jù)中獲取的華盛頓州輪廓。關(guān)鍵是要有一個清晰、一致的顏色條并標注比例尺。風險地圖使用漸變色如viridis表示風險等級并疊加重要的地理要素如主要城市、農(nóng)田分布。在論文中我們將風險地圖與建議的監(jiān)測點直接標注在同一張圖上使建議一目了然。心得圖表的美觀與信息密度平衡。切忌花里胡哨。美賽論文通常是黑白打印所以顏色映射要選擇在灰度下也能區(qū)分的如viridis, plasma。確保所有坐標軸都有標簽單位清晰。多圖排列時對齊、尺寸一致。在提交前一定要將論文導出為PDF檢查圖表是否清晰字體是否過小。我們吃過虧初版圖表在PDF里模糊不清最后時刻緊急重調(diào)了DPI至少300和字體大小。5. 團隊協(xié)作、時間管理與常見避坑指南5.1 四天時間的節(jié)奏把控美賽總共四天時間管理就是生命線。第一天破題與規(guī)劃上午全體成員深入讀題、討論禁止立即查文獻。下午確定初步思路和模型框架并開始撰寫論文的“引言”和“問題重述”部分。晚上分工一人負責文獻調(diào)研和參數(shù)搜集一人開始搭建核心模型代碼框架一人繼續(xù)完善論文前半部分和假設(shè)列表。第二天模型實現(xiàn)與初稿上午模型實現(xiàn)者應產(chǎn)出初步結(jié)果哪怕很粗糙。下午全體會議根據(jù)初步結(jié)果調(diào)整模型細節(jié)。晚上論文撰寫者必須完成模型的數(shù)學描述部分并將初步結(jié)果做成圖表。這一天結(jié)束前論文應有一個完整的草稿包括所有章節(jié)的標題和大部分文字。第三天深入分析與打磨全天進行模型調(diào)試、敏感性分析、結(jié)果深化。論文撰寫者整合所有新結(jié)果開始寫“結(jié)果分析”和“靈敏度分析”。編程者負責生成最終版的所有圖表。晚上必須完成論文初稿的90%包括摘要的草稿。第四天終稿與提交上午集中精力寫摘要、潤色結(jié)論、檢查全文邏輯。摘要反復修改直到精煉準確。下午進行最終排版、交叉檢查、錯別字和語法修正推薦使用Grammarly等工具。至少預留3小時用于最終PDF生成、上傳和提交網(wǎng)絡(luò)擁堵和最后一刻的修改是常態(tài)。5.2 那些我們踩過的“坑”與應對策略坑追求模型復雜度忽視可解釋性。初期我們想過用復雜的神經(jīng)網(wǎng)絡(luò)來分類目擊報告但很快發(fā)現(xiàn)數(shù)據(jù)量太小且模型像個黑箱難以在論文中闡述清楚。對策回歸到貝葉斯模型它的每個參數(shù)都有明確的概率解釋評委容易理解也更容易與題目背景結(jié)合。坑參數(shù)估計憑感覺。最初給擴散系數(shù)D隨便賦了個值。對策我們意識到這不行于是花時間查閱昆蟲擴散的文獻找到了將飛行速度轉(zhuǎn)換為擴散系數(shù)的經(jīng)典生態(tài)學公式并在論文中詳細引證。這極大地提升了模型的可信度?？哟a調(diào)試耗時過長。數(shù)值模擬程序跑不出結(jié)果或結(jié)果怪異。對策從小規(guī)模、簡化模型開始驗證。例如先在一維、參數(shù)恒定的情況下運行Fisher-KPP方程與解析解行波解對比確認代碼正確后再擴展到二維異質(zhì)環(huán)境。設(shè)置斷言assert檢查數(shù)組邊界和數(shù)值范圍?？诱撐膶懽髋c建模脫節(jié)。寫論文的人不知道模型細節(jié)建模的人不關(guān)心文字表達。對策強制同步。每天早晚開短會建模者向?qū)懽髡咧v解當天進展寫作者隨時將寫好的部分給建模者核對。使用共享文檔如Overleaf實時協(xié)作。坑摘要最后才寫倉促完成。摘要是最重要的部分評委可能只看摘要。對策從第二天晚上就開始起草摘要隨著工作推進不斷更新。最后一天花至少1小時集體字斟句酌地修改摘要確保它獨立、完整、準確地概括全文所有亮點。5.3 給未來參賽者的工具箱建議編程語言Python是絕對主流生態(tài)豐富numpy,scipy,pandas,matplotlib,pymc,geopandas。Matlab在微分方程求解和快速畫圖上仍有優(yōu)勢但Python在數(shù)據(jù)預處理、復雜工作流和開源協(xié)作上更勝一籌。文獻與數(shù)據(jù)除了題目附件Google Scholar、USGS美國地質(zhì)調(diào)查局、FAO聯(lián)合國糧農(nóng)組織網(wǎng)站是數(shù)據(jù)寶庫。學會用關(guān)鍵詞組合搜索如 “Vespa mandarinia dispersal rate”, “insect invasion model”。協(xié)作工具OverleafLaTeX在線寫作實時協(xié)作、GitHub代碼版本管理、騰訊會議/釘釘日常溝通、Notion或飛書任務管理。繪圖與可視化matplotlib基礎(chǔ)seaborn美化統(tǒng)計圖plotly可交互圖表可嵌入網(wǎng)頁摘要。對于地圖geopandas結(jié)合contextily可以輕松添加在線底圖。參加美賽是一次高強度、高回報的鍛煉。它逼著你在極短時間內(nèi)將一個模糊的現(xiàn)實問題轉(zhuǎn)化為清晰的數(shù)學問題并給出一個有數(shù)據(jù)、有模型、有見解的解決方案。這個過程里學到的問題分解能力、假設(shè)建模能力和故事講述能力遠比學會幾個算法更重要。希望這份詳細的解題記錄能為你照亮前行的路。記住清晰的思路和嚴謹?shù)谋磉_永遠是贏得比賽的關(guān)鍵。