[論文レビュー] Generating Random Vectors in (Z/pZ)^d Via an Affine Random Process

Martin Hildebrand, Joseph McCollum|ArXiv.org|Jan 20, 2007

Random Matrices and Applications参考文献 7被引用数 4

ひとこと要約

本稿は、有限群 $(\mathbb{Z}/p\mathbb{Z})^d$ 上のアフィン確率過程を研究しており、$\mathbf{X}_{n+1} = T\mathbf{X}_n + \mathbf{B}_n \pmod{p}$ で定義される。ここで $\mathbf{B}_n$ は標準基底ベクトルおよびゼロに一様に分布する。$T$ が絶対値が 1 である複素固有値をもたない場合、$\mathbf{X}_n$ は $O((\log p)^2)$ ステップで一様分布に近づくことが証明されている。一方、$T$ が単位円上にある根の一つである固有値をもつ場合、収束には $\Omega(p^b)$ ステップが必要であり、$b \leq 2$ である。これらの結果は、有限群上のフーリエ解析と上界補題を用いて確立されている。

ABSTRACT

This paper considers some random processes of the form X_{n+1}=TX_n+B_n (mod p) where B_n and X_n are random variables over (Z/pZ)^d and T is a fixed d x d integer matrix which is invertible over the complex numbers. For a particular distribution for B_n, this paper improves results of Asci to show that if T has no complex eigenvalues of length 1, then for integers p relatively prime to det(T), order (log p)^2 steps suffice to make X_n close to uniformly distributed where X_0 is the zero vector. This paper also shows that if T has a complex eigenvalue which is a root of unity, then order p^b steps are needed for X_n to get close to uniform where b is a value which may depend on T and X_0 is the zero vector.

研究の動機と目的

• 固定行列 $T$ と i.i.d. ノイズ $\mathbf{B}_n$ によって駆動される $({\mathbb{Z}}/p{\mathbb{Z}})^d$ 上のアフィン確率過程の混合時間を分析すること。
• 過程 $\mathbf{X}_n$ が $({\mathbb{Z}}/p{\mathbb{Z}})^d$ 上で一様分布に収束する条件を特定すること。
• 特に $T$ の固有値の大きさに注目し、スペクトル的性質を活用することで、アスチの先行研究の境界を改善すること。
• $T$ の固有値が複素数平面の単位円上にあるかどうかに応じて、収束速度の鋭い閾値を確立すること。

提案手法

• 過程は $\mathbf{X}_{n+1} = T\mathbf{X}_n + \mathbf{B}_n \pmod{p}$ でモデル化され、$\mathbf{X}_0 = \mathbf{0}$ かつ $\mathbf{B}_n$ が $\{\mathbf{0}, \mathbf{e}_1, \dots, \mathbf{e}_d\}$ に等確率 $1/(d+1)$ で値をとる。
• $({\mathbb{Z}}/p{\mathbb{Z}})^d$ という有限アーベル群上のフーリエ解析が用いられ、特徴 $\rho_{\mathbf{c}}(\mathbf{b}) = e^{2\pi i \sum b_i c_i / p}$ が定義される。
• $P_n$ の分布のフーリエ変換が分析され、$\hat{P}_n(\mathbf{c}) = \sum_{\mathbf{s}} P_n(\mathbf{s}) \rho_{\mathbf{c}}(\mathbf{s})$ で与えられ、$\mathbf{c} \neq \mathbf{0}$ のときの $|\hat{P}_n(\mathbf{c})|$ の時間的挙動が追跡される。
• 上界補題が適用される：$\|P_n - U\|^2 \leq \frac{1}{4} \sum_{\rho^*} d_\rho \operatorname{Tr}(\hat{P}_n(\rho) \hat{P}_n(\rho)^*)$ であり、$\rho^*$ は非自明な既約表現を走る。
• $T$ の固有値が単位円上にない場合、$|\hat{P}_n(\mathbf{c})|$ の減衰が幾何的であることが示され、$O((\log p)^2)$ の収束が得られる。
• $T$ の固有値が単位根である場合、周期的構造を示す。解析により、適切な線形関数への射影 $p(\mathbf{X}_n)$ が有界なサポートを持つランダムウォークをなすことが示され、混合が遅いことがわかる。

実験結果

リサーチクエスチョン

• RQ1 複素数平面において絶対値が 1 である固有値を持たない $T$ の場合、$\mathbf{X}_{n+1} = T\mathbf{X}_n + \mathbf{B}_n \pmod{p}$ として定義されるアフィン過程の混合時間は何か？
• RQ2 固有値が単位根である場合、一様分布への収束速度にどのような影響を与えるか？
• RQ3 一般の $T$ に対して、$O((\log p)^2)$ の混合時間境界を $O((\log p)\log\log p)$ に改善できるか？
• RQ4 $T$ が単位根の固有値を持つ場合、$\Omega(p^b)$ の下界における $b$ の最大値は何か？
• RQ5 これらの結果は、サイズ $p^d$ の有限体上の確率過程へどの程度拡張可能か？

主な発見

• $T$ が複素数平面において絶対値が 1 である固有値をもたない場合、$p \to \infty$ のとき $n \geq C(\log p)^2$（$C > 0$ は $p$ に依存しない定数）であれば $\|P_n - U\| \to 0$ となる。ただし $p$ と $\det(T)$ は互いに素でなければならない。
• $T$ が単位根である固有値をもつ場合、$n \leq p^b$（$b \leq 2$）であれば $p \to \infty$ のとき $\|P_n - U\| \to 1$ となる。これは混合が遅いことを示している。
• $T$ が単位根の固有値をもつ場合、$\|P_n - U\|$ の下界は $\Omega(p^b)$（$b > 0$）であり、この $b$ は $T$ の構造と初期状態に依存する。
• $T$ が絶対値が 1 である固有値をもたない場合、$n = O((\log p)^2)$ のとき変動距離 $\|P_n - U\|$ は $1/p$ よりも速く 0 に収束する。これは高速混合を示している。
• 証明は、分布のフーリエ変換の評価と、上界補題を用いてスペクトルの減衰を全変動距離に関連付けることに依拠している。
• 結果は鋭く、下界における $b \leq 2$ であり、この境界はタイトである。すなわち、$T$ の構造に応じて $b$ は 2 に限りなく近づけることができる。