1. 從一道題開始標(biāo)準(zhǔn)規(guī)劃問題到底是什么如果你正在準(zhǔn)備數(shù)學(xué)建模競賽或者剛剛開始接觸運(yùn)籌優(yōu)化那么“標(biāo)準(zhǔn)規(guī)劃問題”這個詞一定不陌生。但很多時(shí)候我們只是機(jī)械地套用MATLAB里的linprog或fmincon函數(shù)把系數(shù)矩陣填進(jìn)去然后祈禱得到一個正確的答案。至于為什么這么填、函數(shù)背后在做什么、結(jié)果不理想時(shí)該怎么辦往往是一頭霧水。今天我想結(jié)合自己幾次數(shù)模競賽和實(shí)際項(xiàng)目中的踩坑經(jīng)歷來聊聊標(biāo)準(zhǔn)規(guī)劃問題的MATLAB求解重點(diǎn)不是“怎么用”而是“為什么這么用”以及“用的時(shí)候要注意什么”。所謂“標(biāo)準(zhǔn)規(guī)劃問題”在數(shù)學(xué)建模的語境下通常指那些具有標(biāo)準(zhǔn)數(shù)學(xué)形式的優(yōu)化問題。最常見的就是線性規(guī)劃Linear Programming, LP它的標(biāo)準(zhǔn)形式是求一組決策變量在滿足一系列線性等式或不等式約束的條件下使一個線性目標(biāo)函數(shù)達(dá)到最小或最大。比如經(jīng)典的資源分配、生產(chǎn)計(jì)劃、運(yùn)輸問題都可以歸結(jié)為此類。除此之外還有整數(shù)規(guī)劃IP、**二次規(guī)劃QP**等它們都有各自的標(biāo)準(zhǔn)形式。MATLAB的優(yōu)化工具箱為我們提供了針對這些標(biāo)準(zhǔn)形式的求解器但工具箱不是黑箱理解其輸入輸出的“規(guī)矩”是高效、準(zhǔn)確求解的第一步。很多人拿到一個問題比如“如何安排生產(chǎn)使得利潤最大”會直接去想MATLAB代碼怎么寫。我的建議恰恰相反先忘掉MATLAB拿起筆和紙。第一步也是最重要的一步是數(shù)學(xué)建模即把現(xiàn)實(shí)問題抽象為標(biāo)準(zhǔn)規(guī)劃問題的數(shù)學(xué)形式。這一步?jīng)Q定了你后面所有代碼的骨架。一個清晰的數(shù)學(xué)模型應(yīng)該明確決策變量是什么有幾個分別代表什么目標(biāo)函數(shù)是什么是求最大還是最小表達(dá)式是什么約束條件有哪些是等式還是不等式表達(dá)式是什么決策變量的取值范圍如何是否要求非負(fù)是否是整數(shù)只有把這些都用數(shù)學(xué)符號清晰地表達(dá)出來你才能準(zhǔn)確地將它們“翻譯”成MATLAB求解器能聽懂的語言。2. 線性規(guī)劃linprog的“標(biāo)準(zhǔn)姿勢”與常見陷阱線性規(guī)劃是基礎(chǔ)MATLAB中對應(yīng)的函數(shù)是linprog。它的語法看似簡單但參數(shù)順序和形式有嚴(yán)格規(guī)定這也是新手最容易出錯的地方。2.1linprog的標(biāo)準(zhǔn)形式與參數(shù)映射MATLAB的linprog求解的是如下標(biāo)準(zhǔn)最小化形式min f^T * x subject to: A * x b Aeq * x beq lb x ub其中f是目標(biāo)函數(shù)的系數(shù)列向量x是決策變量向量。A和b對應(yīng)線性不等式約束Aeq和beq對應(yīng)線性等式約束lb和ub是變量的下界和上界。這里第一個關(guān)鍵點(diǎn)就來了你的模型必須轉(zhuǎn)換成這個形式。如果你的原始問題是最大化max那么只需要將目標(biāo)函數(shù)系數(shù)取相反數(shù)轉(zhuǎn)化為最小化問題。例如max 3x1 4x2等價(jià)于min -3x1 -4x2。最終linprog返回的最優(yōu)解x是一樣的但最優(yōu)值fval需要你再取反才能得到原始的最大化目標(biāo)值。第二個關(guān)鍵點(diǎn)是約束的方向。linprog默認(rèn)的不等式是“小于等于”。如果你的約束是“大于等于”比如2x1 x2 10那么需要在不等式兩邊同時(shí)乘以-1轉(zhuǎn)化為-2x1 - x2 -10。這一步轉(zhuǎn)換必須在構(gòu)造矩陣A和向量b之前完成。讓我們看一個簡單的例子某工廠生產(chǎn)兩種產(chǎn)品A和B需要兩道工序。生產(chǎn)一件A需工序一2小時(shí)工序二1小時(shí)生產(chǎn)一件B需工序一1小時(shí)工序二2小時(shí)。工序一每天可用12小時(shí)工序二每天可用9小時(shí)。產(chǎn)品A利潤3元B利潤4元。問如何安排生產(chǎn)使利潤最大建模設(shè)生產(chǎn)A產(chǎn)品x1件B產(chǎn)品x2件。目標(biāo)max z 3*x1 4*x2約束工序一2*x1 x2 12工序二x1 2*x2 9非負(fù)x1 0, x2 0轉(zhuǎn)換為linprog標(biāo)準(zhǔn)形式目標(biāo)由于linprog求最小所以f [-3; -4]。不等式約束恰好是“”所以A [2, 1; 1, 2],b [12; 9]。等式約束無Aeq [],beq []。下界lb [0; 0]上界默認(rèn)為無窮大inf。MATLAB代碼實(shí)現(xiàn)f [-3; -4]; % 目標(biāo)函數(shù)系數(shù)注意負(fù)號 A [2, 1; 1, 2]; b [12; 9]; Aeq []; beq []; lb [0; 0]; [x, fval, exitflag, output] linprog(f, A, b, Aeq, beq, lb); optimal_profit -fval; % 記得把最小化值取反得到最大利潤 disp([最優(yōu)生產(chǎn)計(jì)劃A生產(chǎn) , num2str(x(1)), 件 B生產(chǎn) , num2str(x(2)), 件]); disp([最大利潤為, num2str(optimal_profit), 元]); disp([求解器退出狀態(tài), num2str(exitflag)]); disp(output.message);2.2 解讀輸出exitflag比結(jié)果更重要運(yùn)行上面的代碼你會得到解x和最優(yōu)值fval。但請務(wù)必養(yǎng)成查看exitflag和output信息的習(xí)慣exitflag告訴你求解是否成功以及原因這比單純看一個數(shù)值解重要得多。exitflag 0求解器收斂到一個最優(yōu)解。這是最理想的情況。exitflag 0求解器達(dá)到了最大迭代次數(shù)或函數(shù)計(jì)算次數(shù)限制可能還沒找到最優(yōu)解。這時(shí)你需要檢查結(jié)果是否合理或者通過options參數(shù)增加迭代次數(shù)options optimoptions(linprog, MaxIterations, 10000)。exitflag 0問題無解不可行或無界。這是建?；驍?shù)據(jù)錯誤的高發(fā)區(qū)。-2問題不可行No feasible point found。意味著你給出的約束條件互相矛盾沒有任何一個點(diǎn)能同時(shí)滿足所有約束。比如你要求x1 x2 5同時(shí)又要求x1 x2 10。這時(shí)你需要回頭檢查模型和數(shù)據(jù)的邏輯。-3問題無界Unbounded。在最小化問題中目標(biāo)函數(shù)值可以趨向負(fù)無窮在最大化問題中可以趨向正無窮。通常是因?yàn)榧s束不夠允許決策變量無限增大/減小而不違反約束。例如求min -x1 - x2約束只有x1 0, x2 0那么x1和x2可以無限大目標(biāo)函數(shù)值就無限小。踩坑實(shí)錄在一次比賽中我們模型跑出來結(jié)果好得離譜利潤高到不可思議。當(dāng)時(shí)只顧著高興沒看exitflag。直到最后檢查時(shí)才發(fā)現(xiàn)exitflag -3問題無界原因是我們在轉(zhuǎn)化一個資源約束時(shí)不小心把“”寫成了“”導(dǎo)致約束方向反了相當(dāng)于資源可以無限使用。這個教訓(xùn)讓我銘記永遠(yuǎn)不要相信沒有經(jīng)過exitflag驗(yàn)證的“好結(jié)果”。2.3 處理無可行解與不可行診斷當(dāng)exitflag -2時(shí)如何快速定位是哪個或哪組約束導(dǎo)致了不可行MATLAB沒有內(nèi)置的直接工具但我們可以用一些技巧來診斷。一種實(shí)用的方法是逐步放松約束法。如果你的模型有m個不等式約束可以嘗試每次注釋掉一個或一組約束然后重新求解。如果注釋掉某個約束后問題變得可行了那么這個約束很可能就是導(dǎo)致沖突的“元兇”之一。你需要仔細(xì)檢查這個約束的數(shù)學(xué)表達(dá)式和數(shù)據(jù)是否準(zhǔn)確。另一種思路是引入松弛變量Slack Variables或使用不可行性最小化。但這通常更復(fù)雜。對于競賽或初級應(yīng)用逐步放松法是最直觀的調(diào)試手段。這本質(zhì)上是在模擬“如果這個條件不那么嚴(yán)格是不是就有解了”的過程能幫你快速理解約束之間的沖突關(guān)系。3. 整數(shù)規(guī)劃當(dāng)決策變量不能“分割”時(shí)現(xiàn)實(shí)中的很多問題決策變量必須是整數(shù)。比如生產(chǎn)多少臺設(shè)備不能是半臺、派遣多少輛卡車、某個地點(diǎn)是否建廠0-1決策。這就是整數(shù)規(guī)劃IP特別是0-1規(guī)劃。MATLAB中使用intlinprog函數(shù)求解混合整數(shù)線性規(guī)劃MILP。3.1intlinprog的核心指定整數(shù)變量索引intlinprog的語法和linprog非常相似多了一個關(guān)鍵參數(shù)intcon用于指定哪些決策變量必須是整數(shù)。intcon是一個向量包含整數(shù)變量的索引。例如在之前的工廠問題中如果我們要求生產(chǎn)的產(chǎn)品數(shù)量必須是整數(shù)件很合理那么x1和x2都必須是整數(shù)。假設(shè)x [x1; x2]那么intcon [1; 2]。代碼修改如下f [-3; -4]; A [2, 1; 1, 2]; b [12; 9]; lb [0; 0]; intcon [1, 2]; % x1和x2都是整數(shù)變量 % 注意intlinprog的參數(shù)順序f, intcon, A, b, Aeq, beq, lb, ub [x, fval, exitflag] intlinprog(f, intcon, A, b, [], [], lb); optimal_profit -fval;你會發(fā)現(xiàn)最優(yōu)解從原來的(5, 2)非整數(shù)解變成了一個整數(shù)解比如(4, 2)或(5, 1)具體取決于算法分支總利潤也會相應(yīng)變化。整數(shù)規(guī)劃的最優(yōu)值通常不會優(yōu)于對于最小化問題是不會低于對應(yīng)的線性規(guī)劃松弛問題即去掉整數(shù)限制后的問題的最優(yōu)值。這是理解整數(shù)規(guī)劃性質(zhì)的一個要點(diǎn)。3.2 0-1規(guī)劃建模的巧妙之處0-1變量是整數(shù)變量的特例只能取0或1常用于表示“是/否”、“開/關(guān)”、“選擇/不選擇”這類邏輯決策。intlinprog同樣可以處理只需將變量的上界ub設(shè)為1下界lb設(shè)為0并將其索引加入intcon即可。0-1規(guī)劃的難點(diǎn)和魅力在于建模。如何用線性約束來表達(dá)復(fù)雜的邏輯關(guān)系這里分享幾個經(jīng)典技巧互斥選擇從N個項(xiàng)目中至多選擇K個。設(shè)x_i為0-1變量表示是否選擇項(xiàng)目i。約束可寫為sum(x_i) K。依賴關(guān)系如果選擇項(xiàng)目B則必須選擇項(xiàng)目A。約束為x_B x_A。這意味著當(dāng)x_B1時(shí)x_A也必須為1但當(dāng)x_A1時(shí)x_B可以為0。打包關(guān)系項(xiàng)目A和項(xiàng)目B必須同時(shí)選擇或同時(shí)不選。約束為x_A x_B。固定成本生產(chǎn)某種產(chǎn)品如果生產(chǎn)x0則除了可變成本外還需支付一筆固定成本F。這需要用到一個輔助0-1變量y。設(shè)x為產(chǎn)量M為一個足夠大的數(shù)上界。x M * y如果y0則x必須為0如果y1x可以大于0但受M限制目標(biāo)函數(shù)中加入F * y。 這樣只要x0y就會被“激活”為1從而在目標(biāo)函數(shù)中計(jì)入固定成本F。這個技巧稱為“大M法”是混合整數(shù)規(guī)劃建模的核心技巧之一。選擇M的值需要小心既要足夠大以保證約束有效當(dāng)y1時(shí)x不受此約束限制又不能太大否則會導(dǎo)致數(shù)值計(jì)算困難影響求解速度和穩(wěn)定性。通常取變量x的一個合理的上界即可。經(jīng)驗(yàn)之談處理0-1規(guī)劃或一般整數(shù)規(guī)劃時(shí)求解時(shí)間可能遠(yuǎn)超線性規(guī)劃。intlinprog提供了options參數(shù)來調(diào)整求解器行為比如設(shè)置最大求解時(shí)間MaxTime或相對容差RelativeGapTolerance。在數(shù)模競賽中如果問題規(guī)模較大可以在保證結(jié)果合理性的前提下適當(dāng)放寬RelativeGapTolerance例如設(shè)為0.01或0.05讓求解器在找到可行解并證明其與最優(yōu)解的差距在1%或5%以內(nèi)時(shí)就停止以節(jié)省寶貴時(shí)間。4. 非線性規(guī)劃入門fmincon的靈活與復(fù)雜當(dāng)目標(biāo)函數(shù)或約束條件中出現(xiàn)了非線性項(xiàng)如平方、指數(shù)、三角函數(shù)或變量相乘我們就進(jìn)入了非線性規(guī)劃NLP的領(lǐng)域。MATLAB中功能最強(qiáng)大的通用非線性規(guī)劃求解器是fmincon。它的靈活性很高但設(shè)置也更為復(fù)雜。4.1 從線性到非線性思維轉(zhuǎn)換使用linprog時(shí)我們把所有系數(shù)塞進(jìn)矩陣和向量就行了。但fmincon要求我們以函數(shù)句柄Function Handle的形式來提供目標(biāo)函數(shù)和非線性約束。這意味著你需要單獨(dú)編寫一個或多個MATLAB函數(shù)文件或匿名函數(shù)來計(jì)算這些值。fmincon求解的問題形式一般如下min f(x) subject to: c(x) 0 非線性不等式約束 ceq(x) 0 非線性等式約束 A*x b, Aeq*x beq 線性約束 lb x ub4.2 實(shí)戰(zhàn)一個帶非線性約束的簡單例子假設(shè)我們要優(yōu)化一個簡單問題最小化f(x) x1^2 x2^2約束為x1*x2 1且x1 0, x2 0。轉(zhuǎn)換為fmincon標(biāo)準(zhǔn)形式目標(biāo)函數(shù)f (x) x(1)^2 x(2)^2非線性約束x1*x2 1需要寫成c(x) 0的形式-x1*x2 1 0。所以c (x) -x(1)*x(2) 1。線性約束無非負(fù)約束但我們可以用下界lb表示lb [0; 0]。編寫代碼% 定義目標(biāo)函數(shù)使用匿名函數(shù) objective (x) x(1)^2 x(2)^2; % 定義初始點(diǎn)非常重要非線性規(guī)劃求解結(jié)果嚴(yán)重依賴初始點(diǎn) x0 [2; 2]; % 選擇一個可行的初始點(diǎn)例如(2,2)滿足 x1*x241 % 定義線性約束本例沒有用空數(shù)組 A []; b []; Aeq []; beq []; % 定義變量邊界 lb [0; 0]; ub []; % 無上界 % 定義非線性約束單獨(dú)寫一個函數(shù)或者用匿名函數(shù) % 這里c(x) 0, ceq(x) 0 nonlcon (x) deal(-x(1)*x(2) 1, []); % deal函數(shù)返回兩個輸出c和ceq % 調(diào)用fmincon求解 options optimoptions(fmincon, Display, iter); % 顯示迭代過程便于調(diào)試 [x_opt, fval_opt, exitflag, output] fmincon(objective, x0, A, b, Aeq, beq, lb, ub, nonlcon, options); disp(最優(yōu)解); disp(x_opt); disp(最優(yōu)目標(biāo)值); disp(fval_opt); disp(退出狀態(tài)); disp(exitflag);4.3 初始點(diǎn)選擇與算法選項(xiàng)決定成敗的細(xì)節(jié)對于非線性規(guī)劃初始點(diǎn)x0的選擇至關(guān)重要。fmincon使用基于梯度的局部搜索算法如內(nèi)點(diǎn)法、序列二次規(guī)劃SQP等它只能找到從初始點(diǎn)出發(fā)所能到達(dá)的局部最優(yōu)解而不一定是全局最優(yōu)解。不同的初始點(diǎn)可能導(dǎo)致完全不同的結(jié)果。策略1根據(jù)物理意義或經(jīng)驗(yàn)猜測一個可能接近最優(yōu)解的點(diǎn)作為初始點(diǎn)。策略2如果問題可行域不大可以在可行域內(nèi)隨機(jī)生成多個初始點(diǎn)分別求解然后取目標(biāo)函數(shù)值最好的那個解作為最終結(jié)果。這是一種簡單的“多起點(diǎn)”策略有助于避免糟糕的局部最優(yōu)。策略3對于復(fù)雜問題可以考慮使用全局優(yōu)化算法如GlobalSearch或MultiStart它們會在fmincon的基礎(chǔ)上進(jìn)行多次隨機(jī)起始點(diǎn)的搜索。但這會消耗更多計(jì)算時(shí)間。fmincon的options參數(shù)也非常豐富常用的有Algorithm: 選擇求解算法如interior-point內(nèi)點(diǎn)法默認(rèn)、sqp序列二次規(guī)劃、active-set等。對于不同的問題算法效率可能不同。如果不確定保持默認(rèn)或嘗試sqp。Display: 控制輸出信息iter顯示每次迭代信息調(diào)試用final只顯示最終結(jié)果off不顯示。MaxIterations,MaxFunctionEvaluations: 設(shè)置最大迭代次數(shù)和函數(shù)計(jì)算次數(shù)防止程序在復(fù)雜問題上無休止運(yùn)行。OptimalityTolerance,StepTolerance,ConstraintTolerance: 設(shè)置優(yōu)化的終止容差。通常默認(rèn)值即可如果求解器提前終止或結(jié)果精度不夠可以適當(dāng)調(diào)小這些值例如1e-8。踩坑實(shí)錄曾經(jīng)求解一個工程優(yōu)化問題目標(biāo)函數(shù)有多個“山谷”。第一次隨便設(shè)了個初始點(diǎn)[0,0]結(jié)果收斂到了一個很差的局部最優(yōu)。后來分析了問題背景知道最優(yōu)解大概在某個范圍將初始點(diǎn)改為[5,5]立刻得到了一個好得多的解。所以對于非線性問題永遠(yuǎn)不要忽視初始點(diǎn)的選擇它甚至比調(diào)參更重要。如果結(jié)果不理想換個初始點(diǎn)再試試是最簡單有效的排查方法之一。5. 求解器“報(bào)錯”怎么辦典型問題排查指南在使用MATLAB優(yōu)化工具箱時(shí)你肯定會遇到各種錯誤信息或警告。下面是一些常見問題的排查思路。5.1 “Solver stopped prematurely” 或 “No feasible solution found”這通常意味著求解器在給定的迭代次數(shù)或時(shí)間內(nèi)沒有找到可行解或收斂。檢查模型可行性首先用2.3節(jié)的方法檢查約束是否可能互相矛盾。嘗試放松一些約束比如增大b值或減小A值看問題是否變得可行。調(diào)整求解器選項(xiàng)增加MaxIterations和MaxFunctionEvaluations。對于fmincon還可以嘗試調(diào)整算法Algorithm。縮放問題如果決策變量的數(shù)量級相差巨大例如x1在0~1之間x2在0~100000之間可能會導(dǎo)致數(shù)值計(jì)算困難。盡量對變量進(jìn)行縮放使它們處于相近的數(shù)量級比如0~10或0~100。檢查初始點(diǎn)僅對fmincon確保初始點(diǎn)x0滿足所有約束或至少滿足線性約束和邊界。fmincon對于初始點(diǎn)的可行性有一定要求特別是使用某些算法時(shí)。可以嘗試多個不同的初始點(diǎn)。5.2 結(jié)果不理想或違反直覺求解器給出了一個解但你覺得這個解很奇怪或者目標(biāo)函數(shù)值比你預(yù)想的差很多。驗(yàn)證exitflag確認(rèn)exitflag 0表明求解器是正常收斂的。檢查解是否滿足約束手動將求得的解x_opt代入你的約束條件中計(jì)算看看是否真的滿足所有約束在容差范圍內(nèi)。有時(shí)候數(shù)值計(jì)算會引入微小誤差但大的違反一定有問題。檢查模型是否正確這是最根本的。重新審視你的數(shù)學(xué)建模過程檢查目標(biāo)函數(shù)系數(shù)、約束矩陣A、Aeq、向量b、beq的每一個元素是否填寫正確。一個常見的錯誤是矩陣的維度不匹配或者系數(shù)正負(fù)號弄反。對于非線性問題嘗試多個不同的初始點(diǎn)看是否能得到更好的解?？赡苣愕暨M(jìn)了一個局部最優(yōu)的“坑”里。5.3 性能問題求解太慢對于整數(shù)規(guī)劃或大規(guī)模非線性規(guī)劃求解時(shí)間可能很長。整數(shù)規(guī)劃利用intlinprog的options設(shè)置RelativeGapTolerance。默認(rèn)是1e-4你可以設(shè)為1e-3或5e-3來加速。設(shè)置MaxTime限制最長運(yùn)行時(shí)間。提供初始解對于intlinprog你可以通過x0參數(shù)提供一個可行的整數(shù)初始解這能顯著加快分支定界法的求解過程。簡化模型審視你的模型是否有一些不必要的變量或約束能否通過問題本身的特性進(jìn)行簡化線性化如果可能將非線性部分近似為線性用線性規(guī)劃求解會快得多。使用更高效的算法或工具對于特定類型的問題如二次規(guī)劃quadprog使用專用求解器比通用的fmincon更快。對于超大規(guī)模問題可能需要考慮商業(yè)求解器如Gurobi、CPLEX或者利用問題結(jié)構(gòu)設(shè)計(jì)分解算法。6. 從求解到應(yīng)用結(jié)果分析與模型檢驗(yàn)?zāi)玫角蠼馄鞯妮敵鰔和fval工作只完成了一半。一個負(fù)責(zé)任的建模者必須對結(jié)果進(jìn)行分析和檢驗(yàn)。敏感性分析對于線性規(guī)劃尤其重要優(yōu)化解在多大程度上依賴于模型參數(shù)如果某個資源約束右端項(xiàng)b增加一個單位最優(yōu)目標(biāo)值能改善多少這個“改善率”就是該資源的影子價(jià)格Shadow Price。在linprog中可以通過輸出參數(shù)lambda拉格朗日乘子的下界部分獲得。lambda.ineqlin對應(yīng)不等式約束A*x b的影子價(jià)格。影子價(jià)格高的資源是瓶頸資源增加其供給能帶來較大效益。解的解釋與呈現(xiàn)將數(shù)學(xué)解x翻譯回實(shí)際問題語言。比如x(1)3.5在實(shí)際中可能意味著3.5小時(shí)、3.5噸或者是需要四舍五入為4如果是整數(shù)規(guī)劃則不存在此問題。你需要根據(jù)問題的實(shí)際背景來解釋這個解。模型穩(wěn)健性檢驗(yàn)稍微改變一下模型參數(shù)比如目標(biāo)函數(shù)系數(shù)c或約束右端項(xiàng)b在±10%范圍內(nèi)波動重新求解觀察最優(yōu)解的變化是否劇烈。如果最優(yōu)解變化很大說明模型對參數(shù)很敏感你需要謹(jǐn)慎對待這些參數(shù)取值的準(zhǔn)確性或者在報(bào)告中說明這種敏感性。最后我想說的是MATLAB的優(yōu)化工具箱是一個強(qiáng)大的武器但武器本身不會思考。真正的核心能力在于你將一個模糊的實(shí)際問題清晰、準(zhǔn)確地抽象為一個標(biāo)準(zhǔn)規(guī)劃問題的數(shù)學(xué)模型的能力以及當(dāng)求解器“不聽話”時(shí)你能像偵探一樣根據(jù)exitflag、輸出信息和問題背景一步步排查、調(diào)試、修正模型和代碼的能力。這個過程充滿挑戰(zhàn)但每一次成功的求解都是對邏輯思維和工程實(shí)踐能力的一次扎實(shí)提升。多練、多思考、多踩坑你自然就能形成自己的“數(shù)感”和“碼感”在數(shù)模競賽或?qū)嶋H項(xiàng)目中更加游刃有余。