Графы-экспандеры — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
м (Вершинное расширение)
м (Спектральное расширение)
Строка 154: Строка 154:
  
 
Если <tex>G</tex> является [[Основные определения теории графов#Часто используемые графы|d-регулярным]], возможно определение в терминах [[Линейная алгебра 1 курс|линейной алгебры]] на основе собственных значений [[Матрица смежности графа|матрицы смежности]] <tex>A = A(G)</tex> графа <tex>G</tex>, где <tex>A_{ij}</tex> {{---}} число дуг между вершинами <tex>i</tex> и <tex>j</tex>. Поскольку <tex>A</tex> является симметричной, согласно спектральной теореме, <tex>A</tex> имеет <tex>n</tex> действительных собственных значений <tex>\lambda_{1} \ge \lambda_{2} \ge \cdots \ge \lambda_{n}</tex>. Известно, что эти значения лежат в промежутке <tex>[-d, d]</tex>.
 
Если <tex>G</tex> является [[Основные определения теории графов#Часто используемые графы|d-регулярным]], возможно определение в терминах [[Линейная алгебра 1 курс|линейной алгебры]] на основе собственных значений [[Матрица смежности графа|матрицы смежности]] <tex>A = A(G)</tex> графа <tex>G</tex>, где <tex>A_{ij}</tex> {{---}} число дуг между вершинами <tex>i</tex> и <tex>j</tex>. Поскольку <tex>A</tex> является симметричной, согласно спектральной теореме, <tex>A</tex> имеет <tex>n</tex> действительных собственных значений <tex>\lambda_{1} \ge \lambda_{2} \ge \cdots \ge \lambda_{n}</tex>. Известно, что эти значения лежат в промежутке <tex>[-d, d]</tex>.
Граф регулярен тогда и только тогда, когда вектор <tex>u\in \mathbb {R} ^{n} с u_{i}=1</tex> для всех <tex>i = 1 \ldots n</tex> является собственным вектором матрицы <tex>A</tex>, а его собственное число будет постоянной степенью графа. Таким образом, мы имеем <tex>Au = du</tex>, и <tex>u</tex> — собственный вектор матрицы <tex>A</tex> с собственными значениями <tex>\lambda_{1} = d</tex>, где <tex>d</tex> — степень вершин графа <tex>G</tex>. Спектральный зазор графа <tex>G</tex> определяется как <tex>d-\lambda_{2}</tex> и является мерилом спектрального расширения графа <tex>G</tex>.
+
Граф регулярен тогда и только тогда, когда вектор <tex>u\in \mathbb {R} ^{n}</tex> с <tex>u_{i}=1</tex> для всех <tex>i = 1 \ldots n</tex> является собственным вектором матрицы <tex>A</tex>, а его собственное число будет постоянной степенью графа. Таким образом, мы имеем <tex>Au = du</tex>, и <tex>u</tex> — собственный вектор матрицы <tex>A</tex> с собственными значениями <tex>\lambda_{1} = d</tex>, где <tex>d</tex> — степень вершин графа <tex>G</tex>. Спектральный зазор графа <tex>G</tex> определяется как <tex>d-\lambda_{2}</tex> и является мерилом спектрального расширения графа <tex>G</tex>.
  
 
Известно, что <tex>\lambda_{n} = -d</tex> тогда и только тогда, когда <tex>G</tex> {{---}} [[Двудольные графы|двудольный]]. Во многих случаях, например в лемме о перемешивании, необходимо ограничить снизу не только зазор между <tex>\lambda_{1}</tex> и <tex>\lambda_{2}</tex>, но и зазор между <tex>\lambda_{1}</tex> и вторым максимальным по модулю собственным значением:
 
Известно, что <tex>\lambda_{n} = -d</tex> тогда и только тогда, когда <tex>G</tex> {{---}} [[Двудольные графы|двудольный]]. Во многих случаях, например в лемме о перемешивании, необходимо ограничить снизу не только зазор между <tex>\lambda_{1}</tex> и <tex>\lambda_{2}</tex>, но и зазор между <tex>\lambda_{1}</tex> и вторым максимальным по модулю собственным значением:

Версия 15:56, 20 декабря 2017

Граф-экспандер (или расширяющийся граф, англ. expander graph) — в комбинаторике сильно разреженный граф, при этом связность определяется по вершинам, дугам или спектру. Это конечный ненаправленный мультиграф, в котором любое подмножество вершин, не являясь «слишком большим», имеет «сильную» связность. Различные формализации этих понятий дают различные типы экспандеров: рёберный расширитель, вершинный расширитель, и спектральный расширитель.

Любой связный граф является экспандером, однако различные связные графы имеют различные параметры расширителя. Полный граф имеет лучшие параметры расширителя, но имеет наибольшую возможную степень. Неформально говоря, граф является хорошим экспандером, если он имеет низкую степень и высокий параметр расширителя.

Основные определения

Определение:
Множество соседей для множества вершин [math]S[/math] называется [math]\Gamma(S) = \{v: \exists w \in S:\ (v, w) \in E\}[/math].

Замечание: в графе могут присутствовать кратные ребра и петли. Если у каждой вершины есть петли, то [math]\forall S \supseteq \Gamma(S)[/math].

Однородный (комбинаторный) экспандер

Определение:
Однородный (комбинаторный) экспандер называется граф [math]G = (V, E)[/math] с параметрами [math](n, d, \epsilon) \ (n = |V|[/math], d — степень каждой вершины, константа [math]\epsilon \lt 1)[/math], если выполняется условие: для [math]\ \forall S \subseteq V, \ |S| \leq n/2 \ \ \exists \Gamma(S)[/math], которое достаточно велико, то есть [math]|\Gamma(S)| \gt (1 + \epsilon)|S|[/math].

Замечание: чем больше [math]\epsilon[/math], тем сильнее свойство расширения.

Теорема:
Пусть [math]\epsilon[/math] — некоторое положительное число меньшее 1. Тогда для всех достаточно больших четных [math]d = d(\epsilon)[/math] и всех [math]n[/math] существует однородный (комбинаторный) экспандер с параметрами [math](n, d, \epsilon)[/math].
Доказательство:
[math]\triangleright[/math]

Мы выберем граф случайно и покажем, что с положительной (и даже довольно близкой к 1) вероятностью такой граф оказывается экспандером. Отсюда будет следовать, что экспандеры существуют.

Оценим вероятность того, что полученный в результате граф не является экспандером, если найдется множество вершин [math]S[/math] (состоящее из не более, чем [math]n/2[/math] вершин), все соседи которого лежат в некотором множестве [math]T[/math], состоящем из [math]\lfloor (1 + \epsilon)|S|\rfloor[/math] вершин.

Зафиксируем некоторые множества вершин [math]S[/math] и [math]T[/math]. Зафиксируем номер перестановки [math]\pi_{i}[/math]. Вероятность того, что для каждой вершины [math]v \in S[/math] второй конец ребра {[math]v, pi_{i}(v)[/math]} попадёт в [math]T[/math], равна

[math]\frac{|T|}{n} \cdot \frac{|T| - 1}{n - 1} \cdot \ldots \cdot \frac{|T| - |S| + 1}{n - |S| + 1} \leq (\frac{|T|}{n})^{|S|}[/math].

Поскольку мы выбираем [math]d/2[/math] перестановок независимо, вероятность того, что данное событие произойдёт для всех [math]i[/math] не превосходит [math](\frac{|T|}{n})^{(d/2)|S|}[/math]. Таким образом,

Prob[свойство экспандерности графа нарушено] [math]\leq \sum_{S, T} (\frac{|T|}{n})^{(d/2)|S|}[/math], где суммирование происходит по всем множествам вершин [math]S[/math] размера не более [math]n/2[/math] и по всем множествам [math]T[/math] размера [math]\lfloor(1 + \epsilon)|S|\rfloor[/math].

Оценим сумму:

[math]\sum_{S, T}(\frac{|T|}{n})^{(d/2)|S|} \leq \sum_{s = 1}{n/2}\binom{s}{n} \cdot \binom{(1 + \epsilon)s}{n} \cdot (\frac{(1 + \epsilon)s}{n})^{sd/2}[/math]. (1)

Каждый биномиальный коэффициент [math]\binom{k}{n}[/math] можно оценить сверху величиной [math](\frac{ne}{k})^{k}[/math]. Тогда

[math] \sum_{s = 1}^{n/2} (\frac{ne}{s})^{s} \cdot (\frac{ne}{(1 + \epsilon)s})^{(1 + \epsilon)s} \cdot (\frac{(1 + \epsilon)s}{n})^{sd/2} = \sum_{s = 1}^{n/2} [(s / n)^{d/2 - 2 - \epsilon} \cdot (1 + \epsilon)^{d/2} \cdot \frac{e^{1 + \epsilon} }{(1 + \epsilon)^{1 + \epsilon} }]^s[/math]. (2)

Заметим, что [math]s \leq n/2[/math], а [math]1 + \epsilon \lt 2[/math]. Таким образом, можно подобрать такое

[math]d = d(\epsilon)[/math], чтобы выражение в квадратных скобках в правой части (2) было меньше [math]1/2[/math] при всех значениях [math]\epsilon[/math]. Следовательно, сумма (1) меньше единицы. Что означает, что с положительной вероятностью случайный граф является экспандером с параметрами [math](n, d, \epsilon)[/math].
[math]\triangleleft[/math]

Двудольный экспандер

Определение:
Двудольный экспандер называется двудольный граф [math]G = (L,\ R,\ E)\ (L[/math] и [math]R[/math] — вершины левой и правой доли соответственно, [math]E[/math] — ребра графа [math]G[/math]) с параметрами [math](n,\ m,\ d,\ e)\ (n = |L|,\ m = |R|,\ d[/math] - степень всех вершин в левой доле), если выполняется условие: для [math]\forall S \subseteq L \ |S| \leq k \ \ |\Gamma(S)| \gt (1 - \epsilon)d|S|[/math].

Замечание 1: чем меньше значение [math]\epsilon[/math] в данном определении, тем сильнее свойство расширения.

Замечание 2: в приложениях как правило используют двудольные экспандеры с [math]\epsilon \lt 1/2[/math], а для применения в теории кодирования (для построения экспандерных кодов) часто требуются двудольные экспандеры с ещё меньшими значениями [math]\epsilon[/math].

Теорема:
Пусть [math]\epsilon[/math] — некоторое положительное число. Тогда для [math]\forall \ n[/math] и [math]k \leq n[/math] найдется [math]d = O(\log n)[/math] и [math]m = O(dk)[/math] такие, что [math]\exists[/math] двудольный экспандер с параметрами [math](n, m, d, k, \epsilon)[/math].
Доказательство:
[math]\triangleright[/math]

Выберем случайный граф, то есть для [math]\forall \ v \in L[/math] случайно и независимо выбираем [math]d[/math] соседей в [math]R[/math] (разрешаются кратные ребра). Покажем с большой вероятностью такой граф оказывается экспандером.

Граф не является экспандером, если [math]\exists T \subset R: \ |T| = (1 - \epsilon)d|S|[/math] и [math]\exists S \subset L: \ \Gamma(S) = T[/math].

Поскольку при случайном выборе графа мы приводим все [math]nd[/math] рёбер случайно и независимо, то для [math]\forall[/math] ребра вероятность того, что его правый конец окажется в фиксированном множестве [math]T[/math], равна [math]|T|/m[/math]. Следовательно,

[math]Prob[[/math]свойство экспандерности графа нарушено[math]] \leq \ \sum_{S, T}^{}(\frac{|T|}{m})^{sd}[/math],

где суммирование происходит по всем множествам [math]S \subset L \ |S| \leq k \ \forall T \subset R \ |T| = (1 - \epsilon)d|S| [/math]. Оценим данную сумму сверху:

[math]\sum_{s = 1}^{k} \binom{s}{n} \* \binom{(1 - \epsilon)sd}{m} \* (\frac{(1 - \epsilon)sd}{m})^{sd}[/math].

Оценивая биномиальные коэффициенты, получаем, что сумма не превосходит

[math] \sum_{s = 1}^{k} (\frac{ne}{s})^{s} \cdot (\frac{me}{(1 - \epsilon)sd})^{(1 - \epsilon)sd} \cdot (\frac{(1 - \epsilon)sd}{m})^{sd} \leq \sum_{s = 1}^{k} [\frac{ne}{s} \cdot (\frac{e^{(1 - \epsilon)/\epsilon} \cdot (1 - \epsilon)sd}{m})^{\epsilon d}]^{s} [/math].

Положим [math]m = Const * kd[/math] (с достаточно большим значением [math]Const[/math]), чтобы для всех возможных s выполнялось неравенство [math]\frac{e^{(1 - \epsilon)/\epsilon} \cdot (1 - \epsilon)sd}{m} \leq 1/2[/math]. Тогда выражение в квадратных скобках не превосходит [math]ne \cdot (1/2) ^{\epsilon d}[/math]. Остаётся выбрать [math]d \gt \frac{1}{\epsilon}\log(2en)[/math], и мы получаем

[math]ne \cdot (1/2)^{\epsilon d} \lt 1/2[/math].

Таким образом, для выбранных значений параметров сумм не превосходят 1. Это означает, что с положительной вероятностью случайный двудольный граф является экспандером c параметрами [math](n, \ m, \ k, \ d, \ \epsilon)[/math].
[math]\triangleleft[/math]

Замечание: константы в [math]O[/math] зависят от [math]\epsilon[/math].

теорема о паросочетаниях

Теорема:
Пусть [math]G[/math]двудольный граф. [math]G[/math] имеет совершенное паросочетание тогда и только тогда, когда [math]\forall S \supseteq L |S| \leq \Gamma(S)[/math].
Доказательство:
[math]\triangleright[/math]

Будем доказывать по индукции

Для [math]|L| = 1[/math] очевидно

Предположим, что верно для [math]L: |L| \lt m[/math]

Рассмотрим два случая:

1. [math]\forall S \supseteq L, |S| \neq 0[/math] строго расширяется, то есть [math]|\Gamma(S)| \gt |S|[/math].

Выберем [math]\forall x \supset V(L) : (x, y) \in E[/math]. Тогда рассмотрим двудольный граф [math]G^{*}[/math], у которого левая доля [math]L^{*} = L \diagdown {x}[/math], а правая [math]R^{*} = R \diagdown {y}[/math]. Так как [math]\forall S \supset L[/math] удовлетворяет строгому неравенству теоремы, то каждое подмножество [math]L^{*}[/math] удовлетворяет неравенству, поскольку только одна вершина [math]y[/math] была удалена из [math]R[/math]. Следовательно, по предположению индукции меньший граф [math]G^{*}[/math] имеет паросочетание. К этому паросочетанию добавляем ребро (x, y) что дает совершенное паросочетание.

2. Пусть [math]\exists\ T \in L, \ |T| \neq 0[/math] такое, что [math]|\Gamma(T)| = |T|[/math].

[math]\measuredangle[/math] порожденные графы [math]G_{1} = T \cup \Gamma(T)[/math] и [math]G_{2} = L \diagdown T \cup R \diagdown \Gamma(T)[/math].

По предположению индукции [math]G_{1}[/math] имеет совершенное паросочетание (Заметим, что предположение индукции не может использовано непосредственно на [math]G_{2}[/math]). Пусть [math]S \supseteq L \diagdown T[/math], тогда:

[math] \Gamma_{G_{2}}(S) = \Gamma_{G}(S \cup T) \diagdown \Gamma_{G}(T) \ \Rightarrow[/math] [math]|\Gamma_{G_{2}}| \geqslant |S \cup T| - |T| = |S|[/math]

Неравенство верно, поскольку [math]S \cup T[/math] удовлетворяет [math]|\Gamma_{G}(S \cup T)| \geqslant |S \cup T|[/math] и по предположению [math]|\Gamma_{G}(T)| = |T|[/math].

Следовательно, [math]G_{2}[/math] также удовлетворяет неравенству теоремы и по предположению индукции имеет паросочетание. Объединение совершенных паросочетаний [math]G_{1}[/math] и [math]G_{2}[/math] — паросочетание для [math]G[/math].
[math]\triangleleft[/math]

Определения расширений

Реберное расширение

Рёберное расширение (также изопериметрическое число или константа Чигера) [math]h(G)[/math] графа [math]G[/math] для [math]n[/math] вершин определяется как

[math]h(G) = \min\limits_{0 \lt |S|\leqslant \frac{n}{2}} \frac{|\partial_{\text{out}}(S)|}{|S|}[/math],

где минимум берётся по всем непустым множествам [math]S[/math] не более чем [math]n/2[/math] вершин и [math]\partial(S)[/math]граничные дуги множества [math]S[/math], то есть, множество дуг с единственной вершиной в [math]S[/math].

Вершинное расширение

Вершинное изопериметрическое число [math]h_{out}(G)[/math] (также вершинное раширение) графа [math]G[/math] определяется как

[math]h_{out}(G) = \min\limits_{0 \lt |S|\leqslant \frac{n}{2}} \frac{|\Gamma_{\text{out}}(S)|}{|S|}[/math],

где [math]\Gamma_{\text{out}}(S)[/math]внешняя граница [math]S[/math], то есть множество вершин из [math]V(G)\setminus S[/math], имеющих как минимум одного соседа в [math]S[/math].

[math]\Gamma_{\text{out}}(S) = \{v \in V(G) \diagdown S: \exists w \in S (v, w) \in E\}[/math].

Вершинное изопериметрическое число [math]h_{in}(G)[/math] графа [math]G[/math] определяется как

[math]h_{in}(G) = \min\limits_{0 \lt |S|\leqslant \frac{n}{2}} \frac{|\Gamma_{\text{in}}(S)|}{|S|}[/math],


где [math]\Gamma_{in}(S)[/math]внутренняя граница [math]S[/math], то есть множество вершин [math]S[/math], имеющих как минимум одного соседа в [math]V(G)\setminus S[/math].

[math]\Gamma_{in}(S) = \Gamma(S) \diagdown \Gamma_{in}(S)[/math]

Спектральное расширение

Если [math]G[/math] является d-регулярным, возможно определение в терминах линейной алгебры на основе собственных значений матрицы смежности [math]A = A(G)[/math] графа [math]G[/math], где [math]A_{ij}[/math] — число дуг между вершинами [math]i[/math] и [math]j[/math]. Поскольку [math]A[/math] является симметричной, согласно спектральной теореме, [math]A[/math] имеет [math]n[/math] действительных собственных значений [math]\lambda_{1} \ge \lambda_{2} \ge \cdots \ge \lambda_{n}[/math]. Известно, что эти значения лежат в промежутке [math][-d, d][/math]. Граф регулярен тогда и только тогда, когда вектор [math]u\in \mathbb {R} ^{n}[/math] с [math]u_{i}=1[/math] для всех [math]i = 1 \ldots n[/math] является собственным вектором матрицы [math]A[/math], а его собственное число будет постоянной степенью графа. Таким образом, мы имеем [math]Au = du[/math], и [math]u[/math] — собственный вектор матрицы [math]A[/math] с собственными значениями [math]\lambda_{1} = d[/math], где [math]d[/math] — степень вершин графа [math]G[/math]. Спектральный зазор графа [math]G[/math] определяется как [math]d-\lambda_{2}[/math] и является мерилом спектрального расширения графа [math]G[/math].

Известно, что [math]\lambda_{n} = -d[/math] тогда и только тогда, когда [math]G[/math]двудольный. Во многих случаях, например в лемме о перемешивании, необходимо ограничить снизу не только зазор между [math]\lambda_{1}[/math] и [math]\lambda_{2}[/math], но и зазор между [math]\lambda_{1}[/math] и вторым максимальным по модулю собственным значением:

[math]\lambda=\max\{|\lambda_{2}|, |\lambda_{n}|\}[/math]

Поскольку это собственное значение соответствует некоторому собственному вектору, ортогональному [math]u[/math], его можно определить, используя отношение Рэлея: [math]\lambda=\max\limits_{0\neq v\perp u} \frac{\|Av\|_2}{\|v\|_2},[/math]

где

[math]\|v\|_2=\left(\sum_{i=1}^n v_i^2\right)^{1/2}[/math]

— евклидова норма вектора [math]v\in \mathbb {R} ^{n}[/math].

Нормализованная версия этого определения также широко используется и более удобна для получения определённых результатов. В таком случае рассматривается матрица [math]{\tfrac {1}{d}}A[/math], являющаяся матрицей переходов графа G. Все её собственные значения лежат между [math]-1[/math] и [math]1[/math]. Если граф не регулярен, спектр графа может быть определён аналогичным образом, используя собственные значения матрицы Кирхгофа. Для направленного графа используются сингулярные значения матрицы сопряжения A, которые равны квадратным корням из собственных значений симметричной матрицы [math]A^TA[/math].

Конструирование

Существуют три основные стратегии создания семейств графов расширений. Первая стратегия — алгебраическая и теоретико-групповая, вторая — аналитическая, использующая аддитивную комбинаторику, и третья — комбинаторная, использующая зигзаг-произведение и связанные комбинаторные произведения.

Маргулис-Габбер-Галил

Алгебраическое конструирование, основанное на графах Кэли, известно для различных вариантов экспандеров. Следующее конструирование принадлежит Маргулису и было проанализировано Габбером (Gabber) и Галилом (Galil). Для любого натурального [math]n[/math] строим граф, [math]G_{n}[/math] со множеством вершин [math]\mathbb Z _{n}\times \mathbb {Z} _{n}[/math], где [math]\mathbb {Z} _{n}=\mathbb Z /n \mathbb Z[/math] . Для любой вершины [math](x,y)\in \mathbb {Z} _{n}\times \mathbb {Z} _{n}[/math], её восемь соседей будут

[math](x \pm 2y,y), (x \pm (2y+1),y), (x,y \pm 2x), (x,y \pm (2x+1)).[/math]

Выполняется следующая теорема:

Теорема:
Для всех [math]n[/math] граф [math]Gn[/math] второе по величине собственное число [math]\lambda(G)\leq 5 \sqrt{2}[/math].

Граф Рамануджана

[math]d[/math]-регулярный граф называется графом Рамануджана, если его второе по модулю собственное число не превосходит [math]2{\sqrt {d-1}}[/math]. Любоцкий, Сарнак, Филлипс и Маргулис указали явную конструкцию графов Кэли, являющихся графами Рамануджана. Опишем эту конструкцию. Пусть [math]p[/math] и [math]q[/math] простые числа, [math]p = 1 \bmod 4[/math] и [math]q = 1 \bmod 4[/math]. В качестве группы [math]G[/math] возьмём [math]PGL(2, \mathbb Z/q{\mathbb Z})[/math], т.е. невырожденные матрицы 2 × 2 над полем вычетов по модулю [math]q[/math], профакторизованные по отношению пропорциональности (с обычной операцией матричного умножения). Далее мы зададим в этой группе симметричное множество [math]S[/math]. Выберем такое целое [math]i[/math], что [math]i^2 = −1 \bmod q[/math]. Можно доказать, что тогда имеется ровно [math](p + 1)[/math] целочисленное решение уравнения

[math]a_0^2 + a_1^2 + a_3^2 = p[/math]

такое, что [math]a_0[/math] положительно и нечётно, а [math]a_1[/math], [math]a_2[/math], [math]a_3[/math] чётны. Каждой такой четвёрке сопоставим матрицу

[math]A = \begin{pmatrix} a_0 + ia_1 & a_2 + ia_3\\ -a_2 + ia_3 & a_0 - ia_1 \end{pmatrix}[/math]

Эти матрицы образуют множество S. Нетрудно понять, что граф Кэли [math](G, S)[/math] состоит из [math]\Theta(q^3)[/math] вершин, и степень каждой вершины равна [math](p + 1)[/math]. Свойства данного графа зависят от соотношения [math]p[/math] и [math]q[/math]. Рассмотрим случай, когда p является квадратичным вычетом по модулю [math]q[/math]. Тогда полученный граф Кэли состоит из двух связных компонент (поскольку все матрицы из [math]S[/math] лежат в подгруппе [math]G[/math] индекса два — подгруппе матриц, определитель которых является квадратичным вычетом). Обозначим [math]X^{p,q}[/math] связную компоненту полученного графа. Можно доказать, что у [math]X^{p,q}[/math] второе по абсолютной величине собственное число не превосходит [math]2{\sqrt{p}}[/math], т.е. мы получили граф Рамануджана.

Примеры применения экспандеров

Коды, исправляющие ошибки

С помощью расширяющего графа можно посторить линейный код, позволяющий исправлять ошибки в доле [math]\delta = 1/(2000d)[/math] битов. Чтобы задать линейный код с длиной кодового слова n, достаточно описать его проверочную матрицу [math]H[/math] [math]([/math]слово [math]x \in \{0, 1\}^n[/math] является кодовым словом если и только если [math]Hx^t = 0)[/math]. Другими словами, нужно задать систему линейных уравнений для переменных [math]x_1, . . . , x_n;[/math] решения этой системы и будут кодовыми словами.

Увеличение вероятности успеха в алгоритмах с датчиком случайных чисел

Определение:
Язык [math]L[/math] принадлежит сложностному классу [math]RP[/math], если существует полиномиальный алгоритм [math]A[/math] такой что

1. для [math]x \in L[/math] для всех [math]r \in \{0, 1\}^{poly(n)}[/math] [math]A(x, r) = 1[/math]

2. для [math]x \notin L[/math] не более чем для 1/2000 всех [math]r \in \{0, 1\}^{poly(n)}[/math] может выполняться [math]A(x, r) = 1[/math]


Покажем, что для любого [math]\epsilon \gt 0[/math] полиномиальный вероятностный алгоритм A можно модифицировать таким образом, чтобы вероятность ошибки уменьшилась до [math]\epsilon[/math], а число используемых случайных битов не изменится.

Пусть исходный алгоритм использует [math]k = k(n)[/math] случайных битов для вычислений на входах длины [math]n[/math]. Зафиксируем [math](2^k, 2^k, d)[/math] - экспандер [math]G[/math], где [math]d \gt 1/\epsilon[/math]. Новый алгоритм действует следующим образом: выбирается случайная вершина [math]v[/math] из левой доли графа (для этого требуется [math]k[/math] случайных битов); затем исходный алгоритм [math]A[/math] последовательно запускается на всех [math]d[/math] наборах случайных битов, соответствующих соседям вершины [math]v[/math]. Если все полученные ответы равны [math]1[/math], новый алгоритм также возвращает единицу; в противном случае возвращается ноль. Покажем, что у нового алгоритма вероятность ошибки не превосходит [math]1/d[/math]. В самом деле, обозначим [math]B = B(x)[/math] множество таких вершин [math]w[/math] из правой доли графа, которые соответствуют неверному ответу старого алгоритма на входе [math]x[/math]; аналогично, обозначим [math]S = S(x)[/math] множество таких вершин v из левой доли графа, которые для которых новый алгоритм даёт неверный ответ на входе [math]x[/math]. Очевидно, [math]S[/math] состоит из вершин, все соседи которых лежат в B. Предположим, что [math]S[/math] содержит не менее [math]n/(1000d)[/math] вершин. Выберем среди них ровно [math]n/(1000d)[/math] вершин и назовём это множество [math]S'[/math]. По свойству экспандера, имеем

[math]|\Gamma(S')| \geqslant \frac{7}{8}d\frac{n}{1000d}=7n/8000 \gt n/2000[/math]

Это противоречит тому, что все соседи [math]S'[/math] лежат в [math]B[/math]. В данном случае нам нужна явная в более сильном (чем в первом примере) смысле конструкция экспандера. Размер графа экспоненциално растёт с увеличенем [math]k[/math], и нам необходим алгоритм, который по заданному номеру вершины [math]v[/math] (из левой доли графа) за время [math]poly(k)[/math] находит список номеров всех соседей этой вершины (в правой доле графа).

Хранение множества со сверхбыстрым запросом элементов

Мы организуем хранение m-элементного множества [math]S \subset \{1, . . . , n\}[/math] в виде описания [math]X[/math], состоящего из [math]O(m log n)[/math] битов. При этом проверка принадлежности [math]a \in S[/math] будет производиться чрезвычайно быстро. А именно, мы построим такой вероятностный алгоритм, который по любому входу [math]a[/math] запрашивает из [math]X[/math] один бит; если этот бит оказывается равным единице, то алгоритм отвечает, что [math]a[/math] является элементом [math]S[/math]; в противном случае алгоритм говорит, что [math]a[/math] множеству не принадлежит. При этом для каждого [math]a \in \{1, . . . , n\}[/math] алгоритм ошибается с вероятностью не более [math]1/3[/math].

Чтобы построить нужное нам хранилище [math]X[/math], мы сначала зафиксируем некоторый экспандер, у которого левая доля [math]L[/math] состоит из [math]n[/math] вершин, правая [math]R[/math] из [math]k = O(m log n)[/math] вершин, степень всех вершин левой доли одинакова и равна некоторому [math]d[/math], и для каждого [math]A \subset L[/math] размера не более [math]2m[/math]

[math]|\Gamma(A)| \geqslant \frac{7}{8}d|A|[/math]

[math]X[/math] будет состоять в разметке вершин правой доли нулями и единицами. Эту разметку нужно выбрать таким образом, чтобы у каждой вершины из [math]S[/math] не менее [math]2/3[/math] соседей были помечены единицей, а у каждой вершины не из [math]S[/math] не менее [math]2/3[/math] соседей были помечены нулями.[math]X[/math] будет состоять в разметке вершин правой доли нулями и единицами. Эту разметку нужно выбрать таким образом, чтобы у каждой вершины из [math]S[/math] не менее [math]2/3[/math] соседей были помечены единицей, а у каждой вершины не из [math]S[/math] не менее [math]2/3[/math] соседей были помечены нулями.

Остаётся объяснить, как построить нужную нам разметку правой доли графа. Будем строить её последовательными приближениями. Сначала пометим всех соседей всех вершин из [math]S[/math] единицами, а все остальные вершины – нулями. На данной разметке наш алгоритм с вероятностью [math]1[/math] возвращает правильный ответ для всех [math]a \in S[/math]. Однако для [math]a[/math] не из [math]S[/math] проверочный алгоритм может работать неверно. Обозначим T множество всех таких вершин вне [math]S[/math], у которых более [math]d/3[/math] соседей помечены единицей. Поменяем разметку: пометим всех соседей [math]T[/math] нулём. Теперь разметка может стать плохой для части вершин из [math]S[/math]. Обозначим [math]S'[/math] множество всех таких вершин из [math]S[/math], у которых более [math]d/3[/math] соседей помечены нулями. Далее, поменяем разметку у всех соседей [math]S'[/math] на единицы. После этого может вновь возникнуть множество ‘неправильных’ вершин вне [math]S[/math], и т.д.

Чтобы доказать, что данный процесс в конце концов сойдётся, нужно показать, что на каждом шаге число ‘проблемных’ вершин уменьшается в константу раз. Поскольку все шаги аналогичны, достаточно разобрать самый первый: докажем, что [math]T[/math] в константу раз меньше, чем [math]S[/math]. Мы воспользуемся тем, что для [math]S \cup T[/math] выполнено свойство расширения:

[math](7/8)d(|S|+|T|) \leqslant |\Gamma(S \cup T)| \leqslant d|S| + (2/3)d|T|[/math]

Откуда получаем [math]|T| \leqslant 3/5|S|[/math].

Приложения и полезные свойства

Первоначально интерес к экспандерам возник с целью построения устойчивой сети (телефонов или компьютеров) — число дуг графов расширения с ограниченным значением регулярности растет линейно по отношению к числу вершин.

Экспандеры нашли широкое применение в теории вычислительных машин и систем, для построения алгоритмов, в корректирующих кодах, экстракторах, генераторах псевдослучайных чисел, сетях сортировки и компьютерных сетях. Они также используются для доказательства многих важных результатов в теории вычислительной сложности, таких как [math]SL=L[/math] и Теорема PCP. В криптографии экспандеры используются для создания хеш-функций.

Ниже приведены некоторые свойства экспандеров, считающиеся полезными во многих областях.

Лемма о перемешивании

Лемма о перемешивании утверждает, что для любых двух подмножеств вершин [math]S,T\subseteq V[/math] число рёбер между [math]S[/math] и [math]T[/math] примерно равно числу рёбер в случайном [math]d[/math]-регулярном графе. Приближение тем лучше, чем меньше [math]\lambda =\max\{|\lambda _{2}|[/math],[math]|\lambda _{n}|\}[/math]. В случайном [math]d[/math]-регулярном графе, также как и в случайном графе Эрдёша — Реньи с вероятностью [math]{\tfrac {d}{n}}[/math] выбора ребра, ожидается [math]{\tfrac {d}{n}}\cdot |S|\cdot |T|[/math] рёбер между [math]S[/math] и [math]T[/math].

Более формально, пусть [math]E(S, T)[/math] обозначает число рёбер между [math]S[/math] и [math]T[/math]. Если эти два множества пересекаются, дуги в пересечении считаются дважды, так что

[math]E(S,T)=2|E(G[S\cap T])|+E(S\setminus T,T)+E(S,T\setminus S)[/math]. Лемма о перемешивании утверждает, что

[math]\left|E(S,T)-{\frac {d\cdot |S|\cdot |T|}{n}}\right|\leq d\lambda {\sqrt {|S|\cdot |T|}}[/math], где [math]\lambda [/math] — абсолютное значение нормализованного второго по величине собственного значения.

Недавно Билу (Bilu) и Линайл (Linial) показали, что обратное тоже верно, то есть, при условии выполнения неравенства из леммы второе по величине собственное значение равно [math]O(d\lambda \cdot (1+\log(1/\lambda )))[/math].

Блуждания по экспандеру

Согласно границе Чернова, если выбирать много независимых случайных значений из интервала [math][−1, 1][/math], с большой степенью вероятности среднее выбранных значений будет близко к математическому ожиданию случайной переменной. Лемма о блуждании по экспандеру, согласно статьям Аджтари, Комлоша и Семереди и Гилмана, утверждает, что то же самое верно и для блужданий по экспандеру. Это полезно в теории дерандомизации, поскольку блуждание по экспандеру использует много меньше случайных бит, чем случайная независимая выборка.

Источники информации