【新卒研修資料】数理最適化 / Mathematical Optimization

by BrainPad

Slide 1

Slide 1 text

数理最適化研修

Slide 2

Slide 2 text

2 ©BrainPad Inc. Strictly Confidential 本日のスケジュールタイムライン所要時間(分)* 内容 9:45-10:15 30 ①数理最適化イントロダクション 10:15-11:10 55 ②線形計画定式化 11:10-12:00 50 ③線形計画実装 12:00-13:00 60 昼休憩 13:00-14:00 60 ④整数計画定式化&実装 14:00-15:00 60 ⑤混合整数計画定式化&実装 15:00-16:40 100 ⑥実践演習 16:40:17:40 60 ⑦案件紹介 17:40-18:00 20 ⑧今後学ぶと良いこと *適宜休憩時間を設けます

Slide 3

Slide 3 text

3 ©BrainPad Inc. Strictly Confidential この講義では、数理最適化の基礎知識を習得していただくことが目標です。本研修のゴールビジネス職のゴール最適化でできることや最適化PJの進め方を知る  「最適化」、「最適解」というビジネス用語を使う時に目的関数、制約条件は何なのかまでを意識して”いいかんじ”の中身を言語化できるようにする  数理最適化でできることを理解し、クライアント/PJメンバーとのコミュニケーションを円滑に進められるようにする技術職のゴール簡単な最適化問題であれば自力で定式化・実装できる  最適化問題の定式化のコツと実装の流れを問題演習を通して理解する  現実問題に対しても定式化をして最適化問題に落とし込むことができるようになる

Slide 4

Slide 4 text

©BrainPad Inc. Strictly Confidential ①数理最適化イントロダクション

Slide 5

Slide 5 text

©BrainPad Inc. Strictly Confidential 5 1. 最適化問題・数理最適化とは 2. なぜ最適化を学ぶのか 3. ブレインパッドにおける最適化案件 4. 最適化案件の進め方数理最適化イントロダクション

Slide 6

Slide 6 text

6 ©BrainPad Inc. Strictly Confidential 最適化問題: 与えられた制約条件のもとで目的関数の値を最大(最小)にする解を求める問題数理最適化: 最適化問題を解くことを通じて現実における意思決定や問題解決を実現する技術最適化問題・数理最適化とは最適化問題1 最適化問題2 配送計画問題機械学習における”学習” 目的関数: • 配送トラックの最小化 • 移動距離最小化、. . . 制約条件: • 全ての店舗を訪問 • トラックの積載量の上限、. . . 目的関数: 損失関数の最小化制約条件: 𝜃に対する制約小売A 小売B 小売C 倉庫機械学習モデルパラメータ𝜃 入力データ 𝑥 正解データ 𝑦 予測データ ො 𝑦 学習データ(𝑥, 𝑦) 損失関数 (異なり具合を評価)

Slide 7

Slide 7 text

7 ©BrainPad Inc. Strictly Confidential 実は今日のチーム分けも最適化問題で行っています！今日の講義が終わる頃には下の式の意味がわかるようになっているはずです最適化問題・数理最適化とは集合受講生の集合グループの集合問題定数：受講生 i の最適化のスキル：受講生 i のPythonのスキル：全員のスキルの総和をグループ数で割った値

Slide 8

Slide 8 text

8 ©BrainPad Inc. Strictly Confidential 機械学習と数理最適化はどう違うのかと思われた方もいるかと思います機械学習は予測を行うのに対して、数理最適化は実際に何をするのかの意思決定を行うという違いがあります例として天気予報をもとに傘を持っていくか決める状況を考えます最適化問題・数理最適化とは機械学習（予測）数理最適化（意思決定） *余談ですが講師は折り畳み傘を常に持ち歩くが最適解だと考えています *前述の機械学習における”学習”のように様々な分野の基盤技術としての側面もあります • 気温 • 湿度 • 気象衛生 • etc. . . 降水確率90% 説明変数予測モデル目的変数雨が降りそうなので傘を持っていく機械学習においては説明変数を入力として予測モデルを用いて目的変数を予測します数理モデル(最適化問題) 求解目的関数 • 雨に濡れた時の”嫌さ”と傘を持っていく時の”煩わしさ”の総和を最小化制約条件 • 手荷物の多さ(両手が塞がっていないか？) • etc . . . 最適化問題を解いた結果傘を持っていくのが最適化である数理最適化では予測モデルの結果をもとに数理モデルを用いて意思決定しますこの意思決定を数理モデルに基づいて解決！

Slide 9

Slide 9 text

9 ©BrainPad Inc. Strictly Confidential 実案件においてもデータ分析を通してデータドリブンな意思決定を行いたいという場面はたくさんあります意思決定を行う状況において最適化の活躍の余地があります。以下に広告予算配分を例に説明しますなぜ最適化を学ぶのか予測(機械学習) 意思決定(数理最適化) 可視化まずは実績データの可視化を通して現状の確認を行います例) • 広告費用と利益の散布図 • 広告費用とCVの散布図 • 広告費用の時系列次にどのくらいの広告費用をかけるとどのくらいの利益やCVが得られるかの予測モデルを構築します最後に予測モデルをもとに予算上限を超えない下で利益やCVを最大化する最適化問題を解きます得られた最適な予算配分可視化・予測は「現状どうか？」や「今後どうなる？」には答えてくれるが「どうしたら良い？」には明示的に答えてくれない *そうは言ってもROIやデータが得られないなどの観点から最適化まで行わずに可視化や予測で十分という状況も往々にしてあります広告費用 CV数広告A 広告B 広告費用 CV数広告A 広告B 2025/04 広告費用(万円) 2025/05 2025/06

Slide 10

Slide 10 text

10 ©BrainPad Inc. Strictly Confidential 最適化案件は以下の一連の手続きからなります定式化がアルゴリズムが一発で決まることは稀でお客さんのフィードバックをもらいながら何回もサイクルを回します最適化案件の進め方現実問題最適化問題計算結果問題解決定式化アルゴリズム or ソルバーによる求解得られた解の可視化・確認最適化問題とアルゴリズムの再検討 https://zenn.dev/umepon/articles/50a74c75f32096875d94 を参考に作成

Slide 11

Slide 11 text

11 ©BrainPad Inc. Strictly Confidential 本日の研修では典型問題を通して定式化→ソルバーによる求解の流れを学んでいただき、最後に実践演習として一連のサイクルを手を動かしながら学んでもらいます最適化案件の進め方 https://zenn.dev/umepon/articles/50a74c75f32096875d94 を参考に作成実践演習線形計画定式化&実装整数計画定式化&実装混合整数計画定式化&実装

Slide 12

Slide 12 text

12 ©BrainPad Inc. Strictly Confidential 最適化案件の進め方 | 定式化現実問題最適化問題定式化アルゴリズム or ソルバーによる求解得られた解の可視化・確認最適化問題とアルゴリズムの再検討まず初めに現実問題を数式を使って最適化問題として記述する必要があります目的関数：最大化・最小化したい指標決定変数：最適化により決定したい値制約式：指標の最適化にあたっての制約生産スケジュールの作成を自動化したい。製品の数は〇〇で各商品に対して△△の条件を満たす必要があり…、××のあとは□時間インターバルを空ける必要があり… 定式化 ※定式化イメージ

Slide 13

Slide 13 text

13 ©BrainPad Inc. Strictly Confidential 最適化案件の進め方 | 最適化問題を解く現実問題現実問題を最適化問題として記述できれば、その問題をどう解くかに専念することができます最適化問題を解くアプローチとしては大きく以下の２通りに大別されます (1)汎用ソルバーによる解法幅広い問題に対して適用可能なアルゴリズムを予め実装したソフトウェアを利用する方法制約条件の変更が生じても比較的変更が容易 (2)フルスクラッチアルゴリズムによる解法自分で1から求解アルゴリズムを実装する方法問題固有の性質を使うことができるので求解性能は高い一方で保守性が低くなりやすい求解性能汎用性 (1)汎用解法 (2)フルスクラッチそれぞれのアプローチには性能・汎用性についてのトレードオフが存在する案件ではソルバーが使われることが多い最適化問題計算結果定式化アルゴリズム or ソルバーによる求解得られた解の可視化・確認最適化問題とアルゴリズムの再検討

Slide 14

Slide 14 text

14 ©BrainPad Inc. Strictly Confidential 最適化案件の進め方 | 解の確認・定式化の修正考えた定式化が一発でうまくいくことはまずあり得ません最適化問題の定式化やアルゴリズムを見直しながら改善を繰り返しますソルバーやアルゴリズムによって解を得た後は結果を可視化を行います生産計画における生産順番を表すガントチャート可視化の例工場の設備上、製品Aの次に製品B を生産することはできません定式化の修正 ⋮ 部分的な修正で良いこともあればモデルの大部分を変更しなければならない破滅的変更を伴う場合もあります現実問題最適化問題計算結果問題解決定式化アルゴリズム or ソルバーによる求解得られた解の可視化・確認最適化問題とアルゴリズムの再検討 2024/12/1 12:00 2024/12/4 12:00 2024/12/2 12:00 2024/12/3 12:00 休み製品A 製品B 製品E 製品C 製品D 休み

Slide 15

Slide 15 text

©BrainPad Inc. Strictly Confidential ②線形計画定式化

Slide 16

Slide 16 text

©BrainPad Inc. Strictly Confidential 16 1. 最適化問題の基礎知識 2. 線形計画問題とは 3. 例題 4. 類題演習 5. 質疑応答線形計画定式化

Slide 17

Slide 17 text

17 ©BrainPad Inc. Strictly Confidential ひとくちに最適化問題と言ってもさまざまな最適化問題が存在します決定変数、目的関数、制約条件の性質によって以下のように分類されることがあります最適化問題の基礎知識 | 代表的な最適化問題引用元：梅谷俊治「しっかり学ぶ数理最適化」連続最適化決定変数が実数値のような連続的な値を取る最適化問題離散最適化（組合せ最適化）決定変数が整数値や{0, 1}の2値のような離散的な値をとる最適化問題や、最適解を含む解の集合が順列やネットワークなど組合せ的な構造を持つ最適化問題非線形計画問題(NLP) 2次計画問題(QP) 線形計画問題(LP) 整数計画問題(IP) 混合整数計画問題(MIP) ネットワーク最適化問題

Slide 18

Slide 18 text

18 ©BrainPad Inc. Strictly Confidential 本講義では線形計画問題、整数計画問題、混合整数計画問題について解説します最適化問題の基礎知識 | 代表的な最適化問題引用元：梅谷俊治「しっかり学ぶ数理最適化」線形計画問題(LP) 整数計画問題(IP) 混合整数計画問題(MIP)

Slide 19

Slide 19 text

19 ©BrainPad Inc. Strictly Confidential 最適化問題は以上のような目的関数と制約条件から構成されます 𝑓(𝑥) 𝑥 最小化局所最適解最適解実行可能領域続いて基本的な最適化における用語について確認していきます最適化問題の基礎知識 | 用語最適値引用元：梅谷俊治「しっかり学ぶ数理最適化」 𝑆 用語意味決定変数最適化対象の変数 (上の問題の𝑥) 解決定変数に割り当てられた値実行可能解制約条件を満たす解実行可能領域実行可能解全体の集合 (大域)最適解目的関数の値が最大or最小となる解局所最適解近傍内の任意の解よりも目的関数値が大きいor小さい解最適値最適解における目的関数の値

Slide 20

Slide 20 text

20 ©BrainPad Inc. Strictly Confidential 最適化問題は必ず最適解が存在するとは限りませんそもそも実行可能解が存在しなかったり、いくらでも目的関数を改善できてしまう状況が存在します最適化問題の基礎知識 | 用語 𝑓(𝑥) 𝑥 最小化最適解実行可能領域 𝑓(𝑥) 𝑥 最大化実行可能領域 𝑥を大きくすることでいくらでも目的関数値 𝑓(𝑥) の値を大きくできてしまう最適解が存在実行不能非有界 𝑓(𝑥) 𝑥 最小化実行可能領域が空集合制約条件 𝑥 ≤ 1 かつ 3 ≤ 𝑥 引用元：梅谷俊治「しっかり学ぶ数理最適化」

Slide 21

Slide 21 text

21 ©BrainPad Inc. Strictly Confidential 線形計画問題とは「いくつかの一次不等式や等式を制約として、ある一次関数を最大化・最小化する問題」のことです例えば以下のような問題があります線形計画問題(LP: Linear Programming)とは例1 例2

Slide 22

Slide 22 text

22 ©BrainPad Inc. Strictly Confidential 何やら難しく聞こえるかもしれませんが、実は線形計画問題は高校数学でも登場しています線形計画問題(LP: Linear Programming)とは引用元：高校数学の美しい物語

Slide 23

Slide 23 text

23 ©BrainPad Inc. Strictly Confidential 線形計画問題の基本的な性質線形計画問題(LP: Linear Programming)とは ①最大化と最小化は等価 (−1) × (−1) = 1のようなイメージです最尤推定において”尤度の最大化”を行う代わりに”負の対数尤度の最小化”を行うのもこの性質に基づいています* もちろん実務においては目的関数によって最大化or最小化が自然という場合はあります等価等価 *単調増加関数の性質も使っています *慣習的に連続最適化においては最小化問題として扱うことが多いです

Slide 24

Slide 24 text

24 ©BrainPad Inc. Strictly Confidential 線形計画問題の基本的な性質線形計画問題(LP: Linear Programming)とは ②不等式制約と等式制約は等価変形可能等式は不等式に書き換えられ、不等式も等式に書き換えることができるので全体として不等式制約と等式制約は等価となります変形可能

Slide 25

Slide 25 text

25 ©BrainPad Inc. Strictly Confidential 簡易的なビジネス課題を最適化問題として定式化に取り組みます例題＆類題演習 | 最適化問題の定式化例題 & 類題演習現実問題最適化問題定式化

Slide 26

Slide 26 text

26 ©BrainPad Inc. Strictly Confidential 定式化を行う際は以下の順番を意識すると見通しが良くなります例題 | 最適化問題を定式化する流れ 1. 決定変数を決める 2. 目的関数を決定変数を用いて表現する 3. 制約条件を決定変数を用いて表現する目的関数：最大化・最小化したい指標決定変数：最適化により決定したい値制約式：指標の最適化にあたっての制約最終的に作りたい定式化

Slide 27

Slide 27 text

27 ©BrainPad Inc. Strictly Confidential 実際に簡易的なビジネス課題を最適化問題として定式化してみましょう例題 | 飲料の生産量最適化￥3,000 ￥2,000 あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり2,000円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーを何kgずつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫の通り問題文フルーツジュースフルーツティー材料在庫量リンゴ 2 5 500 オレンジ 4 2 300 材料商品 2kg 4kg 5kg 2kg

Slide 28

Slide 28 text

28 ©BrainPad Inc. Strictly Confidential まず決定変数を定義します例題 | 解説(1/3) 決定変数を決めよう問題文考え方 • 本問で決めたいことは「フルーツジュースとフルーツティーを何kgずつ作るか」 • 従ってひとまず決定変数は以下のようにおいてみる • フルーツジュースを作る量 𝑥1 kg • フルーツティーを作る量 𝑥2 kg ￥3,000 ￥2,000 あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり2,000円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーを何kgずつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫の通りフルーツジュースフルーツティー材料在庫量リンゴ 2 5 500 オレンジ 4 2 300 材料商品 2kg 4kg 5kg 2kg

Slide 29

Slide 29 text

29 ©BrainPad Inc. Strictly Confidential 次に目的関数を決定変数を用いて表します例題 | 解説(2/3) 目的関数を決めよう問題文考え方 • 本問では「売上を最大化したい」 • 売上はジュースの売上とティーの売上の合計 • 式で書くと、 • 売上 = ジュースの売上+ ティーの売上 = 3000 × 𝑥1 + 2000 × 𝑥2 • ジュースを作る量 𝑥1 kg • ティーを作る量 𝑥2 kg あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり2,000円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーを何kgずつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫の通りフルーツジュースフルーツティー材料在庫量リンゴ 2 5 500 オレンジ 4 2 300 材料商品

Slide 30

Slide 30 text

30 ©BrainPad Inc. Strictly Confidential フルーツジュースフルーツティー材料在庫量リンゴ 2 5 500 オレンジ 4 2 300 材料商品最後に制約条件を決定変数を用いて表します例題 | 解説(3/3) 制約条件を決めよう考え方 • 本問での制約条件は材料の在庫数 • 材料としてはリンゴ、オレンジの二つがあるのでそれぞれについて在庫量の制約を考えれば良い • リンゴの材料の制約は、ジュースで使うリンゴ＋ティーで使うリンゴ≦リンゴの在庫 ⇔ 1kgのジュースに使うリンゴの量 × ジュース生産量 + 1kgのティーに使うリンゴの量 × ティー生産量 ≦リンゴの在庫 • 式で書くと • 同様にオレンジの材料の制約は問題文あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり2,000円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーを何kgずつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫の通り

Slide 31

Slide 31 text

31 ©BrainPad Inc. Strictly Confidential 得られた決定変数、目的関数、制約式を一つにすると下記の通りとなります例題 | 解答例題の定式化

Slide 32

Slide 32 text

32 ©BrainPad Inc. Strictly Confidential 下記の問題を演習として定式化してみましょう類題演習 | 飲料の生産量最適化あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり2,000円、リンゴジュースは1kgあたり500円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーとリンゴジュースをいくつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫数の通り • リンゴジュースはフルーツジュース以上の量を作ることとする問題文 • この制約が必要になるのはどんなときでしょうか？ • ⇒ 例えば、リンゴジュースは、 • フルーツジュースとセット販売 • 単品で販売の2通りで売る時などフルーツジュースフルーツティーリンゴジュース材料在庫量リンゴ 2 5 1 500 オレンジ 4 2 0 300 材料商品

Slide 33

Slide 33 text

33 ©BrainPad Inc. Strictly Confidential 例題を参考にしながら、類題の定式化を手元のノートに書いてください。グループで相談してOKです。質問や疑問点があれば近くのTAまで聞いてください！類題演習 | 飲料の生産量最適化作業時間（10分）

Slide 34

Slide 34 text

34 ©BrainPad Inc. Strictly Confidential 類題演習 | 解説(1/3) 決定変数を決めよう考え方 • 基本的な考え方は先ほどと同様です • 本問で決めたいことは「フルーツジュース、フルーツティー、リンゴジュースをどのくらいの量作るか」 • 決定変数をとりあえず以下のようにおいてみる • フルーツジュースを作る量 𝑥1 kg • フルーツティーを作る量 𝑥2 kg • リンゴジュースを作る量 𝑥3 kg まず決定変数を決めます問題文あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり 2,000円、リンゴジュースは1kgあたり500円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーとリンゴジュースをいくつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫数の通り • リンゴジュースはフルーツジュース以上の量を作ることとするフルーツジュースフルーツティーリンゴジュース材料在庫量リンゴ 2 5 1 500 オレンジ 4 2 0 300 材料商品

Slide 35

Slide 35 text

35 ©BrainPad Inc. Strictly Confidential 類題演習 | 解説(2/3) 目的関数を決めよう問題文考え方 • 本問では「売上を最大化したい」 • 売上は、フルーツジュースとフルーツティーとリンゴジュースの売上の合計 • フルーツジュースの売上 = 1kgあたりのフルーツジュースの売値 × フルーツジュースの生産量 • フルーツティーの売上 = 1kgあたりのフルーツティーの売値 × フルーツティーの生産量 • リンゴジュースの売上 = 1kgあたりのリンゴジュースの売値 × リンゴジュースの生産量 • 式で書くと次に決定変数を用いて目的関数を表しますあなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり 2,000円、リンゴジュースは1kgあたり500円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーとリンゴジュースをいくつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫数の通り • リンゴジュースはフルーツジュース以上の量を作ることとするフルーツジュースフルーツティーリンゴジュース材料在庫量リンゴ 2 5 1 500 オレンジ 4 2 0 300 材料商品

Slide 36

Slide 36 text

36 ©BrainPad Inc. Strictly Confidential 類題演習 | 解説(3/3) 制約条件を決めよう考え方 • 本問での制約条件は2種類の制約を考える必要がある • 材料 (リンゴ、オレンジ)の在庫量制約 • リンゴジュースはフルーツジュースよりも多く作る • 在庫数制約は例題と同様に考えることができる • リンゴの在庫制約 • オレンジの在庫制約 • リンゴジュースはフルーツジュースよりも多く作る最後に決定変数を用いて制約条件を表します問題文あなたは飲料を製造・販売している業者です。フルーツジュースは1kgあたり3,000円、フルーツティーは1kgあたり 2,000円、リンゴジュースは1kgあたり500円で販売しています。売上を最大化するには、フルーツジュースとフルーツティーとリンゴジュースをいくつ作ればよいでしょうか？ただし、製造に関して、下記の条件がある • 各商品は、リンゴとオレンジを材料とする • 各商品を1kg作るのに必要な材料の量は、下表の通り • 使える材料の量は、下表の在庫数の通り • リンゴジュースはフルーツジュース以上の量を作ることとするフルーツジュースフルーツティーリンゴジュース材料在庫量リンゴ 2 5 1 500 オレンジ 4 2 0 300 材料商品

Slide 37

Slide 37 text

37 ©BrainPad Inc. Strictly Confidential 得られた決定変数、目的関数、制約式を一つにすると下記の通りとなります類題演習 | 解答類題の定式化

Slide 38

Slide 38 text

©BrainPad Inc. Strictly Confidential ③線形計画実装

Slide 39

Slide 39 text

©BrainPad Inc. Strictly Confidential 39 1. 最適化問題の解き方 2. 例題 3. 類題演習 4. 質疑応答線形計画実装

Slide 40

Slide 40 text

40 ©BrainPad Inc. Strictly Confidential 先程までで最適化問題に定式化する方法を学びました。これ以降はどうやって解くかについて説明します本講義では特に汎用ソルバーを用いて最適化問題を解く方法について学びます最適化問題の解き方モデラーソルバー再掲定式化した問題をソルバーが解釈可能な形にするための変換器線形計画問題ではLPファイルを通してソルバーに問題を渡します* MAXIMIZE 2*x_0 + 4*x_1 + 1*x_2 + 3*x_3 + 0 SUBJECT TO _C1: 3 x_0 + 5 x_1 + x_2 + 4 x_3 <= 8 VARIABLES 0 <= x_0 <= 1 Integer 0 <= x_1 <= 1 Integer 0 <= x_2 <= 1 Integer 0 <= x_3 <= 1 Integer model.py model.lp 予めアルゴリズムが実装されたパッケージ線形計画問題では • Gurobi (有償) • SCIP (無償) • Cbc (無償) などのソルバーが存在します SCIP Cbc (1)汎用ソルバーによる解法幅広い問題に対して適用可能なアルゴリズムを予め実装したソフトウェアを利用する方法制約条件の変更が生じても比較的変更が容易 (2)フルスクラッチアルゴリズムによる解法自分で1から求解アルゴリズムを実装する方法問題固有の性質を使うことができるので求解性能は高い一方で保守性が低くなりやすい求解性能汎用性 (1)汎用解法 (2)フルスクラッチそれぞれのアプローチには性能・汎用性についてのトレードオフが存在する案件ではソルバーが使われることが多いモデラー

Slide 41

Slide 41 text

41 ©BrainPad Inc. Strictly Confidential Pythonでよく使われるモデラーやソルバーとしては以下のようなものがありますモデル作成の流れはどのモデラーでも共通する部分が多いです最適化問題の解き方モデラーソルバー PuLP CP-SAT pywraplp SCIP CPLEX Cbc CP-SAT *使用可能なソルバー、モデラーの組合せが存在する 1. 問題の定義 model = Model(name=‘sample’) 2. 決定変数の実装 x[i] = IntVar(name=f‘x{i}’, LB=0, UB=10) 3. 目的関数の実装 model += 2*x[0]+3*x[1]+x[3] 4. 制約条件の実装 model += 3*x[0]+6*x[1]+2*x[2]<=21 5. 求解 model.optimize() 6. 求めた解の値の取得 x[i].value() モデル作成の流れ機械学習におけるmodel.fit() → model.predict()のように最適化モデラーにおいてもモデリングの流れが存在しますモデラーに応じて細かい書き方は異なりますが、使い方の流れの大枠は共通しています gurobipy 引用：https://www.coin-or.org/、 https://www.gurobi.com/、https://github.com/scipopt、 https://crescointl.com/product/third-party-products/ibm-products/ibm-ilog-cplex-optimization-studio/、https://developers.google.com/optimization/routing/routing_options?hl=ja PySCIPOpt

Slide 42

Slide 42 text

42 ©BrainPad Inc. Strictly Confidential 簡易的なビジネス課題を想定し、ソルバーによる求解に取り組みます例題＆類題演習 | ソルバーによる求解例題 & 類題演習現実問題最適化問題計算結果定式化（アルゴリズム or）ソルバーによる求解

Slide 43

Slide 43 text

43 ©BrainPad Inc. Strictly Confidential 先程までの例題を実際にソルバーに解かせてみましょう今回はPuLPというモデラーとCbcというソルバーを用いて最適化問題を解きます例題（実装）例題の定式化

Slide 44

Slide 44 text

44 ©BrainPad Inc. Strictly Confidential 最初に、どんな問題（最大化/最小化）を解きたいのかを定義します例題 | 解説（1/5）問題の定義をしよう問題文考え方 • LpProblemで問題を定義する。 • name：最適化問題の名前 • sense：目的関数を最大化 or 最小化するか import pulp model = pulp.LpProblem(name= ='例題', sense=pulp.LpMaximize)

Slide 45

Slide 45 text

45 ©BrainPad Inc. Strictly Confidential 次に決定変数を実装します例題 | 解説（2/5）決定変数を実装しよう問題文考え方 x1 = pulp.LpVariable(name='x1', lowBound=0, upBound=None, cat='Continuous') x2 = pulp.LpVariable(name='x2', lowBound=0, upBound=None, cat='Continuous') ● LpVariableで決定変数を定義します。 ○ name：変数の名前を表す ○ lowBound：変数の取る値の最小値 ○ upBound：変数の取る値の最大値 ○ cat：変数の型 ■ continuousは実数

Slide 46

Slide 46 text

46 ©BrainPad Inc. Strictly Confidential 次に目的関数を実装します例題 | 解説（3/5）目的関数を実装しよう問題文考え方 ● 目的関数を以下のように実装します ● 間違えている例 ● 多くのモデラーは複合代入演算子で目的関数、制約式の定義を行うので最初は違和感を感じるかもしれません model += 3000 * x1 + 2000 * x2 model = 3000 * x1 + 2000 * x2

Slide 47

Slide 47 text

47 ©BrainPad Inc. Strictly Confidential 次に制約条件を実装します例題 | 解説（4/5）制約条件を実装しよう問題文考え方 ● 制約条件を以下のように実装します model += 2 * x1 + 5 * x2 <= 500 model += 4 * x1 + 2 * x2 <= 300

Slide 48

Slide 48 text

48 ©BrainPad Inc. Strictly Confidential 最適化問題の実装が完了したので実際にソルバーに解かせてみます例題 | 解説（5/5）解を求めてみよう model.solve() print('Status:', pulp.LpStatus[model.status]) print('x1=', x1.value(), 'x2=', x2.value()) Solveメソッドを呼んで最適化問題を解きます Status: Optimal x1= 31.25 x2=87.5 print文の結果はとなって最適解を得られました変数.value() は忘れやすいので注意

Slide 49

Slide 49 text

49 ©BrainPad Inc. Strictly Confidential 同様の問題を10分で自分で実装してみましょう類題演習類題の定式化

Slide 50

Slide 50 text

50 ©BrainPad Inc. Strictly Confidential 類題演習作業時間（10分）分からない所などあったら、遠慮なく声をかけてください。 5分経過したらヒントを出します。

Slide 51

Slide 51 text

©BrainPad Inc. Strictly Confidential ④整数計画問題定式化&実装

Slide 52

Slide 52 text

©BrainPad Inc. Strictly Confidential 52 1. 整数計画問題とは 2. 例題 3. 類題演習 4. 質疑応答整数計画問題定式化&実装

Slide 53

Slide 53 text

53 ©BrainPad Inc. Strictly Confidential 決定変数が全て整数であるような最適化問題を整数計画問題と呼びます目的関数や決定変数の取り得る値によって問題が分類されます整数計画問題(IP: Integer Programming)とは  0-1整数計画問題 (0-1 integer programming) 全ての変数が0か1かという制約が課される問題です。ある選択肢を「選択する・選択しない」、あるルートを「通る・通らない」といったon/offで表せるものはこの問題になります。  線形整数計画問題決定変数が全て整数で、目的関数と制約条件が線形で表される問題です。  非線形整数計画問題目的関数もしくは制約条件に線形ではないものを含んだより広い問題です。目的関数をとして二乗誤差などを設定した場合はこの問題となります。この問題は本講義で解説しているPuLP+Cbcではモデリング・求解できないので他のモデラー・ソルバーを使う必要があります

Slide 54

Slide 54 text

54 ©BrainPad Inc. Strictly Confidential 下記のビジネス問題を定式化してみますポイントは①決定変数、②目的関数、③制約式の順に作ることです例題 | ナップサック問題問題文あなたは、商品をナップサック（カバン）に詰めて販売している行商です。ナップサックに、どの商品を詰めたら、総収入を最大化できるでしょうか？ただし、商品とナップサックには、次の条件があります。 • 商品は、４種類あり、それぞれ重さと価値が違う（下表参照）。 • ナップサックには、合計で８ｋｇまでしか詰められない。商品 1 2 3 4 価値（万円） 2 4 1 3 重さ（kg） 3 5 1 4 最もお金を稼ぐにはどう詰めるか？

Slide 55

Slide 55 text

55 ©BrainPad Inc. Strictly Confidential 最初に、最適化で決定する量となる決定変数を決めます例題 | 定式化解説（1/3）決定変数を決めようあなたは、商品をナップサック（カバン）に詰めて販売している行商です。ナップサックに、どの商品を詰めたら、総収入を最大化できるでしょうか？ただし、商品とナップサックには、次の条件がある。 • 商品は、４種類あり、それぞれ重さと価値が違う（下表）。 • ナップサックには、合計で８ｋｇまでしか詰められない。問題文考え方 • 本問で決めたいことは、どの商品をナップサックに詰めるか • 決定変数は商品 1 2 3 4 価値（万円） 2 4 1 3 重さ（kg） 3 5 1 4 商品iをナップサックに詰める商品iをナップサックに詰めない

Slide 56

Slide 56 text

56 ©BrainPad Inc. Strictly Confidential 次に、最大化（or 最小化）をしたい事象を目的関数として定義します例題 | 定式化解説（2/3）目的関数を決めようあなたは、商品をナップサック（カバン）に詰めて販売している行商です。ナップサックに、どの商品を詰めたら、総収入を最大化できるでしょうか？ただし、商品とナップサックには、次の条件がある。 • 商品は、４種類あり、それぞれ重さと価値が違う（下表）。 • ナップサックには、合計で８ｋｇまでしか詰められない。問題文考え方 • 本問で決めたいことは、どの商品を詰めると、総収入が最大化できるか • （再掲）決定変数は • 目的関数を決定変数を使って表現すると商品 1 2 3 4 価値（万円） 2 4 1 3 重さ（kg） 3 5 1 4 商品iをナップサックに詰める商品iをナップサックに詰めない

Slide 57

Slide 57 text

57 ©BrainPad Inc. Strictly Confidential 最後に、決定変数や目的関数が満たすべき制約を表した制約条件を決めます例題 | 定式化解説（3/3）制約条件を決めようあなたは、商品をナップサック（カバン）に詰めて販売している行商です。ナップサックに、どの商品を詰めたら、総収入を最大化できるでしょうか？ただし、商品とナップサックには、次の条件がある。 • 商品は、４種類あり、それぞれ重さと価値が違う（下表）。 • ナップサックには、合計で８ｋｇまでしか詰められない。問題文考え方 • ナップサックに8kgまでしか詰められないが制約条件です。 • （再掲）決定変数は • 決定変数を用いると、制約条件は商品 1 2 3 4 価値（万円） 2 4 1 3 重さ（kg） 3 5 1 4 商品iをナップサックに詰める商品iをナップサックに詰めない

Slide 58

Slide 58 text

58 ©BrainPad Inc. Strictly Confidential 得られた決定変数、目的関数、制約式を1つにすると下記の通りになります例題 | ナップサック問題解説例題の定式化

Slide 59

Slide 59 text

59 ©BrainPad Inc. Strictly Confidential 今までは価値や重さ、容量の値をベタ打ちしていましたが、実際は色々な入力データが入ってきますそのような状況に対応できるように集合と総和記号を使って定式化を一般化しましょう例題 | ナップサック問題解説例題の定式化例題の定式化 (一般化) :商品の集合 :商品𝑖の価値 :商品𝑖の重さ :ナップサックの容量集合と定数の定義集合と総和記号(∑)の使い方については次ページ参照

Slide 60

Slide 60 text

60 ©BrainPad Inc. Strictly Confidential さらに複数の制約式をまとめる制約式に現れる足し算をまとめる最適化問題の定式化では集合や総和記号( ∑ )がよく出てきます補足集合と総和記号 ※行列を使って制約をまとめることもあります (例: 𝐴𝑥 ≤ 𝑏, 𝐴 ∈ ℝ𝑚×𝑛, 𝑥 ∈ ℝ𝑛, 𝑏 ∈ ℝ𝑚)

Slide 61

Slide 61 text

61 ©BrainPad Inc. Strictly Confidential ここまでの話を踏まえて定式化を行う時は以下の順番を意識すると見通しが良くなります基本的な流れは先程と同じですが、まず集合を定義してから決定変数を決めるようにしましょう最適化問題を定式化する流れ 1. 集合,定数を定義する 2. 決定変数を決める 3. 目的関数を決定変数を用いて表現する 4. 制約条件を決定変数を用いて表現する :商品の集合 :商品𝑖の価値 :商品𝑖の重さ :ナップサックの容量 :商品𝑖を選択するか否か

Slide 62

Slide 62 text

62 ©BrainPad Inc. Strictly Confidential 先ほどのナップサック問題をPythonとpulpを用いて、実装しましょうコードを書く際も汎用的なコーディングを行うために一般化した定式化を意識しましょう例題 | 実装ナップサック問題例題の定式化例題の定式化 (一般化) :商品の集合 :商品𝑖の価値 :商品𝑖の重さ :ナップサックの容量集合と定数の定義入力の受け取り問題のインスタンス化

Slide 63

Slide 63 text

63 ©BrainPad Inc. Strictly Confidential 最初に、入力データの集合、定数を定義します実際の案件ではエクセルファイル経由でdataclassを用いたりして入力データを受け取ります例題 | 実装解説（1/6）集合、定数の定義をしよう問題文考え方 • 今回は入力データの定義を行う • 商品の集合をリスト、価値を辞書で定義する I = [1,2,3,4] v = {1:2, 2:4, 3:1, 4:3} w = {1:3, 2:5, 3:1, 4:4} C = 8

Slide 64

Slide 64 text

64 ©BrainPad Inc. Strictly Confidential 次に、どんな問題（最大化/最小化）を解きたいのかを定義します例題 | 実装解説（2/6）問題の定義をしよう問題文考え方 • LpProblemで問題を定義する。 • name：最適化問題の名前 • sense：目的関数を最大化 or 最小化するか import pulp model = pulp.LpProblem(name='knapsack', sense=pulp.LpMaximize)

Slide 65

Slide 65 text

65 ©BrainPad Inc. Strictly Confidential 次に、決定変数を実装します例題 | 実装解説（3/6）決定変数を実装しよう問題文考え方 • LpVariableで決定変数を定義する。 • name：変数の名前 • cat：変数の型 • {0, 1}変数は’Binary’と指定する • 商品集合Iを用いた内包表記で1行で定義を行う • pulpでは以下のように書いても同様の辞書に格納された決定変数を定義することができる x = {i: pulp.Lpvariable(name=f’x{i}’, cat=‘Binary’) for i in I} x = pulp.Lpvariable.dicts(’x’, I, cat=’Binary’)

Slide 66

Slide 66 text

66 ©BrainPad Inc. Strictly Confidential 次に、目的関数を実装します例題 | 実装解説（4/6）目的関数を実装しよう問題文考え方 • 目的変数を下記のように実装する • pulp.lpsum(リスト)でリストに格納された変数の総和を表現できる model += pulp.lpSum(v[i]*x[i] for i in I)

Slide 67

Slide 67 text

67 ©BrainPad Inc. Strictly Confidential 次に、制約条件を実装します例題 | 実装解説（5/6）制約条件を実装しよう問題文考え方 • 制約条件を下記の用に実装する。 model += pulp.lpSum(w[i]*x[i] for i in I) <= C

Slide 68

Slide 68 text

68 ©BrainPad Inc. Strictly Confidential 最後に、モデルを実行することで解が求まります例題 | 実装解説（6/6）解を求めてみようモデル実行の実装 model.solve() print('Status:', pulp.LpStatus[model.status]) print('x1=', x1.value(), 'x2=', x2.value(),'x3=', x3.value(), 'x4=', x4.value()) Print文の結果 Status: Optimal x1= 1.0 x2= 0.0 x3= 1.0 x4= 1.0

Slide 69

Slide 69 text

69 ©BrainPad Inc. Strictly Confidential 下記の問題を制限時間10分以内に定式化してください（ご自身のノートに書いてください）類題演習 | ビンパッキング問題問題文あなたは、１０個の商品をいくつかの箱に詰めて、発送しようとしています。どのように商品を組み合わせて箱詰めしたら、箱の使用数を最小化できるでしょうか？ただし、商品と箱には、次の条件がある。 • 商品は、１０種類（１種類につき１個）あり、それぞれ重さが違う（下表参照）。 • 箱１つには、合計１０ｋｇの商品しか詰められない。 • 箱は、合計１０個ある。 • 商品は、余すことなく箱詰めする。商品 1 2 3 4 5 6 7 8 9 10 重さ（kg） 3 5 7 4 8 8 6 9 6 5 ちなみに”ビン” はbottleでなく大箱のこと。

Slide 70

Slide 70 text

70 ©BrainPad Inc. Strictly Confidential ビンパッキング問題を制限時間10分以内に定式化してください（ご自身で手元のノートに書いて下さい）類題演習 | ビンパッキング問題作業時間（10分）分からない所などあったら、遠慮なく声をかけてください。 5分したらヒントを出します。

Slide 71

Slide 71 text

71 ©BrainPad Inc. Strictly Confidential 定式化を行う時は以下の順番を意識すると見通しが良くなります最適化問題を定式化する流れ 1. 集合,定数を定義する 2. 決定変数を決める 3. 目的関数を決定変数を用いて表現する 4. 制約条件を決定変数を用いて表現する :商品の集合 :商品𝑖の価値 :商品𝑖の重さ :ナップサックの容量 :商品𝑖を選択するか否か再掲

Slide 72

Slide 72 text

72 ©BrainPad Inc. Strictly Confidential 最初に、最適化で決定する量となる決定変数を決めます類題演習 | 定式化解説（1/3）決定変数を決めようあなたは、１０個の商品をいくつかの箱に詰めて、発送しようとしています。どのように商品を組み合わせ箱詰めしたら、箱の使用数を最小化できるでしょうか？ただし、商品と箱には、次の条件がある。 • 商品は、１０種類（１種類につき１個）あり、それぞれ重さが違う（下表参照）。 • 箱１つには、合計１０ｋｇの商品しか詰められない。 • 箱は、合計１０個ある。 • 商品は、余すことなく箱詰めする。問題文考え方 • 本問で決めたいことは、箱にどの商品を詰めるか • 決定変数は、 • 箱iを使うかどうか • 箱iに商品jを詰めるかどうか商品 1 2 3 4 5 6 7 8 9 10 重さ (KG) 3 5 7 4 8 8 6 9 6 5 箱iに商品を詰める場合上記以外の場合箱iに商品jを詰める場合上記以外の場合

Slide 73

Slide 73 text

73 ©BrainPad Inc. Strictly Confidential 次に、最大化（or 最小化）をしたい事象を目的関数として定義します類題演習 | 定式化解説（2/3）目的関数を決めようあなたは、１０個の商品をいくつかの箱に詰めて、発送しようとしています。どのように商品を組み合わせ箱詰めしたら、箱の使用数を最小化できるでしょうか？ただし、商品と箱には、次の条件がある。 • 商品は、１０種類（１種類につき１個）あり、それぞれ重さが違う（下表参照）。 • 箱１つには、合計１０ｋｇの商品しか詰められない。 • 箱は、合計１０個ある。 • 商品は、余すことなく箱詰めする。問題文考え方 • 目的関数は、使う箱の数の最小化 • （再掲）決定変数は、 • 箱iを使うかどうか • そのため、目的関数は、商品 1 2 3 4 5 6 7 8 9 10 重さ (KG) 3 5 7 4 8 8 6 9 6 5 箱iに商品を詰める場合上記以外の場合

Slide 74

Slide 74 text

74 ©BrainPad Inc. Strictly Confidential 最後に、決定変数や目的関数が満たすべき制約を表した制約条件を決めます類題演習 | 定式化解説（3/3）制約条件を決めようあなたは、１０個の商品をいくつかの箱に詰めて、発送しようとしています。どのように商品を組み合わせ箱詰めしたら、箱の使用数を最小化できるでしょうか？ただし、商品と箱には、次の条件がある。 • 商品は、１０種類（１種類につき１個）あり、それぞれ重さが違う（下表参照）。 • 箱１つには、合計１０ｋｇの商品しか詰められない。 • 箱は、合計１０個ある。 • 商品は、余すことなく箱詰めする。問題文考え方 • 満たすべき制約は2つある。 • 1つの箱には10Kgしか詰められない i番目の箱に注目すると • 全ての商品を箱詰めする必要がある j番目の商品に着目すると商品 1 2 3 4 5 6 7 8 9 10 重さ (KG) 3 5 7 4 8 8 6 9 6 5 箱iに入れる商品の重さの合計箱iに入る重さ y_i = 1→10 y_i = 0 →0

Slide 75

Slide 75 text

Slide 76

Slide 76 text

Slide 77

Slide 77 text

77 ©BrainPad Inc. Strictly Confidential ビンパッキング問題を実際に実装してみましょう作業時間（10分）分からない所などあったら、遠慮なく声をかけてください。 5分経過したらヒントを出します。類題演習 | 実装

Slide 78

Slide 78 text

Slide 79

Slide 79 text

Slide 80

Slide 80 text

80 ©BrainPad Inc. Strictly Confidential 実数をとる決定変数と整数をとる決定変数が両方含まれるような問題のことを混合整数計画問題(MIP)と呼びます混合整数計画問題も目的関数によってさらに分類できます混合整数計画問題(MIP: Mixed Integer Programming)とは  混合整数線形計画問題問題 (MILP: Mixed Integer Linear Programming) 目的関数が線形な場合の混合整数計画問題です古くから研究が行われており、さまざまな汎用ソルバーが開発されています (Cbcソルバーもそのひとつ)  混合整数非線形計画問題問題 (MINLP: Mixed Integer NonLinear Programming) 目的関数が非線形な場合の混合整数計画問題で、かなり解きづらい問題です一方で現実問題を愚直に定式化するとしばしばMINLPが現れますその場合は最適解を出すことを諦めてヒューリスティックアルゴリズムで解くか、頑張って混合線形計画に定式化して MILPソルバーで解くなどの選択肢があります

Slide 81

Slide 81 text

81 ©BrainPad Inc. Strictly Confidential 例題 | 施設配置問題問題文あなたはメーカーの物流担当者で、倉庫から製品を店舗に配送しようとしています契約可能な倉庫は5つあり、契約すると倉庫費用がかかりますどの倉庫を契約し、どれだけ配送すると、倉庫費用と輸送費用を最小化できるでしょうか？ただし、倉庫と店舗に次の条件がある 1. 倉庫から配送する製品は倉庫容量以下となる 2. 店舗の需要を満たすように配送する 3. 配送費用は、倉庫と店舗間で製品を1単位運ぶのに発生するコストとする倉庫契約費用: 倉庫容量: 店舗需要

Slide 82

Slide 82 text

82 ©BrainPad Inc. Strictly Confidential 例題 | 施設配置問題データ - 店舗倉庫 S1 S2 S3 S4 F1 10 6 7 4 F2 5 10 3 5 F3 4 2 8 7 F4 5 7 5 5 F5 2 1 1 3 倉庫倉庫費用倉庫容量 F1 1000 500 F2 1000 500 F3 1000 400 F4 1000 500 F5 1000 200 店舗年間需要量 S1 200 S2 100 S3 200 S4 100 店舗別需要量倉庫別の費用と容量倉庫店舗間の輸送費用 : 倉庫𝑖から店舗𝑗への輸送コスト

Slide 83

Slide 83 text

83 ©BrainPad Inc. Strictly Confidential 最初に、最適化で決定する量となる決定変数を決めます例題 | 定式化解説（1/3）決定変数を決めようあなたはメーカーの物流担当者で、倉庫から製品を店舗に配送しようとしています。契約可能な倉庫は5つあり、契約すると倉庫費用がかかります。どの倉庫を契約し、どれだけ配送すると、倉庫費用と輸送費用を最小化できるでしょうか？ただし、倉庫と店舗に次の条件がある。 1. 倉庫から配送する製品は倉庫容量以下となる。 2. 店舗の需要を満たすように配送する。 3. 配送費用は工場と店舗間で製品を1単位運ぶと発生する。問題文考え方 • 本問で決めたいことは以下の２つ • どの倉庫と契約するか • どの倉庫からどの店舗へ配送するか • 決定変数は各倉庫𝑖に対して倉庫iを契約する倉庫iを契約しない倉庫iから店舗jに運ぶ製品の量を

Slide 84

Slide 84 text

84 ©BrainPad Inc. Strictly Confidential 次に、最大化（or 最小化）をしたい事象を目的関数として定義します例題 | 定式化解説（2/3）目的関数を決めようあなたはメーカーの物流担当者で、倉庫から製品を店舗に配送しようとしています。契約可能な倉庫は5つあり、契約すると倉庫費用がかかります。どの倉庫を契約し、どれだけ配送すると、倉庫費用と輸送費用を最小化できるでしょうか？ただし、倉庫と店舗に次の条件がある。 1. 倉庫から配送する製品は倉庫容量以下となる。 2. 店舗の需要を満たすように配送する。 3. 配送費用は工場と店舗間で製品を1単位運ぶと発生する。問題文考え方 • 今回の問題は、倉庫費用と輸送費用の合計の最小化 • 倉庫費用は、 • 配送費用は、 :倉庫の集合 :倉庫iの倉庫費用 :倉庫iから店舗jの輸送費用

Slide 85

Slide 85 text

85 ©BrainPad Inc. Strictly Confidential 最後に、決定変数や目的関数が満たすべき制約を表した制約条件を決めます例題 | 定式化解説（3/3）制約条件を決めようあなたはメーカーの物流担当者で、倉庫から製品を店舗に配送しようとしています。契約可能な倉庫は5つあり、契約すると倉庫費用がかかります。どの倉庫を契約し、どれだけ配送すると、倉庫費用と輸送費用を最小化できるでしょうか？ただし、倉庫と店舗に次の条件がある。 1. 倉庫から配送する製品は倉庫容量以下となる。 2. 店舗の需要を満たすように配送する。 3. 配送費用は工場と店舗間で製品を1単位運ぶと発生する。問題文考え方 • 今回の問題で満たすべき制約は2つある。 • 倉庫から配送する製品は倉庫容量以下 • 店舗の需要を満たすように配送する 𝑀𝑖 は倉庫iの容量。倉庫iを使わない場合は、右辺が0となり倉庫iから運ぶ量の合計が0となる

Slide 86

Slide 86 text

Slide 87

Slide 87 text

Slide 88

Slide 88 text

Slide 89

Slide 89 text

Slide 90

Slide 90 text

90 ©BrainPad Inc. Strictly Confidential 例題 | 最適化研修のチーム分け問題文あなたは最適化研修の講師を担当することになりました各受講生33名には予め、最適化の経験(4段階)とPythonの経験(5段階)を答えてもらっています研修のグループワークを行う上で、8グループにバランスよく受講生を割り当てるにはどうすれば良いでしょうか？与えられているデータ id name 最適化スキル Pythonスキル 1 A 3 2 2 B 2 2 … 受講生の集合グループの集合グループごとにスキルポイントを均等に分けられた時の理想的な値

Slide 91

Slide 91 text

91 ©BrainPad Inc. Strictly Confidential まず最低限の割り当てを行うことを考える、とりあえず目的関数はなしで制約条件を考えます例題 | 定式化解説(1/3) あなたは最適化研修の講師を担当することになりました各受講生33名には予め、最適化の経験(4段階)とPythonの経験(5 段階)を答えてもらっています研修のグループワークを行う上で、8グループにバランスよく受講生を割り当てるにはどうすれば良いでしょうか？問題文考え方 • 今回は決定変数を人iをグループgに割り当てるか否かとする • 今回39人を9グループに割り当てるので、各グループは4人か 5人となる • さらに各人に対していずれかのグループに割り当てられる制約 id name 最適化スキル Pythonスキル 1 A 3 2 2 B 2 2 受講生の集合グループの集合

Slide 92

Slide 92 text

92 ©BrainPad Inc. Strictly Confidential 次にバランスの良い割り当てになるためにはどうすれば良いかを考えるます例題 | 定式化解説(2/3) 考え方 • 今回は各グループの経験値の総和が均等に割り振った時の値(AvgPoints)に近くなるようにします • ↑だと実行不能になる可能性があるため条件を緩めます AvgPoints あなたは最適化研修の講師を担当することになりました各受講生33名には予め、最適化の経験(4段階)とPythonの経験(5 段階)を答えてもらっています研修のグループワークを行う上で、8グループにバランスよく受講生を割り当てるにはどうすれば良いでしょうか？問題文 id name 最適化スキル Pythonスキル 1 A 3 2 2 B 2 2 受講生の集合グループの集合

Slide 93

Slide 93 text

93 ©BrainPad Inc. Strictly Confidential 次にバランスの良い割り当てになるためにはどうすれば良いかを考えます例題 | 定式化解説(3/3) 考え方 • この条件のもとで超過分+不足分の最小化を行いますハードな制約条件をソフトにするあなたは最適化研修の講師を担当することになりました各受講生33名には予め、最適化の経験(4段階)とPythonの経験(5 段階)を答えてもらっています研修のグループワークを行う上で、8グループにバランスよく受講生を割り当てるにはどうすれば良いでしょうか？問題文 id name 最適化スキル Pythonスキル 1 A 3 2 2 B 2 2 受講生の集合グループの集合

Slide 94

Slide 94 text

94 ©BrainPad Inc. Strictly Confidential 例題 | 解答ここまでの話を踏まえると以下のような定式化になります研修を通して理解できるようになったでしょうか？集合受講生の集合グループの集合問題定数：受講生 i の最適化のスキル：受講生 i のPythonのスキル：全員のスキルの総和をグループ数で割った値

Slide 95

Slide 95 text

Slide 96

Slide 96 text

96 ©BrainPad Inc. Strictly Confidential これまでで学んだ知識を総動員して実践演習を行います最適化案件で必要なすべての部分を考える必要があり、かなり難易度が高いですが頑張りましょう【実践演習】現実問題の定式化～問題解決まで実践演習の範囲現実問題最適化問題計算結果問題解決定式化アルゴリズム or ソルバーによる求解得られた解の可視化・確認最適化問題とアルゴリズムの再検討

Slide 97

Slide 97 text

97 ©BrainPad Inc. Strictly Confidential コンビニ食の献立最適化をテーマに実践演習を行います実践演習課題かつて関口さんの一週間のコンビニ食の献立を最適化するという社内PJが存在しました本日はこの社内PJの模擬体験を90分で行ってもらいます DOORS記事関口さんnote 演習の進め方 1. 定式化の記述実際のプロジェクトではヒアリングが発生しますが、今回は自作自演でチームでお題を決めて定式化を行ってくださいホワイトボードや手元の紙に定式化を書いてチーム外の人にも定式化を説明できるようにしてください 2. 最適化モデルの実装演習で行ったようにPuLPを用いて実際にモデルの作成を行ってください実装担当の人を決めても、全員で実装しても構いません 3. 解の可視化サンプルコードを参考に必要があれば改変して得られた解の可視化を行ってくださいチームで考えた目的関数や制約条件がしっかり満たされているかの確認などを行ってください 4. モデルの修正・改良得られた解を可視化して解に違和感がある場合はモデルの改善を行いましょう納得のいく解が出た場合は1週間分の献立に拡張やまだ使ってない情報を使うように改良などを考えてみてください

Slide 98

Slide 98 text

Slide 99

Slide 99 text

99 ©BrainPad Inc. Strictly Confidential ここまでお疲れ様でした！一日中講義、演習を行なってものすごく疲れたのではないかと思います。今日は数理最適化入門として色々な話題について勉強してきましたが、”最適化を使う”ハードルは意外と低かったのではないでしょうか？汎用ソルバーを使うことで「裏側のどうやって解くか」を気にせずにすぐにビジネス課題をどう最適化問題に落とし込むかに集中することができます。もちろん解けなかった場合の定式化のデバッグやアルゴリズムの設計が行いやすくなるので後々理論も勉強する必要がありますが徐々に勉強していけば良いです。最適化に対する敷居を高く思わずに、今日の実践演習のように身の回りの問題が最適化問題として定式化できないか？という視点を持っていただけるようになれば幸いです。はじめに

Slide 100

Slide 100 text

100 ©BrainPad Inc. Strictly Confidential 典型問題は「青チャート」みたいなもので、典型問題を組み合わせることで多くの問題は定式化できます以下に参考書籍をいくつか紹介します今後勉強すると良いと思うこと : MIPで典型問題を学ぶ DS新卒勉強会の指定書籍。実務において役立つモデリングや考え方が多く載っている典型問題の実装方法と定式化に関するさまざまなコツを紹介してくれる本混合整数計画問題へのモデリングがかなり詳しく載っている Pythonではじめる数理最適化あたらしい数理最適化 Python言語とGuroiで解く Model Building in Mathematical Programming

Slide 101

Slide 101 text

101 ©BrainPad Inc. Strictly Confidential  MIPでの定式化に慣れてくると、なんでも定式化・最適化できるのでは？という気持ちになってきますただそれは気のせいです. . .  実務で最適化を行う場合には「定式化はできるのだが大規模すぎて汎用ソルバーで問題が全然解けない」といった状況はよく起こります  その場合に「汎用ソルバーに頼らずに自分でアルゴリズムを書く」という選択肢を最終手段として持っておけると非常に強いです今後勉強すると良いと思うこと : ヒューリスティクスアルゴリズムを学ぶ Atcoder Heuristic Contest https://www.terry-u16.net/entry/ahc001 より引用広告配置を題材とした長方形詰め込み問題など、実務に近い問題が出題されています BrainPadプログラミングコンテスト2025（AtCoder Heuristic Contest 046）より引用

Slide 102

Slide 102 text

102 ©BrainPad Inc. Strictly Confidential ここまでPuLP+Cbcを使ってMILPを解く話をしてきましたが、MILPソルバー以外にも最適化ソルバーは存在します例としてCbcなどをはじめとするMILPソルバーとは異なるアルゴリズムで実装されたOR-ToolsのCP-SAT Solverなどがありますこちらのソルバーは非線形な目的関数も扱える代わりに決定変数は整数変数のみとなっていますが、論理式に関して柔軟なモデリングが可能です今後勉強すると良いと思うこと : 他のモデラー・ソルバーの使い方を学ぶ https://developers.google.com/optimization/cp/cp_solver のサンプルコード →not equalのような論理制約をモデリングすることができます

Slide 103

Slide 103 text

103 ©BrainPad Inc. Strictly Confidential 今日はnotebook上で最適化モデルを作って計算を行いましたが、実際の案件で最適化モデルを使うのはプログラミングを知らないお客さんであることが多いですさらに最適化条件を色々変えて最適化問題を解くことを考えると、最適化のモデルをWebアプリケーションにのせて提供することが望ましい場合があります今後勉強すると良いと思うこと : Webアプリケーションを学ぶ BrainPad FASTというソリューションを展開しています https://www.brainpad.co.jp/services/professionals/bp-fast.html

Slide 104

Slide 104 text

株式会社ブレインパッド 106-0032 東京都港区六本木三丁目1番1号六本木ティーキューブ TEL：03-6721-7002 FAX：03-6721-7010 www.brainpad.co.jp [email protected] 本資料は、未刊行文書として日本及び各国の著作権法に基づき保護されております。本資料には、株式会社ブレインパッド所有の特定情報が含まれており、これら情報に基づく本資料の内容は、貴社以外の第三者に開示されること、また、本資料を評価する以外の目的で、その一部または全文を複製、使用、公開することは、禁止されています。また、株式会社ブレインパッドによる書面での許可なく、それら情報の一部または全文を使用または公開することは、いかなる場合も禁じられております。 ©BrainPad Inc.

Slide 105

Slide 105 text

105 ©BrainPad Inc. Strictly Confidential [外部公開用_202508_数理最適化研修のコード（例題・類題のみ） - Google ドライブ] https://drive.google.com/drive/u/1/folders/1fQJ_3irJA8eMTEUWAls-mj1ApX_Bntmh Appendix.