О типичной структуре конфигурационного интернет-графа с известным числом связей Текст научной статьи по специальности «Математика»

Труды Карельского научного центра РАН № 5. 2011. С. 86-96

УДК 519.2

О ТИПИЧНОЙ СТРУКТУРЕ КОНФИГУРАЦИОННОГО ИНТЕРНЕТ-ГРАФА С ИЗВЕСТНЫМ ЧИСЛОМ СВЯЗЕЙ

Ю. Л. Павлов

Институт прикладных математических исследовании Карельского научного центра РАН

Рассматривается конфигурационная модель случайного графа с N вершинами. Степени вершин выбираются независимо из степенного закона распределения с показателем т G (1, 2) при условии, что сумма степеней вершин равна n. Показано, что если N, n ^ œ так, что п — Z(т)N = O(N1/т), где Z(x) - дзета-функция Римана, то предельная структура графа не отличается от типичной структуры случайного графа без условия на сумму степеней.

Ключевые слова: случайные графы, Интернет, конфигурационная модель, случайная структура.

Yu. L. Pavlov. ON THE TYPICAL STRUCTURE OF CONFIGURATION INTERNET GRAPH WITH KNOWN NUMBER OF LINKS

We consider the configuration model of random graph consisting of N vertices. The degrees of vertices are drawn independently from power-law distribution with the exponent т G (1,2) under the condition that the sum of vertex degree is equal to n. We show that if N, n ^ œ in such a way that n — Z(т)N = O(N 1/t), where Z(x) is the Rimann‘s zeta-function, then the limit structure of the graph is the same typical structure as the random graph without the condition on the sum of degrees.

Key words: random graph, Internet, configuration model, random structure.

Введение

За последние 10-15 лет специалисты по случайным графам уделяют значительное внимание разработке моделей сложных телекоммуникационных и социальных сетей, в частности, сети Интернет (см. например, [Durrett, 2007]). В настоящей работе мы рассмотрим одну из наиболее известных таких моделей, получившую название конфигурационной модели со случайными степенями (см. [Reittu, Norros, 2004; Durrett, 2007; Hofstad,

2011]). Пусть граф содержит N вершин, занумерованных числами от 1 до N. Степени вершин являются независимыми одинаково распределенными случайными величинами Пі,..., Пм такими, что

Р{Пі ^ к} = к-Т, (1)

где т - положительный параметр. Если сумма ^м = Пі + • • • + Пм нечетна, то в граф вводится дополнительная вершина с единичной степенью. Формирование графа происходит в два этапа. На первом этапе определяются степени

0

вершин в соответствии с распределением (1), при этом считается, что из вершин выходят «полуребра» [Ив^и, N01X08, 2004], т. е. ребра, для которых смежные вершины еще не определены (все полуребра графа различны). На втором этапе полуребра соединяются равновероятно для образования ребер.

В ряде работ (см. [Ра10^808 et а1., 1999; Ив^и, N01X08, 2004; Биггвй, 2007]) показано, что типичные для многих реальных сетей значения параметра т распределения (1) принадлежат интервалу (1,2). Поскольку такие значения характерны и для сети Интернет, соответствующие графы иногда называются графами Интернет-типа или Интернет-графами.

Предельное поведение структуры и различных числовых характеристик Интернет-графов изучалось во многих публикациях (см. [Ив^и, N01X08, 2004; Биггвй, 2007] и ссылки в них). Обнаружено, что граф имеет единственную гигантскую компоненту связности, число вершины которой пропорционально N, а объемы других компонент бесконечно малы по сравнению с N. Найдена асимптотика математического ожидания объема гигантской компоненты и получена оценка сверху ее диаметра. Предельное распределение объема гигантской компоненты получено в [Павлов, 2007].

В графах Интернет-типа сумма степеней вершин является случайной величиной. Однако представляют интерес и условные графы такого вида при условии, что величина известна и равна п. Такие графы полезны для моделирования сетей, в которых число связей можно оценить. Кроме того, результаты об условных графах можно использовать для нахождения предельных распределений числовых характеристик Интернет-графов без ограничений на путем усреднения условных распределений этих характеристик по распределению общего числа ребер графа. Изучение свойств условных случайных графов Интернет-типа начато в статьях [Павлов, Чеплюкова, 2008; Павлов, 2009]. Основным методом исследования в этих работах была обобщенная схема размещения частиц по ячейкам, введенная и подробно изученная В. Ф. Колчи-ным (см., например, [Колчин,, 2000]).

Согласно (1),

Рк = Р(Ш = к} = — т,—, к = 1,2,...

{ 11 } кт (к + 1)т’

(2)

Отсюда следует, как легко видеть, что

Е П1 = С(т ^ (3)

где Z (т) — значение дзета-функции Римана в точке т. Поэтому при изучении динамики структуры Интернет-графов типичным является случай N, n ^ ж так, что N/n ^ Z (т). Естественно ожидать, что в этом случае конфигурация условного графа асимптотически эквивалентна структуре графа без ограничений на сумму vn и представляет интерес вопрос о том, при каких значениях n такая ситуация сохраняется. В данной статье показано, что это имеет место при условии n — Z(т)N = O(N1/т). По-видимому, такое утверждение сохраняет силу и в более широкой зоне изменения параметров, однако этот вопрос требует дальнейшего исследования. Идеи, на которых основаны приведенные ниже результаты, близки к изложенным в работе [Reittu, Norros, 2004], однако необходимость учета условия vn = n потребовала изменения техники доказательств. Это стало возможным после получения предельных распределений максимальной степени и числа вершин заданной степени в условных Интернет-графах [Павлов, Чеплюкова, 2008; Павлов, 2009].

Основной РЕЗУЛЬТАТ

Введем следующее распределение

qk = (k + 1)pfc+iC-1(т ), k = 0,1,2,... (4)

Пусть <^(z) и F(z) означают производящие функции распределений (2) и (4) соответственно, следовательно,

<p(z) = ^ Pkzk, F(z) = ^ qkzk (5)

k=l

Обозначим также:

k=0

a =

lnln N — ln(r — 1)

exp

2w 2 - t

(6)

где [ж] означает целую часть числа ж, а ад = ад^) - медленно растущая функция [Ибрагимов, Линник, 1965], такая, что и> > 0 и

ад = о(1п1п1п N), 1п1п1п N = о(ад). (7)

Обозначим АN некоторое событие, зависящее от N. Будем говорить, что событие AN асимптотически достоверно (а. д.), если P{AN} ^ 1 при N ^ то. Автор выражает благодарность А. М. Зубкову за обсуждение этого определе-

ния.

Справедливо следующее утверждение.

Теорема. Пусть N п ^ то так, что п — ((т)Ж = 0(Ж1/т). Тогда граф содержит компоненту связности объема V, для которой

87

Е У/Ж ^ 1 - р(д), где д - вероятность вырождения ветвящегося процесса Гальтона-Ватсона с распределением (4) числа прямых потомков каждой частицы. Диаметр этой компоненты а.д. не превосходит 2(а + в).

Ниже приводятся леммы 1-6, с помощью которых в конце статьи будет доказана сформулированная теорема.

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

Связь между подмножествами вершин

Пусть ] - произвольная вершина графа, имеющая степень й = й(Ж, п). Рассмотрим некоторое множество вершин Ш, сумма степеней которых равна О = О(Ж, п). Пусть ^ Ш, й ^ О. Обозначим А событие, состоящее в том, что вершина ] соединена ребром с какой-нибудь вершиной множества Ш.

Лемма 1. Пусть N п ^ то п — £(т)Ж = 0(Ж1/т). Тогда справедливы следующие утверждения.

1. Если йО/п ^ то, то Р{А} ^ 1.

2. Если йО/п ^ 0, то Р{А} ^ 0.

Доказательство. Мы используем идею доказательства леммы 3.8 статьи [Ие1йи, Коггоз, 2004]. Пусть выполнено условие 1. Нетрудно видеть, что из множества й ребер, выходящих из вершины _/, можно выбрать подмножество из (И ребер, удовлетворяющую в условиях леммы следующим соотношениям:

в! ^ й, ^О/п ^ то, й/2/п ^ 0. (8)

Обозначим А1 событие, состоящее в том, что все й/ ребер не идут в Ш и не образуют петель, а через А2 событие, состоящее в наличии петель, образованных какими-то из й/ ребер. Нетрудно видеть, что

Р{А} < Р{А1} + Р{А2} <

(п — й/ — О)(п — ё — О — 1)...(п — 2й/ — О + 1) (п — 1)(п — 3)...(п — 2й/ + 1)

. (п — й/ )(п — й! — 1)...(п — 2й/ + 1)

+ (п — 1)(п — 3)... (п — 2й/ + 1) <

( п — й/ — О У' (п — 2й/ + 1У _п ( й/2 \

\п — 2й/ + 1у + у п — 1 ) \ п /

Отсюда и из (8) следует первое утверждение. Второе доказывается аналогично:

Р{А} =

(п — й — Б)(п — й — Б — 1)...(п — 2й — Б + 1) (п — 1)(п — 3)...(п — 2й + 1)

п — 2й — Б + 1\гі п — 2й + 1 у ^ .

□

Алгоритм построения графа

Рассмотрим один из возможных алгоритмов построения графа на втором этапе его формирования. Заметим, что для безусловных графов аналогичный алгоритм описан в [Иеійи, Копте, 2004] (см. также [Павлов, 2007]). Идея алгоритма следующая. Граф строится с начальной вершины, например, с вершины 1. Далее для всех полуребер начальной вершины по очереди и равновероятно ищутся пары и, таким образом, образуются ребра, инцидентные первой вершине. Вершины, соединенные ребром с начальной (и не совпадающие с ней), называются вершинами первого поколения. Далее для всех свободных по-луребер вершин первого поколения опять последовательно и равновероятно ищутся пары и образуются новые ребра. В результате вершины, соединенные путями длины 2 с начальной вершиной, образуют второе поколение. Этот процесс продолжается до тех пор, пока у вершин очередного поколения не окажется свободных полуребер. Обозначим ^ число вершин ¿-го поколения, £ = 0,1, 2,.... Очевидно, что = 1. Если = 0, а построенный граф еще не содержит все N вершин, то завершено построение только одной компоненты связности и для продолжения формирования графа среди всех не рассмотренных вершин выбирается вершина с наименьшим номером в качестве начальной для построения следующей компоненты связности, и далее процедура осуществляется аналогично. Повторив, если нужно, эти действия несколько раз до исчерпания всех N вершин, мы построим одну из возможных реализаций графа. Нетрудно видеть, что любая такая реализация возникает с вероятностью, равной соответствующей величине, индуцируемой вероятностным пространством, заданным на множестве рассматриваемых графов. Ниже приведем алгоритм построения одной компоненты связности. По ходу работы этого алгоритма составляется последовательность номеров вершин, которые одна за другой присоединяются к графу по мере об. разования ребер: ^,^2,...,І«, где ^ -номер

начальной частицы, ^ - номер первой вершины, соединенной ребром с ^1, а js-номер последней вершины последнего поколения данной компоненты связности. Процесс образования такой последовательности номеров подобен естественной нумерации частиц в ветвящихся процессах Гальтона - Ватсона (так называемая «история семейств», см. [Ватутин, 1993]). Отличие состоит в том, что номера частиц процесса не повторяются, а в последовательности ^1,^2, ...,^8 номера повторяться могут в силу возможности образования петель, циклов и кратных ребер.

В тексте алгоритма мы будем рассматривать вершины как N различных ящиков, содержащих, соответственно, kl,...,kN различных шаров (^ + ... + kN = п). Понятно, что эти шары обозначают полуребра вершин. Удаление двух шаров подряд подразумевает образование очередного ребра графа путем соединения соответствующих полуребер. В этом алгоритме и означает последовательность номеров присоединяемых вершин.

Алгоритм

1. и := {л}, г = 1, k = 2.

2. Удаляется первый шар из ящика ji.

3. Выбирается равновероятно шар из оставшихся и удаляется, номер его ящика обозначается .

4. и := и и {л}.

5. k = k + 1.

6. Если в ji есть еще шары, то переход на шаг 2.

7. г = г + 1.

8. Если в ящиках из и есть еще шары, то переход на шаг 6.

9. Стоп.

ДрЕвовидность компоненты связности

Далее мы покажем, что если Ж, п ^ то, то при больших Ь (и даже при Ь ^ то, если это стремление происходит достаточно медленно) построенный по приведенному алгоритму подграф до поколения Ь включительно является деревом, степени вершин которого сравнительно невелики. Следующие леммы 2 и 3 являются развитием соответствующих результатов работ [Ие^и, Коггоз, 2004] и [Павлов, 2007]. Обозначим А* = А*(Н), Ь = 0,1,..., событие, состоящее в том, что такой подграф является деревом и в нем нет вершин, степень которых больше или равна Н, где Н = Н(Ж) -медленно растущая функция такая, что Н = 0(1п Ж).

Лемма 2. Пусть N, n ^ то так, что n — Z(t)N = O(N1/т), t = o(lnlnh). Тогда P{At} = 1 + O(h—(т-1)t).

Доказательство. Нетрудно видеть, что P{Ao} = P{ni < h|vN = n}

= 1 — (P{vn = n})-1 ^pk P{vn-i = n — k}.

k^h

(9)

Согласно лемме 3 из [Павлов, 2007] при N ^ то

sup |N1/т P{vn = k} —f ((k—NZ(t))/N1/т)| ^ 0, k

(10)

где f(x) - плотность устойчивого распределения с показателем т и характеристической функцией

exp

Г(2 — т) т- 1 1

t пт\ пт

1 +1Щ ^т]с<кТ г

Учитывая одновершинность устойчивых распределений [Ибрагимов, Линник, 1965] и (1), из (9) находим, что

Р{Ао} = 1 + 0(Н-т),

следовательно, для Ь = 0 утверждение леммы 2 доказано.

Обозначим Ф(Ь) характеристическую

функцию случайной величины пь Из (2) нетрудно получить, что

Ф(Ь) = 1 + (ей — 1)Ф(ей, т, 1), (11)

где Ф(,г,з,а) - транцендентная функция Лер-ча (см., например, [Иа]о1е1, Sedgewick, 2009]), задаваемая соотношениями

Ф(,г, s, а) =

ГО i

¿0 (1 + a)S

r(s)

xs— 1e—ax

1 — ze—x

dx,

(12)

Г(з) - значение гамма-функции в точке 8. Используя (11) и формулу обращения для решетчатых распределений, можно уточнить (10) следующим образом:

Р{^ = ^

= Ж-1/т (/((k — NZ (т ))/Ж1/т) + 0(Ж-1/т)),

(13)

причем это равенство выполняется равномерно относительно k таких, что величина ^ — NZ(т))Ж-1/т лежит в любом фиксированном конечном интервале.

89

т

1

Рассмотрим Р{А1}. Ясно, что

Р{Аі} = £ Р{Аі,пі = к}. (14)

к<Ь

Обозначим £(1) степень первой вершины первого поколения. Тогда

Р{А1,П1 = 1} = ^ Р{А1,П1 = 1,С(1) = 8}.

(15)

Пусть означает число вершин степени г в графе. В теореме 6 статьи [Павлов, 2009] показано, что при Ж, п ^ то так, что п—Z(т)Ж = о(Ж1/т) для целых неотрицательных k

Р{^г = ^

1 + о(1) ( (к — )2

л/^РГ)

ехр

2^г (1 — рг)

(16) (к-№рг)

равномерно относительно , 4 7 в лю-

(і-Рг )

бом фиксированном конечном интервале. Повторив доказательство этой теоремы при условии

0 < Сі ^ (п — £(т)N)/Nі/т ^ С2 < то, здесь и далее символы ,... означают некоторые

положительные постоянные, легко получить, что для таких же к

Р{^г = к}

(1 + о(1))/((п — С (т )N )/^ (1 — Рг ))і/т)

(1 — Рг )і/тУ2ЩТ—РГ) / ((п — С (т )Nі/т)

х ехр < —

(к — )2

2^г (1 — рг)

(17)

= Рі^і

Р{Аі,Пі = 1,Є(і) = 1}

(п — 3)(п — 5)...3 ■ 1

Отсюда и из (15) следует:

Р{Аі ,пі = 1} = Рі(1 —Е &-і)(1+0^і/т-і)).

(19)

Из (2) при к ^ то получаем :

Рк ~ т/кт+і,

поэтому с учетом (4)

£

(20)

(21)

Отсюда и из (19) вытекает :

Р{А1,П1 = 1} = Р1(1 + 0(Н1-т)). (22)

Подобным же образом можно рассмотреть поведение вероятности Р{А1, П1 = 2}. Тогда

Р{Аі, пі = 2} = £ Р{Аі, пі = 2, Є(і) = «}.

(23)

Заметим, что

Р{Аі,пі = 2,Є(і) = 5}

= £ Р{Аі,пі =2,Є(і) = 1,Є(і) = і}, (24)

г<Ь,

(1)

где Ц2 - степень второй вершины первого по-

коления (такая вершина, очевидно, существует, поскольку событие Аі подразумевает, что у начальной частицы нет петель). Легко видеть, опять учитывая (4), (16), (17) и условие ^ = п, что

Р{Аі,пі =2,Є(і) = 1,Є(і) = 1}

Npl N(п — 5)(п — 7)...3 ■ 1

(1)

Р2

2

(п — 1)(п — 3)...3 ■ 1

Учитывая соотношения N і/2 = о^і/т і), п— с (т )N = і/т), процедуру построения гра-

фа и свойства функции Л, из (4), (16) и (17) получаем, что

К1 + Тм 0 = Р2^о (1 + o(Nі/т-і))

и, аналогично, для 2 ^ і < Л

Р{Аі,пі = 2,Є(і) = 1,Є(і) = і}

Р2<?0<?г

(п — 1)(п — 3)...3 ■ 1

= Р2С-і(т )(1+0^і/т-і)) = Рідо(1+0^і/т-і))

Рассуждая аналогично для 2 ^ в ^ Л, находим, что

-і (1 + o(Nі/т-і))

1 + -і/2)^ Отсюда и из (21), (24) находим, что

Р{Аі, пі = 1, Є(і) = в} = Рі^-і (1+0^і/т-і)).

(18)

Р{Аі,пі = 2,Є(і) = 1}

= Р29о (1 + 0(Лі-т)) (1 + o(Nі/т-і))2 .

Продолжая этот процесс, получаем, что при 2 ^ в < Л

Р{Аі,пі = 2,Є(і) = в}

2

= Р298-1 (1 + 0(Н1-т)) (1 + о(Ж1/т-1))2 .

Поэтому из (21) и (23) следует, как легко видеть, что

Р{А1,П1 = 2} = Р2 (1 + 0(Н1-т))2 .

Теперь становится понятным, что для k = 1, 2, ...

Р{А1,П1 = k} = (1 + 0(Н1-т))к .

Отсюда и из (1), (3), (14) и (32) следует: Р{А1} = ^Рк (1 + О^Н1-т)) = 1 + 0(Н1-т)

к<Ь

(25)

Далее рассмотрим Р{А1}. Проводя выкладки, подобные приведенным выше, находим, что

Р{А2,Пі = 1} = Рі ^ &-і( 1 &)

¿-і

¿<Ь

¿-і

(26)

¿<ь

Рассмотрим выражение

¿і

¿<ь

с і(тН2^¿(1 — ^9*)

¿=і

— £ іРІ(1 — £Ы

¿і

) ¿ і)

(27)

¿>ь

Е £

Ф(А, т, 1) — (1 — А)Ф(А, т, 1)

1 — (1 — А)Ф(А, т, 1)

Ф(Л,т, 1) = Л-і(С (т)+ 0((1 — А)т-і)).

ті

Отсюда и из явного вида А находим:

Е Є = С (т)+ 0(Л-(т-і)2).

Заметим, что в силу (20)

£ ірі(1 — £ < С5 £ і-т = 0(Л

іт

¿>ь

¿>ь

Обозначим А = 1 — ^ и рассмотрим вспомогательную случайную величину £, имеющую распределение

Р{£ = k} = АкРк(1 — (1 — А)Ф(А, т, 1))-1, (28)

где k = 1, 2,..., а Ф(А, т, 1) определено в (12). Нетрудно найти, что

поэтому из (26) - (29) следует:

Р{А2 ,П1 = 1} = Р1(1 + 0(Н(т-1)2)). Таким же образом находим, что для г = 1, 2, ...

Р{А2,П1 = г} = р^1 + 0(Н-(т-1)2)), следовательно

Р{А2} = 1 + 0(Н-(т-1)2).

Учитывая (25) и продолжая этот процесс по индукции, получаем утверждение леммы. □

Связь с вЕтвящимися процессами

Рассмотрим начинающийся с одной частицы ветвящийся процесс Гальтона - Ватсона С, в котором число прямых потомков начальной частицы имеет распределение (2), а число прямых потомков остальных частиц имеет распределение (4). Обозначим ^(С) число частиц Ь-го поколения. Связь между рассмотренными случайными графами и процессом С можно описать следующим образом.

Лемма 3. Пусть Ж, п ^ то так, что п — Z(т)Ж = 0(Ж1/т), Ь = о(1п Ж/ 1п1п Ж). Тогда для т = 0,1, 2,...

Р{£* = т} = Р{Я*(С) = т}(1+0(Н-(т-1)‘-1)). Доказательство. Ясно, что

Р{^о = 1} = Р{^о(С) = 1} = 0.

Из леммы 2 и (13) следует :

Р{£1 = т} = Р{п1 = т|^ = п}

= рт(1 + -і/т)),

(30)

поэтому, в силу равенства Р{^1(С) = т} = рт, получаем утверждение леммы для Ь = 1. Пусть Ь = 2. Согласно лемме 2,

и в наших условиях А ^ 1, Е £ ^ Z (т). Для того, чтобы получить более точную оценку Е £, воспользуемся асимптотическим разложением Ф(А, т, 1) при А ^ 1 (см. [Иа^е^ Sedgewick, 2009]):

Р{^2 = т} = Р{^2 = т, Аі} + 0(Лі-т). (31) Рассмотрим

Р{^2 = т, Аі} = £ Р{пі = к, ^2 = т, Аі}.

к<Ь

(32)

91

Из (18) и (21) находим, что

Р {пі = 1,^2 = т, Аі}

= Рі«т-1 (1 + o(Nі/т-і))(1 + 0(Лі-т)). Учитывая (16), (17) и условие п — ((т)N = о^і/т ), получаем :

^Р{пі = 2,^2 = 2,£(і) =3

к=і

Р{пі = 2,^2 = 0,Аі} = Р22

Npl

х(1 + 0(№-і/2))іп — 5>!п — 7І-3 • 1

1 1 ;;(п — 1)(п — 3)...3 ■ 1

= Р2«0(1 + o(Nі/т-і)). (33)

Нетрудно видеть далее, что

Р {пі = 2,^2 = 1,Аі} =

2 £ Р{пі = 2,^2 = 1,е(і) = к,е(і) = 1,Аі},

і<к<Ь

(п — 1)(п — 3)...3 ■ 1

х 1 + О

1

) = Р29і9о2(1 + o(Nі/т-і))3

Подобным же образом для 2 ^ к ^ Л приходим к равенствам:

Р {пі = 2,^2 = 1,Є(і) = к,Є(і) = 1,Аі}

= Р29і9о9к-і(1 + o(Nі/т-і))3, поэтому из леммы 2 и (34) находим:

Р {пі = 2,^2 = 1,Аі}

= 2р25о9і(1 + o(Nі/т-і))3(1 + 0(Лі-т)). (35)

Дальнейшие рассуждения проводятся по той же схеме. Ясно, что

Р {пі = 2,^2 = 2,Аі}

= 2 Р{пі = 2, ^2 = 2, Є(і) = 3, Є(і) = 1, Аі}

+ Р{пі = 2, ^2 = 2, Єіі) = 2, Є2і) = 2, Аі}.

Обозначим £(2) и первую и вторую вершины второго поколения. Нетрудно видеть, что

Р{П1 = 2,^2 = 2, £(1) =3,£21) = 1,А1}

Є(і) = 1,Й2) = к,Аі}. (37)

Очевидно также, что

Р{пі = 2, ^2 = 2, Є(і) = 3, Є(і) = 1, Є(2) = 1, Аі}

ГО

(і) (і) (2)

^Р{пі = 2,^2 = 2,Є(і) =3,Є(і) = 1,£

«=і

= 1, е22) = в,Аі}.

(38)

Далее получаем:

Р{пі = 2, ^2 = 2, е(і) = 3, е(і) = 1, е(2) = 1,

{(2) = в, Аі} = Р2<?о<?29«-і(1 + о^і/т-і))4, поэтому из (38) следует:

Р{пі = 2, ^2 = 2, Є(і) = 3, Є(і) = 1, Є(2) = 1, Аі}

= Р2 «2«2(1 + o(Nі/т-і))4. (39)

Теперь рассмотрим

Р{п і = 2, ^2 = 2, і ) = 3, і ) = 1, Є(2) = 2, А і}

(34)

где обозначения 1) и £(1) имеют тот же смысл, что и в доказательстве леммы 2. Используя (16), (17), получаем:

Р{П1 = 2,^2 = 1, £(1) =2,£21) = 1,А1}

= Р2 ( ^ ) Жр24-(п — 7)(п — 9).”3 ' 1

£ Р{пі = 2, ^2 = 2, е(і) = 3, е(і) = 1, Є(2)

«=і

= 1. е22) = в}.

(40)

Как обычно, находим:

Р{П1 = 2, ^2 = 2, £(1) = 3, £(1) = 1, £(2) = 2,

£(2) = 8, А1} = Р24о91929.8-1(1 + о(Ж1/т-1))4 и из (40) получаем:

Р{П1 = 2, ^2 = 2, £(1) = 3, £(1) = 1, £(2) = 2,

А1} = Р2?о9192(1 + о(Ж1/т-1))4. (41)

Продолжая этот процесс и учитывая (39), (41), видим, что

Р{П1 = 2,^2 = 2, £(1) = 3, £(1) = 1, £(2) = ^ А1}

= Р29о929к:

-і(1 + о^і/т-і))4

(36) поэтому из (37) следует равенство:

Р {пі = 2,^2 = 2,е(і) =3,е2і) = 1,Аі}

= Р29о«2(1 + o(Nі/т і))4.

(42)

2

Рассматривая следующую вероятность по той же схеме, находим:

Р {пі = 2,2 = 2,е(і) =2,е2і) =2,Аі}

= Р2«2(1 + o(Nі/т-і))4.

Отсюда и из (36), (42) вытекает:

Р {пі = 2,2 = 2,Аі}

= (2р2«о92 + Р2«2)(1 + o(Nі/т-і))4. (43)

Подобно тому, как были получены равенства

(33), (35) и (42), из леммы (2) находим:

Р{пі = 2,2 = т, Аі}

= Р2 ^ «кі «к2 (1 + о(Лт-і)).

кі, к2 ^ 0, к1 + к2 = т

Аналогичным образом можно вывести, что Р{пі = к, 22 = т, Аі}

= Рк ^ «кі ^ ...9кк (1 + 0(Лі-т)).

кі, ...,кк ^ 0, кі + ... + кк = т

Отсюда и из (1), (31), (32) следует:

Р{2 = т}

ГО

= ЕРк X] «кі«к2 ...9кк (1+ 0(Лі-т)).

к=і кі,..., кк ^ 0, кі + ... + кк = т

(44)

Из определения ветвящегося процесса С и (30), (44) ясно, что утверждение леммы верно при Ь ^ 2. Как и при доказательстве леммы 2, продолжим рассмотрение вероятностей Р{2 = т} по индукциии и придем к утверждению леммы 3. □

Ядро графа

В работе [Ивійи, Копте, 2004] показано, что в Интернет-графе (без ограничений на число связей) вершины, имеющие достаточно большие степени, концентрируются вокруг вершины с максимальной степенью і * (если таких вершин несколько, то в качестве і* можно взять любую из них) и соединены с ней путями небольшой (порядка 1п1п N) длины. Такой подграф получил название ядра. Не удивительно, что подобная структура возникает и в рассматриваемых условных графах.

Обозначим $1 = (т — 1)/т, П(м) - степень вершины ' *, е = е(Ж) = ад/ 1п Ж, где ад определено в (7). Пусть и означает множество вершин следующего вида:

и = {' : п- ^ ^1+£}. (45)

Согласно теореме 4 из [Павлов, 2009], в наших условиях при г = (Ж/7)1/т(1 + о(1)), где 7 -некоторая положительная постоянная, справедливо соотношение

Р{П(М) < г} У в-7, (46)

откуда ясно, что ' * € ^1. Нетрудно видеть, что

У то. (47)

Учитывая (45) - (47), с помощью леммы 1 получаем, что любая вершина из и1 а.д. соединена ребром с '*. Обозначим

г-1

$г = (т — 1)г/т + е £(т — 1)к, г = 2, 3,..., (48) к=о

и пусть

= {' : ЖЛ+£ ^ пг < Ж*-1+£}, г = 2, 3,...,.

(49)

Множество вершин графа

а

С = и иг (50)

г=1

назовем ядром.

Лемма 4. Справедливы следующие утверждения.

1. Любая вершина ядра, принадлежащая множеству иг+1, г = 1, 2,..., а.д. соединена ребром с какой-нибудь вершиной множества иг.

2. Любая вершина ядра, не принадлежащая

«+2

множеству {] иг, в = 1, 2,..., а.д. не

г=1

имеет ребер, соединяющих ее с вершиной множества из.

Доказательство. Из (1), (45), (49) нетрудно показать, что при Ж У то справедливо соотношение Р{' € иг} Ж (*+е)т. Используя известное приближение биномиального распределения нормальным, находим, что число вершин в иг асимптотически эквивалентно Ж 1-(л+е)т. Отсюда и из (49) получаем, что сумма степеней вершин из и а. д. пропорциональна Ж1-(т-1)(* +е). Из (48) следует, что

$г+1 = (т — 1)$г + е, (51)

поэтому из (47), (49) с помощью леммы 1 и соотношения n/N ^ Z (т) приходим к первому утверждению леммы 4.

Как мы видим, число степеней вершин в Us пропорционально N1-(т-1)(gs+e). Из (49) следует, что степень любой вершины, не принад-

s+2

лежащей U Ui, меньше, чем Nös+2+£. Отсюда

i=1

с помощью (51) и леммы 1 получаем второе утверждение леммы 4. □

Из леммы 4 очевидным образом вытекают такие следствия.

Следствие 1. Ядро а.д. состоит из непере-секающихся слоев U1 , U2,..., концентрически расположенных вокруг вершины j*. Кратчайший путь из вершины j* в любую вершину ядра, принадлежащую слою Ui, i = 2, 3,..., проходит последовательно через вершины, принадлежащие слоям U1, U2,..., Ui-1, и он не содержит ребер, соединяющих непосредственно вершины из слоев, номера которых отличаются больше, чем на 2.

Следствие 2. Длина кратчайшего пути, соединяющего любые две вершины ядра, а.д. не превосходит 2а.

Хорошо известно [Севастьянов, 1971], что определенная в теореме вероятность вырождения q удовлетворяет уравнению

q = F(q), 0 < q < 1. (52)

Выберем произвольную вершину j* и обозначим Dj событие, состоящее в том, что эта вершина принадлежит компоненте связности, содержащей ядро.

Лемма 5. Пусть N ^ то. Тогда P{Dj} ^

1 - ^(q).

Доказательство. Найдем сначала вероятность вырождения процесса G. Следуя [Севастьянов, 1971], введем обозначение

ГО

Ft(z, G) = ^ P{Zt(G) = k}zk, t = 0,1,2,... fc=0

Известно [Севастьянов, 1971], что Ft(z,G) = ^>(Ft(z,G)). Полагая z = 0, отсюда получаем:

P{Zi+1(G)=0} = ^(P{Zt (G) = 0}).

Переходя в этом равенстве к пределу при t ^ то и учитывая (5), (52), находим вероятность вырождения G:

lim P{Zt(a) = 0} = <^(q). (53)

Обозначим Dj,t событие, состоящее в том, что к моменту t закончилось построение первой связной компоненты графа, не содержащей вершин из ядра, с помощью описанного выше алгоритма, работа которого началась с вершины j. Ясно, что в этом случае

P{Dj,t} < P{Z = 0}. (54)

Положим h = ln N, t ^ ln ln ln ln N. Из леммы

2 следует :

P{Dj,t} ^ P{Dj,t, At}

= P{Zt = 0} + O((ln N )-(T-1)2). (55)

Утверждение леммы 5 следует из леммы 3 и соотношений (53) - (55). □

В лемме 3 установлено соответствие между рассматриваемыми случайными графами и ветвящимся процессом G. С помощью этого соответствия ниже, в лемме 6, будет показано, что асимптотически эквивалентны события, состоящие в том, что процесс G выродился и что соответствующий подграф, построенный с помощью алгоритма, не содержит вершин ядра.

Лемма 6. Пусть N, n — то так, что n — Z(т)N = O(N1/т). Если вершина графа принадлежит компоненте связности, содержащей ядро, то расстояние от нее до ядра, а. д. не превосходит в.

Доказательство. Мы будем следовать схеме рассуждений, приведенных в [Reittu, Norros, 2004] при доказательстве аналогичного утверждения о безусловных случайных графах.

Для последовательности номеров вершин ji,j2,..., полученной с помощью алгоритма, введем обозначения:

^ = ^(ji,j2,..., ) = inf {k : jfc = j, 1<k}, k k

A = inf{k : ^(nji — 1) ^ Sk = ^ nji,

¿=1 ¿=1

m = [4в/ P{n1 > 1}].

Докажем, что

P{A ^ m|vN = n} —— 1. (56)

Действительно, нетрудно видеть, что в силу леммы 3 [Павлов, 2007] при n — то

P{nJl > 1|vn = n} — P{n1 > 1}. (57)

Обозначим I(A) индикатор события A. Для случайной величины X, имеющей биномиальное распределение с параметрами n,p, справедлива оценка [Феллер, 1984]:

P{X ^ np/2} ^ 4(1 — p/2)/np, (58)

С94)

поэтому с помощью (57), (58) находим, что P{A ^ m|vN = n}

{Ш

^ I(j > 1) > e|vN = n

¿=1

{Ш

^ I(nj, > 1) > m P{n1 > 1}/4|vn = n ^

Ш

^ I(nji > 1) > m P{^1 > 1|vn = n}/2|vN

¿=1

= n j — P{P{n1 > 1|vn = n} > P{n1 > 1}/2|vn

= n} ^ 1 — 8/m P{n1 > 1|vn = n} + o(1) — 1,

откуда и следует (56).

По определению A для любого и эквивалентны неравенства u ^ A и ^ и + в. Отсюда нетрудно видеть, что неравенства

и ^ A и Su ^ A + в (59)

также эквивалентны. Для всех l = o(N)

P{SCT-1 < 1|vn = n}

i

= ^ P{ct = i, Si-1 < 1|vn = m}. (60)

¿=1

Используя лемму 3 [Павлов, 2007], находим, что

P{ct = i, Sj-1 < 1|vn = n} ^ C11/n, поэтому из (60) следует:

P{SCT-1 < 1|vn = n} ^ C112/n.

Отсюда и из (56), (59) получаем:

P{ct > A|vn = n} = P{£o-_1 ^ A + e|vN = n} ^ P{SCT-1 ^ m+в|vn = n} —P{A > m|vN = n} ^ 1 —(m+e)2/n—P{A > m|vN = n} — 1. (61)

Обозначим L = inf{i : j G G}. Оценим эту величину сверху. Для того, чтобы вершина принадлежала ядру, необходимо (см. (6), (48), (50), (51)), чтобы ее степень была не меньше, чем Nga+e. Как легко видеть из (51), это значит, что при достаточно больших N степень такой вершины должна превосходить N(3-т)е/(2-т). Отсюда и из (1) следует, как нетрудно проверить, что для любой вершины j справедливо соотношение P{j G G} ~

Ж-(3-т)те/(2-т). Учитывая нормальное приближение биномиального распределения, для случайной величины V, равной объему ядра, получаем неравенство V ^ С2Ж1-(3-т)те/(2-т), откуда следует, что число степеней вершин ядра больше, чем С3Ж1-(т-1)(3-т)е/(2-т). Тогда из (6), (49) и условия п/Ж У ^(т) находим:

Р{'1,..., 'в) П С = 01^м = п}

< (1 — С4Ж(1-т)(3-т)е/(2-т))в у 0, откуда

Р{£ < в} У 1. (62)

□

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

Предположим, что построенный с помощью алгоритма и начинающийся с рассматриваемой вершины подграф содержит 5 < то поколений вершин. Положим Н = о(Ж(3-т)е/(2-т)). В силу лемм 2 и 3 это значит, что ветвящийся процесс О к моменту 5 выродился, а построенный подграф не содержит вершин ядра. Если же, наоборот, процесс О не выродился до рождения в частиц, то из (61) и определения а следует, что мы можем считать идентичными траекторию процесса до этого момента и соответствующую реализацию подграфа и, в силу (62), это значит, что ядро достигается.

Теперь нетрудно доказать теорему. Первое утверждение очевидным образом следует из леммы 5 и свойств биномиального распределения, а второе - из следствия 2 и леммы 6.

Литература

Ватутин В. А. Распределение расстояния до корня минимального поддерева, содержащего все вершины данной высоты // Теория вероятностей и ее применения. 1993. Т. 38, № 2. С. 273-287.

Ибрагимов И. А., Линник Ю. В. Независимые и стационарно связанные величины. М.: Наука, 1965. 424 с.

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

Павлов Ю. Л. Предельное распределение объема гигантской компоненты в случайном графе Интернет-типа // Дискретная математика. 2007. Т. 19, № 3. С. 22-34.

Павлов Ю. Л. О предельных распределениях степеней вершин в условных Интернет-графах // Дискретная математика. 2009. Т. 21, № 3. С. 1423.

Павлов Ю. Л., Чеплюкова И. А. Случайные графы Интернет-типа и обобщенная схема размещения // Дискретная математика. 2008. Т. 20, №

3. С. 3-18.

Севастьянов Б. А. Ветвящиеся процессы. М.: Наука, 1971. 436 с.

Феллер В. Введение в теорию вероятностей и ее применения. Т. 1. М.: Мир, 1984. 527 с.

Durrett R. Random Graph Dynamics. N. Y.: Cambridge University Press, 2007. 221 p.

Faloutsos C., Faloutsos P., Faloutsos M. On power-law relationship of the Internet topology // Computer Communications Rev. 1999. Vol. 29. P. 251-262.

Flajolet P., Sedgewick R. Analytic Combinatorics. Cambridge: Cambridge University Press, 2009. 824 p.

Hofstad R. Random Graph and Complex Networks. Department of Mathematics and Computer Sciences. Eindhoven University of Thechnology. 2011.

URL:http://www.uin.tue.nl/ rhofstad/ Notes

RGGN.pdf (дата обращения: 11.04.2011)

Reittu H., Norros I. On the power-low

random graph model of massive data networks // Performance Evaluation. 2004. Vol. 55, N. 1. P. 3-23.

СВЕДЕНИЯ ОБ АВТОРБ:

Павлов Юрий Леонидович

зав. лаб. теории вероятностей и компьютерной статистики, д. ф.-м. н.

Институт прикладных математических исследований КарНЦ РАН

ул. Пушкинская, 11, Петрозаводск, Республика

Карелия, Россия, 185910

эл. почта: pavlov@krc.karelia.ru

тел.: (8142) 761218

Pavlov, Yury

Institute of Applied Mathematical Research, Karelian Research Centre, Russian Academy of Science 11 Pushkinskaya St., 185910 Petrozavodsk, Karelia, Russia

e-mail: pavlov@krc.karelia.ru tel.: (8142) 761218