第0章 技術解説書インデックスおよび理論概要¶

本インデックスは、BMS Part Tuner のコアロジックを構成する音響信号処理、離散数学、および時間軸シミュレーションの理論的背景について、大学レベルの応用数学・信号処理理論の枠組みに則り体系化した技術解説書の全体像を示すものである。

0.1 ターゲット読者および前提知識¶

本ドキュメント群は、単なるアプリケーションの使用方法の提示にとどまらず、「極限の処理性能がいかにして数理的・アルゴリズム的に導出されるか」に関心を持つアルゴリズムエンジニア、コンピュータサイエンスの研究者、およびハードウェア最適化のスペシャリストを対象とする。

0.2 本書で使用する数学記号（Notation）¶

本書の各章にわたり共通して使用される主要な数学記号とその定義を以下に示す。

• \(N\) : 音声波形のサンプル長（時間軸解像度）
• \(M\) : 入力となる音声ファイル数、あるいは処理対象のノーツイベント数
• \(L\) : システム全体における総ノーツ数（通常 \(L \ge M\)）
• \(W\) : 短時間エネルギー計算等に用いる移動窓（ウィンドウ）の幅
• \(k\) : 1小節内の細分化された音符数（解像度）
• \(\alpha(x)\) : アッカーマン関数の逆関数（極めてゆっくりと成長する関数）
• \(x(t), y(t)\) : 離散時間 \(t\) における音響信号（波形）の振幅
• \(\mathcal{O}(\cdot)\) : アルゴリズムの漸近的な計算量オーダー

0.3 最最適化がもたらすシステム全体のインパクトと次元の呪い¶

数万規模の音声ファイル群から「聴覚上のクオリティを一切損なわずに」重複を検出し統合する処理、および高解像度フォーマット（bmson）から表現能力に制限のあるBMSフォーマットへのダウンコンバート（逆写像）処理は、素朴に実装した場合、以下の致命的な技術的障壁に直面し破綻する。

• 次元の呪い: ペア間の総当り比較による指数関数的な時間計算量の増大。
• 浮動小数点演算の丸め誤差（桁落ち）: 極小振幅区間における精度崩壊。
• 頻繁なヒープアロケーション: 大量のオブジェクトアロケーションに伴うGCストール。

本システムでは、これらの課題を「数理的定理による計算空間の次元削減」と「DODやSIMD命令によるハードウェア駆動最適化」の双方を組み合わせることで克服し、数時間を要していた処理時間を数秒〜数十秒へと劇的に（\(99.9\%\)以上）短縮した。その理論的アプローチを各章にて詳述する。

0.3.1 各最適化フェーズにおける計算量削減効果¶

本システムにおいて適用されている数学的変形によるオーダー削減の全貌を以下に示す。

※表内の各変数の定義：

• \(N\) ：波形サンプル数
• \(M\) ：音声ファイル数・ノーツイベント数
• \(L\) ：総ノーツ数
• \(W\) ：窓幅
• \(k\) ：小節内音符数
• \(\alpha\) ：アッカーマン関数の逆関数
• \(\mathcal{C}\) ：BMS最大チャンネル数（1296）

0.3.2 パイプライン統合時におけるトータル計算量の削減¶

各フェーズにおける独立した最適化を統合したパイプラインにおいて、システム全体の時間および空間の複雑性は以下のように劇的に圧縮される。

1. 素朴な実装（総当り・時間領域シフト全探索・動的アロケーション彩色）：

• 破綻の理由: すべての音声ペア（\(M^2\)）に対して時間領域のズレ探索（\(N^2\)）を実行し、さらにプレイアビリティ検証でも総当り（\(L^2\)）を行うため、入力ファイル数やノーツ数が微増するだけで演算回数が指数関数的に跳ね上がり、処理が完全に膠着する。また、ダウンコンバート時にオブジェクト生成を乱発（\(\mathcal{O}(M)\)）することで、ヒープの肥大化とGCストールによる重大な遅延を引き起こす。

2. 本システム（LSH ＋ FFT ＋ Union-Find ＋ 走査線 ＋ DOD彩色）：

• 削減のメカニズム:
  • LSH および SimHash ハック: 音声の総当り比較空間を \(\mathcal{O}(M^2) \to \mathcal{O}(M)\) の線形オーダーへと確率的に圧縮する。ハミング距離は XOR および POPCNT 命令のわずか 2 クロックで算出される（第6章）。
  • FFT: 2波形間のアライメント積和演算を \(\mathcal{O}(N^2) \to \mathcal{O}(N \log N)\) へと削減する（第5章）。
  • Union-Find: 抽出された重複ペアから同値群をまとめる処理を、\(\mathcal{O}(M^2) \to \mathcal{O}(M \alpha(M))\) の実質的な線形時間で完了させる（第6章）。
  • Sweep-Line: プレイフィールを完全保存する区間衝突検証、および逆変換時の多重発音彩色パッキングを、それぞれ \(\mathcal{O}(L \log L)\) または \(\mathcal{O}(M \log M)\) に抑制する（第7章、第1章）。
  • DOD / Array-Backed Pool: チャンネル割当の状態管理を、動的アロケーションを廃した定数サイズのフラットな配列ポインタ操作に射影し、空間計算量を実質定数 \(\mathcal{O}(1)\) まで圧縮。GCストールを物理的に排除する（第1章）。

この直交するすべての数理の組み合わせ（ \(\mathcal{O}(M^2 \cdot N^2 + L^2) \Longrightarrow \mathcal{O}(M \cdot N \log N + M \alpha(M) + L \log L)\) ）こそが、数時間におよぶ力技のバッチ処理をわずか数十秒のバックグラウンドタスクへと昇華させる「異次元のパフォーマンス」の真の論理的根拠である。

0.4 処理アーキテクチャと数理理論のマッピング¶

巨大なBMSファイル群が入力されてから、最適化またはダウンコンバートされて出力されるまでのデータフローと、各フェーズを支える数学理論（各章）の対応関係を以下の図に示す。

flowchart TB
    %% 1. グローバル起点
    A["入力: 未最適化のBMS / bmson<br>と音声ファイル"]

    %% 2. 前処理フェーズ (bmsonの場合のみ)
    subgraph Phase4 ["前処理: 逆写像と格子・チャンネルパッキング"]
        direction TB
        I("第1章: <br>ユークリッド互除法<br>による最適格子量子化")
        J("第1章: <br>区間グラフ彩色 <br>＆<br>Array-Backed Pool <br>によるGCフリーパッキング")
        I --> J
    end

    %% 3. 並列実行・メモリプール基盤
    subgraph Platform ["共通実行基盤: 並列処理とゼロ・アロケーションメモリ空間"]
        direction TB
        K("第8章: <br>アムダールの法則による<br>タスク並列化 (Parallel.ForEach)")
        L("第8章: <br>ArrayPool <br>＆ Span による<br>LOH断片化防止と<br>空間 O(1) 圧縮")
        K <--> L

        %% 4. 最適化パイプラインの定義 (基盤上で実行される)
        subgraph Phase1 ["フェーズ1: 信号の正規化と特徴抽出"]
            direction TB
            B("第2章:<br> 無音トリミングと<br>短時間振幅包絡線")
            C("第3章:<br> LSH /<br> SimHashによる<br>特徴パッキング")
            B --> C
        end

        subgraph Phase2 ["フェーズ2: 並列比較と同値関係の抽出"]
            direction TB
            D("第4章: <br>SIMD超並列による<br>ピアソン相関係数")
            E("第5章: <br>FFTアライメントによる<br>位相補正")
            D --> E
        end

        subgraph Phase3 ["フェーズ3: グラフ縮約と論理検証"]
            direction TB
            F("第6章: <br>Union-Findによる<br>同値類グラフの構築")
            G("第7章: <br>Sweep-Lineによる<br>プレイアビリティ<br>等価性証明")
            F --> G
        end

        %% 内部ノードを直接繋ぐことでメインの縦軸（体幹）を形成
        L -->|"バッファ借用 & <br>並列実行"| B
        C -->|"ハッシュバケツ分割 & POPCNT比較"| D
        E -->|"高相関ペアの抽出"| F
    end

    %% 5. フローの結合と合流ノード
    MergeNode{"フォーマット共通化<br>(BMS統一)"}

    A -->|"bmsonの場合"| I
    A -->|"BMSの場合"| MergeNode
    J -->|"ダウンコンバート済BMS"| MergeNode

    MergeNode -->|"統合済ファイルリスト"| K

    H["出力: 最適化済BMS と<br>統合音声 (代表元)"]
    G -->|"プレイアビリティ等価性確認済"| H

    %% 6. 非類似・除外フロー (横へ逸れる分岐)
    Z1["出力: 統合除外 (WAV保持)"]
    Z2["出力: 統合除外 (WAV保持)"]
    Z3["出力: 統合除外 (WAV保持)"]
    Z4["出力: 統合除外 (WAV保持)"]

    C -.->|"ハミング距離超過"| Z1
    D -.->|"相関係数未満"| Z2
    E -.->|"アライメント後相関未満"| Z3
    G -.->|"同時発音干渉"| Z4

    style Z1 fill:#fafafa,stroke:#9e9e9e,stroke-dasharray: 5 5
    style Z2 fill:#fafafa,stroke:#9e9e9e,stroke-dasharray: 5 5
    style Z3 fill:#fafafa,stroke:#9e9e9e,stroke-dasharray: 5 5
    style Z4 fill:#fafafa,stroke:#9e9e9e,stroke-dasharray: 5 5

    %% 7. スタイリング
    style B fill:#e1f5fe,stroke:#01579b,stroke-width:2px
    style C fill:#e1f5fe,stroke:#01579b,stroke-width:2px
    style D fill:#e8f5e9,stroke:#1b5e20,stroke-width:2px
    style E fill:#e8f5e9,stroke:#1b5e20,stroke-width:2px
    style F fill:#fff3e0,stroke:#e65100,stroke-width:2px
    style G fill:#fff3e0,stroke:#e65100,stroke-width:2px
    style I fill:#f3e5f5,stroke:#4a148c,stroke-width:2px
    style J fill:#f3e5f5,stroke:#4a148c,stroke-width:2px
    style K fill:#eceff1,stroke:#37474f,stroke-width:2px
    style L fill:#eceff1,stroke:#37474f,stroke-width:2px

0.5 理論と実装の対応マップ¶

抽象的な数学理論が、本プロジェクトの実際のコードベースにおいてどこに実装されているかを示す。ドキュメントを読み進めながら実際のソースコードを参照することで、理論と実践のギャップを埋めることができる。

0.6 各章の構成と解説対象の数理¶

アーキテクチャの各コンポーネントを支える厳密な数学的証明と実装への射影について、詳細は以下の各章を参照されたい。

• 第1章：bmsonからBMSへの逆変換における離散格子量子化と彩色パッキング数理¶

  • 数理トピック: 連続時間座標をユークリッドの互除法により最適な離散格子へ射影する量子化数理、および多重発音の競合を区間グラフ彩色問題として解決する数理。
  • 実装への射影: C#のGCオーバーヘッドを完全に無力化する、空間 \(\mathcal{O}(1)\) のアロケーションフリーな Array-Backed Pool / Handle パターン（DOD）の実装。

• 第2章：短時間振幅包絡線と無音トリミングの数理¶

  • 数理トピック: ゼロクロス点での誤検知やポップノイズを防ぎつつ、二乗エネルギー累積和系列を用いて任意の窓幅を \(\mathcal{O}(1)\) で評価する定数時間局所エネルギー計算。
  • 実装への射影: AudioFeatureExtractor における平方根回避の二乗比較最適化。

• 第3章：LSHと超並列SIMD演算による高速化数理（実装極限）¶

  • 数理トピック: SimHash定理に基づき、数千次元の浮動小数点特徴量を単一の 64-bit 整数に詰め込み、近似相関距離を XOR ＋ POPCNT のハードウェア命令（\(\mathcal{O}(1)\)）に射影する離散代数。
  • 実装への射影: FastWaveCompare での Vector256 命令と、パイプラインハザードを防止する自動コンパイルループアンロール。

• 第4章：音響信号の等価性判定とピアソン相関係数¶

  • 数理トピック: 浮動小数点演算の桁落ちを防ぐ Welford のアルゴリズム（1パス漸化式計算）と、相関係数の数理的定義。
  • 実装への射影: PearsonAccumulator における定数空間 \(\mathcal{O}(1)\) での逐次オンライン累積演算。

• 第5章：周波数領域における時間アライメント¶

  • 数理トピック: 時間領域の \(\mathcal{O}(N^2)\) 探索を、畳み込み定理を用いて周波数領域の \(\mathcal{O}(N \log N)\) へと射影する数学的証明。
  • 実装への射影: FftAlignmentEngine でのゼロパディング境界条件処理と高速相互相関。

• 第6章：Union-Find グラフ理論による音声定義の同値類分割¶

  • 数理トピック: 類似関係の推移律の破れを「連結成分」として再定義し、経路圧縮とランク結合を用いてアッカーマン関数の逆関数 \(\alpha(M)\) で同値類を構築する証明。
  • 実装への射影: UnionFind クラスでのオブジェクトインスタンス生成を排除した、高速フラットポインタ配列。

• 第7章：時間軸シミュレーションエンジンとプレイアビリティ保存¶

  • 数理トピック: 単なる静的置換が招くコムフィルタ効果や同時発音数の限界を証明。Sweep-Lineを用いて、複雑な発音干渉を \(\mathcal{O}(L \log L)\) で精密にシミュレート。
  • 実装への射影: SimulationEngine による最適化安全性検証。

• 第8章：アムダールの法則とGCフリーなメモリ空間理論¶

  • 数理トピック: アムダールの法則によるタスク並列化（Thread-Level Parallelism）のスループット限界と、Large Object Heap (LOH) 断片化を回避するためのメモリプール理論。
  • 実装への射影: ParallelAudioComparisonEngine における ArrayPool<T> と Span<T> を用いた、完全並列かつアロケーションゼロのI/O・解析パイプライン。

0.7 理論の実証およびベンチマーク評価¶

本解説書で展開された数学的証明およびオーダー削減は、机上の空論ではない。実際のハードウェア上で実行された公式ベンチマーク（BenchmarkDotNet）において、その圧倒的なパフォーマンスが実証されている。

0.7.1 フルパイプライン実行における実測パフォーマンス¶

以下は、最適化パイプライン全体（音声波形の解析、LSHバケッティング、SIMD波形比較、FFTアライメント、Union-Find縮約、Sweep-Line検証まで）を統合実行した際の実測データである。

0.7.2 統計データに基づくアーキテクチャの堅牢性の論理的検証¶

• 平均実行時間 (Mean: 3.018 s)¶

  数万規模の計算空間を持つ巨大な譜面データのフル最適化が、わずか約3秒で完了している。これは、\(\mathcal{O}(M^2 \cdot N^2)\) の素朴な実装であれば数時間を要する処理が、直交する数学的次元圧縮によって実用的なリアルタイム処理へと昇華された決定的な証拠である。

• 標準偏差と誤差 (StdDev: 0.052 s / Error: 0.059 s)¶

  平均実行時間 3秒 に対して、標準偏差がわずか \(1.7\%\) 程度（約52ミリ秒）に収束している。この極めて低い分散は、GCの予測不能な停止やOSレベルのキャッシュミスが、アルゴリズムの力によって数理的に抑え込まれ、「常に一定の最大スループットで決定論的に実行されていること」を統計学的に証明している。

• 初回実行ペナルティと定常状態¶

  C# (.NET) ランタイムのアーキテクチャ上、初回実行時には JIT コンパイラによる機械語への翻訳、SIMD命令の Tier 1（最適化）への昇格、およびCPUキャッシュのウォームアップが発生する。そのため、アプリケーション起動直後の1回目の処理は上記より時間を要する。 本ベンチマークはこれらのウォームアップフェーズを厳格に分離・除外し、「ハードウェアとエンジンが完全に温まった定常状態における純粋な限界性能」を抽出している。

• メモリ空間 (Allocated: 331.06 MB)¶

  音声波形のオンメモリ展開と数万回の区間検証を伴う巨大な処理でありながら、総アロケーションが 331 MB に抑え込まれている。これは、第7章等で解説した Array-Backed Pool やインデックスポインタを用いた DOD により、実行時のヒープへの無駄なオブジェクト生成（GC爆発のトリガー）が物理的に排除されている結果である。