歴史 巡回セールスマン問題の起源は不明である。1832年の巡回セールスマン向けハンドブックにはこの問題が言及されており、ドイツ とスイス を巡る例が掲載されているが、数学的な扱いは含まれていない。[ 3 ]
ウィリアム・ローワン・ハミルトン、1850年頃 TSPは、19世紀にアイルランドの数学者ウィリアム・ローワン・ハミルトン とイギリスの数学者トーマス・カークマン によって数学的に定式化されました。ハミルトンのイコシアンゲームは、 ハミルトン閉路を 見つけることに基づく娯楽パズルでした。[ 4 ] TSPの一般形は、1930年代にウィーンとハーバード の数学者によって初めて研究されたようで、特にカール・メンガー は、問題を定義し、明らかな総当たりアルゴリズムを検討し、 最近傍 ヒューリスティックの非最適性を指摘しました。
我々は、有限個の点間の距離が既知である点について、それらの点を結ぶ最短経路を見つけるタスクを、メッセンジャー問題 と呼ぶ(実際にはこの問題は各郵便配達員、いずれにせよ多くの旅行者によって解決されるべきである)。もちろん、この問題は有限回の試行で解決できる。試行回数を与えられた点の順列の数以下に抑えるルールは知られていない。出発点から最も近い点へ、次にその最も近い点へ、といったように進むというルールは、一般に最短経路を与えない。[ 5 ]
この問題が数学的に初めて検討されたのは、1930年代にスクールバスの経路問題を解決しようとしていたメリル・M・フラッドでした。 [ 6 ] プリンストン大学 のハスラー・ホイットニーはこ の問題に関心を持ち、「48州問題」と名付けました。「巡回セールスマン問題」というフレーズを使用した最初の出版物は、1949年のランド研究所の ジュリア・ロビンソン による報告書「ハミルトンゲーム(巡回セールスマン問題)について」です。[ 7 ] [ 8 ]
1950年代と1960年代には、サンタモニカ のランド研究所が 問題解決のステップに対して賞を提供したことをきっかけに、この問題はヨーロッパとアメリカの科学界でますます人気が高まった。 [ 6 ] ランド研究所のジョージ・ダンツィヒ 、デルバート・レイ・フルカーソン 、セルマー・M・ジョンソン は、この問題を整数線形計画 問題として表現し、その解法として切除平面 法を開発した。彼らは、この主題に関する先駆的な論文を執筆し、その中で、これらの新しい方法を用いて、巡回路を構築し、他の巡回路がそれより短くならないことを証明することで、49都市のインスタンスを最適解に導いた。しかし、ダンツィヒ、フルカーソン、ジョンソンは、ほぼ最適な解が与えられた場合、少数の追加の不等式(カット)を追加することで最適解を見つけたり、最適解を証明したりできるのではないかと推測した。彼らはこの考えを用いて、ストリングモデルを使用して最初の49都市の問題を解決した。彼らは、49都市問題の解決には26回のカットしか必要ないことを発見した。この論文はTSP問題に対するアルゴリズム的アプローチを示したものではないが、その中に含まれるアイデアは、後にTSPの正確な解法を作成する上で不可欠であった。ただし、これらのカットを作成するためのアルゴリズム的アプローチを見つけるには15年かかることになる。[ 6 ] ダンツィヒ、フルカーソン、ジョンソンは、おそらく初めて、カット平面法に加えて分岐限定法アルゴリズムを使用した。 [ 6 ]
1959年、ジリアン・ビアードウッド 、JH・ハルトン、ジョン・ハマーズリーは、 ケンブリッジ哲学協会 誌に「多数の点を通る最短経路」というタイトルの記事を発表しました。ビアードウッド・ハルトン・ハマーズリーの定理は、巡回セールスマン問題に対する実用的な解決策を提供します。著者らは、自宅やオフィスから出発し、一定数の場所を訪れてから出発地点に戻るセールスマンの最短経路の長さを決定するための漸近式を導き出しました。
その後数十年間、この問題は数学 、コンピュータ科学 、化学 、物理学 、その他の科学分野の多くの研究者によって研究されました。しかし、1960年代には、最適な解を求めるのではなく、最適な長さの倍数で長さが制限されることが証明できる解を生成する新しいアプローチが考案され、そうすることで問題の下限が作成されました。これらの下限は、分岐限定法のアプローチで使用されます。これを行う方法の1つは、グラフの最小全域木 を作成し、そのすべてのエッジを2倍にすることで、最適なツアーの長さが最小全域木の重みの最大2倍であるという境界が得られることです。[ 6 ]
1976 年、Christofides と Serdyukov は (互いに独立して) この方向で大きな進歩を遂げました。[ 10 ] Christofides –Serdyukov アルゴリズムは 、最悪の場合でも最適解より最大 1.5 倍長い解を生成します。このアルゴリズムはシンプルで高速であったため、多くの人がほぼ最適な解法に取って代わることを期待しました。しかし、この改善への期待はすぐには実現せず、Christofides–Serdyukov アルゴリズムは、2011 年に「グラフィカル」 TSP のサブセットに対して (ごくわずかに) 改善された近似アルゴリズムが開発されるまで、最悪のシナリオで最良の方法であり続けました。[ 11 ] 2020 年に、このわずかな改善が完全な (メトリック) TSP に拡張されました。[ 12 ] [ 13 ]
リチャード・M・カープは 1972年にハミルトン閉路 問題がNP完全で あることを示し、これは巡回セールスマン問題(TSP)のNP困難性 を意味する。これにより、最適経路を見つける際の計算上の困難さに対する数学的な説明が得られた。
1970年代後半から1980年代にかけて大きな進歩があり、グロッチェル 、パドバーグ、リナルディらは、切除平面法と 分岐限定法 を用いて、最大2,392都市のインスタンスを正確に解くことに成功した。
1990年代には、Applegate 、Bixby 、Chvátal 、Cookらが、 最近の多くの記録的解法で使用されているプログラムConcorde を開発しました。Gerhard Reineltは1991年に、さまざまな難易度のベンチマークインスタンスのコレクションであるTSPLIBを公開し、多くの研究グループが結果を比較するために使用してきました。2006年、Cookらは、マイクロチップレイアウト問題で与えられた85,900都市のインスタンスを通る最適なツアーを計算しました。これは現在、解決されたTSPLIBインスタンスの中で最大です。数百万の都市を持つ他の多くのインスタンスでは、最適なツアーから2~3%以内の解が確実に見つかります。[ 14 ]
説明
グラフ問題として 4都市の対称型TSP TSPは、都市を頂点、 経路を辺 、経路の距離を辺の重みとする無向重み付きグラフ としてモデル化できます。これは、指定された頂点を出発点とし、他のすべての頂点を ちょうど一度ずつ訪問した後に、その頂点 で終了するという最小化問題です。多くの場合、モデルは完全グラフ (つまり、すべての頂点のペアが辺で接続されている)です。2つの都市間に経路が存在しない場合、十分な長さの辺を追加することで、最適な経路に影響を与えることなくグラフを完成させることができます。
非対称と対称 対称型TSP では、2つの都市間の距離はどちらの方向でも同じで、無向グラフ を形成します。この対称性により、可能な解の数は半分になります。非対称型TSP では、両方向に経路が存在しない場合や、距離が異なる場合があり、有向グラフ を形成します。交通渋滞、一方通行の道路、出発地と到着地で料金が異なる都市間の航空運賃などは、非対称型のTSP問題を生み出す可能性のある現実的な要因です。
グラフ理論 の観点から同等の定式化を行うと、次のようになります。完全重み付きグラフ (頂点は都市、辺は道路、重みはその道路のコストまたは距離を表す)が与えられたとき、重みが最小のハミルトン閉路 を見つけます。これは、完全でない重みなしグラフにハミルトン路(または閉路)が存在するかどうかを問うハミルトン路問題 よりも一般的です。出発都市に戻るという要件は、問題の計算複雑性を変化させません。 ハミルトン経路問題を 参照してください。 関連するもう 1 つの問題は、ボトルネック巡回セールスマン問題です。これは、最も重みのある エッジ の重みが最小となる重み付きグラフ のハミルトン閉路を見つける問題です。現実世界の例としては、大型バスで狭い道路を避けることが挙げられます。[ 15 ] この問題は、明らかな輸送および物流分野以外にも、実用上非常に重要な問題です。典型的な例は、プリント基板製造です。PCB に穴を開けるための ドリル マシンのルートのスケジューリングです。ロボットによる機械加工や穴あけアプリケーションでは、「都市」は機械加工する部品または穴あけする (さまざまなサイズの) 穴であり、「移動コスト」にはロボットの再調整時間が含まれます (単一機械ジョブシーケンス問題)。[ 16 ] 一般化された巡回セールスマン問題 (「巡回政治家問題」とも呼ばれる)は、(1 つ以上の)「都市」を持つ「州」を扱い、セールスマンは各州からちょうど 1 つの都市を訪問しなければならない。 1 つの応用例は、刃の交換を最小限に抑えるために、切削ストック問題の解を順序付ける際に見られる。もう 1 つは、 半導体 製造における穴あけに関するものである。たとえば、米国特許第 7,054,798 号 を参照。 Noon と Bean は、一般化された巡回セールスマン問題を、都市の数は同じだが距離行列 が変更された標準的な TSP に変換できることを示した。 順序付け問題とは、都市間に優先順位関係が存在する一連の都市を訪問する問題を扱うものである。 Google の面接でよく聞かれる質問の一つに、データ処理ノード間でデータをどのようにルーティングするかというものがあります。ルーティングはデータ転送にかかる時間によって異なりますが、ノードの計算能力やストレージ容量も異なるため、データの送信先を決定する問題が複雑化します。巡回購買者問題 とは、一連の製品を購入する責任を負う購買者に関する問題です。購買者は複数の都市でこれらの製品を購入できますが、価格は都市によって異なり、またすべての都市で同じ製品が販売されているわけではありません。目的は、総コスト(移動コスト+購入コスト)を最小限に抑えつつ、必要なすべての製品を購入できるような、都市間の経路を見つけることです。
TSPは整数線形計画 問題として定式化できます。[ 17 ] [ 18 ] [ 19 ] いくつかの定式化が知られています。注目すべき定式化は、Miller–Tucker–Zemlin (MTZ) 定式化と Dantzig–Fulkerson–Johnson (DFJ) 定式化です。DFJ 定式化の方が強力ですが、MTZ 定式化も特定の状況では依然として有用です。[ 20 ] [ 21 ]
これら2つの定式化に共通するのは、都市に番号を付けることである。1 、 … 、 n {\displaystyle 1,\ldots ,n} そして、c 私 j > 0 {\displaystyle c_{ij}>0} 都市からのコスト(距離)私 {\displaystyle i} 市内へj {\displaystyle j} 式における主な変数は以下のとおりです。
x 私 j = { 1 道は都市から 私 市内へ j 0 さもないと。 {\displaystyle x_{ij}={\begin{cases}1&{\text{経路は都市 }}i{\text{から都市 }}jへ向かう。\\0&{\text{それ以外の場合。}}\end{cases}}} これらの変数が0/1であるため、定式化は整数計画問題となります。その他の制約はすべて純粋に線形です。特に、このプログラムの目的は、巡回経路の長さを最小化することです。
∑ 私 = 1 n ∑ j ≠ 私 、 j = 1 n c 私 j x 私 j 。 {\displaystyle \sum _{i=1}^{n}\sum _{j\neq i,j=1}^{n}c_{ij}x_{ij}.} さらなる制約がなければ、{ x 私 j } 私 、 j {\displaystyle \{x_{ij}\}_{i,j}} は、エッジの集合のすべての部分集合を効果的に網羅しますが、これはツアーのエッジの集合とは非常に異なり、すべてのエッジが網羅される自明な最小値を可能にします。x 私 j = 0 {\displaystyle x_{ij}=0} したがって、どちらの定式化にも、各頂点には正確に 1 つの入力エッジと 1 つの出力エッジが存在するという制約があり、これは次のように表現できます。2 n {\displaystyle 2n} 線形方程式
∑ 私 = 1 、 私 ≠ j n x 私 j = 1 {\displaystyle \sum _{i=1,i\neq j}^{n}x_{ij}=1} のためにj = 1 、 … 、 n {\displaystyle j=1,\ldots ,n} そして∑ j = 1 、 j ≠ 私 n x 私 j = 1 {\displaystyle \sum _{j=1,j\neq i}^{n}x_{ij}=1} のために私 = 1 、 … 、 n 。 {\displaystyle i=1,\ldots ,n.} これらは、選択されたエッジの集合が局所的にツアーのように見えることを保証する一方で、選択されたエッジが複数のツアーを構成し、それぞれが頂点のサブセットのみを訪れる可能性があるため、すべての頂点を訪れるツアーが1つ 存在するという全体的な要件に違反する解も許容します。おそらく、この全体的な要件こそがTSPを難しい問題にしているのです。MTZとDFJの定式化は、この最後の要件を線形制約としてどのように表現するかが異なります。
それに加えてx 私 j {\displaystyle x_{ij}} 上記の変数にはそれぞれ私 = 1 、 … 、 n {\displaystyle i=1,\ldots ,n} ダミー変数u 私 u_i 都市訪問の順序を記録し、都市から数える1 {\displaystyle 1} 解釈 はu 私 < u j u_{i}<u_{j}} 都市を意味する私 {\displaystyle i} 都市を訪れる前にj 。 {\displaystyle j.} 特定のツアー(値にエンコードされたもの)x 私 j {\displaystyle x_{ij}} 変数)に対して満足のいく値を見つけることができるu 私 u_i 変数を作成するu 私 u_i 都市から出発する場合、そのツアーに沿ったエッジの数に等しい。1 {\displaystyle 1} 市内へ私 。 {\displaystyle i.} [ 22 ]
線形計画法は非厳密な不等式を好むため(≥ {\displaystyle \geq } ) 厳格すぎる(> {\displaystyle >} ) 我々は、次のような制約を課したいと考えている。
u j ≥ u 私 + 1 {\displaystyle u_{j}\geq u_{i}+1} もしx 私 j = 1. {\displaystyle x_{ij}=1.} 単に要求するだけでu j ≥ u 私 + x 私 j {\displaystyle u_{j}\geq u_{i}+x_{ij}} それでは達成できない 、なぜならそれにはu j ≥ u 私 {\displaystyle u_{j}\geq u_{i}} いつx 私 j = 0 、 {\displaystyle x_{ij}=0,} これは正しくありません。代わりにMTZはn ( n − 1 ) {\displaystyle n(n-1)} 線形制約
u 私 − u j + 1 ≤ ( n − 1 ) ( 1 − x 私 j ) {\displaystyle u_{i}-u_{j}+1\leq (n-1)(1-x_{ij})} すべての異なる私 、 j ∈ { 2 、 … 、 n } 、 {\displaystyle i,j\in \{2,\dotsc ,n\},} ここで定数項 n − 1 {\displaystyle n-1} 十分な余裕を持たせることでx 私 j = 0 {\displaystyle x_{ij}=0} 関係を課さないu j {\displaystyle u_{j}} そしてu 私 。 {\displaystyle u_{i}.}
そのやり方u 私 {\displaystyle u_{i}} 変数は、単一のツアーがすべての都市を訪れることを強制し、少なくとも増加します1 {\displaystyle 1} ツアーの各ステップごとに、ツアーが都市を通過する場合にのみ減額が許可される 1. {\displaystyle 1.} その制約は、都市を通らないすべてのツアーで違反されるだろう。 1 、 {\displaystyle 1,} だからそれを満たす唯一の方法は、ツアーが通過する都市 1 {\displaystyle 1} また、他のすべての都市も通過します。
したがって、MTZによるTSPの定式化は、以下の整数線形計画問題となる。
ミニ ∑ 私 = 1 n ∑ j ≠ 私 、 j = 1 n c 私 j x 私 j : x 私 j ∈ { 0 、 1 } 私 、 j = 1 、 … 、 n ; ∑ 私 = 1 、 私 ≠ j n x 私 j = 1 j = 1 、 … 、 n ; ∑ j = 1 、 j ≠ 私 n x 私 j = 1 私 = 1 、 … 、 n ; u 私 − u j + 1 ≤ ( n − 1 ) ( 1 − x 私 j ) 2 ≤ 私 ≠ j ≤ n ; 2 ≤ u 私 ≤ n 2 ≤ 私 ≤ n 。 {\displaystyle {\begin{aligned}\min \sum _{i=1}^{n}\sum _{j\neq i,j=1}^{n}c_{ij}x_{ij}&\colon &&\\x_{ij}\in {}&\{0,1\}&&i,j=1,\ldots ,n;\\\sum _{i=1,i\neq j}^{n}x_{ij}={}&1&&j=1,\ldots ,n;\\\sum _{j=1,j\neq i}^{n}x_{ij}={}&1&&i=1,\ldots ,n;\\u_{i}-u_{j}+1\leq {}&(n-1)(1-x_{ij})&&2\leq i\neq j\leq n;\\2\leq u_{i}\leq {}&n&&2\leq i\leq n.\end{aligned}}} 最初の等式は、各都市へは他の都市からちょうど1つだけ到達できることを要求し、2番目の等式は、各都市から他の都市へはちょうど1つだけ出発できることを要求する。最後の制約は、すべての都市を網羅するツアーは1つだけであり、すべての都市をまとめて網羅する2つ以上の互いに素なツアーは存在しないことを保証する。
都市に1、…、nの 番号を付け、以下を定義してください。
x 私 j = { 1 道は都市から 私 市内へ j 0 さもないと。 {\displaystyle x_{ij}={\begin{cases}1&{\text{the path goes from city }}i{\text{ to city }}j\\0&{\text{otherwise.}}\end{cases}}} 取るc 私 j > 0 {\displaystyle c_{ij}>0} を都市i から都市j までの距離とする。すると、TSP は次の整数線形計画問題として記述できる。
ミニ ∑ 私 = 1 n ∑ j ≠ 私 、 j = 1 n c 私 j x 私 j : ∑ 私 = 1 、 私 ≠ j n x 私 j = 1 j = 1 、 … 、 n ; ∑ j = 1 、 j ≠ 私 n x 私 j = 1 私 = 1 、 … 、 n ; ∑ 私 ∈ Q ∑ j ≠ 私 、 j ∈ Q x 私 j ≤ | Q | − 1 ∀ Q ⊊ { 1 、 … 、 n } 、 | Q | ≥ 2. {\displaystyle {\begin{aligned}\min &\sum _{i=1}^{n}\sum _{j\neq i,j=1}^{n}c_{ij}x_{ij}\colon &&\\&\sum _{i=1,i\neq j}^{n}x_{ij}=1&&j=1,\ldots ,n;\\&\sum _{j=1,j\neq i}^{n}x_{ij}=1&&i=1,\ldots ,n;\\&\sum _{i\in Q}{\sum _{j\neq i,j\in Q}{x_{ij}}}\leq |Q|-1&&\forall Q\subsetneq \{1,\ldots ,n\},|Q|\geq 2.\\\end{aligned}}} DFJ 定式化の最後の制約 (部分巡回路除去 制約と呼ばれる) は、どの真部分集合 Q も部分巡回路を形成できないことを保証するため、返される解は単一の巡回路であり、より小さな巡回路の和集合ではありません。直感的には、都市の各真部分集合 Q に対して、この制約は Q 内の都市の数よりもエッジの数が少ないことを要求します。Q 内のエッジの数と Q 内の都市の数が同じであれば、それは Q 内の都市の部分巡回路を表すことになります。これは可能な制約の数を指数関数的に増加させるため、実際には行生成 によって解決されます。[ 23 ]
解を計算する NP困難問題に対する従来のアプローチは以下のとおりです。
正確なアルゴリズム を考案するが、それは問題の規模が小さい場合にのみ、比較的高速に動作する。「準最適」またはヒューリスティックなアルゴリズム 、つまり妥当な時間内に近似解を提供するアルゴリズムを考案する。 問題の特殊なケース(「部分問題」)を見つけ出し、それに対してより優れた、あるいは正確なヒューリスティックが可能であることを示す。
正確なアルゴリズム 最も直接的な解決策は、すべての順列 (順序付き組み合わせ)を試して、どれが最も安価かを確認することです(総当たり探索 を使用)。このアプローチの実行時間は、多項式係数の範囲内です。O ( n ! ) {\displaystyle O(n!)} 、これは都市数の階乗なので、この解決策は20都市だけでも非現実的になります。
動的計画法 の初期の応用例の一つに、問題を時間で解くヘルド・カープアルゴリズム がある。O ( n 2 2 n ) {\displaystyle O(n^{2}2^{n})} [ 24 ]
総当たり探索による7都市の対称TSPの解法。注:順列の数:(7 − 1)!/2 = 360 これらの時間制限を改善するのは難しいようです。たとえば、TSP の古典的な厳密アルゴリズム が時間で実行できるかどうかはまだ決定されていません。O ( 1.9999 n ) {\displaystyle O(1.9999^{n})} 存在する。 AmbainisらによるTSPに対する現在の最良の量子厳密アルゴリズムは時間で実行される。 O ( 1.728 n ) {\displaystyle O(1.728^{n})} [ 26 ]
その他のアプローチとしては、以下のようなものがある。
7都市の巡回セールスマン問題(TSP)を単純な分岐限定法で解く。注:順列の数は総当たり探索よりもはるかに少ない。 TSPLIB の 15,112 のドイツの町に対する正確な解は、1954 年にGeorge Dantzig 、Ray Fulkerson 、およびSelmer M. Johnson によって提案された線形計画法に基づく 切断平面法を使用して 2001 年に発見されました。計算は、 ライス大学 とプリンストン大学 にある 110 個のプロセッサのネットワークで実行されました。総計算時間は、単一の 500 MHz Alpha プロセッサ で 22.6 年に相当するものでした。2004 年 5 月には、スウェーデンのすべての 24,978 の町を訪問する巡回セールスマン問題が解決されました。約 72,500 キロメートルの長さのツアーが発見され、より短いツアーは存在しないことが証明されました。[ 29 ] 2005 年 3 月、回路基板上の 33,810 点すべてを訪問する巡回セールスマン問題がConcorde TSP Solver を 使用して解決されました。長さ 66,048,945 ユニットの巡回経路が見つかり、より短い巡回経路は存在しないことが証明されました。計算には約 15.7 CPU 年かかりました(Cook et al. 2006)。2006 年 4 月には、85,900 点のインスタンスがConcorde TSP Solver を使用して解決され、136 CPU 年以上かかりました。Applegate et al. (2006) を参照してください。
ヒューリスティックアルゴリズムと近似アルゴリズム さまざまなヒューリスティック と近似アルゴリズムが 考案されており、これらは迅速に良好な解をもたらします。これらには、マルチフラグメントアルゴリズム が含まれます。最新の手法では、非常に大規模な問題(数百万の都市)に対して、妥当な時間内に、最適解からわずか2~3%しか離れていない解を高い確率で見つけることができます。[ 14 ]
ヒューリスティックにはいくつかのカテゴリーが認められている。
構成的ヒューリスティクス 7つの都市からなるTSPに対する最近傍アルゴリズム。開始地点を変更すると解も変化する。 最近傍法(NN)アルゴリズム (貪欲法の 一種)では、セールスマンは次に訪れるべき最も近い未訪問都市を選択します。このアルゴリズムは、効果的に短いルートを迅速に生成します。平面上にランダムに分布するN個の 都市の場合、このアルゴリズムは平均して最短経路よりも25%長い経路を生成します。[ 30 ] ただし、NNアルゴリズムが最悪のルートを生成するような、特別な都市分布が多数存在します。[ 31 ] これは、非対称TSPと対称TSPの両方に当てはまります。[ 32 ] Rosenkrantzら[ 33 ] は、NNアルゴリズムの近似係数がΘ ( ログ | V | ) {\displaystyle \Theta (\log |V|)} 三角不等式を満たすインスタンスの場合。NN アルゴリズムの変形である最近傍フラグメント (NF) 演算子は、最も近い未訪問都市のグループ (フラグメント) を接続し、連続する反復でより短いルートを見つけることができます。[ 34 ] NF 演算子は、より良い解のみが受け入れられるエリート主義モデルでさらに改善するために、NN アルゴリズムによって得られた初期解にも適用できます。
点の集合のビトニックツアー とは、それらの点を頂点とする最小周長の単調多角形のことであり、 動的計画法 を用いて効率的に計算することができる。
もう1つの構成的ヒューリスティックで あるMatch Twice and Stitch (MTS)は、2つの連続したマッチング を実行します。2番目のマッチングは、最初のマッチングのすべてのエッジを削除した後に実行され、サイクルのセットを生成します。次に、サイクルをステッチして最終的なツアーを生成します。[ 35 ]
クリストフィデスとセルジュコフのアルゴリズム マッチングを作成する 上記のマッチングによって作成されたグラフに対してショートカットヒューリスティックを使用する ChristofidesとSerdyukovのアルゴリズムは 同様の概要に従いますが、最小全域木と別の問題である最小重み完全マッチング の解法を組み合わせています。これにより、最適解の最大1.5倍のTSPツアーが得られます。これは最初の近似アルゴリズムの1つであり、 扱いにくい問題 に対する実用的なアプローチとして近似アルゴリズムに注目を集めるきっかけの一つとなりました。実際、「アルゴリズム」という用語は、後になるまで近似アルゴリズムに一般的に拡張されることはありませんでした。Christofidesアルゴリズムは当初、Christofidesヒューリスティックと呼ばれていました。[ 10 ]
このアルゴリズムは、グラフ理論の結果を利用することで物事を異なる視点から捉え、最小全域木のコストを2倍にすることで得られたTSPの下限を改善します。オイラーグラフが与えられた場合、 オイラーツアー を見つけることができます。 O ( n ) {\displaystyle O(n)} 時間 、 [ 6 ] したがって、TSP の都市を頂点とするオイラーグラフがあれば、オイラー巡回路を見つけるためのこのような方法を使用して TSP の解を見つけることができることが容易にわかります。三角不等式 により、TSP 巡回路はオイラー巡回路より長くはならないことがわかっているので、TSP の下限が得られます。このような方法は以下に説明します。
この問題に対する最小全域木を求めなさい。 オイラーグラフを作成するために、すべての辺の複製を作成します。 このグラフのオイラー路を求めなさい。 TSPに変換する:ある都市が2回訪問された場合、ツアー内でその都市の前の都市から次の都市へのショートカットを作成する。 下限値を改善するには、オイラーグラフを作成するより良い方法が必要です。三角不等式により、最適なオイラーグラフは最適な巡回セールスマンツアーと同じコストを持つ必要があります。したがって、最適なオイラーグラフを見つけることは、TSP と同程度に困難です。これを行う 1 つの方法は、複雑度が のアルゴリズムを使用して最小重みマッチングを行うことです。 O ( n 3 ) {\displaystyle O(n^{3})} [ 6 ]
グラフをオイラーグラフにするには、まず最小全域木から始めます。次に、奇数次の頂点をすべて偶数にする必要があります。そのため、奇数次の頂点のマッチングを追加し、すべての奇数次の頂点の次数を 1 増やします。[ 6 ] これにより、すべての頂点の次数が偶数になるグラフが得られ、これはオイラーグラフとなります。上記の方法を適用すると、Christofides と Serdyukov のアルゴリズムが得られます。
この問題に対する最小全域木を求めなさい。 奇数個の都市の集合を用いて、この問題のマッチング問題を作成してください。 このグラフのオイラー路を求めなさい。 ショートカットを使用してTSPに変換します。
ペアワイズ交換 2-opt反復の例 ペアワイズ交換法、または2-opt 法は、2つのエッジを繰り返し削除し、それらを2つの異なるエッジで置き換えて、エッジ削除によって生じた断片を再接続し、より短い新しい巡回路を形成する手法です。同様に、3-opt 法は3つのエッジを削除し、それらを再接続してより短い巡回路を形成します。これらはk -opt法の特殊なケースです。Lin -Kernighan という名称は、2-opt法の誤称としてよく耳にしますが、実際にはLin-Kernighan法はより一般的なk -opt法です。
ユークリッド空間のインスタンスの場合、2-optヒューリスティックは平均してクリストフィデスのアルゴリズムよりも約5%優れた解を与えます。貪欲アルゴリズム で得られた初期解から始めると、平均移動回数は再び大幅に減少し、 O ( n ) {\displaystyle O(n)} ただし 、ランダムスタートの場合、平均移動回数は O ( n ログ ( n ) ) {\displaystyle O(n\log(n))} 。これはサイズの増加としては小さいものの、小さな問題の場合、ランダムスタートの初期移動回数は、貪欲ヒューリスティックから開始した場合と比較して 10 倍になります。これは、このような 2-opt ヒューリスティックが、交差などのソリューションの「悪い」部分を利用するためです。このようなタイプのヒューリスティックは、車両経路問題の ヒューリスティック内で、経路ソリューションを再最適化するためによく使用されます。 [ 30 ]
k -optヒューリスティック、またはリン・カーニハンヒューリスティックリン・カーニハン・ヒューリスティックは、 V -opt(可変最適化)手法の特殊なケースです。その手順は以下のとおりです。
与えられたツアーから、互いに素なk 個の辺を削除する。 残りの断片を再構成してツアーを作成します。ただし、互いに素な部分ツアーは残さないでください(つまり、断片の端点同士を接続しないでください)。これにより、検討対象のTSPは、より単純な問題に簡略化されます。 各フラグメントのエンドポイントは、2k − 2個の 他の可能性に接続できます。つまり、利用可能なフラグメントのエンドポイントの総数2k個 のうち、検討対象のフラグメントの2つのエンドポイントは使用できません。このような制約付き2k都市 TSPは、総当たり法を用いて元のフラグメントの最小コストの組み合わせを見つけることで解決できます。k -opt法の中で最もよく知られているのは、1965年にベル研究所 のシェン・リンによって導入された3-opt法です。3 -opt法の特殊なケースとして、辺が互いに分離していない場合(2つの辺が隣接している場合)があります。実際には、削除する辺のうち2つが隣接しているこの特殊な部分集合に3-変更を限定することで、一般的な3-opt法のような組み合わせコストをかけずに、2-opt法よりも大幅に改善できることがよくあります。このいわゆる2.5-opt法は、得られるツアーの質とツアーの達成に必要な時間の両方において、2-opt法と3-opt法のほぼ中間に位置します。
V -optヒューリスティック可変オプト法は、k- オプト法と関連があり、その一般化です。k-オプト法は元の巡回路から固定数(k)のエッジを削除しますが、 可変オプト 法は削除するエッジセットのサイズを固定しません。代わりに、探索プロセスが続くにつれてセットが大きくなります。このファミリーで最もよく知られている方法は、Lin-Kernighan法(上記では2-オプトの誤称として言及)です。Shen Lin とBrian Kernighanは 1972年にこの方法を初めて発表し、20年近くにわたって巡回セールスマン問題を解くための最も信頼できるヒューリスティックでした。より高度な可変オプト法は、1980年代後半にベル研究所でDavid Johnsonとその研究チームによって開発されました。これらの方法(Lin-Kernighan-Johnson と呼ばれることもあります)は、Lin-Kernighan法をベースに、タブー探索 と進化的計算の アイデアを追加したものです。基本的な Lin–Kernighan 手法では、少なくとも 3-opt であることが保証された結果が得られます。Lin–Kernighan–Johnson 法では、Lin–Kernighan ツアーを計算し、次に、少なくとも 4 つのエッジを削除してツアーを別の方法で再接続する突然変異と呼ばれる操作によってツアーを摂動し、新しいツアーをV -opt します。この突然変異は、 Lin–Kernighan 法で特定された局所最小値 からツアーを移動させるのに十分な場合が多いです。V - opt 法は、この問題に対する最も強力なヒューリスティックとして広く認識されており、ハミルトン サイクル問題や他の非計量 TSP など、他のヒューリスティックでは対応できない特殊なケースにも対応できます。長年にわたり、Lin–Kernighan–Johnson 法は、最適解が既知であるすべての TSP の最適解を特定し、この手法が試された他のすべての TSP の既知の最良解を特定してきました。
アリコロニー最適化 人工知能 研究者のマルコ・ドリゴは1993年に、 ACS (アリコロニーシステム )と呼ばれるアリコロニーのシミュレーション を使用してTSPの「良い解」をヒューリスティックに生成する方法を説明した。 [ 37 ] これは、実際のアリが食料源と巣の間の最短経路を見つける行動をモデル化しており、各アリが他のアリによって残された道しるべフェロモン に従うことを好むことから生じる創発的な 行動である。
ACSは、マップ上の多数の仮想アリエージェントを送り出し、様々な経路を探索させます。各アリは、都市までの距離と、その都市への境界線上に堆積された仮想フェロモンの量を組み合わせたヒューリスティックに基づいて、次に訪れる都市を確率的に選択します。アリは探索を続け、通過する境界線ごとにフェロモンを堆積させ、すべてのアリがツアーを完了するまで続けます。この時点で、最短のツアーを完了したアリが、完了したツアー経路に沿って仮想フェロモンを堆積させます(グローバルトレイル更新 )。堆積されるフェロモンの量はツアーの長さに反比例します。つまり、ツアーが短いほど、堆積されるフェロモンの量が多くなります。
1) アリは可能な経路の中から1つを選び、その経路上にフェロモンの痕跡を残します。2) すべてのアリは異なる経路を移動し、溶液の質に比例したフェロモンの痕跡を残します。3) 最良の経路の各辺は、他の辺よりも強化されています。4) 蒸発によって、質の悪い溶液は消滅します。この地図はイヴ・オーブリーの作品です。。 7都市のTSPに対するアリコロニー最適化アルゴリズム:フェロモンマップ上の赤と太線は、より多くのフェロモンが存在することを示しています。
特別なケース
メトリック メトリックTSP( デルタTSP またはΔ-TSPとも呼ばれる)では、都市間の距離は三角不等式 を満たします。
TSPの非常に自然な制約は、都市間の距離が三角不等式を 満たすような距離 を形成することを要求することです。つまり、Aから B への直接接続は、中間地点C を経由する経路よりも決して遠くならないということです。
d A B ≤ d A C + d C B {\displaystyle d_{AB}\leq d_{AC}+d_{CB}} 。すると、辺は頂点の集合上に距離関数 を構築する。都市を平面上の点とみなすと、多くの自然な距離関数が 距離関数となるため、多くの自然なTSP(巡回セールスマン問題)の事例がこの制約を満たす。
以下に、さまざまな指標に対するメトリックTSPの例をいくつか示します。
ユークリッドTSP(下記参照)では、2つの都市間の距離は、対応する点間のユークリッド距離です。 直線型TSPでは、2つの都市間の距離は、それぞれのx 座標とy 座標の差の絶対値の合計で表されます。この距離は、マンハッタン距離 またはシティブロック距離と呼ばれることがよくあります。 最大距離法 では、2点間の距離は、それらのx 座標とy 座標の差の絶対値の最大値となる。最後の2つの指標は、例えば、プリント基板 に所定の穴をあける機械の経路設計などに用いられます。マンハッタン距離は、まず一方の座標を調整し、次に他方の座標を調整する機械に対応し、新しい点への移動時間は両方の移動時間の合計となります。最大距離は、両方の座標を同時に調整する機械に対応し、新しい点への移動時間は、2つの移動のうち遅い方となります。
TSPの定義では、都市を2回訪問することは許可されていませんが、多くのアプリケーションではこの制約は必要ありません。そのような場合、対称的で非メトリックなインスタンスをメトリックなインスタンスに縮小できます。これにより、元のグラフが都市間の距離が完全なグラフに置き換えられます。d A B {\displaystyle d_{AB}} は、元のグラフにおけるA とB間の 最短経路 長に置き換えられます。
ユークリッド ユークリッド平面 上の点の場合、巡回セールスマン問題の最適解は、すべての点を通る単純な多角形 、つまり点の多角形化を形成します。 [ 38 ] 交差のある非最適解は、局所最適化によって交差のないより短い解にすることができます。ユークリッド距離は 三角不等式に従うため、ユークリッドTSPはメトリックTSPの特殊なケースを形成します。しかし、入力点が整数座標を持つ場合でも、それらの距離は一般に平方根 の形をとり、ツアーの長さは根号の和 となるため、異なるツアーの長さを正確に比較するために必要な記号計算を 実行することは困難です。
一般的な TSP と同様に、正確なユークリッド TSP は NP 困難ですが、根号の和の問題が、決定バージョンが NP に属し、したがって NP 完全であることを証明する際の障害となっています。距離を整数に丸めた問題の離散化バージョンは NP 完全です。有理座標と実際のユークリッド距離では、ユークリッド TSP はPSPACE のサブクラスである計数階層に属することが知られています。任意の実座標では、可能な入力が非可算個あるため、ユークリッド TSP はそのようなクラスには属せません。これらの複雑さにもかかわらず、ユークリッド TSP は、一般的な距離の場合よりも近似がはるかに容易です。例えば、ユークリッドTSPのインスタンスに関連付けられたグラフの最小全域木はユークリッド最小全域木であり、 n 個の点(エッジの数よりかなり少ない)に対して期待されるO ( n log n )の時間で計算できます。これにより、上記の三角不等式を持つTSPの単純な2近似アルゴリズムがより高速に動作します。
一般に、任意のc > 0 に対して、d はユークリッド空間の次元数であるが、TSP の幾何学的インスタンスに対して、最適値の最大 (1 + 1/ c ) 倍の長さのツアーを見つける多項式時間アルゴリズムが存在する。
O ( n ( ログ n ) O ( c d ) d − 1 ) {\displaystyle O{\left(n(\log n)^{O(c{\sqrt {d}})^{d-1}}\right)}} 時間。これは多項式時間近似スキーム (PTAS)と呼ばれます。 Sanjeev Arora とJoseph SB Mitchellは、 ユークリッドTSPのPTASを同時に発見した功績により、2010年にゲーデル賞 を受賞しました。
実際には、より単純な、保証の弱いヒューリスティックが引き続き使用されている。
非対称 ほとんどの場合、TSPネットワーク内の2つのノード間の距離は両方向で同じです。AからBまでの距離とBからAまでの距離が等しくない場合を非対称 TSPと 呼び ます。非対称TSPの実用的な応用例としては、 道路レベルのルーティング(一方通行、スリップロード、高速道路などによって非対称になる)を用いた経路最適化が挙げられます。
スタッカークレーン問題は 、非対称TSPの特殊なケースと見なすことができます。この問題では、入力は距離空間内の点の順序付きペアで構成され、ツアーはこれらのペアを順番に訪問する必要があります。これらの点のペアは、非対称TSPのノードと見なすことができ、非対称な距離は、ペアの最初の点から2番目の点まで、そして2番目の点から次のペアの最初の点まで移動する際のコストの合計を反映しています。
対称への変換 非対称TSPグラフを解くのはやや複雑になる場合があります。以下は、ノードA 、B 、C間のすべての可能なパス重みを含む3×3行列です。1つのオプションは、サイズ N の非対称行列をサイズ2Nの対称 行列 に変換することです。[ 43 ]
サイズを2倍にするには、グラフ内の各ノードを複製して、2番目のゴーストノード を作成します。このゴーストノードは、非常に低い(場合によっては負の)重みを持つ「ゴースト」エッジで元のノードに接続され、ここでは − w と表記されます。(あるいは、ゴーストエッジの重みは0で、他のすべてのエッジに重み w が追加されます。)上に示した元の 3×3 行列は左下に表示され、元の行列の転置行列は右上に表示されます。行列のどちらのコピーも、対角線が − w で表される低コストのホップパスに置き換えられています。新しいグラフでは、元のノードを直接リンクするエッジはなく、ゴーストノードを直接リンクするエッジもありません。
ゴーストノードと対応する元のノードを結ぶ「ゴースト」エッジの重み − w は、すべてのゴーストエッジが新しいグラフ上の任意の最適な対称 TSP 解に属することを保証するために十分に小さくなければなりません ( w = 0 が常に十分に小さくなるわけではありません)。結果として、最適な対称ツアーでは、各元のノードがゴーストノードの隣に現れます (たとえば、可能なパスは A → A ′ → C → C ′ → B → B ′ → A です)。そして、元のノードとゴーストノードを再びマージすることで、元の非対称問題の (最適な) 解が得られます (この例では、A → C → B → A です)。
計算複雑性 この問題はNP 困難 であることが示されており(より正確には、複雑性クラス FP NPに対して完全です。 関数問題 を参照)、決定問題 バージョン (「コストと数値xが与えられたとき、 x より安い往復ルートがあるかどうかを判定する」) はNP 完全 です。ボトルネック巡回セールスマン問題 も NP 困難です。この問題は、都市がユークリッド距離 で平面上にある場合や、その他多くの制約のある場合でも NP 困難のままです。各都市を「一度だけ」訪問するという条件を削除しても NP 困難性は解消されません。平面の場合、各都市を一度だけ訪問する最適なツアーが存在するためです (そうでなければ、三角不等式 により、繰り返し訪問をスキップするショートカットはツアーの長さを増加させません)。
TSP、特にそのユークリッド 版は、認知心理学 の研究者の注目を集めています。人間は、ほぼ線形的な方法で、ほぼ最適な解を迅速に生成できることが観察されており、そのパフォーマンスは、10~20ノードのグラフでは1%効率が低く、120ノードのグラフでは11%効率が低い範囲に及びます。[ 64 ] [ 65 ] 人間がこの問題に対してほぼ最適な解を正確に生成する明らかな容易さから、研究者は、人間が1つ以上のヒューリスティックを使用しているという仮説を立てており、おそらく最も有力な2つの理論は、凸包仮説と交差回避ヒューリスティックです。[ 66 ] [ 67 ] [ 68 ] しかし、追加の証拠は、人間のパフォーマンスはかなり多様であり、個人差とグラフの幾何学がタスクのパフォーマンスに影響を与えることを示唆しています。[ 69 ] [ 70 ] [ 71 ] それにもかかわらず、結果は、TSP におけるコンピュータのパフォーマンスは、人間がこれらの問題に対して使用する方法を理解し、模倣することによって改善される可能性があることを示唆しており、[ 72 ] また、人間の思考のメカニズムに関する新たな洞察にもつながっています。[ 73 ] Journal of Problem Solving の最初の号は、TSP における人間のパフォーマンスというトピックに特化しており、[ 74 ] 2011 年のレビューでは、この主題に関する数十の論文がリストされています。[ 73 ]
2011年に行われた動物認知 に関する研究「ハトにバスの運転をさせよう」(児童書『ハトにバスの運転をさせてはいけない! 』にちなんで名付けられた)では、巡回セールスマン問題に関連付けて、実験室で複数の餌場の間を移動するハトの飛行パターンを研究することで、ハトの空間認知能力 が調べられました。最初の実験では、ハトを実験室の隅に置き、近くのエンドウ豆の入った餌場まで自由に飛ばせました。研究者たちは、ハトが次にどの餌場を選ぶかを決める際に、主に近さを利用していることを発見しました。2番目の実験では、ハトがすべての餌場を訪れる必要がある場合、あらゆる機会に最も近い餌場に飛ぶことは非常に非効率的になるように餌場を配置しました。2番目の実験の結果は、ハトは依然として近さに基づく解決策を好むものの、「近さに基づく効率的なルートと非効率的なルートの移動コストの差が大きくなると、ルートに沿って数歩先まで計画を立てることができる」ことを示しています。[ 75 ] これらの結果は、霊長類以外の動物で行われた他の実験と一致しており、霊長類以外の動物の中には複雑な移動経路を計画できるものもいることを証明している。これは、霊長類以外の動物が比較的高度な空間認知能力を持っている可能性を示唆している。
ベンチマーク TSPアルゴリズムのベンチマークには、TSPLIB [ 79 ] はTSPおよび関連問題のサンプルインスタンスのライブラリです。その多くは実際の都市のリストと実際のプリント回路 のレイアウトです。[ 80 ]
注記 ↑ ラベ、マルティーヌ。ラポルト、ギルバート。ロドリゲス・マルティン、インマクラダ。サラザール・ゴンサレス、フアン・ホセ(2004年5月)。 「リングスター問題: 多面体解析と正確なアルゴリズム」。ネットワーク 。43 (3): 177–189 .土井 : 10.1002/net.10114。ISSN 0028-3045。 ↑ 既に最適解から0.05%以内の精度で解決されているTSP世界一周問題を参照してください。 ↑ 「Der Handlungsreisende – wie er sein soll und was er zu tun hat, um Aufträge zu erhalten und eines glücklichen Erfolgs in seinen Geschäften gewiß zu sein – von einem alten Commis-Voyageur」コミッションを獲得し、彼のビジネスが幸せに成功することを確信してください – 古い委員会航海者 による) ↑ ハミルトンとカークマンの初期の研究についての議論は、ビッグス、ロイド、ウィルソン著『グラフ理論、1736~1936年』 (クラレンドン・プレス、1986年)に掲載されている。 ↑ Schrijver (2005) で引用および英語訳。ドイツ語原文:「Wir bezeichnen als Botenproblem (weil diese Frage in der Praxis von jedem Postboten, übrigens auch von vielen Reisenden zu lösen ist) die Aufgabe, für endlich viele Punkte, deren paarweise Abstände bekannt sind, den kürzesten die Punkte verbindenden」 Weg zu finden. Dieses は、Regeln の自然な状態を示し、Anzahl der Versuche unter die Anzahl der Permutationen der gegebenen Punkte herunterdrücken würden、sind nicht bekannt です。ツム・ネクストゲレゲネンPunkt, dann zu dem dieem nächstgelegenen Punkt gehen usw., liefert im allgemeinen nicht den kürzesten Weg.」 1 2 3 4 5 6 7 8 Lawler, EL (1985). The Travelling Salesman Problem: A Guided Tour of Combinatorial Optimization (Repr. with corrections. ed.). John Wiley & Sons. ISBN 978-0-471-90413-7 。↑ ロビンソン、ジュリア(1949年12月5日)。 ハミルトンゲーム(巡回セールスマン問題)について (PDF) (技術報告書)。カリフォルニア州サンタモニカ:ランド研究所。RM-303 。 2020年 5月2日 取得 – 国防技術情報センター経由。 ↑ メンガーとホイットニーの関係、およびTSPの研究の発展に関する詳細な説明は、 Schrijver (2005) に記載されています。 1 2 3 van Bevern, René; Slugina, Viktoriia A. (2020). "メトリック巡回セールスマン問題に対する 3/2 近似アルゴリズムに関する歴史的注記". Historia Mathematica . 53 : 118– 127. arXiv : 2004.02437 . doi : 10.1016/j.hm.2020.04.003 . ↑ Klarreich, Erica (2013年1月30日). 「コンピュータ科学者が悪名高い巡回セールスマン問題の新たな近道を発見」 . WIRED . 2015年 6月14日 閲覧 . ↑ Klarreich, Erica (2020年10月8日). 「コンピュータ科学者が巡回セールスマンの記録を破る」 . Quanta Magazine . 2020年 10月13日 閲覧 . ↑ Karlin, Anna R. ; Klein, Nathan; Gharan, Shayan Oveis (2021), "メトリックTSPのための(わずかに)改善された近似アルゴリズム", Khuller, Samir ; Williams, Virginia Vassilevska (eds.), STOC '21: 53rd Annual ACM SIGACT Symposium on Theory of Computing, Virtual Event, Italy, June 21-25, 2021 , pp. 32– 45, arXiv : 2007.01409 , doi : 10.1145/3406325.3451009 , ISBN 978-1-4503-8053-9 1 2 Rego, César; Gamboa, Dorabela; Glover, Fred; Osterman, Colin (2011), "Traveling salesman problem heuristics: leading methods, implementations and latest advances", European Journal of Operational Research , 211 (3): 427– 441, doi : 10.1016/j.ejor.2010.09.010 , MR 2774420 。↑ マクギンティ、ジョー・クレイブン(2017年8月12日~13日)。 「スクールバスのルートをどう修正するか?MITに電話しよう」 (PDF) 。 ウォール・ストリート・ジャーナル 。p. A2。2018年4月12日に オリジナル (PDF) からアーカイブ済み 。 ↑ Behzad, Arash; Modarres, Mohammad (2002)、「一般化巡回セールスマン問題を巡回セールスマン問題に効率的に変換する新しい方法」、 第15回国際システム工学会議議事録(ラスベガス) ↑ Papadimitriou, CH; Steiglitz, K. (1998), Combinatorial optimization: algorithms and complexity , Mineola, NY: Dover 、308-309ページ。↑ タッカー、AW(1960)、「有向グラフと整数計画法について」、IBM数学研究プロジェクト(プリンストン大学) ↑ ダンツィヒ、ジョージ・B. (1963)、線形計画法とその拡張 、プリンストン、ニュージャージー州:プリンストン大学出版局、545–7頁、 ISBN 0-691-08000-3 第6刷、1974年。 ↑ Velednitsky, Mark (2017). "非対称巡回セールスマン問題における DFJ 多面体が MTZ 多面体に含まれることを示す簡潔な組み合わせ論的証明". Operations Research Letters . 45 (4): 323–324 . arXiv : 1805.06997 . doi : 10.1016/j.orl.2017.04.010 . ↑ Bektaş, Tolga; Gouveia, Luis (2014). "Miller–Tucker–Zemlin 部分巡回路除去制約へのレクイエム?". European Journal of Operational Research . 236 (3): 820– 832. doi : 10.1016/j.ejor.2013.07.038 . ↑ CE Miller、AW Tucker、RA Zemlin。1960年。「巡回セールスマン問題の整数計画法による定式化」。J . ACM 7、4(1960年10月)、326-329。DOI: https://doi.org/10.1145/321043.321046 ↑ Dantzig, G.; Fulkerson, R.; Johnson, S. (1954年11月). 「大規模巡回セールスマン問題の解法」. Journal of the Operations Research Society of America . 2 (4): 393– 410. doi : 10.1287/opre.2.4.393 . ↑ ベルマン(1960) 、ベルマン(1962) 、ヘルド& カープ(1962) ↑ アンバイニス、アンドリス。バロディス、カスパール。イレイド、ヤーニス。コカイニス、マーティンズ。プルーシス、クリシュヤーニス。ヴィフロフス、イェヴゲニス(2019)。 「指数時間動的計画法アルゴリズムの量子高速化」 。 離散アルゴリズムに関する第 30 回 ACM-SIAM シンポジウムの議事録 。 pp. 1783–1793 。 土井 : 10.1137/1.9781611975482.107 。 ISBN 978-1-61197-548-2 。↑ 巡回セールスマン問題 - 分岐限定法 ( YouTube) 。 ハンガリー行列アルゴリズム のように、行と列を減らして実を結ばない枝を切る方法。↑ アップルゲート、デヴィッド;ビクスビー、ロバート。フヴァータル、ヴァシェク。クック、ウィリアム。ケルド州ヘルスガウン(2004 年 6 月)。 「スウェーデンの最適なツアー」 。 2020 年 11 月 11 日 に取得 。 1 2 Johnson, DS ; McGeoch, LA (1997). "巡回セールスマン問題: 局所最適化のケーススタディ" (PDF) . In Aarts, EHL; Lenstra, JK (eds.). Local Search in Combinatorial Optimisation . London: John Wiley and Sons Ltd. pp. 215– 310. ↑ Gutina, Gregory; Yeob, Anders; Zverovich, Alexey (2002年3月15日). "Traveling salesman should not be greedy: domination analysis of greedy-type heuristics for the TSP" . Discrete Applied Mathematics . 117 ( 1–3 ): 81–86 . doi : 10.1016/S0166-218X(01)00195-0 . >↑ Zverovitch, Alexei; Zhang, Weixiong; Yeo, Anders; McGeoch, Lyle A.; Gutin, Gregory; Johnson, David S. (2007), "Experimental Analysis of Heuristics for the ATSP", The Traveling Salesman Problem and Its Variations , Combinatorial Optimization, Springer, Boston, MA, pp. 445– 487, CiteSeerX 10.1.1.24.2386 , doi : 10.1007/0-306-48213-4_10 , ISBN 978-0-387-44459-8 ↑ Rosenkrantz, DJ; Stearns, RE; Lewis, PM (1974年10月14日~16日). 巡回セールスマン問題に対する近似アルゴリズム . 第15回スイッチングおよびオートマタ理論に関する年次シンポジウム (SWAT 1974). doi : 10.1109/SWAT.1974.4 . ↑ Ray, SS; Bandyopadhyay, S.; Pal, SK (2007). "TSPとマイクロアレイ遺伝子順序付けにおける組み合わせ最適化のための遺伝的演算子". Applied Intelligence . 26 (3): 183– 195. CiteSeerX 10.1.1.151.132 . doi : 10.1007/s10489-006-0018-y . ↑ Kahng, AB; Reda, S. (2004). "Match Twice and Stitch: A New TSP Tour Construction Heuristic". Operations Research Letters . 32 (6): 499– 509. doi : 10.1016/j.orl.2004.04.001 . ↑ Alatartsev, Sergey; Augustine, Marcus; Ortmeier, Frank (2013年6月2日). 「近傍を持つ巡回セールスマン問題に対する挿入ヒューリスティックの制約」 (PDF) . 自動計画およびスケジューリングに関する国際会議議事録 . 23 : 2– 10. doi : 10.1609/icaps.v23i1.13539 . ↑ Dorigo, Marco; Gambardella, Luca Maria (1997). "巡回セールスマン問題に対するアリのコロニー". Biosystems . 43 (2): 73– 81. Bibcode : 1997BiSys..43...73D . CiteSeerX 10.1.1.54.7734 . doi : 10.1016/S0303-2647(97)01708-5 . PMID 9231906 . ↑ Quintas, LV; Supnick, Fred (1965). "最短ハミルトン回路のいくつかの性質について". The American Mathematical Monthly . 72 (9): 977– 980. doi : 10.2307/2313333 . JSTOR 2313333 . MR 0188872 . ↑ Jonker, Roy; Volgenant, Ton (1983). "非対称な巡回セールスマン問題を対称な巡回セールスマン問題に変換する". Operations Research Letters . 2 ( 161–163 ): 1983. doi : 10.1016/0167-6377(83)90048-2 . ↑ Arlotto, Alessandro; Steele, J. Michael (2016), "Beardwood–Halton–Hammersley theorem for stationary ergodic sequences: a counterexample", The Annals of Applied Probability , 26 (4): 2141– 2168, arXiv : 1307.0221 , doi : 10.1214/15-AAP1142 ↑ Few, L. (1955). "n 点を通る最短経路と最短道路". Mathematika . 2 (2): 141– 144. doi : 10.1112/s0025579300000784 . ↑ Fiechter, C.-N. (1994). "大規模巡回セールスマン問題に対する並列タブー探索アルゴリズム" . Disc. Applied Math . 51 (3): 243– 267. doi : 10.1016/0166-218X(92)00033-I . ↑ Held, M.; Karp, RM (1970). "巡回セールスマン問題と最小全域木". Operations Research . 18 (6): 1138– 1162. Bibcode : 1970OpRes..18.1138H . doi : 10.1287/opre.18.6.1138 . ↑ Goemans, Michel X. ; Bertsimas, Dimitris J. (1991). "ユークリッド巡回セールスマン問題に対するHeldとKarpの下限の確率的解析". Mathematics of Operations Research . 16 (1): 72– 89. doi : 10.1287/moor.16.1.72 . ↑ Johnson, DS; McGeoch, LA; Rothberg, EE (1996). "Held-Karp巡回セールスマン限界の漸近実験解析" (PDF) . In Tardos, Éva (ed.). Proceedings of the 7th Annual ACM-SIAM Symposium on Discrete Algorithms . Philadelphia: Society for Industrial and Applied Mathematics. pp. 341–350 . ISBN 978-0-89871-366-4 2013年6月16日にオリジナル(PDF) からアーカイブされました。↑ Christine L. Valenzuela と Antonia J. Jones による記事(2007年10月25日、 Wayback Machine に アーカイブ済み) ↑ Orponen, P.; Mannila, H. (1987). 近似保存還元について: 完全問題とロバスト尺度' (レポート). ヘルシンキ大学コンピュータサイエンス学部. 技術レポート C-1987–28. ↑ Serdyukov, Anatoliy I. (1978)、 「О некоторых экстремальных обходах в графах」 [ グラフにおけるいくつかの極端な散歩について ] ( PDF) 、 Upravlyaemye Sistemy (ロシア語)、 17 : 76–79 ↑ Svensson, Ola; Tarnawski, Jakub; Végh, László A. (2018). "非対称巡回セールスマン問題に対する定数係数近似アルゴリズム" . 第50回ACM SIGACT理論計算機科学シンポジウム論文集 . Stoc 2018. ロサンゼルス: ACM Press. pp. 204–213 . doi : 10.1145/3188745.3188824 . ISBN 978-1-4503-5559-9 。↑ Traub, Vera ; Vygen, Jens (2020年6月8日). 「ATSPのための改良された近似アルゴリズム」 . 第52回ACM SIGACT理論計算機科学シンポジウム議事録 . Stoc 2020. シカゴ、イリノイ州: ACM. pp. 1–13 . arXiv : 1912.00670 . doi : 10.1145/3357713.3384233 . ISBN 978-1-4503-6979-4 。↑ Traub, Vera; Vygen, Jens (2024). 巡回セールスマン問題の近似アルゴリズム . Cambridge University Press. ISBN 9781009445436 。↑ Macgregor, JN; Ormerod, T. (1996年6月)、「巡回セールスマン問題における人間のパフォーマンス」、 Perception & Psychophysics 、 58 (4): 527–539 、 doi : 10.3758/BF03213088 、 PMID 8934685 。↑ Dry, Matthew; Lee, Michael D.; Vickers, Douglas; Hughes, Peter (2006). "ノード数が異なる視覚的に提示された巡回セールスマン問題における人間のパフォーマンス". The Journal of Problem Solving . 1 (1). CiteSeerX 10.1.1.360.9763 . doi : 10.7771/1932-6246.1004 . ↑ Rooij, Iris Van; Stege, Ulrike; Schactman, Alissa (2003年3月1日). "ユークリッド巡回セールスマン問題における凸包と巡回交差:人間のパフォーマンス研究への示唆". Memory & Cognition . 31 (2): 215–220 . CiteSeerX 10.1.1.12.6117 . doi : 10.3758/bf03194380 . PMID 12749463 . ↑ MacGregor, James N.; Chu, Yun (2011). "巡回セールスマン問題および関連問題における人間のパフォーマンス:レビュー" . The Journal of Problem Solving . 3 (2). doi : 10.7771/1932-6246.1090 . ↑ MacGregor, James N.; Chronicle, Edward P.; Ormerod, Thomas C. (2004年3月1日). "凸包か交差回避か?巡回セールスマン問題における解法ヒューリスティクス" . Memory & Cognition . 32 (2): 260– 270. doi : 10.3758/bf03196857 . PMID 15190718 . ↑ Vickers, Douglas; Mayo, Therese; Heitmann, Megan; Lee, Michael D; Hughes, Peter (2004). "3種類の視覚的に提示された最適化問題における知能とパフォーマンスの個人差". Personality and Individual Differences . 36 (5): 1059– 1071. doi : 10.1016/s0191-8869(03)00200-9 . ↑ Kyritsis, Markos; Gulliver, Stephen R.; Feredoes, Eva (2017年6月12日). "ユークリッド巡回セールスマン問題を解く際の交差回避ヒューリスティック違反の認識". Psychological Research . 82 (5): 997–1009 . doi : 10.1007 /s00426-017-0881-7 . PMID 28608230 . ↑ Kyritsis, Markos; Blathras, George; Gulliver, Stephen; Varela, Vasiliki-Alexia (2017年1月11日) 「ユークリッド巡回セールスマン問題におけるパフォーマンスの予測 因子 としての方向感覚と誠実さ」 Heliyon . 3 (11) e00461. Bibcode : 2017Heliy...300461K . doi : 10.1016/j.heliyon.2017.e00461 . PMC 5727545. PMID 29264418 . ↑ Kyritsis, Markos; Gulliver, Stephen R.; Feredoes, Eva; Din, Shahab Ud (2018 年 12 月). "ユークリッド巡回セールスマン問題における人間の行動: ヒューリスティクスと図形効果の計算モデル化". Cognitive Systems Research . 52 : 387– 399. doi : 10.1016/j.cogsys.2018.07.027 . 1 2 MacGregor, James N.; Chu, Yun (2011), "巡回セールスマン問題および関連問題における人間のパフォーマンス:レビュー" , Journal of Problem Solving , 3 (2), doi : 10.7771/1932-6246.1090 。↑ Journal of Problem Solving 1(1)、2006年、2014年6月6日取得。↑ Gibson, Brett; Wilkinson, Matthew; Kelly, Debbie (2012年5月1日). 「ハトにバスを運転させよう:ハトは部屋の中で将来のルートを計画できる」。Animal Cognition . 15 (3): 379–391 . doi : 10.1007/s10071-011-0463-9 . PMID 21965161 . ↑ Jones, Jeff; Adamatzky, Andrew (2014), "縮小するブロブによる巡回セールスマン問題の計算" (PDF) , Natural Computing : 2, 13, arXiv : 1303.4969 , 2017年6月4日に オリジナル (PDF)からアーカイブ済み、 2016年 1月26日 取得 ↑ Morell, Virginia (2012年9月21日). 「空飛ぶ数学:ミツバチが巡回セールスマン問題を解決する」 . Wired . ISSN 1059-1028 . 2025年 11月30日 取得 。 ↑ 「マルハナバチは移動セールスマン問題をその場で解決する」 。 ニュー・サイエンティスト 。2017年12月11日。 2025年 11月30日 閲覧 。 ↑ "TSPLIB" . GitHub . 2025年 12月28日 取得 . ↑ Reinelt, Gerhard (1991年11月)「TSPLIB – 巡回セールスマン問題ライブラリ」 ORSA Journal on Computing 3 (4) . Institute for Operations Research and the Management Sciences (INFORMS): 376–384 . doi : 10.1287/ijoc.3.4.376 . ↑ ギア、ダンカン (2012 年 4 月 26 日)。 」 「『旅するセールスマン』の映画は、PがNPに等しい場合の影響を考察している」。Wired UK 。 2012年 4月26日 取得 。↑ モナ・リザがNP困難であるとき エヴリン・ラム著、サイエンティフィック・アメリカン、2015年4月31日
参考文献 Applegate, DL; Bixby, RM; Chvátal, V.; Cook, WJ (2006) 『巡回セールスマン問題 』プリンストン大学出版局、ISBN 978-0-691-12993-8 。アレンダー、エリック。ピーター・ビュルギッサー;ケルトゴー・ペダーセン、ヨハン。 Mitersen、Peter Bro (2007)、「数値解析の複雑さについて」(PDF) 、SIAM J. Comput. 、38 (5): 1987–2006 、CiteSeerX 10.1.1.167.5495 、doi : 10.1137/070697926 。Arora, Sanjeev (1998)、「ユークリッド巡回セールスマン問題およびその他の幾何学的問題に対する多項式時間近似スキーム」(PDF) 、Journal of the ACM 、45 (5): 753–782 、doi : 10.1145/290179.290180、MR 1668147 。Beardwood, J.; Halton, JH; Hammersley, JM (1959年10月)、「多数の点を通る最短経路」、ケンブリッジ哲学協会紀要 、55 (4): 299–327 、Bibcode : 1959PCPS...55..299B、doi : 10.1017/s0305004100034095 。Bellman, R. ( 1960)、「組合せ過程と動的計画法」、Bellman, R.、Hall, M. Jr. (編)、『組合せ解析、応用数学シンポジウム論文集 10』 、アメリカ数学会、pp. 217–249 。Bellman, R. (1962)、「巡回セールスマン問題の動的計画法による処理」、Journal of the Association for Computing Machinery 、9 :61–63 、doi :10.1145/321105.321111 。Berman, Piotr; Karpinski, Marek (2006)、「(1,2)-TSP に対する 8/7 近似アルゴリズム」、第 17 回 ACM-SIAM 離散アルゴリズムシンポジウム (SODA '06) 論文集 、pp. 641–648 、CiteSeerX 10.1.1.430.2224 、doi : 10.1145/1109557.1109627、ISBN 978-0-89871-605-4 ECCC TR05-069 。Christofides, N. (1976)、「巡回セールスマン問題に対する新しいヒューリスティックの最悪ケース分析」 、技術報告書388、カーネギーメロン大学産業経営大学院、ピッツバーグ 。Hassin, R.; Rubinstein, S. (2000), "Better approximations for max TSP", Information Processing Letters , 75 (4): 181–186 , CiteSeerX 10.1.1.35.7209 , doi : 10.1016/S0020-0190(00)00097-1 。Held, M. ; Karp, RM (1962)、「シーケンス問題への動的計画法アプローチ」、Journal of the Society for Industrial and Applied Mathematics 、10 (1): 196–210 、doi : 10.1137/0110015 。Kaplan, H.; Lewenstein, L.; Shafrir, N.; Sviridenko, M. (2004)、「有向正則多重グラフの分解による非対称TSPの近似アルゴリズム」、第44回IEEEコンピュータサイエンス基礎シンポジウム論文集 、pp. 56–65 。Karpinski, M.; Lampis, M.; Schmied, R. (2015)、「TSP の新しい近似不可能性の限界」、Journal of Computer and System Sciences 、81 (8): 1665–1677 、arXiv : 1303.6437 、doi : 10.1016/j.jcss.2015.06.003 Kosaraju, SR; Park, JK; Stein, C. (1994)、「ロングツアーとショートスーパーストリングス」 「、第35回IEEEコンピュータサイエンス基礎シンポジウム論文集 、IEEEコンピュータソサエティ、pp. 166–177 」 。Larson, Richard C.; Odoni, Amedeo R. (1981)、「6.4.7: ネットワークモデルの応用 § ルーティング問題 § § ユークリッドTSP」、都市オペレーションズリサーチ 、Prentice-Hall、ISBN 978-0-13-939447-8 OCLC 6331426 。Padberg, M.; Rinaldi, G. (1991), "大規模対称巡回セールスマン問題の解決のための分岐限定法アルゴリズム", SIAM Review , 33 (1): 60–100 , Bibcode : 1991SIAMR..33...60P , doi : 10.1137/1033004 。Papadimitriou, Christos H. (1977)、「ユークリッド巡回セールスマン問題はNP完全である」、Theoretical Computer Science 、4 (3): 237–244 、Bibcode : 1977TComS...4..237P、doi : 10.1016/0304-3975(77)90012-3 、MR 0455550 。Papadimitriou, Christos H.; Yannakakis, Mihalis (1993)、「距離が1と2の巡回セールスマン問題」、Mathematics of Operations Research 、18 : 1–11 、doi : 10.1287/moor.18.1.1 。Schrijver, Alexander (2005). 「組合せ最適化の歴史について(1960年まで)」。K. Aardal 、GL Nemhauser 、R. Weismantel編『離散最適化ハンドブック (PDF)』 所収。アムステルダム:Elsevier。pp. 1–68 。 Serdyukov, AI (1984)、「最大値の巡回セールスマン問題の推定値を持つアルゴリズム」 」、ウプラヴリャエミエ・システミー 、25 : 80–86 。Steinerberger, Stefan (2015)、「巡回セールスマン定数の新しい境界」、Advances in Applied Probability 、47 (1): 27–36 、arXiv : 1311.6338 、doi : 10.1239/aap/1427814579 。Woeginger, GJ (2003)、「NP困難問題に対する厳密なアルゴリズム:概説」、Combinatorial Optimization – Eureka, You Shrink! Lecture notes in computer science、vol. 2570 、Springer、pp . 185–207 。
さらに読む Adleman, Leonard (1994)、「組み合わせ問題の解の分子計算」(PDF) 、Science 、266 (5187):1021–4 、Bibcode :1994Sci...266.1021A、CiteSeerX 10.1.1.54.2565 、doi :10.1126/science.7973651、PMID 7973651 、 2005年2月6日にオリジナル(PDF) からアーカイブ済み Babin, Gilbert; Deneault, Stéphanie; Laportey, Gilbert (2005)、「対称巡回セールスマン問題に対するOr-optヒューリスティックの改善」、オペレーションズ・リサーチ学会誌 、Cahiers du GERAD、G-2005-02 (3)、モントリオール:意思決定分析研究グループ:402–407 、CiteSeerX 10.1.1.89.9953 、JSTOR 4622707 クック、ウィリアム (2012)。『巡回セールスマンを追い求めて:計算の限界における数学 』プリンストン大学出版局。ISBN 978-0-691-15270-7 。Cook, William ; Espinoza, Daniel; Goycoolea, Marcos (2007)、「TSP に対するドミノパリティ不等式を用いた計算」、INFORMS Journal on Computing 、19 (3): 356–365 、doi : 10.1287/ijoc.1060.0204Cormen, Thomas H. ; Leiserson, Charles E. ; Rivest, Ronald L. ; Stein, Clifford (2009年7月31日). 「35.2: 巡回セールスマン問題」 .アルゴリズム入門 (第2 版). MIT Press. pp. 1027–1033 . ISBN 978-0-262-03384-8 。Dantzig, GB ; Fulkerson, R. ; Johnson, SM (1954)、「大規模巡回セールスマン問題の解法」、Operations Research 、2 (4): 393–410 、doi : 10.1287/opre.2.4.393、JSTOR 166695、S2CID 311786 Garey, Michael R.; Johnson, David S. (1979). "A2.3: ND22–24". Computers and Intractability: A Guide to the Theory of NP-completeness . WH Freeman. pp. 211–212 . ISBN 978-0-7167-1044-8 。 Goldberg, DE (1989)、Genetic Algorithms in Search, Optimization & Machine Learning 、Reading: Addison-Wesley、Bibcode : 1989gaso.book.....G、ISBN 978-0-201-15767-3 Gutin, G.; Yeo, A.; Zverovich, A. (2002年3月15日). "巡回セールスマン問題は貪欲であってはならない: TSPに対する貪欲型ヒューリスティックの支配分析" . Discrete Applied Mathematics . 117 ( 1–3 ): 81–86 . doi : 10.1016/S0166-218X(01)00195-0 . Gutin, G.; Punnen, AP (2007年5月18日).巡回セールスマン問題とその変種 . Springer US. ISBN 978-0-387-44459-8 。 Johnson, DS ; McGeoch, LA (1997)、「巡回セールスマン問題:局所最適化の事例研究」、Aarts, EHL; Lenstra, JK (編)、『組合せ最適化における局所探索』 (PDF) 、John Wiley and Sons Ltd.、pp. 215–310 Lawler, EL; Shmoys, DB; Kan, AHG Rinnooy; Lenstra, JK (1985). 『巡回セールスマン問題 』 John Wiley & Sons, Incorporated. ISBN 978-0-471-90413-7 。 MacGregor, JN; Ormerod, T. (1996)、「巡回セールスマン問題における人間のパフォーマンス」、Perception & Psychophysics 、58 (4): 527–539 、doi : 10.3758/BF03213088 、PMID 8934685 Medvedev, Andrei; Lee, Michael; Butavicius, Marcus; Vickers, Douglas (2001年2月1日). 「視覚的に提示された巡回セールスマン問題における人間のパフォーマンス」. Psychological Research . 65 (1): 34–45 . doi : 10.1007/s004260000031 . PMID 11505612 . Mitchell, JSB (1999)、「ギロチン分割は多角形分割を近似する:幾何学的TSP、 k -MST、および関連問題に対する単純な多項式時間近似スキーム」、 SIAM Journal on Computing 、28 (4): 1298–1309 、doi : 10.1137/S0097539796309764Rao, S.; Smith, W. (1998). 「スパナーとバニヤンによる幾何学的グラフの近似」 「. STOC '98: Proceedings of the thirtieth annual ACM symposium on Theory of computing . pp. 540– 550. CiteSeerX 10.1.1.51.8676 . Rosenkrantz, Daniel J.; Stearns, Richard E.; Lewis, Philip M. II (1977). 「巡回セールスマン問題に対するいくつかのヒューリスティックの分析」. SIAM Journal on Computing . 6 (5). SIAM (Society for Industrial and Applied Mathematics): 563–581 . doi : 10.1137/0206041 . ウォルショー、クリス(2000)『巡回セールスマン問題への多層的アプローチ 』CMS Press Walshaw, Chris (2001), 『巡回セールスマン問題に対する多段階Lin-Kernighan-Helsgaunアルゴリズム』 、CMS Press