まだ PuLPのdocumentに載っていない(執筆時)関数ですが sequentialSolveという多段階最適化を行う関数が実装されています import pulp from pulp import LpVariable # 問題の宣言 problem = pulp.LpProblem() # 変数の宣言(連続変数, 上限2) x = LpVariable('x', upBound=2, cat='Continuous') y = LpVariable('y', upBound=2, cat='Continuous') z = LpVariable('z', upBound=2, cat=…