1. 項(xiàng)目概述從一篇論文到一套方法論每年數(shù)學(xué)建模競(jìng)賽結(jié)束后總有一批“神作”在圈內(nèi)流傳2023年C題的優(yōu)秀論文就是其中之一。我花了將近一周時(shí)間把能找到的高分論文都捋了一遍不是為了膜拜而是想搞清楚一件事在題目看似常規(guī)的“蔬菜類商品自動(dòng)定價(jià)與補(bǔ)貨決策”背后這些頂尖隊(duì)伍到底做對(duì)了什么才讓他們的論文從海量作品中脫穎而出這不僅僅是看他們用了什么模型更是要拆解他們構(gòu)建問題、分析數(shù)據(jù)、設(shè)計(jì)方案的完整邏輯鏈條。對(duì)于準(zhǔn)備參賽的同學(xué)或者任何需要處理復(fù)雜優(yōu)化決策問題的人來說這種“逆向工程”式的分析價(jià)值遠(yuǎn)超單純學(xué)習(xí)一個(gè)算法。它關(guān)乎如何將模糊的實(shí)際問題轉(zhuǎn)化為清晰、可解、且能打動(dòng)評(píng)委的數(shù)學(xué)模型。簡(jiǎn)單來說2023年C題要求參賽者為一個(gè)生鮮超市的蔬菜品類設(shè)計(jì)自動(dòng)化的定價(jià)與補(bǔ)貨策略。題目給出了過去四年多個(gè)蔬菜品類的銷售數(shù)據(jù)包括成本、售價(jià)、銷量、損耗等。核心目標(biāo)就兩個(gè)第一基于歷史數(shù)據(jù)建立定價(jià)模型實(shí)現(xiàn)利潤(rùn)最大化第二設(shè)計(jì)補(bǔ)貨決策模型在滿足需求的同時(shí)最小化損耗。聽起來像是經(jīng)典的“報(bào)童問題”升級(jí)版但實(shí)際做起來坑多得能絆倒一個(gè)排。很多隊(duì)伍折戟沉沙要么是模型過于理想化脫離實(shí)際要么是求解復(fù)雜到算不出來要么是論文寫得像實(shí)驗(yàn)報(bào)告缺乏洞察。而我分析的這些優(yōu)秀論文它們共同的成功秘訣在于沒有把題目當(dāng)成一道數(shù)學(xué)題而是當(dāng)成一個(gè)真實(shí)的商業(yè)決策模擬。它們深入挖掘了數(shù)據(jù)背后的故事比如蔬菜銷售的強(qiáng)季節(jié)性、不同品類間的關(guān)聯(lián)性、損耗的非線性特征并將這些洞察巧妙地融入了模型假設(shè)和約束條件中。接下來我就把自己從這些“學(xué)霸筆記”里提煉出的核心思路、關(guān)鍵技術(shù)選型、實(shí)操細(xì)節(jié)以及那些容易踩坑的地方系統(tǒng)地分享給你。無論你是備戰(zhàn)下一次數(shù)模競(jìng)賽還是在工作中面臨類似的資源優(yōu)化問題相信這些來自實(shí)戰(zhàn)頂端的經(jīng)驗(yàn)都能給你帶來實(shí)實(shí)在在的啟發(fā)。2. 核心思路拆解優(yōu)秀論文的四大共性特征看完幾十篇論文后我發(fā)現(xiàn)高分作品在思路上有明顯的共性。它們沒有追求算法的炫技而是在問題理解的深度、模型結(jié)合的巧度、以及方案表述的清晰度上做到了極致。下面我把它拆解為四個(gè)關(guān)鍵特征。2.1 特征一問題重構(gòu)與多階段建模題目直接問“定價(jià)”和“補(bǔ)貨”但高手絕不會(huì)直接套用現(xiàn)成的定價(jià)模型或庫存公式。他們首先做的是問題重構(gòu)。他們會(huì)問定價(jià)和補(bǔ)貨真的是兩個(gè)獨(dú)立問題嗎顯然不是。補(bǔ)貨量取決于對(duì)需求的預(yù)測(cè)而需求又受價(jià)格直接影響。因此這本質(zhì)上是一個(gè)聯(lián)合決策優(yōu)化問題。大部分優(yōu)秀論文都采用了多階段建模的框架數(shù)據(jù)預(yù)處理與特征工程階段這不是簡(jiǎn)單清洗而是深度挖掘。例如除了常規(guī)的日期、銷量、成本他們會(huì)構(gòu)建衍生特征周銷量增長(zhǎng)率、品類間銷量相關(guān)系數(shù)、歷史損耗率與銷售量的比值等。有一篇論文特別指出葉菜類如菠菜的損耗對(duì)前一日剩余庫存量極為敏感而非簡(jiǎn)單的線性關(guān)系這個(gè)洞察直接影響了后續(xù)損耗模型的建立。需求預(yù)測(cè)階段這是整個(gè)模型的基石。幾乎所有高分論文都摒棄了單一的預(yù)測(cè)模型而是采用了集成或分層的預(yù)測(cè)思路。對(duì)于整體趨勢(shì)如年度、季度可能使用時(shí)間序列模型如SARIMA對(duì)于節(jié)假日、周末等季節(jié)性因素采用季節(jié)性分解或LightGBM/XGBoost這類樹模型對(duì)于新品或數(shù)據(jù)稀疏的品類則采用類似品類的遷移學(xué)習(xí)或貝葉斯平滑方法。關(guān)鍵點(diǎn)在于他們不僅預(yù)測(cè)了“點(diǎn)估計(jì)”平均銷量更預(yù)測(cè)了“概率分布”或“置信區(qū)間”為后續(xù)的魯棒優(yōu)化打下了基礎(chǔ)。聯(lián)合優(yōu)化階段在此階段將預(yù)測(cè)的需求作為隨機(jī)變量或帶參數(shù)的分布嵌入到一個(gè)優(yōu)化模型中。目標(biāo)函數(shù)通常是最大化期望利潤(rùn)而決策變量就是價(jià)格和補(bǔ)貨量。約束條件則復(fù)雜得多包括庫存容量約束、資金約束、最小展示量約束超市貨架不能空、以及最重要的——損耗與庫存和時(shí)間的動(dòng)態(tài)關(guān)系約束。很多論文在這里引入了條件風(fēng)險(xiǎn)價(jià)值CVaR的概念不再只追求平均利潤(rùn)最高而是要求“在95%的情況下利潤(rùn)不低于某個(gè)值”這更符合企業(yè)穩(wěn)健經(jīng)營(yíng)的需求。策略輸出與模擬驗(yàn)證階段優(yōu)化模型求解出的可能是一組復(fù)雜的函數(shù)或規(guī)則。優(yōu)秀論文會(huì)將其轉(zhuǎn)化為可執(zhí)行的“策略表”或“決策樹”。例如“當(dāng)A蔬菜昨日剩余庫存大于X且預(yù)測(cè)明日氣溫高于Y度時(shí)采取降價(jià)Z%的策略并補(bǔ)貨W公斤”。最后他們一定會(huì)用歷史數(shù)據(jù)或模擬數(shù)據(jù)對(duì)策略進(jìn)行回測(cè)并與簡(jiǎn)單的基準(zhǔn)策略如固定價(jià)格、按昨日銷量補(bǔ)貨對(duì)比用清晰的圖表展示利潤(rùn)提升、損耗降低的具體幅度。注意切忌將多階段建模做成流水賬。每個(gè)階段之間的邏輯傳遞必須清晰。例如需求預(yù)測(cè)階段的輸出誤差如何在優(yōu)化階段被考慮是通過隨機(jī)規(guī)劃還是在目標(biāo)函數(shù)中加入懲罰項(xiàng)這一點(diǎn)必須在論文中明確闡述這是評(píng)委考察建模連貫性的重點(diǎn)。2.2 特征二成本與損耗模型的精細(xì)化處理這是區(qū)分普通論文和優(yōu)秀論文的最大分水嶺。很多隊(duì)伍對(duì)損耗的處理過于粗糙簡(jiǎn)單地用“銷量*固定損耗率”來計(jì)算這完全不符合生鮮零售的現(xiàn)實(shí)。高分論文普遍將損耗分為兩大類自然損耗隨時(shí)間推移必然發(fā)生的腐敗、脫水。這部分通常建模為關(guān)于庫存時(shí)間Age的函數(shù)例如指數(shù)衰減模型損耗量 初始庫存 * (1 - exp(-k * t))其中k是腐敗速率系數(shù)不同蔬菜葉菜vs根莖的k值差異巨大。機(jī)會(huì)損耗因未能及時(shí)售出而被迫丟棄的部分。這部分與需求不確定性和定價(jià)策略強(qiáng)相關(guān)。如果定價(jià)過高導(dǎo)致銷量低于預(yù)期剩余庫存就會(huì)轉(zhuǎn)化為機(jī)會(huì)損耗。更精彩的在于對(duì)成本的分解。除了顯性的進(jìn)貨成本優(yōu)秀論文還考慮了損耗成本不僅是貨物成本還包括處理損耗的人工、垃圾清運(yùn)等隱形成本。缺貨成本貨架空空導(dǎo)致的客戶滿意度下降和潛在銷售損失。這部分難以量化但有的論文巧妙地用“歷史同期同類產(chǎn)品銷量”作為機(jī)會(huì)成本的參照。資金占用成本庫存占用的資金本可以用于其他投資這部分成本在長(zhǎng)期運(yùn)營(yíng)中不可忽視。有一篇讓我印象深刻的論文它甚至為不同的蔬菜品類賦予了不同的“新鮮度價(jià)值衰減曲線”。例如西蘭花一旦花蕾變黃價(jià)值急劇下跌適合采用“階梯式折扣”定價(jià)策略而土豆的貶值曲線則平緩得多可以采用“緩慢降價(jià)”策略。這種基于商品特性的差異化建模極大地提升了模型的現(xiàn)實(shí)解釋力。2.3 特征三算法選型與求解的務(wù)實(shí)性面對(duì)這樣一個(gè)包含隨機(jī)變量、非線性約束的優(yōu)化問題算法求解是另一大挑戰(zhàn)。我觀察到頂尖隊(duì)伍在算法選型上非常務(wù)實(shí)遵循“能用簡(jiǎn)單的不用復(fù)雜的能結(jié)合使用的不單一依賴”的原則。預(yù)測(cè)算法如前所述融合模型是主流。具體工具上Python的statsmodels庫用于傳統(tǒng)時(shí)間序列scikit-learn和LightGBM用于機(jī)器學(xué)習(xí)回歸。他們通常會(huì)做特征重要性分析來驗(yàn)證構(gòu)建的特征是否有效比如發(fā)現(xiàn)“前一周平均銷量”比“上月同期銷量”預(yù)測(cè)力更強(qiáng)這本身就是一個(gè)有價(jià)值的發(fā)現(xiàn)。優(yōu)化算法對(duì)于簡(jiǎn)化后的確定性模型如將需求取為預(yù)測(cè)均值線性規(guī)劃LP或混合整數(shù)線性規(guī)劃MILP可以通過PuLP、Gurobi等求解器高效解決。這是很多論文的起點(diǎn)。對(duì)于考慮不確定性的隨機(jī)規(guī)劃或魯棒優(yōu)化模型直接求解可能非常困難。常用的方法是樣本平均近似SAA從需求預(yù)測(cè)的概率分布中抽取大量場(chǎng)景例如1000個(gè)將隨機(jī)優(yōu)化問題轉(zhuǎn)化為一個(gè)大規(guī)模但確定性的優(yōu)化問題再用求解器處理。雖然計(jì)算量大但借助云計(jì)算平臺(tái)或高效求解器是可實(shí)現(xiàn)的。對(duì)于模型過于復(fù)雜、無法用傳統(tǒng)優(yōu)化方法求解的情況不少論文轉(zhuǎn)向了啟發(fā)式算法如遺傳算法GA、模擬退火SA。這里的關(guān)鍵不是證明算法全局最優(yōu)而是證明該算法找到的解顯著且穩(wěn)定地優(yōu)于人工經(jīng)驗(yàn)策略或其他基準(zhǔn)算法。論文中需要展示算法參數(shù)的調(diào)優(yōu)過程如交叉率、變異率和收斂曲線。仿真驗(yàn)證無論用什么優(yōu)化算法得到策略都必須通過仿真來評(píng)估。他們會(huì)編寫一個(gè)模擬器模擬長(zhǎng)達(dá)數(shù)月的每日運(yùn)營(yíng)隨機(jī)生成符合歷史規(guī)律的需求應(yīng)用自己的定價(jià)補(bǔ)貨策略計(jì)算每日利潤(rùn)和損耗。這個(gè)過程通常重復(fù)成千上萬次蒙特卡洛模擬以得到策略性能的統(tǒng)計(jì)分布平均利潤(rùn)、利潤(rùn)方差、最大回撤等。這種“用仿真檢驗(yàn)優(yōu)化”的思路極大地增加了論文結(jié)論的可信度。2.4 特征四論文表述與可視化呈現(xiàn)再好的模型如果表達(dá)不清也難獲高分。優(yōu)秀論文在寫作和呈現(xiàn)上堪稱范本。摘要采用“問題-方法-結(jié)果-結(jié)論”的黃金結(jié)構(gòu)。開門見山指出研究問題簡(jiǎn)述核心方法如“我們建立了基于隨機(jī)規(guī)劃的多階段聯(lián)合決策模型”用數(shù)據(jù)說話展示核心結(jié)果如“相比基準(zhǔn)策略平均利潤(rùn)提升15.2%損耗降低22.7%”最后總結(jié)模型優(yōu)點(diǎn)。模型假設(shè)假設(shè)不是隨便寫的。每一條假設(shè)都有其目的要么是為了簡(jiǎn)化問題如“忽略運(yùn)輸時(shí)間”要么是基于數(shù)據(jù)分析的結(jié)論如“假設(shè)單個(gè)品類的需求相互獨(dú)立因?yàn)橄嚓P(guān)系數(shù)分析顯示大部分品類間相關(guān)性小于0.3”。后者顯示了從數(shù)據(jù)到模型的推導(dǎo)過程。圖表摒棄了軟件默認(rèn)的丑陋圖表。所有圖表都精心設(shè)計(jì)顏色簡(jiǎn)潔、標(biāo)注清晰。典型的高分圖表包括?；鶊D展示資金流、貨物流采購(gòu)、銷售、損耗的走向和比例。熱力圖展示不同蔬菜品類在不同季節(jié)的利潤(rùn)率或損耗率。時(shí)間序列對(duì)比圖將模型推薦的補(bǔ)貨量、定價(jià)與實(shí)際歷史數(shù)據(jù)、基準(zhǔn)策略進(jìn)行多曲線對(duì)比一目了然。參數(shù)敏感性分析圖展示關(guān)鍵參數(shù)如損耗率系數(shù)、需求波動(dòng)率在±20%范圍內(nèi)變動(dòng)時(shí)對(duì)總利潤(rùn)的影響。這證明了模型的魯棒性。代碼與附錄雖然正文不展示代碼但他們會(huì)提及核心算法的實(shí)現(xiàn)工具如PythonPandasPuLP并在附錄中給出關(guān)鍵的算法流程圖或偽代碼。清晰的流程圖能讓評(píng)委快速理解你的求解邏輯。3. 從零到一的實(shí)操要點(diǎn)與避坑指南分析了這么多優(yōu)秀特征我們來看看如果自己動(dòng)手該如何一步步實(shí)現(xiàn)以及過程中有哪些“坑”需要提前避開。3.1 數(shù)據(jù)預(yù)處理不僅僅是清洗題目給的數(shù)據(jù)通?！芭K”且包含信息。預(yù)處理的目標(biāo)是將其轉(zhuǎn)化為適合建模的“干凈”特征同時(shí)盡可能保留和增強(qiáng)有用信息。異常值處理不要武斷地刪除。先分析異常值產(chǎn)生的原因。例如銷量為0的日子是門店關(guān)門節(jié)假日、還是缺貨如果是節(jié)假日這本身就是一個(gè)重要的特征應(yīng)該用“是否節(jié)假日”的啞變量來標(biāo)記而不是刪除。對(duì)于明顯是錄入錯(cuò)誤的極端值如成本為負(fù)才用箱線圖或3σ原則進(jìn)行剔除或修正。缺失值填補(bǔ)時(shí)間序列的缺失值不宜用整體均值填補(bǔ)。應(yīng)采用前向填充、線性插值或者基于同類產(chǎn)品同期數(shù)據(jù)進(jìn)行填補(bǔ)。例如菠菜某天數(shù)據(jù)缺失可以用同期的油菜、生菜的數(shù)據(jù)趨勢(shì)來估算。特征工程這是拉開差距的地方。除了從日期衍生出年、月、日、周幾、是否節(jié)假日、是否周末等還可以構(gòu)建滯后特征前1天、前7天、前30天的銷量、價(jià)格。滾動(dòng)統(tǒng)計(jì)特征過去7天的平均銷量、銷量標(biāo)準(zhǔn)差、平均價(jià)格。交互特征自身價(jià)格與競(jìng)品替代品平均價(jià)格的比值。外部特征如果允許雖然題目沒給但可以合理推斷并作為模型假設(shè)。例如可以引入公開的天氣數(shù)據(jù)溫度、濕度作為預(yù)測(cè)特征并在論文中說明其合理性。數(shù)據(jù)標(biāo)準(zhǔn)化/歸一化對(duì)于涉及距離計(jì)算的聚類或神經(jīng)網(wǎng)絡(luò)模型是必須的。對(duì)于樹模型如LightGBM則不一定需要。但價(jià)格和銷量單位差異大進(jìn)行適當(dāng)?shù)目s放如Min-Max縮放有助于優(yōu)化求解時(shí)的數(shù)值穩(wěn)定性。實(shí)操心得建立一個(gè)可復(fù)現(xiàn)的數(shù)據(jù)預(yù)處理流水線Pipeline。使用Python的sklearn.pipeline或自定義函數(shù)確保從原始數(shù)據(jù)到建模數(shù)據(jù)的每一步都是透明的。這不僅能保證一致性在調(diào)整模型時(shí)也能快速重新生成特征。3.2 需求預(yù)測(cè)模型構(gòu)建實(shí)戰(zhàn)我們以構(gòu)建一個(gè)融合模型為例演示具體步驟。步驟1數(shù)據(jù)劃分不要隨機(jī)劃分時(shí)間序列數(shù)據(jù)必須按時(shí)間順序劃分。例如用前3年的數(shù)據(jù)作為訓(xùn)練集最后1年的數(shù)據(jù)作為測(cè)試集來模擬“用歷史預(yù)測(cè)未來”的真實(shí)場(chǎng)景。步驟2基準(zhǔn)模型建立先建立幾個(gè)簡(jiǎn)單的基準(zhǔn)模型用于后續(xù)對(duì)比樸素模型預(yù)測(cè)值等于前一天的值。季節(jié)性樸素模型預(yù)測(cè)值等于上一周同一天的值。簡(jiǎn)單移動(dòng)平均過去7天的平均值。 計(jì)算這些基準(zhǔn)模型在測(cè)試集上的均方根誤差RMSE和平均絕對(duì)百分比誤差MAPE。步驟3高級(jí)模型構(gòu)建與融合模型1SARIMA。使用statsmodels庫通過ACF/PACF圖確定參數(shù)范圍用網(wǎng)格搜索尋找最優(yōu)的(p,d,q)(P,D,Q,s)參數(shù)組合。SARIMA能很好地捕捉固定趨勢(shì)和季節(jié)模式。模型2LightGBM。將上一步構(gòu)建的所有特征滯后、滾動(dòng)、日期特征等作為輸入。需要注意防止過擬合合理設(shè)置num_leaves,learning_rate, 并啟用早停法early_stopping。模型融合簡(jiǎn)單有效的融合方法是加權(quán)平均。權(quán)重如何確定可以在測(cè)試集或單獨(dú)劃出一個(gè)驗(yàn)證集上以最小化預(yù)測(cè)誤差為目標(biāo)求解最優(yōu)權(quán)重。更高級(jí)的方法可以用線性回歸以各模型預(yù)測(cè)值為特征真實(shí)值為目標(biāo)來學(xué)習(xí)融合權(quán)重。步驟4不確定性量化對(duì)于優(yōu)化而言點(diǎn)預(yù)測(cè)不夠我們需要知道預(yù)測(cè)的分布。有兩種常用方法分位數(shù)回歸使用LightGBM的objectivequantile分別訓(xùn)練多個(gè)分位數(shù)如0.05, 0.5, 0.95的模型從而得到預(yù)測(cè)區(qū)間。Bootstrap法對(duì)訓(xùn)練數(shù)據(jù)進(jìn)行有放回抽樣構(gòu)建多個(gè)不同的預(yù)測(cè)模型用這些模型的預(yù)測(cè)結(jié)果集合來近似分布。# 示例使用LightGBM進(jìn)行分位數(shù)回歸預(yù)測(cè)核心代碼片段 import lightgbm as lgb # 準(zhǔn)備數(shù)據(jù) train_data lgb.Dataset(X_train, labely_train) # 訓(xùn)練中位數(shù)回歸模型0.5分位數(shù) params { objective: quantile, alpha: 0.5, # 分位數(shù) metric: quantile, boosting_type: gbdt, num_leaves: 31, learning_rate: 0.05, feature_fraction: 0.9 } model_median lgb.train(params, train_data) # 訓(xùn)練上界回歸模型0.95分位數(shù) params[alpha] 0.95 model_upper lgb.train(params, train_data) # 訓(xùn)練下界回歸模型0.05分位數(shù) params[alpha] 0.05 model_lower lgb.train(params, train_data) # 預(yù)測(cè) pred_median model_median.predict(X_test) pred_upper model_upper.predict(X_test) pred_lower model_lower.predict(X_test) # 此時(shí)對(duì)于每個(gè)預(yù)測(cè)點(diǎn)我們得到了一個(gè)區(qū)間 [pred_lower, pred_upper]3.3 聯(lián)合優(yōu)化模型建模與求解假設(shè)我們采用樣本平均近似SAA方法處理不確定性。模型建立設(shè)決策變量為每日對(duì)每個(gè)品類i的補(bǔ)貨量 (Q_i)定價(jià) (P_i)。 目標(biāo)最大化期望利潤(rùn)。 約束庫存平衡約束、容量約束、非負(fù)約束等。由于需求 (D_i) 不確定我們采用SAA。假設(shè)我們從預(yù)測(cè)分布中抽取了N個(gè)需求場(chǎng)景 (D_i^1, D_i^2, ..., D_i^N)。則確定性等價(jià)模型可以寫成 最大化( \frac{1}{N} \sum_{n1}^{N} \sum_{i} [ P_i * min(Q_i I_{i,prev}, D_i^n) - C_i*Q_i - h * I_i^n - v * L_i^n ] ) 其中(I_{i,prev}) 是上日庫存。(min(Q_i I_{i,prev}, D_i^n)) 表示在場(chǎng)景n下的實(shí)際銷售量不能超過需求。(C_i) 是進(jìn)貨成本。(h) 是單位庫存持有成本。(I_i^n) 是場(chǎng)景n下的期末庫存。(v) 是單位損耗成本。(L_i^n) 是場(chǎng)景n下的損耗量這里可以建模為 (L_i^n \theta * I_i^n)簡(jiǎn)化線性損耗或更復(fù)雜的函數(shù)。還需要加上庫存平衡方程(I_i^n max(Q_i I_{i,prev} - D_i^n, 0))。模型求解這是一個(gè)大規(guī)模的線性/非線性規(guī)劃問題取決于損耗函數(shù)的形式?？梢允褂脤I(yè)的優(yōu)化求解器如Gurobi、CPLEX學(xué)術(shù)免費(fèi)或開源的PuLP調(diào)用CBC求解器來求解。# 示例使用PuLP定義和求解一個(gè)簡(jiǎn)化的確定性補(bǔ)貨模型單品類需求已知 import pulp # 創(chuàng)建問題 prob pulp.LpProblem(Single_Product_Inventory, pulp.LpMaximize) # 參數(shù) demand 100 # 預(yù)測(cè)需求 cost 2 # 單位進(jìn)貨成本 price 5 # 售價(jià) holding_cost 0.1 # 單位持有成本 initial_inventory 20 # 決策變量 order_qty pulp.LpVariable(OrderQuantity, lowBound0, catContinuous) sales pulp.LpVariable(Sales, lowBound0, catContinuous) ending_inventory pulp.LpVariable(EndingInventory, lowBound0, catContinuous) # 目標(biāo)函數(shù)利潤(rùn) 銷售收入 - 進(jìn)貨成本 - 持有成本 prob price * sales - cost * order_qty - holding_cost * ending_inventory # 約束條件 prob sales demand, SalesCannotExceedDemand prob sales initial_inventory order_qty, SalesCannotExceedAvailable prob ending_inventory initial_inventory order_qty - sales, InventoryBalance # 求解 prob.solve(pulp.PULP_CBC_CMD(msgFalse)) print(f狀態(tài): {pulp.LpStatus[prob.status]}) print(f最優(yōu)補(bǔ)貨量: {order_qty.varValue:.2f}) print(f預(yù)期利潤(rùn): {pulp.value(prob.objective):.2f})對(duì)于多品類、多場(chǎng)景的復(fù)雜模型決策變量和約束會(huì)呈指數(shù)級(jí)增長(zhǎng)需要仔細(xì)設(shè)計(jì)模型結(jié)構(gòu)并可能需要借助更強(qiáng)大的求解器或分解算法。4. 常見問題、調(diào)試技巧與進(jìn)階思考在實(shí)際操作中你會(huì)遇到各種各樣的問題。下面是我總結(jié)的一些典型問題及其解決思路。4.1 模型求解速度慢或無法求解這是最常遇到的問題。問題根源模型規(guī)模太大場(chǎng)景太多、品類太多、非線性約束導(dǎo)致非凸、整數(shù)變量導(dǎo)致組合爆炸。排查與解決簡(jiǎn)化模型先從最簡(jiǎn)化的版本開始。例如先考慮單品類、需求為確定值、損耗為線性函數(shù)的模型。確保這個(gè)簡(jiǎn)單模型能快速求解。減少場(chǎng)景數(shù)SAA中場(chǎng)景數(shù)N不需要極大。可以通過場(chǎng)景縮減技術(shù)用少量具有代表性的場(chǎng)景如通過聚類生成來近似原分布。線性化/凸化嘗試將非線性約束如包含max/min的損耗函數(shù)用線性約束和大M法進(jìn)行等價(jià)轉(zhuǎn)換。如果目標(biāo)函數(shù)或約束是非凸的求解會(huì)非常困難需要考慮凸近似或改用啟發(fā)式算法。松弛整數(shù)變量如果補(bǔ)貨量被建模為整數(shù)變量箱/包可以先松弛為連續(xù)變量求解得到最優(yōu)解后再對(duì)結(jié)果進(jìn)行取整并評(píng)估取整帶來的性能損失。如果損失可接受則可行。使用啟發(fā)式算法當(dāng)精確算法失效時(shí)轉(zhuǎn)向遺傳算法、模擬退火等。雖然不能保證全局最優(yōu)但一個(gè)好的啟發(fā)式算法通常能找到高質(zhì)量的解。關(guān)鍵是要設(shè)計(jì)好的編碼方式、適應(yīng)度函數(shù)和鄰域搜索操作。4.2 模型結(jié)果不理想或違反常識(shí)優(yōu)化結(jié)果出來了但補(bǔ)貨量是負(fù)數(shù)或者定價(jià)高得離譜這顯然是模型或數(shù)據(jù)有問題。問題根源約束條件不完整或錯(cuò)誤目標(biāo)函數(shù)權(quán)重設(shè)置不合理輸入數(shù)據(jù)存在量綱不一致或異常值。排查與解決檢查約束逐條檢查約束條件是否反映了現(xiàn)實(shí)限制。例如是否遺漏了“補(bǔ)貨量非負(fù)”、“價(jià)格不低于成本”這樣的基本約束庫存平衡約束的符號(hào)是否正確分析對(duì)偶變量/影子價(jià)格在LP/MILP求解后求解器會(huì)給出每個(gè)約束的影子價(jià)格。它表示該約束右邊常數(shù)每增加一個(gè)單位目標(biāo)函數(shù)值的變化量。如果某個(gè)約束的影子價(jià)格非常高說明該約束是“緊”的是當(dāng)前利潤(rùn)的主要瓶頸。這能幫你理解模型的限制在哪里。進(jìn)行敏感性分析系統(tǒng)地改變關(guān)鍵參數(shù)如需求預(yù)測(cè)的均值、損耗率系數(shù)觀察最優(yōu)解如何變化。如果最優(yōu)解對(duì)某個(gè)參數(shù)極其敏感那么你需要花更多精力去校準(zhǔn)這個(gè)參數(shù)或者說明你的策略在該參數(shù)不確定時(shí)風(fēng)險(xiǎn)很高?；厮莸筋A(yù)測(cè)階段檢查輸入給優(yōu)化模型的需求預(yù)測(cè)值是否合理。用歷史數(shù)據(jù)對(duì)比看預(yù)測(cè)是否存在系統(tǒng)性偏差。4.3 論文寫作與呈現(xiàn)中的痛點(diǎn)痛點(diǎn)一模型描述過于晦澀。技巧用“比喻公式文字”三重描述。例如在介紹庫存平衡約束時(shí)可以寫“這就像是一個(gè)水池庫存每天有進(jìn)水補(bǔ)貨也有出水銷售剩下的就是池子里的水期末庫存。用公式表示就是期末庫存 期初庫存 補(bǔ)貨量 - 銷售量。”痛點(diǎn)二結(jié)果分析流于表面。技巧不要只說“模型提升了利潤(rùn)”。要深入分析為什么提升了。是減少了特定品類在特定季節(jié)的損耗還是通過動(dòng)態(tài)定價(jià)在需求高時(shí)多賺了錢結(jié)合具體的圖表數(shù)據(jù)給出有因果關(guān)系的解釋。例如“如圖5所示我們的模型在夏季對(duì)綠葉蔬菜采用了更激進(jìn)的降價(jià)策略相比基準(zhǔn)策略平均低8%這雖然降低了單價(jià)但使得銷量提升了25%最終該品類在夏季的總利潤(rùn)提升了15%同時(shí)損耗降低了30%?！蓖袋c(diǎn)三代碼、模型、論文脫節(jié)。技巧建立一條清晰的“數(shù)據(jù)流水線”。從原始數(shù)據(jù)到最終結(jié)果的每一個(gè)關(guān)鍵步驟都用腳本自動(dòng)化。確保論文中的每一個(gè)數(shù)字、每一張圖都能由你的代碼重新生成。這不僅能避免低級(jí)錯(cuò)誤在回答評(píng)委質(zhì)疑時(shí)也更有底氣。4.4 從優(yōu)秀到卓越可以深入的進(jìn)階方向如果你已經(jīng)掌握了上述基本框架還想讓論文更出彩可以考慮以下進(jìn)階方向考慮競(jìng)爭(zhēng)因素題目假設(shè)是單家超市。但在現(xiàn)實(shí)中定價(jià)會(huì)受到周邊競(jìng)爭(zhēng)對(duì)手的影響?？梢砸氩┺恼摰乃枷爰僭O(shè)競(jìng)爭(zhēng)對(duì)手也采用類似策略建立動(dòng)態(tài)定價(jià)博弈模型。引入消費(fèi)者行為模型需求不僅僅是價(jià)格的函數(shù)還可能受到促銷、陳列方式的影響?？梢試L試引入基礎(chǔ)的消費(fèi)者選擇模型如Logit模型讓需求預(yù)測(cè)更精細(xì)。多目標(biāo)優(yōu)化除了利潤(rùn)企業(yè)可能還關(guān)心客戶滿意度缺貨率低、供應(yīng)鏈穩(wěn)定性補(bǔ)貨量波動(dòng)小?？梢詫栴}建模為多目標(biāo)優(yōu)化使用帕累托前沿來展示不同目標(biāo)之間的權(quán)衡。強(qiáng)化學(xué)習(xí)這是一個(gè)非常前沿的嘗試。將定價(jià)和補(bǔ)貨決策視為一個(gè)序列決策問題超市是智能體利潤(rùn)是獎(jiǎng)勵(lì)使用深度強(qiáng)化學(xué)習(xí)算法來學(xué)習(xí)最優(yōu)策略。這種方法對(duì)數(shù)據(jù)量和算力要求高但一旦實(shí)現(xiàn)會(huì)非常亮眼。在論文中你需要詳細(xì)說明狀態(tài)空間、動(dòng)作空間、獎(jiǎng)勵(lì)函數(shù)的設(shè)計(jì)以及訓(xùn)練過程。最后我想強(qiáng)調(diào)的是數(shù)學(xué)建模競(jìng)賽和解決真實(shí)商業(yè)問題的核心是相通的理解問題本質(zhì)做出合理假設(shè)利用數(shù)學(xué)工具構(gòu)建簡(jiǎn)化但有效的模型并用數(shù)據(jù)驅(qū)動(dòng)決策。2023年C題的這些優(yōu)秀論文正是這一理念的完美體現(xiàn)。它們的故事告訴我們成功不在于用了多么高深的算法而在于是否構(gòu)建了一個(gè)自洽、合理、且能產(chǎn)生實(shí)際價(jià)值的解決方案閉環(huán)。當(dāng)你下次面對(duì)一個(gè)復(fù)雜問題時(shí)不妨也試試這種“先拆解、再建模、后驗(yàn)證”的思路或許會(huì)有意想不到的收獲。