JP / EN

PHCpack Benchmark Test 3

【MAL-Seeker 実験記録】主要80問題、その最初の一巡

Target Benchmark Suite:

PHCpack Main Database — Final 45 Polynomial Systems

Range: cyclic7 — rose
Result: 45 / 45 systems reached the official isolated-solution count
Milestone: One complete pass through the 80-system main collection

1. roseへ到達するまで

第一報ではboonからnoon5までの25問題、第二報ではredcyc7からcaprasseまでの10問題を記録しました。そして今回、caprasseの次に並ぶcyclic7から、PHCpack公式データベースの主要コレクション末尾に置かれたroseまで、残る45問題の検証が完了しました。

ここで重要なのは、同じ解法を45回繰り返したのではないという点です。序盤の多段再開探索、中盤の問題削減と非飽和化、終盤の総次数ホモトピー、多斉次ホモトピー、変数スケーリング。問題ごとに異なる構造を読み取り、必要な道具を一つずつMAL-Seekerへ組み込んできました。

45 / 45今回、公式解数へ到達
5,694表中の捕捉解数合計
約3時間57分採用ログの実行時間合計
80 problems主要コレクションを一巡

2. 巨大巡回系の突破 — cyclic7からkatsura5

最初の10問題では、924解を持つ7次巡回系と、1,152解を持つ8次巡回系が立ちはだかりました。チェックポイントへ解集合を保存し、探索領域を広げながら不足分だけを追うことで、長時間探索を「続きから始められる計算」へ変えています。

ProblemOfficial SolutionsSolutions FoundRun Time (sec)Method / Note
cyclic79249247,476.584*560→658→728→854→924。再開探索+巡回対称性。
extcyc79249242,917.936*868解から再開し、残る56解を補完。
ipp2161649.35011変数の6R逆位置問題を複素探索。
puma16160.073PUMAロボットの位置・姿勢問題。
redeco5880.098縮約された5次元経済モデル。
trinks10100.092PoSSoテスト系を全解捕捉。
cassou161690.610Cassou-Noguès系を複数境界で探索。
cyclic81,1521,152367.001巡回対称性展開を用いて全解捕捉。
extcyc81,1521,1522,907.289*1,120解から極端領域を再探索し32解を補完。
katsura5323212.189磁性モデルの32解を複素領域で捕捉。

* 再開に使用した連続した探索ステージの実行時間を合算しています。

3. 方程式の形を変えて勝つ — quadfor2からrabmo

次の15問題では、力任せの初期値探索よりも定式化が主役になりました。des18_3は縮約と総次数ホモトピー、fbrfive4は非飽和化、kinemaは多斉次構造の利用によって公式解数へ到達しています。

ProblemOfficial SolutionsSolutions FoundRun Time (sec)Method / Note
quadfor2220.0142代表解を捕捉。多重度込みでは4。
redeco616164.793縮約された6次元経済モデル。
virasoro2002001.241Virasoro代数構成系の200解。
chemequ161612.040炭化水素燃焼の化学平衡。
d1484844.51912変数の疎なD1ベンチマーク。
fbrfive12363681.8115点を通る4節リンク機構。
kin1484841.39912変数の機構学問題。
quadgrid5521.462格子上関数の補間求積公式。
redeco7323229.949縮約された7次元経済モデル。
wood991.345Wood関数の最適化から導かれる系。
chemequs16162.262スケール化された化学平衡系。
des18_346460.155縮約+総次数ホモトピー。
fbrfive436360.208非飽和化+多斉次ホモトピー。
kinema404045.1833群の多斉次構造を利用。
rabmo16160.680多次元最適求積公式。

4. スケーリングとホモトピー — redeco8からrose

最後の20問題では、総次数ホモトピーと多斉次ホモトピーが探索の中核へ入りました。特にcohn2、cohn3、discret3では、そのままでは少数の解しか見えなかった系をスケーリングし、経路追跡へ渡すことで一気に公式値へ到達しています。

ProblemOfficial SolutionsSolutions FoundRun Time (sec)Method / Note
redeco864640.289総次数ホモトピー。
wright32320.060総次数ホモトピー。
cohn218180.799正規化スケーリング+総次数ホモトピー。
des22_24424211.94210変数のdessin d'enfant系。
fourbar36360.957多斉次ホモトピーによる4節リンク設計。
ku10220.01310群の多斉次構造を利用。
rbpl404058.363並列ロボット左手問題。公式40解を捕捉。
rediff3770.005反応拡散系を総次数ホモトピーで追跡。
cohn31021023.801正規化スケーリング+多斉次ホモトピー。
discret31281285.262係数スケーリング+多斉次ホモトピー。
gaukwa2222.9792節点ガウス求積公式。
lorentz11110.0394次元Lorentzアトラクタの平衡点。
rbpl2440401.25524実解を含む並列ロボット系。
reimer51441445.848総次数ホモトピー+置換対称性。
comb300016161.3783000度の燃焼化学モデル。
eco5880.0155次元一般均衡問題。
gaukwa36621.232多斉次ホモトピー+置換対称性。
lumped440.011集中定数化学反応モデル。
rbpl24es40400.600スケール化された並列ロボット系。
rose1361360.408一般経済均衡モデル。主要コレクションの終着点。

5. 45問が残した三つの転換点

探索は、一回で終わらなくてよい

cyclic7は560、658、728、854、そして924解へ進みました。途中結果を失わず、次の探索へ渡す再開機構によって、時間のかかる全解探索を積み上げ可能な工程へ変えました。

式の構造は、計算量そのものである

fbrfive4の非飽和化、des18_3の縮約、kinemaの多斉次分割は、同じ数学的対象でも「どう見せるか」で難易度が変わることを示しました。全解探索はソルバーの性能だけでなく、定式化との共同作業です。

ニュートン探索から経路追跡へ

cohn2、cohn3、discret3を境に、総次数・多斉次ホモトピーが本格的に加わりました。局所的に解へ落ちるだけでなく、構造から生成した始点を連続的に終点まで運ぶ。MAL-Seekerは探索器から、多項式系を解き切るための複合エンジンへ一段進みました。

今回の到達点:

今回追加した45問題は、すべてPHCpack公式の孤立解数と一致しました。第一報・第二報を合わせ、boonからroseまでの主要80問題に対する最初の一巡が完了しました。ただし、初期の実験記録には未完の結果もそのまま残しています。これは「すべてを完全制覇した」という宣言ではなく、80種類の地形を実際に歩き、次に直すべき場所まで見えたという一区切りです。

6. 解析環境およびエンジン仕様

  • 解析エンジン: MAL-Seeker (Project Antares Ver 3.1.6)
  • 開発プロジェクト: Project Bellatrix V1.0.4
  • 主な探索機構: 5SS / Polynomial Newton / Resume / Symmetry Closure
  • 経路追跡: Total-Degree / M-Homogeneous Homotopy
  • PC: MacBook Pro (Mac17,2)
  • OS: macOS Tahoe 26.4 (Build 25E246)
  • CPU: Apple M5, 10-core (4P + 6E), arm64
  • メモリ: 32GB

7. 結論 — 完成ではなく、座標を得た

PHCpackの主要80問題は、数値解析、ロボット機構学、化学平衡、経済均衡、求積公式、代数構造まで、まるで異なる世界を一つの多項式系データベースへ押し込んだ地図でした。その端から端まで一度歩いたことで、MAL-Seekerが得意な地形と、まだ装備が不足する地形の輪郭が見えました。

roseへの到達はゴールではありません。公式ページには別枠の新規テスト例や、特異解のデフレーション、代数集合の分解という次の領域が続いています。けれど、ここで一度立ち止まり、主要コレクション一巡の記録を残します。迷宮の出口を見つけたのではなく、ようやく迷宮全体の座標を手に入れたのです。