Модели случайных графов и их применения

УДК 519.175.4

А.М. Райгородский

Московский государственный университет им. М.В. Ломоносова Московский физико-технический институт

ООО «Яндекс»

Модели случайных графов и их применения

В статье дается обзор основных современных направлений в теории случайных графов. Делается акцент на связь моделей случайного графа с транспортной проблематикой.

Ключевые слова: случайный граф, Интернет, транспортные сети.

I. Введение

Теория графов играет огромную роль в фундаментальной и прикладной математике. Ей посвящены сотни монографий и тысячи — если не десятки тысяч — статей. Разумеется, мы не можем ставить перед собою цель дать на этих страницах сколь-нибудь подробное изложение теории графов. Нас будет интересовать лишь одно направление, которое с каждым годом становится все более актуальным. В рамках этого направления графы изучаются с вероятностной точки зрения. Типичная постановка вопроса (говоря не совсем строго) такова: велика ли вероятность того, что граф обладает данным свойством? Вопрос исключительно важный, и мы в этом не раз убедимся ниже. Правда, в нем ни слова не сказано о том, как именно мы понимаем термин «вероятность». Всякий человек, имеющий представление об аксиоматике Колмогорова, хорошо знает, что можно вложить множество разных смыслов в этот термин. И его можно действительно определять по-разному. В зависимости от определения получится та или иная модель случайного графа. С чисто математических позиций любая такая модель имеет право на существование. Однако для приложений — в том числе приложений к транспортной проблематике — некоторые из этих моделей более интересны, некоторые — менее. Соответственно ниже мы расскажем о двух классах моделей, каждый из которых за десятилетия, прошедшие с момента своего появления, зарекомендовал себя плодотворным как в рамках «чистой» математики, так и в рамках ее разнообразных приложений, среди которых надежность транспортной сети, рост Интернета и других социальных и биологических сетей, теория алгоритмов и пр.

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

II. Модель Эрдеша—Реньи

Этот раздел мы посвятим описанию модели случайного графа, которая возникла исторически первой. На рубеже 50-х и 60-х годов ХХ века эту модель предложили классики современной комбинаторики и теории вероятностей П. Эрдеш и А. Реньи (см. [1-3]). Отметим, что Эрдеш — это, пожалуй, одна из самых ярких фигур в математике ХХ века. Ему принадлежат сотни статей и задач, которые оказали огромное влияние на развитие многих математических дисциплин. Реньи также сыграл значительную роль в формировании венгерской вероятностной школы, и его именем назван математический институт в Будапеште.

ІІ.1. Формальное описание модели

Пусть дано множество Уп = {1, ..., п}, элементы которого мы назовем вершинами. Именно на этом множестве мы будем «строить» случайный граф. Понятно, стало быть, что случайным будет множество ребер графа. Мы не хотим сейчас рассматривать графы с кратными ребрами (мультиграфы), графы с петлями (псевдографы) и ориентированные графы (орграфы). Поэтому мы считаем, что потенциальных ребер у графа не больше, чем СП штук. Будем соединять любые две вершины і и і ребром с некоторой вероятностью р Є [0,1] независимо от всех остальных С2 — 1 пар вершин. Иными словами, ребра появляются в соответствии со стандартной схемой Бернулли, в которой С2 испытаний и «вероятность успеха» р. Обозначим через Е случайное множество ребер, которое возникает в результате реализации такой схемы. Положим О = (Уп,Е). Это и есть случайный граф в модели Эрдеша-Реньи.

Если записывать приведенное только что определение в формате аксиоматики Колмогорова, то мы имеем вероятностное пространство

О(п,р) = (Пп,^п,Рп,р);

в котором

Пп = {О = (Уп,Е)}, Тп = 2Пп, Рп,р(О) =

= р1Е1

Здесь через |А| обозначена мощность множества А (количество элементов в нем), а 2А — это совокупность всех подмножеств множества А.

Элемент сигма-алгебры Тп — это набор графов. Если нам хочется найти вероятность, с которой граф на п вершинах обладает данным свойством А, то мы просто берем множество Л € Тп, состоящее из всех графов, для которых выполнено свойство А, и вычисляем

Рпр(Л) = ^ Рпр(О).

Таким образом, вероятность того, например, что случайный граф связен, — это величина, равная сумме вероятностей всех связных графов (на фиксированном множестве вершин). Казалось бы, все совсем просто и мы вряд ли имеем шансы наткнуться здесь на нечто особенно интересное. Однако дело обстоит прямо противоположным образом: специфика вероятностных методов, которые эффективно работают в задачах о случайных графах, позволит нам пронаблюдать весьма нетривиальные и, главное, неожиданные явления, которые возникают даже в этой простой модели и которые влекут приятные следствия для приложений.

Прежде чем двигаться дальше, сделаем еще ряд полезных замечаний. Во-первых, если р = то, как видно из формулы (1), вероятность любого графа равна 2-Сп. Иными словами, в этом специальном случае все графы равновероятны и всякое утверждение о вероятности какого-либо свойства — это, по сути, утверждение о доле графов, данным свойством обладающих.

В действительности, мы не только не обязаны предполагать, что р = ^ (хотя и этот случай очень важен), мы даже можем считать, что с ростом величины п (числа вершин) вероятность р возникновения ребра изменяется. Иначе говоря, р = р(п) — любая функция, значения которой заключены между нулем и единицей. Как правило, в науке о случайных графах важны даже не сами вероятности событий, но их предельные значения. Почему это так, мы скоро увидим.

Скажем, наконец, что свойство выполнено почти всегда, если его вероятность стремится к единице при п ^ж>.

11.2. Транспортная интерпретация модели

Представим себе, что в некоторой стране есть 10 городов, которые попарно соединены дорогами. Это довольно сильное предположение, но пока сохраним его. Допустим, каждая из дорог за определенный срок изнашивается (то есть становится непроезжей) с известной вероятностью д. При этом износ данной дороги никак не зависит от совокупного износа остальных дорог. Спрашивается: какова максимальная вероятность д, при

которой с вероятностью больше 1/2 не исчезнет возможность перемещения между любыми двумя городами? По существу, это вопрос о надежности транспортной сети: чем выше искомая вероятность q, тем, разумеется, сеть надежнее.

Нетрудно видеть, что вопрос о надежности сети — это в свою очередь вопрос о связности случайного графа. В самом деле, сопоставим каждому городу вершину i € V10. Тогда «дорога» между «городами» i и j — это ребро. Износ дороги — это исчезновение ребра. Значит, утверждение «дорога изнашивается с вероятностью q» равносильно утверждению «ребро появляется с вероятностью p = 1 — q». Таким образом, нас интересует, какова минимальная вероятность р, при которой в модели Эрдеша-Реньи G(n,p) вероятность связности графа больше половины (граф, скорее, связен, чем несвязен).

Понятно, что если мы заменим число 10 другим числом, то соответствующее минимальное p может измениться. Этим и обусловлено наше желание рассматривать не только постоянные р, но и нетривиальные функции р = p(n).

В параграфе II.4 мы обсудим ряд строгих утверждений, касающихся сформулированного выше вопроса. Однако есть у нас и более срочное дело: все-таки предположение о том, что города связаны дорогами попарно, чересчур сильное, и в следующем параграфе мы приведем модификацию модели Эрдеша-Реньи, в рамках которой это предположение можно будет адекватно ослабить.

II.3. Обобщения модели Эрдеша-Реньи

Пусть по-прежнему Vn = {1, ..., n}. Однако теперь вероятность ребра между вершинами i и j мы обозначим через pj. Иными словами, мы, как и раньше, проводим ребра независимо друг от друга, но с разными вероятностями. В формате аксиоматики Колмогорова мы получаем вероятностное пространство

G(n,pij ) (^n,Fn,Pn,pij ^

в котором

Пп = {G = (Vn,E)}, Fn = , Pn,pij (G) =

= П pij • П (1 — pij).

(i,j)^E (i,j)^E

Важный частный случай описанного пространства получается, коль скоро мы фиксируем некоторый граф Hn = (Vn,En) и полагаем

p f (i,j) € En,

pij \0, (ijj) € En.

Иначе говоря, ребра графа Hn возникают в случайном графе независимо друг от друга с одной и той же вероятностью p = p(n) € [0,1], а ребра,

которых в графе Нп нет, не возникают в случайном графе вовсе. Этот вариант модели принято обозначать О(Нп,р). В ней

Рп (О) = р 1 Е 1(1 — р) )ЕП - I Е I.

Ясно, что модель О(Нп,р) и есть та самая модель, которая вполне адекватна вопросу о надежности транспортной сети. На сей раз мы не обязаны предполагать, что города попарно соединены дорогами; мы можем с самого начала зафиксировать граф дорог Нп и следить за износом его ребер.

ІІ.4. Некоторые математические результаты о надежности сети

Прежде всего справедлива следующая теорема Эрдеша-Реньи.

Теорема 1. Рассмотрим модель О(п,р). Пусть р = с п . Если с > 1, то почти всегда случайный граф связен. Если с < 1, то почти всегда случайный граф не является связным.

Этот довольно простой с точки зрения доказательства факт мы обоснуем в параграфе ІІ.5. Однако при всей своей формальной простоте теорема 1 несет весьма содержательную и в каком-то смысле неожиданную информацию. Действительно, вернемся к вопросу о надежности сети. Пусть число п городов, попарно соединенных дорогами, растет. Тогда, разумеется, величина р = с п довольно быстро стремится к нулю. Тем не менее теорема 1 утверждает, что вероятность сохранения связности графа при уничтожении его ребер с вероятностью ц = 1 — р стремится к единице. Грубо говоря, если городов 1000, то мы можем позволить дорогам разрушаться с вероятностью « 0,993, так что в результате с вероятностью, близкой к единице, перемещение между любыми двумя городами останется возможным. Поначалу это кажется противоречащим интуиции, но, при здравом размышлении, становится понятно, в чем здесь смысл. Дорог у нас Сп = О(п2), вероятность износа дороги равна 1 — 0(1пп/п) (мы пишем / = О(д) для функций ] и д, если существуют константы сі,с2 > 0, с которыми выполнено сід ^ ^ С2д). Значит, ожидаемое количество

неизношенных дорог имеет порядок п 1п п. Этого хватает для сохранения связности.

При определенной аккуратности в выкладках, которые мы частично проведем в параграфе ІІ.5, можно доказать, например, такой факт.

Теорема V. Рассмотрим модель О(п,р). Пусть р = с п . Если с > 3, то при п > 100

Рп,р{Є) ^ 1---•

п

Этот результат совсем замечателен своей конкретностью. Получается, что при той же тысяче городов и вероятности износа дороги

1 — 3 '"ост00 ~ ,98 вероятность сохранения связности не меньше, чем 0,999!

Теорема 1 любопытна еще и тем, что в ней наблюдается резкий скачок от «почти всегда связности» к «почти всегда несвязности». Функция р(п) = служит своего рода рубежом, преодоление которого означает переход от ненадежности к надежности. Такой переход принято называть фазовым, а соответствующую функцию р(п) — пороговой.

Следующая теорема содержит в себе еще более глубокую информацию о природе связности-надежности. Она была доказана самими Эрдешем и Реньи (см. [1-3]).

Теорема 2. Рассмотрим модель О(п,р). Пусть р = Если с < 1, то найдется такая константа в = в (с), что почти всегда размер каждой связной компоненты случайного графа не превосходит в 1п п. Если же с > 1, то найдется такая константа 7 = 7(с), что почти всегда в случайном графе есть ровно одна компонента размера ^ 7п.

И снова мы имеем фазовый переход — резкое изменение свойств случайного графа при преодолении некоторого порога. В данном случае порогом служит функция р = ^. Оказывается, что если вероятность ребра в с > 1 раз «ниже» порога, то все связные компоненты графа, скорее всего, крошечные — имеющие логарифмический от общего числа вершин размер; если же вероятность ребра в с > 1 раз «выше» порога, то, скорее всего, найдется компонента с числом вершин порядка п. Такая компонента называется гигантской.

Теорема 2 допускает различные уточнения. Например, можно утверждать, что при с > 1, помимо единственной гигантской компоненты, в случайном графе ничего сколь-нибудь крупного почти никогда не возникает: все остальные компоненты снова логарифмические. Можно еще аккуратнее описывать размер гигантской компоненты. В действительности, верно не только неравенство ^ 7п, но и асимптотика ~ 7п. А В.Е. Степанов доказал, что, опять же при р = с > 1, размер гигантской компоненты асимптотически нормален (см. [4-6]).

В целом изменение свойств случайного графа при изменении вероятности ребра р принято трактовать как эволюцию графа. Нам кажется, что несколько правильнее говорить о своего рода истории мира. Сначала (при р «С имеет место «феодализм» — весь граф поделен на несвязанные между собой логарифмические кусочки. Затем (при р > ^) возникает «империя» — гигантская компонента. Наконец, при р импе-

рия уничтожает «окраины» и добивается Мирового господства — связности.

В терминах надежности смысл теоремы 2 также очевиден: можно еще в 1п п раз уменьшить вероятность сохранности отдельной дороги, и, однако же, если не вся страна, то значительная ее

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

Глубокий интерес представляет, конечно, устройство мира «внутри» фазовых переходов, то есть при р ~ 1 и при р ~ . В первом случае все

совсем сложно, и мы отсылаем читателя к книгам [7-9]. Во втором случае можно сформулировать, например, следующий понятный результат.

Теорема 3. Пусть р = 1пп+с+°(1) _ Тогда

Рп,р

.

В частности, при р = вероятность стремится к в-1.

Здесь уже речь не идет о «почти всегда связности» или «почти всегда несвязности». Здесь асимптотическая вероятность связности есть, но она лежит в строгих пределах от нуля до единицы.

Все, о чем мы говорили до сих пор, касалось модели 0(п,р). Естественно, модель 0(Нп,р), будучи более адекватной реальности, является и более сложной для изучения. Главный результат относительно этой модели принадлежит Г.А. Маргу-лису (см. [11]).

Теорема 4. Пусть {Нп} — последовательность графов, реберная связность которых стремится к бесконечности при п ^ж>. Тогда существует пороговая, функция р для свойства связности случайного графа в модели 0(Нп,р). Иными словами, функция р такова,, что для любой функции р1 ^ С1р, где С1 < 1, случайный граф в модели 0(Нп,р\) почти всегда не связен, но для любой функции р2 ^ С2р, где С2 > 1, случайный гра,ф в модели 0(Нп,р2) почти всегда связен. Под реберной связностью графа понимается минимальное количество его ребер, удаление которых влечет потерю связности.

Теорема 4 нетривиальна и ее доказательство (а также массу ссылок на близкие результаты) можно найти в книге [8]. Разумеется, поиск пороговой функции, существование которой доказывается в теореме 4, — это всякий раз сложная задача, повязанная на специфику графов из последовательности {Нп}.

Практический смысл теоремы 4 банален: надо строить дороги так, чтобы связность получающегося графа неуклонно росла. К сожалению, в России положение, как правило, противоположное. Даже в Москве есть много улиц, перекрытие которых означает фактическую потерю связности. Например, таковы улицы, проходящие под железными дорогами: они крайне редки и служат единственными лазейками с одной стороны полотна на

другую.

В следующем параграфе мы докажем теорему

1, а в параграфе 11.6 мы обсудим основные идеи доказательства теоремы 2. Отметим, что дополнительную информацию о поведении случайных гра-

фов в модели Эрдеша-Реньи можно почерпнуть из книг [7-10].

11.5. Доказательство теоремы 1

Сперва обсудим случай с> 1.

Введем случайную величину на пространстве С(п,р):

Хп

Хп

0, если О связен, к, если у О ровно к компонент.

Таким образом, Хп принимает неотрицательные целые значения, причем Хп = 1. Нам нужно показать, что Рп,р(Хп = 0) ^ 1 при п ^ж>. Это равносильно асимптотике Рп,р(Хп ^ 1) ^ 0. По неравенству Чебышёва Рп,р(Хп ^ 1) ^ МХп, и нам остается обосновать стремление к нулю математического ожидания.

Представим Хп в виде суммы

ХП = ХП,1 + ... + Xn,n—1,

где Хп,к = Хп,к(О) — число к-вершинных компонент графа О. Занумеруем все к-элементные подмножества множества вершин Уп случайного графа в некотором (произвольном) порядке: Кі, ..., КСк. Тогда в свою очередь

Хп,к = Хп,к, 1 + ... + Хп,к,СП ,

коль скоро

Хп,к,і — Хп,к,і(О) —

1, если Кі образует компоненту в О,

0, иначе.

В итоге

-1 ск

МХп = ЕЕ МХп,к,і.

к=1і=1

Очевидно,

МХпХі = Рп,р(Кі образует компоненту в О) ^

^ Рпр(из Кі в Уп \ Кі нет ребер в О).

Получая последнее неравенство, мы просто пренебрегли условием связности той части графа О, которая «сидит» на множестве вершин Кі (такую часть принято называть индуцированным подграфом и обозначать О\к€). Далее,

Рпр(из Кі в У \ Кі нет ребер в О) = (1 — р)к(

п-к)

и, значит,

г-1 С к

1

мхп ^ — р)к(п-к) = Т, ск(1 — р)к(

п- к) .

к=1 і=1

к=1

Последняя сумма симметрична в том смысле, что ее слагаемые при к и п — к равны. Рассмотрим к = 1:

п(1 — р)п-1 < пе-р(п-1) = пе-

с

— Є

е

1 Ч с

- =«(!), п

поскольку С> 1.

Оставшаяся часть рассуждения состоит в доказательстве того, что слагаемые с к > 1 и к < п — 1 пренебрежимо малы по сравнению с первым слагаемым. Соответствующую выкладку мы пропустим. Если же поверить в ее справедливость, то получится, что вся сумма доминируется первым и последним слагаемыми, а стало быть, и она стремится к нулю.

Теорема 1 для случая с > 1 доказана.

Теперь рассмотрим случай с < 1. Обозначим через Хп количество изолированных вершин в случайном графе. Запишем

Хп = Хп,1 + ••• + Хп,п.

где

Хп,к = Хп,и(С) =

( 1, если вершина к € Уп изолированная в О, 1 0, иначе.

Тогда

МХп = МХп,1 +... + МХп,п.

В свою очередь

МХп,к = Рпр(кО) = (1 — р)п-\

Таким образом,

МХп = п(1 — р)п-1 = п(1 — р)п(1 + о(1)) =

= (1 + а(1))пв-с 1п п = (1 + а(1))п1-с.

Заметим, что ввиду неравенства с < 1 выполнено МХп ^ ж.

Посчитаем дисперсию случайной величины Хп:

ВХп = мхп — (МХп)2 = М (Хп , 1 +... + Хп , п)2 —

—(МХп)2 = мхп ,1 +... + мхп , п+

+ ^2 МХп,гХп,з — (МХп)2 = МХп,1 + ...

... + МХп,п + ^ МХп,гХп,э — (МХп)2 =

= МХп + ]Т МХп,гХп,з — (МХп)2.

Далее,

МХп^Хп^ = Рп,р(г и ] изолированы в О) = (1 — р)2п-1 =

= (1 + о(1))(1 — р)2п, то есть

]ГМХп,гХп,э = п(п — 1)(1 + о(1))(1 — р)2п =

= (1 + а(1))п2-2с = (1 + о(1))(МХп )2.

В итоге

ВХп = МХп + (1 + о(1))(МХп)2 — (МХп)2 = = о((МХп)2).

По неравенству Чебышёва

Рп,р(О) ^ Рп,р(Хп = 0) = Рп,р(Хп ^ 0) = = Рп,р( — Хп ^ 0) = Рп,р(МХп — Хп ^ МХп) ^ ^ РХп _ п,

. (МХпу о(

и вторая часть теоремы доказана.

11.6. Идеи доказательства теоремы 2

Метод, о котором мы будем здесь говорить, восходит к Р. Карпу (см. [12]), и в таком виде он описан в книге [9]. Мы лишь перечислим ниже основные шаги рассуждения.

Простейший ветвящийся процесс

Пусть Zl, .... Zt, ... — независимые пуассоновские величины с одним и тем же средним Л. Положим

Ус = 1,Уг = Уг-1 + Zt — 1.

Представлять себе описанный только что процесс можно так. В начальный момент времени есть одна частица. Затем она приносит Zl потомков и умирает. Заметим, что она может умереть, даже не принеся потомства, так как величина Zl равна нулю с положительной вероятностью. На следующем шаге все повторяется: какая-то частица (порядок роли не играет) порождает Z2 новых частиц, а сама гибнет. И так далее. Популяция может выродиться, а может и жить вечно. Хорошо известно, что имеют место следующие результаты.

Теорема 5. Пусть Л ^ 1. Тогда с веро-

ятностью 1 процесс Уг вырождается, то есть Р(3 I : У < 0) = 1.

Теорема 6. Пусть Л> 1. Пусть 7 € (0,1) — единственное решение уравнения 1 — 7 = .

Тогда процесс Уг вырождается в вероятностью 1 — 7, то есть Р(3 £ : Уг ^ 0) = 1 — ^.

Доказательства теорем 5 и 6 можно найти, например, в [9]. Впрочем, это стандартные факты теории ветвящихся процессов (см. [13]). Забегая вперед, скажем, что величина 7 в теореме 6 и одноименная величина в теореме 2 суть одно и то же. Вероятность вырождения процесса Уг и размер гигантской компоненты напрямую связаны.

Ветвящийся процесс на случайном графе

Пусть дан граф О = (У,Е): конкретный граф, не случайный. Зафиксируем какую-нибудь его вершину У1 Є У. Назовем ее «живой», а все остальные вершины — «нейтральными». Выберем среди нейтральных вершин всех соседей вершины «1. После этого объявим вершину V1 «мертвой», ее соседей — живыми, а все остальные вершины — нейтральными. Снова зафиксируем какую-нибудь живую вершину «2. Выберем всех ее соседей среди нейтральных. Вершину V2 отправим в царство мертвых, а в живых останутся все, кто был жив, кроме «2, и новорожденные «потомки» вершины V2 . Продолжая этот ветвящийся процесс, мы в конце концов получим кладбище вершин и ни одной живой вершины. Кладбище будет в аккурат компонентой, содержащей «1.

Обозначим число живых вершин в момент времени Ь через Уі, число нейтральных вершин — через N1, а число потомков очередной живой вершины, отправляющейся в последний путь, — через Zі. Тогда, очевидно,

Уо = 1,Уі = Уі-1 + Zt — 1.

Разумеется, все введенные величины зависят от графа О и от последовательности выбираемых вершин V1, ... Если граф О посчитать случайным, то при любом выборе вершин V1, ... получатся случайные величины Уі, N4, Zt на пространстве

О(п,р). На самом деле, ясно, конечно, что распределения этих величин не зависят от V1, ...; поэтому мы нигде зависимость от вершин и не указываем.

Сразу понятно, что Zt не является пуассонов-ской. Тем не менее она похожа на пуассоновскую. Дело в том, что она имеет биномиальное распределение с «числом испытаний» N1 и вероятностью «успеха» р. Правда, число испытаний само случайно. По счастью, это не проблема. Удается доказать, что Уг имеет вид

Уі = Іг + 1 — Ь, Іг - Віпот(п — 1,1 — (1 — р)г).

Подробности можно найти в [9], и мы их здесь не касаемся.

Как известно, биномиальное распределение сходится к пуассоновскому, коль скоро вероятность успеха обратно пропорциональна числу испытаний. Нечто подобное имеет место и у нас (р = с/п), и ровно на этом мы сыграем в итоге.

Случай с < 1

Положим Ьо = [в 1п п], где в = в (с) — константа, которую мы подберем позднее. Нам хочется доказать, что с большой вероятностью каждая из компонент случайного графа имеет размер не больше Ьо. Но размер компоненты — это момент вырождения процесса Уг на случайном графе. Значит, интересующее нас утверждение можно записать

в следующем виде:

Рп,р(3 VI : Уг0 > 0) ^ 0, п ^ ж. Поскольку

Рп,р(3 VI : Уг0 > 0) < пРп,р(Уг0 > 0), достаточно найти такое в, при котором

1

Дальнейшие выкладки будут слегка неаккуратными, но при желании их можно сделать безукоризненно строгими. Мы же хотим максимально прояснить основную суть подхода. Итак,

Рп,р(Уіо > 0) = Рп,р(Іі0 > Ьо) « ~ Рп,р(Віпот(п,1 — (1 — р)і0) > Ьо) «

~ Рп,р(Віпот(п,рго) > Ьо) «

ОС

1 х2

а/27Г

їр-прі0 у/прЬ0

Поскольку с < 1, нижний предел интегрирования имеет порядок л/^о- Стало быть, весь интеграл не превосходит величины е-6г°. Выберем в таким, чтобы е_<5(° оказалось меньше, чем е_21пп = А-,

’ П2 '

и в случае с < 1 теорема доказана.

Случай с > 1

Этот случай гораздо сложнее предыдущего. Здесь ветвящийся процесс на графе нужно «запускать» не один раз, а многократно. Только так удается доказать, что почти наверняка хотя бы в одном запуске возникнет гигантская компонента. Подробности можно найти в [9], мы же лишь поясним, откуда в текущей ситуации появляется константа 7 из формулировки теоремы и почему она совпадает с одноименной константой из теоремы 6.

Грубо говоря, идея следующая. Нам хочется доказать, что есть гигантская компонента. Тогда, как следствие, нам нужно, чтобы ветвящийся процесс на графе не вырождался даже при Ь « 7п. Иными словами, необходимо, чтобы

Рп,р(Уг ^ 0) ^ 0, Ь « 7п, п ^ж.

У нас р = ^. Значит, при і ~ ап выполнено

1 — (1 — р)г - 1 — е-рг - 1 — е-са.

Применим центральную предельную теорему к

Рп,р(Уг < 0) « Рп,р(Віпот(п,1 — е-са) < ап).

Интегрирование пойдет от минус бесконечности до

ап — п

у/п(1 - е-са)е-са'

Если а < 1 — в-са, то мы получим искомое стремление вероятности к нулю. Если а > 1—в-са, то вероятность, напротив, будет стремиться к единице. Таким образом, критическое значение а, вплоть до которого есть именно стремление к нулю, — это решение уравнения а = 1 — в-са или, что равносильно, 1 — а = в-са. А это и есть уравнение из теоремы 6, коль скоро Л мы заменяем на С.

III. Модели Барабаши—Альберт

В этом разделе мы поговорим о самых современных моделях случайных графов, которые призваны описывать рост различных сетей — социальных, биологических, транспортных. Но в первую очередь речь пойдет об Интернете. В 90-е годы ХХ века, когда интернет только зарождался, исследователи уже задались вопросом, каким законам подчиняется рост интернета и какова наиболее адекватная модель для описания свойств этой сети. Одними из первых здесь были А.-Л. Бараба-ши и Р. Альберт. Они нашли ряд важных эмпирических закономерностей в поведении Интернета и на их основе придумали модель, которую впоследствии по-разному формализовывали многие авторы. Соответственно мы построим наше изложение следующим образом. В первом параграфе мы расскажем о результатах Барабаши-Альберт. Во втором параграфе мы опишем модель Б. Боллобаша и О. Риордана, которая весьма неплохо ложится на статистики Барабаши-Альберт. В третьем параграфе мы обсудим возможные уточнения модели Боллобаша-Риордана.

111.1. Наблюдения Барабаши—Альберт

В своих работах [14-16] Барабаши и Альберт, а также Х. Джеонг описали те статистики Интернета, которые легли в основу науки о росте этой сети — науки, имеющей глубокие приложения как в собственно проблематике Интернета, так и в многочисленных близких дисциплинах. В действительности, большинство реальных сетей (социальные, биологические, транспортные и пр.) имеют похожую «топологию».

Итак, сперва договоримся о том, что мы понимаем под сетью Интернет. Это так называемый веб-граф, вершины которого суть какие-либо конкретные структурные единицы в Интернете: речь может идти о страницах, сайтах, хостах, владельцах и пр. Для определенности будем считать, что вершинами веб-графа служат именно сайты. Ребрами же мы будем соединять те вершины, между которыми имеются ссылки. При этом разумно проводить столько ребер между двумя вершинами, сколько есть ссылок между соответствующими сайтами. Более того, ребра естественно считать направленными. Таким образом, веб-граф ориентирован и он может иметь кратные ребра, петли и даже кратные петли (ссылки вполне могут идти с одной страницы данного сайта на дру-

гую его страницу). Это такой «псевдомультиорграф». Сразу понятно, что для подобного «зверя» модель Эрдеша-Реньи вряд ли подходит.

Теперь мы готовы перечислить самые основные моменты исследования Барабаши-Альберт. По существу, этих моментов всего три. Во-первых, веб-граф — это весьма разреженный граф. У него на Ь вершинах примерно Ы ребер, где к ^ 1 — некоторая константа. Для сравнения, у полного графа на Ь вершинах С^ = О(Ь2) ребер. Однако — и это, во-вторых, — диаметр веб-графа исключительно скромен. В 1999 году он имел величину (см. [16]) 5--7. Это хорошо всем известное свойство любой социальной сети, которое принято в обыденной речи характеризовать выражением «мир тесен». Например, говорят о том, что любые два человека в мире «знакомы через 5-- 6 рукопожатий». Точно так же и сайты: «кликая» по ссылкам, можно с любого сайта на любой другой перейти за 5-- 7 нажатий клавиши компьютерной мыши. Конечно, тут есть важная оговорка. Некоторые едва появившиеся сайты могут не быть связаны с внешним по отношению к ним миром. Несколько правильнее сказать, что в веб-графе есть гигантская компонента и уже ее диаметр невелик. Таким образом, веб-граф очень специфичен: будучи разреженным, он тем не менее в известном смысле тесен.

В-третьих, у веб-графа весьма характерное распределение степеней вершин. Эмпирическая вероятность того, что вершина веб-графа имеет степень й, оценивается как с/3^, где Л « 2,1, а С — нормирующий множитель, вычисляемый из условия «сумма вероятностей равна 1». Этот любопытный факт роднит Интернет с очень многими реальными сетями — биологическими, социальными, транспортными. Все они подчиняются степенному закону, только у каждой из них свой показатель Л.

Ввиду перечисленных наблюдений не остается никаких сомнений в том, что модель Эр-деша-Реньи не применима для описания роста Интернета и подобных сетей. Если подбором вероятности р еще можно добиться разреженности и тесноты (хотя и не с теми параметрами), то степенной закон совсем уж не имеет отношения к схеме Бернулли, в рамках которой появляются ребра обычного случайного графа. В модели О(п,р) степень каждой вершины случайного графа биномиальна с параметрами п — 1 и р, и при тех р, которые мало-мальски гарантируют разреженность (то есть при р = О(1/п)), указанное биномиальное распределение аппроксимируется пуассоновским, а вовсе не степенным.

Сами Барабаши и Альберт предложили очень разумный взгляд на процесс формирования интернета. Давайте считать, сказали они, что в каждый момент времени появляется новый сайт, и этот сайт ставит фиксированное количество ссылок на своих предшественников. На кого он предпочтет

сослаться? Наверное, на тех, кто и так уже популярен. Можно допустить, что вероятность, с которой новый сайт поставит ссылку на один из прежних сайтов, пропорциональна числу уже имевшихся на тот сайт ссылок.

Модели случайных графов, основанные на описанной идее, называются моделями предпочтительного присоединения. В своих работах Бараба-ши и Альберт никак не конкретизировали, какую именно из этих моделей они предлагают рассматривать. А эти модели исключительно разнородны по своим свойствам. Ведь можно ставить ссылки независимо друг от друга, а можно еще и зависимости между разными ссылками с одного сайта учитывать. В итоге удается доказать даже такой забавный факт (см. [17]).

Теорема Т. Пусть f (n), n ^ 2, — произвольная целочисленная функция, такая, что f (2) = О, f (n) ^ f (n +1) ^ f (n) + 1 для всех n ^ 2 и f (n) ^ ж при n ^ ж. Тогда существует такая модель типа Барабаши-Альберт, что в ней с вероятностью, стремящейся к единице при n ^ ж, случайный граф содержит в точности f (n) треугольников.

Одну из наиболее правильных спецификаций модели Барабаши-Альберт предложили в начале двухтысячных годов Б. Боллобаш и О. Риордан. В следующем параграфе мы ее обсудим.

III.2. Модель Боллобаша-Риордана

Наиболее полно эта модель описана в книге [В] и обзоре [17]. Также имеется малодоступная книга [18]. Мы представим здесь две основных и, по сути, совпадающих модификации этой модели. В одной дается динамическое, а в другой статическое описание случайности. Интуитивно более понятна динамическая модификация, с нее и начнем.

Динамическая модификация

Сперва построим последовательность (случайных) графов {Gk}, в которой у графа с номером n число вершин и ребер равно n. Затем сделаем из нее последовательность {GП}, в которой у графа с номером n число вершин равно n, а число ребер равно kn, k Є N.

Итак, пусть Gi = ({1},{(1,1)}), то есть в начальный момент времени есть одна вершина и одна петля. Пусть теперь граф GП-і уже построен. У него вершины образуют множество {1, ..., n-1}, а ребер у него тоже n — 1 штука. Добавим вершину n и ребро (n,i), у которого i Є {1, ..., n}. Ребро (n,n) будет появляться с вероятностью 9п1_1 ; ребро (n,i) возникнет с вероятностью 9^;е°'\, где deg* — степень вершины i в графе Gn-1. Очевидно, что распределение вероятностей задано корректно, поскольку

n—1 degi 1 2n — 2 І

> -----5---1------=----------1------= 1.

^ 2n - 1 2n - 1 2n - 1 2n - 1

i=

Случайный граф Gn построен, и он удовлетворяет принципу предпочтительного присоединения.

Осталось перейти к Gn. Берем G^. Это граф с kn вершинами и kn ребрами. Делим множество его вершин на последовательные куски размера k:

{1, ..., k}, {k + 1, ..., 2k}, ..., {k(n — 1) + 1, ..., kn}.

Объявляем каждый кусок «вершиной», а ребра сохраняем, то есть если были ребра внутри куска, то будут кратные петли, а были ребра между двумя различными кусками — будут кратные ребра. Внешне — вполне себе Интернет, как мы его и представляли. Вершин стало n, а ребер — по-прежнему kn. Цель реализована.

Статическая модификация, или LCD-модель

Введем такой объект, который называется линейной хордовой диаграммой. Вообще-то, он возник в топологии и теории узлов (см., например, [19]), но его комбинаторика оказывается напрямую связана с формированием веб-графа.

Итак, зафиксируем на оси абсцисс на плоскости 2n точек: 1, 2, 3, ..., 2n. Разобьем эти точки на пары и элементы каждой пары соединим дугой, лежащей в верхней полуплоскости. Полученный объект назовем линейной хордовой диаграммой (linearized chord diagram или, короче, LCD по-английски). Дуги в нем могут пересекаться, лежать друг под дружкой, но не могут иметь общих вершин. Количество различных LCD легко считается, оно равно

_ (2??.)! 2пп\ '

По каждой LCD построим граф с n вершинами и n ребрами. Действуем так. Идем слева направо по оси абсцисс, пока не встретим впервые правый конец какой-либо дуги. Пусть этот конец имеет номер i . Объявляем набор {І, ... , i } первой вершиной будущего графа. Снова идем от

i +І направо до первого правого конца i2 какой-либо дуги. Обявляем второй вершиной графа набор {iі + 1, ..., i2}. И так далее. Поскольку правых концов у дуг в данной диаграмме n штук, получаем всего n вершин. А ребра порождаем дугами. Иными словами, две вершины соединяем ребром, коль скоро между соответствующими наборами есть дуга. Ребра ориентируем справа налево. Аналогично возникают петли. Дуг n и ребер n.

Теперь считаем LCD случайной, то есть полагаем вероятность каждой диаграммы равной 1/ln. Возникают случайные графы. Можно показать, что такие графы по своим вероятностным характеристикам практически неотличимы от графов Gn.

Графы с n вершинами и kn ребрами получаем тем же способом, что и в предыдущем пункте.

Некоторые результаты

Замечательна модель Боллобаша-Риордана не только тем, что с ее помощью наводится порядок в «каше», которую «заварили» Барабаши и Альберт, но еще и тем, что она полностью адекватна эмпирическим наблюдениям. Прежде всего справедлива

Теорема S. Для любого k ^ 2 и любого е > 0

P

ln n ln ln n

ln n ln ln n

На первый взгляд, утверждение кажется непонятным. Ну, хорошо: диаметр плотно скон-

центрирован (по вероятности) около величины 1пп/ 1п1пп. А у нас ведь какие-то 5-7 были? Так ничего странного. Вершин в Интернете образца 1999 года около 107. Значит,

7ln10

ln 107

lnln107 ln7 + lnln10

6.

Фантастическое попадание. Отметим, что при недавней проверке с другими цифрами эмпирика снова подтвердилась.

Теорема 8 доказана в работе [20] авторами модели. А в работе [21] была внесена ясность и в вопрос о распределении степеней вершин.

Теорема 9. Для любого к ^ 1 и любого ! < п1/1Ъ

M

|{i = 1, ..., n і degGni = d}|

2 k

Поскольку к — константа, выражение в правой части (2) имеет вид const/d3. Да это же в точности степенной закон! Правда, в формулировке теоремы написано математическое ожидание, а не вероятность, но одно из другого получается за счет мартингальных неравенств и соответствующих теорем о плотной концентрации меры около среднего (см. [21]).

У теоремы 9 есть все же два неприятных момента. Первый состоит в том, что степень d в степенном законе, который в ней устанавливается, равна не 2,1, а 3. Второй — это ограничение d ^ , которое ставит крест на практической

применимости теоремы. Даже при n « 1012, чего в природе (пока) не бывает, мы имеем лишь d ^ 104/'\ и это нелепо.

Последний недостаток недавно устранил Е.А. Гречников — исследователь-разработчик в Яндексе, который получил более точный результат практически без ограничений на d. Статья Гречникова еще не опубликована.

Первым же недостатком занимались много и, в частности, предлагали различные альтернативные модели. Одну из таких моделей мы обсудим в параграфе ІІІ.3. Но прежде скажем еще несколько слов о свойствах LCD-модели.

Пусть Н — фиксированный граф. Обозначим через §(Н,ОП) случайную величину, равную количеству подграфов графа ОП, изоморфных графу

Н. Как распределена эта величина? Изучали ее математическое ожидание в разных специальных случаях. Например, в работе [17] приводится громоздкая общая формула и пара ее симпатичных следствий, которые мы выпишем и здесь.

Теорема 10. Пусть к ^ 2. Пусть также Кз — полный граф на трех вершинах. Тогда

m тз ,Gn)) =

k  (k - 1) • (k+ 1) 48

(ln n)3

при n ^ ж.

Теорема ll. Пусть фиксированы k ^ 2 и

l ^ 3. Пусть также Сі — цикл на l вершинах. Тогда

M($(Сі^П)) = (і + o(1))  Ck,i  (lnn)1

при n ^ ж, где Ck,i — это положительная константа. Более того, при k ^ ж имеем Ck,i = ©(kl).

Студенты МФТИ А. Рябченко и Е. Самосват недавно (в несколько иной, но очень близкой модели) установили следующий общий факт.

Теорема l2. Пусть задан граф H, степени вершин которого равны di, ..., ds. Обозначим через §(di = m) число вершин в H, степень каждой из которых равна т. Тогда

m{${h,gd) = 0(j7^ds=o)-(v^:)#^=1)-(inj7O#№=2)).

Зависимость от k занесена в константу ©.

Надо полагать, что нечто подобное было известно и авторам статьи [17], но мы ничего похожего в литературе не встречали. А такая запись результата очень удобна. Скажем, в теореме 10 речь идет про K3 . Ясно, что для K3 выполнено

^(di = 0) = $(di = 1)=0, $(di = 2) = 3.

По теореме 12

M ($(K3 ,Gn)) = e((ln n)3),

и это прекрасно согласуется теоремой 10. Аналогично можно разобраться и с циклами (теорема 11). А если взять K4 — полный граф на четырех вершинах, — то теорема 12 скажет, что средняя его встречаемость в веб-графе постоянна. Иными словами, «тетраэдров» в веб-графах почти не бывает.

Отметим, что в реальном вебе случаются не только тетраэдры, но и клики куда большей мощности. Это связано с деятельностью спамеров, которые искусственно расставляют ссылки, желая повысить рейтинги сайтов, заплативших за раскрутку. Спам в модели Боллобаша-Риордана не учтен, и это тоже минус.

n

III.3. Модель копирования

Здесь мы опишем еще одну очень интересную модель, которая также призвана объяснить феномен степенного закона в реальных сетях. Эта модель возникла практически в одно время с моделью Барабаши-Альберт. Она принадлежит Р. Ку-мару, П. Рагхавану, С. Раджагопалану, Д. Сива-кумару, А. Томкинсу и Э. Упфалу (см. [22]).

Фиксируем a Є (0,1) и d ^ 1, d Є N. Случайный граф будет расти, и это будет похоже на процесс из пункта III.2. Однако здесь процесс будет устроен сильно по-другому.

В качестве начального графа возьмем любой d-регулярный граф (граф, у которого степень каждой вершины равна d). Пусть построен граф с номером t. Обозначим его Gt = (Vt,Et). Здесь Vt = {ui, ..., us}, где s отличается от t на число вершин начального графа, то есть на некоторую константу, выражаемую через d. Добавим к Gt одну новую вершину us+1 и d ребер, выходящих из us+i. Для этого сперва выберем случайную вершину p Є Vt (все вершины в Vt равновероятны). Одно за другим строим ребра из us+i в Vt. На шаге с номером i, i Є {1, ..., d}, разыгрываем случайную величину, которая с вероятностью a принимает значение 1 («монетка падает решкой кверху») и с вероятностью 1 — a принимает значение 0 («монетка падает орлом кверху»). Если вышла единица, то выпускаем ребро из us+i в случайную вершину из Vt (все вершины в Vt равновероятны). Если вышел ноль, то берем i-го по номеру соседа вершины p. Последнее действие всегда возможно, так как по построению у каждой вершины не менее d соседей.

Интуиция за всем этим примерно такая. Появляется новый сайт. Проставляя очередную ссылку, его владелец с некоторой вероятностью будет ориентироваться на кого-то из своих предшественников. Скажем, сайт посвящен автомобилям. Вероятно, владелец возьмет один из уже существовавших сайтов про автомобили и скопирует оттуда ссылку (с точки зрения стороннего наблюдателя, вполне случайную). Это ситуация, когда монетка выпала орлом кверху (p — это сайт, с которого копируются ссылки). Однако при простановке ссылки владелец может и никого не копировать, а случайно (по нашему мнению) цитировать кого-то из предшественников. Это случай выпадения решки. Таким образом, І — a — это вероятность копирования или, если угодно, вероятность выбора, мотивированного тематикой сайта.

Основной результат из [22] — это теорема 13.

Теорема l3. Пусть Nt,r — это математическое ожидание числа вершин степени r в графе Gt. Тогда

г Nt,r

lim ——

t^O t

I — - '

0 Г

Пафос теоремы в том, что в ней мы снова приходим к степенному закону. Более того, если вероятность копирования близка к 1 (а величина а — к нулю), то показатель степени может равняться ожидаемой величине 2,1, чего до сих пор у нас не было.

В целом, распределение степеней вершин в модели копирования очень похоже на распределение степеней вершин в модели Боллобаша-Риор-дана. В остальном, модели сильно разнятся. Например, в модели Боллобаша-Риордана практически отсутствуют плотные двудольные подграфы (см. теорему 12); в модели копирования таких подграфов полно. Это особенно важно ввиду того, что спамерские структуры, о которых мы вскользь говорили в конце пункта III.2, зачастую образуют именно двудольные графы с плотной пе-релинковкой.

Работа выполнена при финансовой поддержке гранта РФФИ 09-01-00294, гранта Президента РФ МД-8390.2010.1, гранта поддержки ведущих научных школ НШ-8784.2010.1.

Литература

1. Erdos P., Renyi A. On random graphs

I // Publ. Math. Debrecen. — 1959. — V. 6. — P. 290-297.

2. Erdos P., Renyi A. On the evolution of random graphs. // Publ. Math. Inst. Hungar. Acad. Sci. — 1960. — V. 5. — P. 17-61.

3. Erdos P., Renyi A. On the evolution of random graphs // Bull. Inst. Int. Statist. — Tokyo, 1961. — V. 38. — P. 343--347.

4. Степанов В. Е О вероятности связности случайного графа gm(t) // Теория вероятностей и ее применения. — 1970. — Т. 15, вып. 1. — С. 55-67.

5. Степанов В.Е. Фазовый переход в случайных графах // Теория вероятностей и ее применения. — 1970. — Т. 15, вып. 2. — C. 187-203.

6. Степанов В.Е. Структура случайных графов gn(x\h) // Теория вероятностей и ее применения. — 1972. — Т. 17, вып. 3. — С. 227-242.

7. Колчин В.Ф. Случайные графы. — М.: Физ-матлит, 2004. — 256 c.

8. Bollobas B. Random Graphs. — Cambridge Univ. Press, 2001. — 520 c.

9. Алон Н., Спенсер Дж. Вероятностный метод. — М: Бином. Лаборатория знаний, 2007. — 320 c.

10. Janson S., Luczak T, Rucinski A. Random graphs. — New York: Wiley, 2000. — 333 с.

11. Маргулис Г.А. Вероятностные характеристики графов с большой связностью // Проблемы передачи информации. — 1974. — Т. 10. — С. 101-108.

12. Karp R. The transitive closure of a random digraph // Random structures and algorithms. — 1990. — V. 1. — P. 73-94.

13. Карлин С. Основы теории случайных процессов. — М: Мир, 1971. — 536 с.

14. BaraЪasi L.-A., Albert R. Emergence of

scaling in random networks jj Science. — 1999. — V. 2В6. — P. 509--512.

15. BaraЪasi L.-A., Albert R., Jeong H. Scale-

free characteristics of random networks: the topology of the world-wide web jj Physica. — 2000. —

V. A281. — P. 69-77.

16. Albert R., Jeong H., BaraЪasi L.A. Diameter

of the world-wide web jj Nature. — 1999. —

V. 401. — P. 130-131.

1Т. Bol^as B., Riordan O. Mathematical

results on scale-free random graphs jj Handbook of graphs and networks. — Weinheim: Wiley-VCH,

2003. — P. 1-34.

1S. Райгородский А.М. Экстремальные задачи теории графов и анализ данных. — М: Регулярная и хаотическая динамика, 2009. — 120 с.

19. Stoimenow A. Enumeration of chord diagrams and an upper bound for Vassiliev invariants // J. Knot Theory Ramifications. — 1998. — V. 7, N 1. — P. 93-114.

20. Bollobas B, Riordan O. The diameter of a scale-free random graph // Combinatorica. —

2004. — V. 24, N 1. — P. 5-34.

21. Bollobas B, Riordan O., Spencer J., Tusnady G. The degree sequence of a scale-free random graph process // Random Structures Algorithms. — 2001. — V. 18, N 3. — P. 279-290.

22. Kumar R., Raghavan P., Rajagopalan S., Sivakumar D., Tomkins A., Upfal E. Stochastic models for the web graph // Proc. 41st Symposium on Foundations of Computer Science. — 2000.

Поступила в редакцию 15.10.2010.