[论文解读] 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$ 和独立同分布噪声 $\mathbf{B}_n$ 驱动的 $({\mathbb{Z}}/p{\mathbb{Z}})^d$ 上仿射随机过程的混合时间。
• 确定过程 $\mathbf{X}_n$ 在 $({\mathbb{Z}}/p{\mathbb{Z}})^d$ 上收敛到均匀分布的条件。
• 通过利用 $T$ 的谱性质，特别是其特征值的模，改进先前的界（尤其是 Asci 的结果）。
• 基于 $T$ 的特征值是否位于 $\mathbb{C}$ 中的单位圆上，建立收敛速度的精确阈值。

提出的方法

• 将过程建模为 $\mathbf{X}_{n+1} = T\mathbf{X}_n + \mathbf{B}_n \pmod{p}$，初始条件 $\mathbf{X}_0 = \mathbf{0}$，且 $\mathbf{B}_n$ 以相等概率 $1/(d+1)$ 取值于 $\{\mathbf{0}, \mathbf{e}_1, \dots, \mathbf{e}_d\}$。
• 使用有限阿贝尔群 $({\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)$ 的收敛速度。
• 当特征值为单位根时，过程表现出周期性结构；分析表明，对某个适当的线性泛函的投影 $p(\mathbf{X}_n)$ 执行具有有界支撑的随机游走，导致混合缓慢。

实验结果

研究问题

• RQ1 当 $T$ 在 $\mathbb{C}$ 中无模为 1 的特征值时，仿射过程 $\mathbf{X}_{n+1} = T\mathbf{X}_n + \mathbf{B}_n \pmod{p}$ 在 $({\mathbb{Z}}/p{\mathbb{Z}})^d$ 上的混合时间是多少？
• RQ2 特征值为单位根的存在如何影响收敛到均匀分布的速度？
• RQ3 对于一般 $T$，能否将 $O((\log p)^2)$ 的混合时间上界改进为 $O((\log p)\log\log p)$？
• RQ4 当 $T$ 有单位根特征值时，$\Omega(p^b)$ 下界中 $b$ 的最大可能值是多少？
• RQ5 这些结果在多大程度上可推广到大小为 $p^d$ 的有限域上的随机过程？

主要发现

• 若 $T$ 无模为 1 的复特征值，则当 $n \geq C(\log p)^2$（$C > 0$ 为与 $p$ 无关的常数）且 $p$ 与 $\det(T)$ 互素时，有 $\|P_n - U\| \to 0$ 当 $p \to \infty$。
• 当 $T$ 有单位根特征值时，若 $n \leq p^b$（$b \leq 2$），则 $\|P_n - U\| \to 1$ 当 $p \to \infty$，表明混合缓慢。
• 当 $T$ 有单位根特征值时，收敛速度的下界为 $\Omega(p^b)$（$b > 0$），且该 $b$ 取决于 $T$ 的结构和初始状态。
• 当 $T$ 无模为 1 的特征值时，若 $n = O((\log p)^2)$，则变差距离 $\|P_n - U\|$ 的衰减快于 $1/p$，表明快速混合。
• 证明依赖于对分布傅里叶变换的有界性估计，并利用上界引理将谱衰减与总变差距离关联。
• 结果是紧确的，即下界中的 $b \leq 2$，且该界是紧的，因为 $b$ 可根据 $T$ 的不同而任意接近 2。