1day internship Fixstars Amplify 20220808/組み合わせ最適化講習会

Fixstars Corporation www.fixstars.com
Copyright © Fixstars Group Copyright © Fixstars Group
【1dayインターン】
組み合わせ最適化講習会

View Slide

Copyright © Fixstars Group
2
本日のAgenda
第一部 (13:00～13:50): 会社紹介
• 会社概要
• 組合せ最適化問題及び事例のご紹介
• Fixstars Amplify のご紹介
第二部 (14:00～15:30): ハンズオンワークショップ
• 数の分割問題
• 搬送経路最適化 (7/25)、シフト最適化 (8/8)、生産計画最適化 (9/9)
• 発展的課題
第三部 (15:30～16:00): 今後について
• ぜひ情報発信を！
• 研究での活用事例
• Amplifyの業務
• 採用情報

View Slide

Q&Aについて
• 質問がある場合はZoomのチャット機能で「全員宛」で発信してくだ
さい
• その場で回答できるものは回答します
• または後ほど時間の許す限りお答えします
3

View Slide

【第一部】
会社紹介

View Slide

5
フィックスターズの概要
会社名株式会社フィックスターズ
本社所在地
東京都港区芝浦3-1-1 msb Tamachi
田町ステーションタワーN 28階
設立 2002年8月
上場区分東証プライム（証券コード：3687）
代表取締役社長三木聡
資本金 5億5,446万円（2021年9月現在）
社員数（連結） 258名（2021年9月現在）
主なお客様
キオクシア株式会社
株式会社ネクスティエレクトロニクス
株式会社日立製作所
キヤノン株式会社
グループ体制
株式会社フィックスターズ
株式会社Fixstars AutonomousTechnologies
株式会社ネクスティエレクトロニクスとのJV
自動運転向けソフトウェア開発
当社完全子会社
米国での営業及び開発
概要
株式会社Sider
株式会社Smart Opinion
開発支援SaaS「Sider」を運営
当社連結子会社
乳がんAI画像診断支援事業を運営
Fixstars Solutions, Inc.
オスカーテクノロジー株式会社
当社連結子会社
ソフトウェア自動並列化サービスを提供
株式会社Fixstars Amplify
量子コンピューティングのクラウド事業を運営
2021/10/1 設立

View Slide

6
ソフトウェア高速化サービス (概要)
• お客様にソースコードをご提供いただき、
• 最適化やアルゴリズムの改良を施して高速化してお返しします
当社お客様
オリジナルソースコードのご提供
高速化したソースコード
コンサルティング高速化サポート
性能評価
ボトルネックの特定
アルゴリズムの改良・開発
ハードウェアへの最適化
レポート作成
レポートやコードへのQ&A
実製品への組込み支援

View Slide

フィックスターズの強み
フィックスターズは、コンピュータの性能を最大限に引き出し大量データの高速処理を実現する、
高速化のエキスパート集団です。
7
低レイヤ
ソフトウェア技術
アルゴリズム
実装力
各産業・研究
分野の知見

View Slide

フィックスターズの量子技術への取り組み
8
2017年
日本で初めて
D-Wave Systems
社と提携
2018年
NEDOのプロジェクト
に採択
「イジングマシン共通
ソフトウェア基盤の研
究開発」
2019年
SIPの研究開発に参画
「光・量子を活用した
Society 5.0実現化技術：
光電子情報処理」
2021年
•2月：量子アニーリングクラウドサービス
「Fixstars Amplify」提供開始
•10月：子会社Fixstars Amplifyを設立
•11月：Q-STAR 量子技術による新産業創出協議
会に特別会員として加入
次世代技術を先取りし
今ある課題の解決を目指す

View Slide

Fixstars Amplifyの量子技術への取り組み
9
一般社団法人量子技術による新産業創出協議会
東洋経済主催シンポジウム
「ビジネスを劇的に変える量子コンピューティングの可
能性」にCEO平岡が登壇しました
モデレーター
九州大学藤澤克樹教授
パネリスト
ロート製薬株式会社
代表取締役会長山田邦雄氏
株式会社東芝
代表執行役社長CEO 島田太郎氏
株式会社Fixstars Amplify
代表取締役社長平岡卓爾
出展：一般社団法人量子技術による新産業創出協議会 https://qstar.jp/
東洋経済新報社主催「量子コンピューティングの実用化に関する」パネルディ
スカッションの基調講演では、実用可能な段階に入った量子コンピューティン
グの現状について、専門家の立場から意見を述べました。
パネルディスカッションでは、九州大学マス・フォア・インダストリ研究所の
藤澤克樹教授、東芝株式会社代表執行役社長 CEO で、量子技術による新産業
創出協議会（Q-STAR）実行委員長を務める島田太郎氏、ロート製薬代表取締
役会長の山田邦雄氏も参加され、専門家や経営者のお立場から、量子コン
ピューティングの現在地を語っています。本セッションは、東洋経済オンライ
ン並びに同社YouTubeで公開されています。
東洋経済オンライン記事：https://toyokeizai.net/articles/-/583900
Fixstars Amplify 代表取締役会長の三木が理事を務めている「量子技術に
よる新産業創出協議会（Q-STAR）」は、量子関連の産業・ビジネスの創出
を目的として設立された協議会です。IT企業で開発が進む量子技術を、実
際のビジネスにどのように役立てていくのかを研究し、量子技術を応用した
新産業を創出する取り組みを加速しています。
2022年5月23日、社団法人化を記念したシンポジウムでは、CEOの平岡が
登壇し「量子技術の産業化加速に向けた課題と取組み」というテーマでパネ
ルディスカッションを行いました。
画像：FixstarsCorporation youtubeチャンネルより
https://www.youtube.com/watch?v=uS-_BZqvqUo
代表取締役会長
三木が参加
代表取締役社長
平岡が参加

View Slide

組合せ最適化・イジングマシン及び
製造業における事例の紹介

View Slide

11
量子アニーリング・イジングマシンと組合せ最適化問題
量子アニーリング・イジングマシン ⇒ 組合せ最適化問題を解くための専用マシン
スケジューリング配送計画スマートシティー集積回路設計
膨大な選択肢から、制約条件を満たし、ベストな選択肢を探索する（組合せ最適化問題）
参考: 慶應義塾大学田中宗准教授「量子コンピュータ最前線とイジングマシンの可能性」

View Slide

12
目的関数:
制約条件:
解を取得
定式化
実装
問題設定
膨大な解候補（組合せ）から最適解を選ぶ
・解候補一つ一つの計算は可能
・候補数が膨大ですべての解候補を計算できない
これを最小化（最大化）する解が最適
解が必ず満たすべき条件
数式で表現
数式をPythonのプログラムで記述
Amplifyが最適解を探索
バイキングで最も安く
必要な栄養が取れる組合せは？
目的関数: 合計金額（最小化）
制約条件:
炭水化物：300g以上
タンパク質：150g以上
脂質：50g以上
(10g単位) 炭水化物タンパク質脂質金額
ごはん 8g 1g 1g 10円
パン 7g 1g 2g 12円
ハンバーグ 1g 5g 4g 50円
焼き魚 1g 8g 1g 35円
組合せ最適化問題の例
ごはん：380g
焼き魚：140g
金額：870円
最適メニュー
組合せ最適化問題は統一的なフレームワークで解くことができる

View Slide

13
製造業における品質管理の4M
本日の事例や
ワークショッ
プの対象
Man
Machine
Method
Material
4Mとは

View Slide

14
最も効率的な作業員のシフトは？
チーム編成とタスク割り当て
解を取得
定式化
実装
問題設定
目的関数：業務要求の満足度最適化
・スキルや業務時間平滑化、個人の希望などの考慮
制約条件：チーム人数や一人1タスクなどの制約
・従来、ホワイトボードを使い、朝・夕30分程度ずつかけて作成
・Amplifyが10秒で推薦配置を提示

View Slide

15
最も効率的な生産スケジュールは？
製品製造工程順序の最適化
October 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
Machine 1
Machine 2
Machine 3
Machine 4
Machine 5
Machine 6
Machine 7
Machine 8
Machine 9
Machine 10
解を取得
定式化
実装
問題設定
目的関数：製造完了時間（最小化）
・段取り時間、納期、稼動率などを考慮
制約条件：機械数や製造可能な製品の制約
・従来技術だと数時間かかることも
・Amplifyがより高速に妥当な解を提示

View Slide

16
装置1
装置2
装置x
A B B
B B B B
C E E D
A
電気機器製造メーカー A社様
複数の製品事業部から様々なプリント基板の注文を受け、生産を行う部門
生産する基板に応じて製造装置の部品や材料を交換する「段取り時間」
が必要。段取り時間を考慮した効率的な生産スケジュールを作成したい
従来は、専任者が、一日数回・毎回数十分かけて経験に基づいてスケジ
ュールを作成。更なる生産性向上やノウハウ継承のため、生産スケジュ
ール作成の自動化に着手
課
題
次期フェーズでは、Amplify
の活用領域の拡大を検討中！
最適化未経験のご担当者様1人がプログラム試作開始
から約1～2ヵ月間取り組んでこの効果を実現
現在は試作段階で、実運用に向けてモデルを改良中！
効
果段取りのための製造装置の停止回数の削減！
(10%以上削減)
段取り時間
段取り時間段取り時間
生産スケジュール作成の時間・コストの大幅な削減！
(一日あたり数時間 → 数分)

View Slide

17
最も効率的な部材や製品の搬送方法は？
AGV（無人搬送車）の待ち時間を最小化するリアルタイム経路選択
解を取得
定式化
実装
問題設定
目的関数：各AGVの次の一定時間の待ち時間（最小化）
制約条件：各AGVが衝突しない、など
・Amplifyが他のAGVの行動を考慮し、次の最適動作をリアルタイムに指示
・最短経路だけでなく、迂回・交差点制御なども実現

View Slide

Fixstars Amplify: メディア掲載
https://amplify.fixstars.com/ja/news/media
Software Design
「はじめての量子プログラミング体験」
（2021年6月号～2022年1月号まで連載）
Interface
「Pythonで体験！量子コンピュータ」
（2022年6月号）
18

View Slide

19
第二部ワークショップでの深掘りテーマ
シフト作成自動化生産計画最適化
人の直観で時間をかけて行っ
ていた生産ラインや物流倉庫
の業務シフト作成を、スキル
や勤務時間などの条件をもと
に最適化します
製造工場の設備の利用割り当て
(ジョブショップスケジューリ
ング) を、納期や段取り時間
などを考慮して最適化します
8/8セミナーのテーマ 9/9セミナーのテーマ
経路指示リアルタイム制御
倉庫を走行する多数の搬送ロ
ボット (AGV) が効率よく動
作するよう、最適経路だけで
なく迂回や交差点での待機な
どリアルタイムに指示します
7/25セミナーのテーマ

View Slide

Fixstars Amplify のご紹介

View Slide

CONFIDENTIAL
チーム
21
代表取締役社長CEO 平岡卓爾
2004年東京大学大学院工学系研究科物理工学専攻修
了（修士）
量子テレポーテーションの研究で有名な古澤研
究室の最初の学部生。ゼロから実験環境を構築
する稀有な機会を経験した。
2004年株式会社エリジオン入社
3次元CADデータ、形状処理ソフトウェアの開発。
工場DXに繋がる3次元点群処理ソフトウェア
「InfiPoints」の立ち上げを手掛けた。
2021年株式会社Fixstars Amplify入社
取締役CTO 松田佳希
2011年東京工業大学大学院物性物理学専攻博士課程
修了
量子アニーリングの理論研究で有名な西森研究
室の出身
2011年東京大学物性研究所助教
2013年株式会社フィックスターズ入社
2017年量子事業立ち上げ、経済産業省NEDO・内閣府
SIPの量子関連プロジェクトに参画。
量子コンピューティング技術を用いたアプリケ
ーション開発プラットフォーム・クラウドサー
ビスの開発責任者として事業開発に取り組む。
2020年早稲田大学グリーンコンピューティングシステ
ム研究機構客員講師（兼務）
2021年株式会社Fixstars Amplify CTO就任

View Slide

22
量子コンピュータとその周辺
1
量子コンピュータ
IBM/Google/Rigetti/IonQ
2
イジングマシン
富士通/日立/東芝/Fixstars
3
量子
アニーリング
D-Wave/NEC
1. 量子コンピュータ
量子ゲート方式
古典汎用コンピュータの上位互換。
量子力学の重ね合わせ状態を制御
する量子ゲートを操作し、特定の
問題を汎用的かつ高速に処理する。
2. イジングマシン
二値二次多項式模型
二次の多変数多項式で表される目
的関数の最適化問題 (QUBO) を扱
う
専用マシン。変数は0,1または±1。
統計物理学におけるイジング模型
(磁性体の性質を表す模型) に由来。
様々な実装により実現されている。
3. 量子アニーリング方式
量子焼きなまし法
イジングマシンの一種であり、量子焼
きなまし法の原理に基づいて動作する。
量子イジング模型を物理的に搭載した
プロセッサで実現する。
自然計算により低エネルギー状態が
出力される。
Amplify AE

View Slide

23
クラウドサービス：Fixstars Amplify
様々な量子コンピュータ・イジングマシンに対応したアプリケーションを開発出来る、クラウドプラットフォームです。
量子コンピューティング時代を見据え、シンプルで効率的な開発環境の提供を目指しています。
・・・
適用分野 (一例)
金融物流ライフサイエンス
様々なマシンに対応
各社が提供する量子アニーリング・イジングマシンを
Fixstars Amplifyから利用することが出来ます。
シンプルで効率的なアプリ開発
複雑で専門性の高いプロセスを自動化し、効率的にマシ
ンを使うための学習コストを、圧倒的に低くします。
PoCから実問題まで対応
大規模問題の入力と高速実行が可能で、PoCや実問題を
視野に入れたアプリケーション開発が行えます。
すぐに開発を開始可能
開発環境と実行環境がセットで提供されるため、すぐに
開発を開始することが出来ます。
主な特長
サービス概要

View Slide

24
二次計画問題
• 最適化問題の分類
• 組合せ最適化問題
• 決定変数が離散値 (整数など)
• 整数計画問題 (決定変数が整数)
• 0-1整数計画問題 (決定変数が二値)
• 連続最適化問題
• 決定変数が連続値 (実数など)
• 量子アニーリング・イジングマシン
Quadratic 二次形式
Unconstrained 制約条件なし
Binary 0-1整数 (二値)
Optimization 計画
𝑓 = ෍
𝑖𝑄𝑖𝑗
𝑞𝑖
𝑞𝑗
+ ෍
𝑖
𝑄𝑖𝑖
𝑞𝑖
𝑞𝑖
∈ 0,1 or 𝑞𝑖
∈ ±1
問題設定
(入力・定数)
決定変数
(出力・変数)
評価値
(目的関数)
• QUBO模型 (0-1整数二次計画問題)

View Slide

25
イジングマシンの実行手順
1. 数理モデル検討解きたい課題の「目的関数」「決定変数」「制約条件」を検討する
2. QUBO定式化 (論理)
「 2値決定変数+二次形式」で「目的関数」と「決定変数」を記述 (変換) する
「制約条件」は直接扱えないので「ペナルティ関数」で表現する
3. QUBO定式化 (物理)
各マシンの仕様や制限に準拠した形式にQUBO模型を変換する
(例: 二次項に制約がある場合は「グラフマイナー埋め込み」問題を解く)
4. 入力データの準備各マシンのSDKやAPI仕様に合わせてQUBO模型 (物理) をデータ化する
5. マシンの実行マシンを実行して出力の変数値やエネルギー値(コスト値)を解析する
上記の逆の手順を辿り解きたい課題の「決定変数」を解釈する
SWによる
支援と自動化

View Slide

26
Fixstars Amplify の特長
– いつでも開発環境と実行環境がセットのため
すぐにプログラミングと実行が出来る
– 誰でもハードウェアや専門的な知識が不要
無料で開発がスタート可能
– 高速に 26万ビットクラスの大規模問題の
高速処理と高速実行が可能
– あらゆる一般に公開されている全てのイジング
マシンを利用可能

View Slide

27
Fixstars Amplify クラウドの構成
Amplify クラウドが量子アニーリングマシン・イジ
ングマシンの大規模実行環境を提供
Amplify SDK で書かれた最適化エンジンは
他社製も含め全ての商用イジングマシンで実行可能
量子アニーリング・
イジングマシン (Amplify AE)
組合せ最適化問題の定式化やマシンを高度に操るた
めの最先端の技術や知見をライブラリ化したシンプ
ルなインターフェースを提供
専門知識が不要で量子アニーリングイジングマシン
を用いた最適化エンジンの開発に取り組める
Amplify SDK
組合せ最適化問題を含むアプリケーションは
Amplify SDK を用いることで効率的かつ簡便に組合
せ最適化の定式化が可能
研究・開発ではAmplifyクラウドは無償提供されるた
め直ぐに開発をスタートできる
アプリケーション

View Slide

28
Amplify の対応マシン
https://news.fixstars.com/2361/ : IBM Quantum対応 (6/14発表)
https://news.fixstars.com/2199/ : Gurobi対応 (5/18発表)

View Slide

29
Fixstars Amplifyによる課題解決
◼ イジングマシンのための革新的な開発環境
簡単多くのマシンに対応始めやすい
✓ SDKをインストールするだけ
ですぐに使える (pip install
amplify)
✓ ハードウェアの専門知識不要
でアプリケーションが開発で
きる
✓ 進化の早いマシンの発展に追
従すべての量子アニーリング/
イジングマシンに対応
✓ 26万ビット級のアニーリング
マシン実行環境が利用可能
✓ 研究・開発用途には開発環境
と実行環境が無償で利用可能
✓ 多くのチュートリアル、サン
プルコードを整備・拡充
多くの人が始めやすい開発プラットフォームを無料で提供
https://amplify.fixstars.com

View Slide

Fixstars Amplify の技術

View Slide

31
Amplify SDKのワークフロー
◼ 従来のプログラミングワークフローとの比較
実行方法

View Slide

32
Amplify SDKによるシンプルプログラミング
数独を解くサンプルアプリ富士通・デジタルアニーラの設定用コード
SDKなし
最適化しても
200行以上
出典: Wikipdia
SDKなし
59行
SDKあり
30行程度
SDKあり
1行
日立CMOSアニーリングマシンの設定用コード
SDKなし
183行
SDKあり
1行

View Slide

33
Amplify による巡回セールスマン問題 (TSP) の実装例

View Slide

34
Amplify Annealing Engine
• NVIDIA GPU V100/A100 で動作
• 独自の並列化シミュレーテッドアニーリングアルゴリズム
• WEB経由で計算機能を提供
• Amplify SDK の実装を直ぐに実行可能
• 社会課題への取り組み・PoC・検証が加速
• 商用マシンでは最大規模かつ最高速レベル
• 128Kビット (全結合) / 256Kビット超 (疎結合)
Amplify Cloud

View Slide

35
Amplify Annealing Engine
◼ 実問題やPoCを視野に入れたアプリケーション開発にも対応
研究・開発利用において無償提供

View Slide

36
Amplify 標準提供マシン
最先端の量子技術
最先端の並列化技術
NVIDIA A100/V100
最大26万ビット以上
D-Wave Advantage
5000量子ビット
出典:https://www.nvidia.com/ja-jp/data-center/dgx-2/

View Slide

【第二部】
ハンズオンワークショップ

View Slide

Q&Aについて
• 質問がある場合はZoomのチャット機能で「全員宛」で発信してくだ
さい
• その場で回答できるものは回答します
• または後ほど時間の許す限りお答えします
38

View Slide

39
ワークショップの準備 (1)
• ご自身のPC (ブラウザ上) でPythonプログラミングを行います。Google
colaboratoryを使うので、事前にログイン出来ることを確認をお願いします
(Googleアカウントが必要です)
• https://colab.research.google.com/
• Fixstars Amplify のトークンを取得済みか確認をお願いします。まだの人は
Fixstars Amplifyホームページのサインアップよりユーザ登録の上、無料トー
クンの取得をお願いします (1分で終わります)
• https://amplify.fixstars.com/ja/register

View Slide

40
• お持ちのAmplifyのトークンを用いて、以下のURLにあるサンプルコードが動くか確認
をお願いします。サンプルコードは閲覧のみ可能ですので、「ドライブにコピー」の上、
ご自身のトークンを入力し、Shift+Enterで実行をしてください (警告が出る場合がありますが、「このまま
実行」を選択下さい)
• https://colab.research.google.com/drive/1evYBKqKfVrEzrQOa-SWwciROfvqjL8qm?usp=sharing
この部分に、ご自身のトークン番号(32桁)を入力の上、Shift+Enterで実行下さい。
ご自身のトークン番号は、Amplify HPよりご確認いただけます
• ご自身のトークンを入力の上、以下の結果が出力されればOKです

View Slide

41
• ワークショップで使うサンプルコードを以下のURLより取得して下さい
• それぞれのサンプルコードにご自身のトークンを入力いただく必要があります。それぞ
れのサンプルコードを「ドライブにコピー」の上、トークンを入力し実行して下さい
Step1 https://colab.research.google.com/drive/1M9_fty7GQ4gPVa87lkLA9jdJR2OEqvvC?usp=sharing
Step2 https://colab.research.google.com/drive/1ZSgLMVgZTjIGODDy_YIJy29zr827OI6Y?usp=sharing
Step3 https://colab.research.google.com/drive/100KyhLplCh9oZ854BEXY0aadeE3IkSWl?usp=sharing
Step4 https://colab.research.google.com/drive/1wzffv95TAr1cx7j_WiYnXEWzU1z3uIZE?usp=sharing
➢ シフト最適化
➢ 数の分割問題
https://colab.research.google.com/drive/1zi_uMgxFfmIvr_ag-fNr6NVJBg5QWs94#scrollTo=xMdHzdYLPPT7

View Slide

・数の分割問題

View Slide

43
数の分割問題
【問題】{2,10,3,8,5,7,9,5,3,2}を、集合内の数の和がもう一方の集合内の数の和と等しくなるよう二
つの集合に分けよ
2,10,3,8,5,7,9,5,3,2
NP完全問題: とても難しい問題として知られている → 全通り試すしか方法は無い (210 = 1,024通り)
𝐴0
𝐴1

View Slide

44
手でやってみよう
まず適当に大きい順に並び替えて交互に一個発見！
分割方法は23通り存在する (対称を除く。重複を除くと9個の解が存在)
並び替え交互に
入れ
替え

View Slide

45
数の分割問題
目的関数の定式化
決定変数を準備
𝐴0
に割当てる場合は1
𝐴1
に割当てる場合は0
実装
𝑓 = ෍
𝑖=0
9
𝑎𝑖
∙ 𝑞𝑖
−
1
2
෍
𝑖=0
9
𝑎𝑖
2
i番目数決定変数
0 2 q_0 (0 or 1)
1 10 q_1 (0 or 1)
2 3 q_2 (0 or 1)
3 8 q_3 (0 or 1)
4 5 q_4 (0 or 1)
5 7 q_5 (0 or 1)
6 9 q_6 (0 or 1)
7 5 q_7 (0 or 1)
8 3 q_8 (0 or 1)
9 2 q_9 (0 or 1)
𝐴0
の合計から集合全
体の半分を引いた数を
最小化

View Slide

46
数の分割問題
求解
結果の取得
Amplify AE

View Slide

・最適シフト作成

View Slide

48
最適シフト作成
【問題】業務で求められる役割・役職・スキルと個人の能力や要求を考慮した最適シフトとは
業務要求
各従業員情報
・
・
・
要求管理者数要求スキル量（各従業員が持つスキル値の合計）
ライン責任者主任技師合計前工程後工程組立て合計
LineA 1 1 2 8 4 5 17
LineB 1 1 2 6 9 3 18
LineC 1 1 2 7 6 5 18
合計 3 3 6 21 19 13 53
要求配置(0=NG, 1=OK, 2=要求)
従業員ID LineA LineB LineC
0 1 0 0
1 1 1 1
2 1 2 0
23 1 1 1
24 1 1 0
担当可能役職スキル値(1=初級, 2=中級, 3=上級)
従業員ID 責任者主任技師技師前工程後工程組立て
0 1 1 1 3 3 3
1 0 1 1 3 3 2
2 0 0 1 2 1 0
・
・
・
・
・
・
23 0 0 1 0 0 2
24 0 0 1 1 0 1
マッチング

View Slide

49
① 従業員は同時に1ラインのみに配置可
② 各ラインが要求する管理職を配置することなど
目的関数:
制約条件:
最適シフト作成
解を取得
定式化
実装
問題設定各ラインの要求を満たすシフトの作成
複数の要求のバランス（充足率の平均の最大化、
分散の最小化、配置要求とのマッチングの最大化）
LineA LineB LineC 充足率

View Slide

50
自動配置
（アニーリング）
手動配置
一部事前配置
自動配置結果の微調整
各種条件を満たすよう
に、未配置のメンバー
を一括割り当て
事例紹介: 人員配置の効率化・自動最適化
関連インタビュー動画(4:00～):
https://www.youtube.com/watch?v=UyvomUej2Xw
充足率の平均の最大化、充足率の分散の最小化を
最適化する問題をイジングマシンで実行

View Slide

51
ワークショップ: 問題設定
15名の従業員を5名ずつ3つの生産ラインへ振り分けます。各従業員は、各ラインに対するスキル値を持
ちます。各ラインに配置された従業員のスキル値 (ラインスキル値) の合計がなるべく高く、また、各ラ
インのラインスキル値のばらつきが少ない、という2つの目的のバランスの取れたシフトの作成を目指し
ます。全てを一度にやるのは難しいので4つのステップに分けてアルゴリズムの完成を目指します
Aライン
(5名)
Bライン
(5名)
Cライン
(5名)
組合せは
約1400万通り!

View Slide

52
ワークショップ: 試しに人の手でやってみましょう

View Slide

53
ワークショップ: 4 Step
まず、2つの制約だけを考慮して配置シフトを求めます
制約➀: 従業員は同時に1つの製造ラインのみに配置が可能
制約②: 各ラインの配置人数が5名ずつになること
Step1に「ラインスキル値の合計 (= 全ラインの合計スキル値) を最大化」という目的➀を
追加し、複数の解の候補から目的を実現するシフトを求めます
解の候補多数あり
Step1
Step2
Step2に「ラインスキル値のばらつきを最小化」という目的②を追加して、2つの目的を
同時に実現するシフトを求めます
Step3
Step4
Step3に目的➀と目的②の重みを調整する「パラメーター」を追加し、最適なバランスの
シフトを作成します
15名の従業員を5名ずつ3つの生産ラインへ振り分けます。各従業員は、各ラインに対するスキル値を持
ちます。各ラインに配置された従業員のスキル値 (ラインスキル値) の合計がなるべく高く、また、各ラ
インのラインスキル値のばらつきが少ない、という2つの目的のバランスの取れたシフトの作成を目指し
ます。全てを一度にやるのは難しいので4つのステップに分けてアルゴリズムの完成を目指します

View Slide

54
Step1
Step1のサンプルコードのレビュー
(尚、本ワークショップでは、最適化のコードにフォー
カスし、下準備や可視化のコードの詳細は割愛します)

View Slide

55
Step1
決定変数の準備
実装
イジングマシン
による計算で最
適な(0,1)の組合
せを探す
従業員3は
line_Cに
配置
BinaryPoly型
(15×3)
= 45 [qbit]
1: 配置
0: 非配置
決定変数

View Slide

56
定式化
෍
𝑙
3
𝑞𝑖,𝑙
= 1
෍
𝑖
15
𝑞𝑖,𝑙
= 5
制約②: 各ラインの配置人数が要求人数（5名）
と一致すること → equal_to制約 (等式制約)
Step1
制約➀: 従業員は同時に1つの製造ラインのみに
配置が可能 → one_hot制約

View Slide

57
Step1
実装
https://amplify.fixstars.com/ja/demo
【補足】
制約条件の取り扱いに関する詳細は、
こちらにあるチュートリアルも合わ
せてご参照下さい

View Slide

58
Step1
求解
無料版は1ジョブ10秒まで設定可有
料版では1分まで設定可能
Amplify AE
• modelに格納してマ
シンに投げます
• 制約条件だけを与え
た場合、制約条件を
満たす解を探してき
てくれます

View Slide

59
Step1
結果の取得
可視化
確
認
同じ人が複数のラインに配置されることなく、各ラインに5人ずつ配置するという二つの制約を満たすシフトを作ることが
できました。但し、解の候補はたくさんあり、最適化の余地も大きいそうです

View Slide

60
Step2 Step1に「ラインスキル値の合計 (= 全ラインの合計スキル値) を最大化」と
いう目的➀を追加し、複数の解の候補から目的を実現するシフトを求めます
ここでは、工場全体のアウトプットの最大化を目指し、各従業員をできるだけ高いスキル値を持つラインへ配置する
ことを目指します
Aライン(5名)
ラインスキル値
Bライン(5名)
Cライン(5名)
全ラインのスキ
ル値の合計
これを最大化
したい！

View Slide

61
Step2のサンプルコードのレビュー
Step2 Step1に「ラインスキル値の合計 (= 全ラインの合計スキル値) を最大化」と

View Slide

62
Step2
定式化
各列の要素同士の掛け算の和 ⇒ 各ラインのラインスキル値
Step1に「ラインスキル値の合計 (= 全ラインの合計スキル値) を最大化」と
𝑠𝑘𝑖𝑙𝑙_𝑠𝑐𝑜𝑟𝑒 = ෍
𝑙
3
෍
𝑖
15
𝑞𝑖,𝑙
∙ 𝑠𝑖,𝑙
決定変数
目的➀: 全ラインの合計スキル値の最大化

View Slide

63
Step2
実装
制約条件には適切な値の重みを設定する必要があります。典型的
なスキル値より大きくしておく必要があるため 150 としました
: 追加コード
イジングマシンは、このobjectiveの値が最小になる組合せを探し
ます (全ラインの合計スキル値は大きいものを選びたいので、
skill_scoreにマイナスをつけたものをobjectiveとしています)

View Slide

64
Step2
結果の取得
可視化
確
認
ばらつき大
二つの制約を満しながら、全ラインの合計スキル値が最大化されたシフトを作ることができました (Step1の全ラインの合
計スキル値は1,265) 。但し、ライン間のスキル値のばらつきが大きいので、更なる最適化をかけたい状況です

View Slide

65
Step3 Step2に「ラインスキル値のばらつきを最小化」という目的②を追加して、
2つの目的を同時に実現するシフトを求めます
目的②: ラインスキル値のばらつき (分散)
を最小化
定式化実装
𝑣𝑎𝑟𝑖𝑎𝑛𝑐𝑒 =
σ𝑙
3 σ
𝑖
15 𝑞𝑖,𝑙
∙ 𝑠𝑖,𝑙
2
3
−
σ𝑙
3 σ
𝑖
15 𝑞𝑖,𝑙
∙ 𝑠𝑖,𝑙
3
2
: 追加コード
ばらつきは小さい方が選ばれるようにしたいのでプラスで足します

View Slide

66
結果の取得
可視化
確
認
Step3 Step2に「ラインスキル値のばらつきを最小化」という目的②を追加して、
2つの目的を同時に実現するシフトを求めます
ばらつきなし！
二つの制約を満しながら、全ラインの合計スキルが高く、各ライン間のばらつきが全くないシフトを作ることができました
（Step1の全ラインの合計スキル値は1,265で、Step2は1,475）。最後に、目的➀と目的②のバランスをチューニングし
て最適なシフトの作成を目指します。

View Slide

67
Step4 Step3に目的➀と目的②の重みを調整する「パラメーター」を追加し、
最適なバランスのシフトを作成します
実装
: 追加コード

View Slide

68
最適なバランスのシフトが完成しました！
結果の取得
可視化
確
認
Step4 Step3に目的➀と目的②の重みを調整する「パラメーター」を追加し、
最適なバランスのシフトを作成します
ばらつき極小！

View Slide

69
Step1
スキル値合計: 1,265
ワークショップ: おさらい
制約のみからスタートして、複数の目的を加え、重みを調整することで最適なバランスのシフトを作りました
スキル値最大化
スキル値ばらつき調整
重みを調整
制約
Step2
Step3
Step4
制約
制約
スキル値ばらつき調整
制約

View Slide

70
類似の実アプリの紹介

View Slide

課題2
最近、従業員番号0の田中さんと従業員番号1の佐藤さんの間
にトラブルが発生して、この二人を同じラインに配属できな
くなりました。この二人を同じチームには入れないという条
件を加えて、ワークショップと同じ目的を実現するための最
適シフトを求めるためには、サンプルコードをどのように変
えればよいでしょうか
課題1
各従業員の各ラインのスキル値を10倍の値とした場合に、ワークショップと同じ目的を実現するための最
適シフトを求めて下さい (スキル値の値を変えるだけだと解は求まりません)
発展課題

View Slide

72
Fixstars Amplify: オンラインデモ & チュートリアル
https://amplify.fixstars.com/ja/demo

View Slide

7月25日(月)
第一部 13:00~13:50
- フィックスターズの紹介
- 組合せ最適化問題・イジングマシン
の紹介
- Fixstars Amplifyの紹介
第二部 14:00~15:30
- Fixstars Amplifyを用いた搬送経
路最適化のワークショップ
- Q&A
第三部 15:30~16:00
- 今後について
- 採用情報
8月８日(月)
第一部 13:00~13:50
の紹介
第二部 14:00~15:30
- Fixstars Amplifyを用いたシフト
最適化のワークショップ
- Q&A
第三部 15:30~16:00
- 採用情報
9月９日(金)
第一部 13:00~13:50
の紹介
第二部 14:00~15:30
- Fixstars Amplifyを用いた生産計
画最適化のワークショップ
- Q&A
第三部 15:30~16:00
- 採用情報
次回のご案内

View Slide

【第三部】
今後について

View Slide

ぜひ情報発信を！

View Slide

課題2
最近、従業員番号0の田中さんと従業員番号1の佐藤さんの間
にトラブルが発生して、この二人を同じラインに配属できな
くなりました。この二人を同じチームには入れないという条
件を加えて、ワークショップと同じ目的を実現するための最
適シフトを求めるためには、サンプルコードをどのように変
えればよいでしょうか
課題1
各従業員の各ラインのスキル値を10倍の値とした場合に、ワークショップと同じ目的を実現するための最
適シフトを求めて下さい (スキル値の値を変えるだけだと解は求まりません)
発展課題

View Slide

SDKの復習や発展課題にチャレンジ
ブログやQiita等で公開された場合はぜひご連絡ください
https://sparse-dense.blogspot.com/2022/06/ibm-quantumfixstars-
amplify.html

View Slide

研究での活用事例

View Slide

研究での活用事例
大学研究室概要参考
慶應義塾大学村松研究室 Phase-fieldモデルの量子アニーリングシミュレー
タ (ジブロックコポリマー相分離構造への適用)
https://www.nature.com/articles/s41598-
022-14735-4
東京大学
/JFLI
Philippe
Codognet
Modeling the Costas Array Problem in QUBO
for Quantum Annealing
https://link.springer.com/chapter/10.1007/
978-3-031-04148-8_10
名古屋大学片桐研究室 Amplifyを用いたCMOSアニーリングマシンの特性
の分析
https://ipsj.ixsq.nii.ac.jp/ej/?action=pages_
view_main&active_action=repository_view_
main_item_detail&item_id=213136&item_n
o=1&page_id=13&block_id=8
東北大学小松研究室組み合わせクラスタリングによるアニーリングマ
シンの評価
https://ipsj.ixsq.nii.ac.jp/ej/?action=pages_
view_main&active_action=repository_view_
main_item_detail&item_id=218973&item_n
o=1&page_id=13&block_id=8
山梨大学鈴木研究室量子アニーリングによる疎行列直接解法向けフィ
ルイン削減オーダリング
https://ieeexplore.ieee.org/document/9691
960
アカデミックの研究においても Amplify が利用されています！

View Slide

Amplifyの業務

View Slide

81
Amplify開発に必要とされる技術
Amplify SDK

View Slide

82
Amplify開発に必要とされる技術
Amplify SDK
• C/C++による Python SDK 開発
• 各種イジングマシンや量子コンピュー
タの活用
• マルチコアCPU・マルチGPU環境にお
ける並列計算（CUDA・C/C++）
• 安定したサービス運用
• ユーザー課題の理解と組合せ最適化問
題への定式化
• PythonやC/C++を用いたアルゴリズム
開発
• Web技術やクラウド技術を活用したア
プリケーション開発

View Slide

採用情報

View Slide

就業体系
84
裁量労働制もしくはフレックスタイム制を適用しています
※入社時のスキルによって1、2のいずれか
所定労働時間 8時間休憩1時間
1 裁量労働制
推奨：始業時間10:00、終業時間19:00
※専門業務型裁量労働制により9時間働いたものとみなします
固定残業：30時間
2 フレックスタイム制
コアタイム 10:00～14:00
固定残業：15時間

View Slide

福利厚生：各種手当
85
子供手当交通費ドリンク補助
業績・成果により変動
扶養義務のあるお子様
一人につき 5,000円/月
全額支給
上限 60,000円/月
奨励金 100％
※本社株式
社内自販機 10円
賞与従業員持株会
SO制度
ストックオプション制度あり

View Slide

福利厚生：休日休暇
86
完全週休2日制
土日祝日 10～20日 3日
夏季休暇
産前産後・介護休暇
慶弔休暇
有給休暇
年末年始休暇

View Slide

社外活動支援
87
※本社規定に準じます
スキルアップ
サポート
セミナーや勉強会、学会発表・参加、プロコンなどの
社外活動を奨励しています。
MBA/ PhD
学費・交通費補助
MBAやPhDの取得希望者に学費と交通費等のサポートを
行っています。
実績：MBA 早稲田大学 2名
PhD 名古屋大学 1名・北海道大学 1名
語学研修補助
社員の外国語コミュニケーション向上のため、
オンラインレッスン（英語・日本語）受講のための
サポートを行っています。
資格取得補助業務に関連のある資格に対して、合格奨励金の支給や
受験料の補助など、各種サポートが受けられます。

View Slide

社内イベント / 制度
88
※本社規定に準じます
社内勉強会
エンジニア主体で週1回以上開かれます。
2021年度：Intel Software Developer Manualを読む会
並列(分散)アルゴリズム勉強会
社内プロコン
年に1度開催します。
プロコン勢が多いため盛り上がります。
2021年：「イジングモデル最適化」
社内交流
おやつタイム・社内新聞・アドベントカレンダーなど
社内の情報を共有したり、メンバー同士の交流の場が
随時あります。
社内大学
知見のある社員が講師となり開講。
業務時間内に受けることができます。
（詳細は次ページへ）

View Slide

成長支援（社内大学2020）
89
テーマや内容は年によって変更となることがあります。
講座名テーマ内容
高速化講座 CPUを使った高速化
x86を例にした現代のCPUアーキテクチャと
最適化技法についての理解と習得
CUDAを使った高速化
CUDAを使った高速なプログラムの書き方、
CUDAプログラムの解析方法、高速化方法の習得
性能モデルに基づく
高速化
性能モデルに基づいたメモリ階層を活用する
高速化手法やCompute-intensiveな処理の
高速化手法の習得
AI講座 AI・機械学習
機械学習アルゴリズムの概念及び深層学習技術の習
得、AI・深層学習技術を利用した問題解決の実装
CV講座コンピュータビジョン古典的な画像処理による手法から、より新しい
機械学習ベースの手法まで、コンピュータービジョ
ン関連技術の理解と習得

View Slide

使用ツール
90
G
Suite
全社員使用ツール
・こだわりのキーボードやマウスが使えます。
・開発に必要なソフトウェア、技術書の購入に特に制限はありません。
・各種オープンソースや開発ツール、テストツールなども自由です。
・申請は必要ですが、自分の機材を持ち込むことも可能です。

View Slide

本棚
91
本棚の本、技術雑誌は自由に読むことができます。
また、本の購入にも制限はなく相談の上、簡単に購入することができます。
技術書以外にも、漫画やボードゲームもあります。

View Slide

開発環境
92
大きな机
高性能な椅子
複数モニター（2枚以上）
リラクゼーションスペース
マッサージチェア
10円ドリンク
ウォーターサーバー

View Slide

選考の流れ

View Slide

新卒選考
94
応募書類選考
オファー
面談
一次面接
適性検査
最終面接
応募
弊社HPの応募フォームよりご応募ください。
書類選考
お送りいただいた内容に基づき書類選考を行います。
一次面接
担当者、役職者との面接を受けていただきます。
※ 研究発表（15分）・ライブコーディング試験を実施します。

View Slide

新卒選考
95
応募書類選考
オファー
面談
一次面接
適性検査
最終面接
最終面接
社員全員との面接を受けていただきます。
オファー面談
最終面接を通過された方は、条件提示の面談をさせていただきます。

View Slide

ライブコーディング試験
96
1 2 3
問題理解コーディングコード解説
問題の読解・
アルゴリズムの理解
Wandboxを使用して
コーディング
記述いただいた
コードの説明
◆ 時間は30分程度 ◆疑問点は随時質問いただけます
▽ネット検索はご遠慮いただいております

View Slide

Wandbox
97

View Slide

Thank You
お問い合わせ窓口 : [email protected]

View Slide

1day internship Fixstars Amplify 20220808/組み合わせ最適化講習会

1day internship Fixstars Amplify 20220808/組み合わせ最適化講習会

fixstars

More Decks by fixstars

Other Decks in Programming

Featured

Transcript