Решения сетевых игр с попарным взаимодействием

2019 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 15. Вып. 1 _ПРИКЛАДНАЯ МАТЕМАТИКА. ИНФОРМАТИКА. ПРОЦЕССЫ УПРАВЛЕНИЯ_

ПРОЦЕССЫ УПРАВЛЕНИЯ

УДК 519.71 МБС 91А12

Решения сетевых игр с попарным взаимодействием*

М. А. Булгакова

Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

Для цитирования: Булгакова М. А. Решение сетевых игр с попарным взаимодействием // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2019. Т. 15. Вып. 1. С. 147-156. https://doi.org/10.21638/11702/ spbu10.2019.112

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

Ключевые слова: кооперативные игры, выпуклость, вектор Шепли, характеристическая функция.

Введение. Теория кооперативных сетевых игр — важная часть современной теории игр, которая будет использоваться для построения решений в играх на сетях со многими участниками. Эта теория включает в себя кооперативную траекторию, стратегии, ее порождающие, выигрыш вдоль кооперативной траектории, а также распределение выигрыша между игроками и анализ динамической устойчивости решений. Сеть иллюстрирует наличие взаимодействия между парами игроков. Впервые попарное взаимодействие на сети, т. е. взаимодействие один на один между игроками, в некооперативной форме было исследовано в [1].

Попарное взаимодействие было рассмотрено и на примере распространения информации и дезинформации в социальных сетях в работе [2]. Эффективность и устойчивость сетей, зависящих от внешних факторов, таких как суммарные предельные издержки, были показаны в [3]. Подходам для нахождения оптимального поведения

* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 17-51-53030).

© Санкт-Петербургский государственный университет, 2019

в многошаговых играх посвящена работа [4]. Условия сильной динамической устойчивости в двухшаговой игре с попарным взаимодействием были найдены в [5]. Динамические свойства кооперативных решений в игре п-лиц были рассмотрены в [6].

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

Модель игры. Пусть N — конечное множество игроков, которые принимают решения в двух состояниях, N\ = п ^ 2. Будем через г обозначать состояние игры. Игра начинается в состоянии г\ — в нем каждый игрок г е N выбирает свое поведение Ь1 = (Ъ1!,... ,Ь'1п) — п-мерный вектор предложений связи другим игрокам. Введем следующие обозначения: Ы. с N \ г — те игроки, которым игрок г е N может предложить связь, значение а. е {0,...,п - 1} равно максимальному числу связей игрока г. Если Ы. = N \ {г}, игрок г может предложить связь всем игрокам и, если а. = п - 1, игрок г может поддерживать любое число связей.

Для каждого поведения Ь1 существует такое подмножество Qi с Ы., для которого выполняется выражение

Ь1 = { 1, если ] е Qi

113 I 0, в противном случае, при условии

<

Условие (1) означает, что число возможных связей ограничено для каждого игрока. Также очевидно, что \Qi \ ^ а..

Связь г] имеет место тогда и только тогда, когда Ь! = Ь1 = 1. Сформированные связи г] образуют ребра сети д, вершинами которой являются игроки, т. е. если Ь! = Ь! = 1, то в сети д появляется ребро с концевыми вершинами г и ].

Обозначим через Щ(д) соседей игрока г в сети д, т. е. Щ(д) = {] е N \ {г} : г] е д}. Далее для краткости иногда вместо ^(д) будем писать N.1. Результатом первого состояния является сеть д(Ь1,... ,ЬП). После ее формирования игроки переходят в состояние (д), которое обусловливается сетью (от сети зависит множество соседей N.1, а следовательно, и правило взаимодействия между игроками). Во втором состоянии г2(д) соседи по сети играют попарно в одновременные биматричные игры, после чего игроки получают выигрыши и игра заканчивается. Другими словами, имеем двухшаговую игру Г21 (д), которая является частным случаем многошаговых неантагонистических игр. Адаптируя к данному случаю определение стратегий, примем, что в рассматриваемом случае стратегия — это правило, которое для каждого игрока определяет множество его соседей в первом состоянии, а именно вектор Ь1, и поведение в каждой биматричной игре во втором состоянии в соответствии с сетью, которая сформирована в первом состоянии, — Ь2. Обозначим через и. = (Ь1,Ь2), г е N, стратегию игрока г в двухшаговой игре Г21 (д). Рассчитаем выигрыш игрока г как Н.(г2), где (г1,г2) — траектория, реализованная в ситуации и = (и1(-),...,ип(-)) в игре Г21 (д). Так как в первом состоянии игроки не получают выигрыши, то функция выигрыша в игре Г01 (д) с начальной позицией г1 определяется следующим образом:

К.(г1; и) = К.(г1; щ(-),...,ип(-)) = Ы(г2).

Характеристическая функция. Во втором состоянии игра представляет собой семейство попарных одновременных биматричных игр } между соседями на сети. А именно, пусть г е N,2 е Ы^. Тогда г играет с 2 в биматричную игру ^г] с матрицами выигрышей А] и С]г игроков г и 2 соответственно:

= / aij a11 ij ai2j1 aij a12 ij a22 • aij \ a1k ij a2k = / j c11 c21 c12 j c22 j \ c1k c2k

ij \a'm1 ij am2 • ij  a , . mkl j \cm1 cm2  j mk /

api Z 0, cpi Z 0, p = l,...,m, l = l,...,k.

Константы m и k одинаковы для всех i и j. Обозначим через Tf2 (g) подыгру игры Г, которая происходит в состоянии х2. Опишем такую игру в кооперативной форме. Найдем характеристическую функцию для каждого подмножества (коалиции) S с N как нижнее (максиминное) значение антагонистической игры двух лиц коалиции S и дополнительной коалиции N \ S, построенной на основе игры Г^2 (g), при этом выигрыш коалиции S рассматривается как сумма выигрышей игроков, входящих в S. Супераддитивность характеристической функции следует из ее определения. Обозначим, что

maxmin aije, p = l,...,m, l = l,...,k,

p i pi

maxmin cj\, p = 1,...,m, l = 1,...,k,

i p pi

и v(z2; S), S с N, — нижнее значение антагонистической игры Tf2 (g).

Теорема 1. Функция v(z2; S) определяется по следующим формулам:

Wij = Wji =

{0})= 0,

К^ЦЕ Е т^а« +<$) + £ £ ьк, в с Ж,

у

2 геЫ Р

Доказательство. Формула (2) очевидна.

Докажем (3). Поскольку игрок г, действуя против коалиции N \ {г}, играет с игроками 2 из N п Ыг в независимые биматричные игры, то в каждой из этих биматричных игр он может себе гарантировать наибольший выигрыш , следовательно, наибольший выигрыш, который игрок г может гарантировать себе во всей игре, — сумма Ы]. Это и есть в данном случае максиминный выигрыш г против N\ {г}.

Формула (3) справедлива.

Докажем формулу (4) для произвольной коалиции Каждый из игроков, входящих в играет в независимые попарные биматричные игры как с игроками в N п 5, так и с игроками в N n{N \ 5}. В первом случае игроки, взаимодействующие внутри 5, всегда могут выбрать стратегию, максимизирующую их суммарный выигрыш,

т. е. Y maxpe(aj + cj). Во втором случае игроки, взаимодействующие с игрока-

jtNinS ' Р Р

ми из N \ S, могут гарантировать себе лишь нижнее значение игры, т. е. величину Y wik. Таким образом, максимальный суммарный выигрыш, который может обес-

kiNi\S

печить коалиция S, будет равен

^Е Е ma„xK><£) + E Е wik-

AiiSjiNinS piiSkiNi\S

Формула (5) следует из определения максимального суммарного выигрыша игроков.

В формулах (4) и (5) коэффициент ^ исключает повторное суммирование выигрышей по одним и тем же ребрам.  Рассмотрим кооперативную форму двухшаговой игры rzi (g). Предположим, что игроки выбирают стратегии ui,i е N, которые максимизируют их суммарный выигрыш в игре rzi (g), т. e.

Y, Ki(z1;ü1,...,un) = max ^ Ki(z1; ui ,...,un).

iiN u iiN

Ситуацию u = (uui,... ,un) будем называть кооперативным поведением, соответствующую траекторию (zi,z2) — кооперативной траекторией, в данном случае состоящую из двух состояний.

Как и раньше, для коалиции S £ N определим характеристическую функцию v(zi; S) как максимин в антагонистической игре двух лиц коалиции S и дополнительной коалиции N \ S, где выигрыш коалиции S равен сумме выигрышей ее участников и стратегия S — элемент декартова произведения множеств стратегий игроков, входящих в S.

Обозначим через v(zi; S), S с N, нижнее значение антагонистической игры

rzi (g).

Теорема 2. Функция v(zi; S) определяется по следующим формулам:

v

v

2 iiS jiNi (g)nS p'

v

2 iiN jiNi p'

Доказательство. Доказательство аналогично доказательству теоремы 1, однако здесь в (7) в отличие от (4) отсутствует второе слагаемое. Оно равно нулю, поскольку минимизирующая стратегия b1 игроков из N \ S в первом состоянии состоит в отказе от связей с игроками из S. Это следует из неотрицательности выигрышей во всех играх.

Выпуклость игры. Характеристическая функция называется супермодулярной, а соответствующая ей игра выпуклой [7], если для любых произвольных коалиций X с N и Y с N выполняется неравенство

v

Теорема 3. В подыгре Г^ (д) характеристическая функция (2)-(5) супермоду-лярна.

Доказательство. Докажем неравенство (9) для характеристической функции (2)-(5). Будем для сокращения записи вместо шах^Дар], + арг£) писать тг3. Имеем

г;

2 гАеХиУ геХиУ

]еМ€ кеМ€\ХиУ г*3

у

2 г] еХ гХ

3 еМ^ \Х

г*]

1

00 = е

геУ 2 гз^у

кеЫ^У

г*3

ю

гхпУ гхпУ

]ХпУ кеЩ\ХпУ

г*3 3N

Вычитая из (10) выражения (11), (12) и прибавляя (13), получим неравенство

Е т3 > 0,

гХ\У \Х

которое следует из неотрицательности выигрышей. Таким образом, неравенство (9) в игре Г^2 (д) выполняется.

Итак, построенная характеристическая функция v(z2; 5) супермодулярна. Это гарантирует непустоту С-ядра и принадлежность вектора Шепли С-ядру.

Вектор Шепли. Поскольку вектор Шепли в подыгре Г^2 (д) всегда принадлежит С-ядру, то его значимость в данном классе задач значительно возрастает. Приведем формулу для вектора Шепли для одной конкретной сетевой игры с попарным взаимодействием.

Определим дележ как вектор = ((1 [«],...,£„[«]), который удовлетворяет условиям коллективной рациональности, т. е. = N), и индивидуальной

рациональности, т. е. ^ «(•; {г}), для всех г е N. Обозначим множество всех

дележей в игре Г21 (д) через I(«). В качестве дележа возьмем вектор Шепли [8] =

(Vl[v],...,Vn[v]), где

= Е

SsN,i£S п-

Вычисление компонент — достаточно трудный процесс для большого числа

игроков. Предположим следующее. Пусть М1 = N\{1}, а1 = п-1 и Мг = {1}, аг = 1 для г Ф 1. Тогда в первом состоянии игры в целях максимизации суммарного выигрыша игроки должны выбрать следующие поведения:

Ь1 = {

Ьг = {

V

Поведение (15) формирует сеть-звезду в первом состоянии с центральным игроком 1 (рис. 1), в которой \М1\ = п - 1 и ¡N1 = 1, г * 1.

Рис. 1. Сеть-звезда

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

Сеть имеет центральную симметрию, которая позволяет предположить, что формула (14) может быть упрощена.

Утверждение 1. Для сети-звезды с центральным игроком 1, компоненты вектора Шепли )], £ = 1,2, представимы в виде [9]

-лК^)] =

{1}) + т,(т1з- Ш)) , г = 1,

3*1 J

{г}) + ти - ти], г * 1.

Доказательство: Сначала докажем выражение (16) для подыгры Г^ (д). Зафиксируем игрока г € Б, г * 1 и рассмотрим его предельный вклад го(х2; Б)-у(г2; Б\ {г}) в коалицию Б. Принимая во внимание формулы (2)—(5) для характеристической функции у(г2; •), получаем, что

)(г2; Б) -у(г2; Б \ {г}) =

ти - ти, 1 € Б, тц, 1 £ Б.

Подставляя (2)—(5) в формулу компонент вектора Шепли (14), для игрока г * 1 имеем

Б£М -МБ

= (тц - т1г) £

^ „ (\Б\-1)!(п-\Б\)!

+ 1"г Ъ -\-

Бс№ ;гёБ;1£Б п!

" (\Б\-1)!(п-\Б\)!( п - 2

= (тц - т\г) £

( П - 2 ) (\Б\ - 2/

(\Б\-1)\(п-\Б\)\(п-2\

= (ти - т\1) £

151-1 + .а V1 п~\в\

=2 п

1 п

1 Г il lin

= 9 Lw» + mii ~ wi \  Следуя свойству эффективности вектора Шепли, сразу находим, что

VI[>'(^2)] = V(Z2;N) - = \ Е Е тп - \ Е N1 + - №11г] =

гф1 2 iaN jaNi 2 гф1

= Е тн - \ Е N1 + mii - wi1'] = ^ Е [и,11г + mii - ч1] •

С учетом того, что v(z2; {1}) = EjN1 = w1 и v(z2;{г}) = EjN wj = W1, « * 1, выражения (16) справедливы для подыгры Г S (g).

Теперь докажем (16) для игры Г21 (g). Аналогично, зафиксируем игрока г € S, г * 1 и рассмотрим его максимальный вклад v(z1; S) - v(z\;S \ {г}) в коалицию S. Используя выражения (6)—(8) для характеристической функции v(z1; •), получим

v(z1; S)-v(z1;S\{«}) = | ^^ 1 € S

Подставляя их в формулу для компонент вектора Шепли (14), для игрока г * 1 имеем Г Г- M V (|S|-!)!("-|£|)! Г сл с rui

Ы'Фи] = L -;-1Фь5)-Фь^4W)J =

SsN ;iaS n-

= mii L -;-= mii L-;- K.| о =

SSN;1,iaS n! |S|=2 n! MS| —2/

n |S| — 1 m1i = mli L —(-77 = ~7T-

Снова, следуя свойству эффективности вектора Шепли, можно легко найти компоненту вектора Шепли для игрока 1:

viK^i)] = v{zi;N) - Е ViK^i)] = \ Е Е тп - \ Е т« =

i*1 2 iaN jaNi 2 i*1

= Emii- ^Emii = ^Emii-i*1 2 i*1 2 i*1

Напомним, что по (6)—(8) v(z\; {г}) = 0 для любого игрока г € N, выражение (16) справедливо для двухшаговой игры Г21 (g).

Пример. Рассмотрим игру 4 лиц, N ={1, 2,3,4}. С учетом принятых выше предположений игрокам в первом состоянии с целью максимизации прибыли выгодно выбрать поведения b1 = (0,1,1,1), b2 = b1 = b1 = (1,0,0,0). Получим сеть-звезду, приведенную на рис. 2.

Рис. 2. Сеть-звезда для 4 игроков Вестник СПбГУ. Прикладная математика. Информатика... 2019. Т. 15. Вып. 1 153

Заданы следующие матрицы:

(3;2) (1;4)) . с ((2;4) (1;3)) . с ((3;2) (1;4)) (2;0) (3; 1) ), ^13°31 = ((2;3) (3;0)) Л14°41 = ((5;0) (2;2)).

Значения т^ и т3: т.12 = 5, т1з = 6, т.14 = 5, т12 = 2, т21 = 1, т1з = 2, тз1 = 3, т>14 = 2, т41 = 2.

Вычислим характеристическую функцию для подыгры ГБ2(д), используя выражения (2)-(5):

у(г2; {0})= 0, у(г2; {1}) = 6, «(¿2; {2}) = 1, «(¿2; {3}) = 3, «(¿2; {4}) = 2, «(¿2; {12}) = 9, «(¿2; {13}) = 10, «(¿2; {14}) = 9, «(¿2; {23}) = 4, «(¿2;{34}) = 5, «(¿2;{24}) = 3, «(¿2;{234}) = 6, «(¿2;{123}) = 13, «(¿2; {124}) = 12, «(¿2; {134}) = 13, «(¿2; {1234}) = 16.

Далее определим характеристическую функцию для двухшаговой игры Г21 (д), применяя выражения (6)-(8):

«(¿1; {0}) = 0, «(¿1; {1}) = «(¿1; {2}) = «(¿1; {3}) = «(¿2; {4}) = 0, «(¿1; {12}) = 5, «(¿1; {13}) = 6, «(¿1; {14}) = 5, «(¿1; {23}) = «(¿1; {24}) = «(хх; {34}) = 0,

«(¿1; {234}) = 0, «(¿1; {123}) = 11, «(л; {124}) = 10, «(хх; {134}) = 11, «(¿1; {1234}) = 16.

Теперь вычислим значение вектора Шепли для подыгры Г Б (д) и двухшаговой игры Г:

—[«(¿2)] = (8, 2,3.5,2.5), -[«(¿1)] = (8, 2.5, 3,2.5).

Для подыгры ГБ2 (д) доказана выпуклость характеристической функции, из этого следует, что вектор Шепли принадлежит С-ядру. Проверим принадлежность вектора Шепли С-ядру для случая двухшаговой игры Г21 (д). Напомним, что С-ядро в игре Гг1 (д) — это множество дележей х = (х1,...,хп), которые удовлетворяют таким условиям:

^ «

геБ

N

^Хг = «

г=1

Подставляя значения характеристической функции (6)-(8) и вектора Шепли —[«(¿1)] в условия (17), (18), получим верные неравенства. Это означает, что в данной задаче вектор Шепли принадлежит С-ядру.

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

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

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

Автор благодарит рецензентов за ценные замечания. Литература

1. Dyer M., Mohanaraj V. Pairwise-interaction games // ICALP. 2011. Vol. 1. P. 159-170.

2. Acemoglua D., Ozdaglarb A., Parandeh Gheibib A. Spread of (mis)information in social networks // Games and Economic Behavior. 2010. Vol. 70, iss. 2. P. 194-227.

3. Konig M. D., Battistonb S., Napoletanoc M., Schweitzerb F. The efficiency and stability of R& D networks // Games and Economic Behavior. 2012. Vol. 75, iss. 2. P. 694-713.

4. Петросян Л. A., Седаков A. A. Многошаговые сетевые игры с полной информацией // Математическая теория игр и ее приложения. 2009. Т. 1, № 2. С. 66-81.

5. Bulgakova M. A., Petrosyan L. A. About strongly time-consistency of core in the network game with pairwise interactions // Proceedings of 2016 Intern. Conference "Stability and Oscillations of Nonlinear Control Systems". 2016. P. 157-160.

6. Kuzyutin D., Nikitina M. Time consistent cooperative solutions for multistage games with vector payoffs // Operations Research Letters. 2017. Vol. 45, iss. 3. P. 269-274.

7. Shapley L. S. Cores of convex games // Intern. Journal of Game Theory. 1971. Vol. 1. P. 11-26.

8. Shapley L. S. A value for га-person games // Contributions to the theory of games. Vol. II / eds by H. W. Kuhn, A. W. Tucker. Princeton: Princeton University Press, 1953. Vol. 28. P. 307-317.

9. Petrosyan L., Bulgakova M., Sedakov A. Time-consistent solutions for two-stage network games with pairwise interactions // Game Theory for Networking Applications. 2018. P. 15-23.

Статья поступила в редакцию 8 мая 2017 г.

^атья принята к печати 18 декабря 2018 г.

Контактная информация:

Булгакова Мария Александровна — аспирант; maribulgakova@mail.ru

Solutions of network games with pairwise interactions *

М. А. Bulgakova

St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation

For citation: Bulgakova M. A. Solutions of network games with pairwise interactions. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2019, vol. 15, iss. 1, pp. 147-156. https://doi.org/10.21638/11702/spbu10.2019.112

This article is devoted to cooperative network games with pairwise interaction. We consider a two-stage game, the first stage of which represents a network-formation stage, and the second is simultaneous bimatrix games, which take place between neighbours over the network. The characteristic function is constructed, its supermodularity is proved for the case of a one-step subgame starting with the second stage. For a special class of networks (star-network),

* The work is supported by Russian Found of Fundamental Research (grant N 17-51-53030).

a simplified formula for the Shapley vector is found, which does not require the calculation

of the values of the characteristic function over all coalitions, but only over coalitions of

dimension no more than two.

Keywords : cooperative games, convexity, Shapley value, characteristic function. References

1. Dyer M., Mohanaraj V. Pairwise-interaction games. ICALP, 2011, vol. 1, pp. 159—170.

2. Acemoglua D., Ozdaglarb A., Parandeh Gheibib A. Spread of (mis)information in social networks. Games and Economic Behavior, 2010, vol. 70, iss. 2, pp. 194—227.

3. Konig M. D., Battistonb S., Napoletanoc M., Schweitzerb F. The efficiency and stability of R& D networks. Games and Economic Behavior, 2012, vol. 75, iss. 2, pp. 694—713.

4. Petrosyan L. A., Sedakov A. A. Mnogoshagovye setevye igry s polnoy informatsiey [Multistage network games with full information]. Matematicheskaia teoria igr i ee prilozhenia [Mathematical game theory and applications], 2009, vol. 1, no. 2, pp. 66—81.

5. Bulgakova M. A., Petrosyan L. A. About strongly time-consistency of core in the network game with pairwise interactions. Proceedings of 2016 Intern. Conference "Stability and Oscillations of Nonlinear Control Systems", 2016, pp. 157—160.

6. Kuzyutin D., Nikitina M. Time consistent cooperative solutions for multistage games with vector payoffs. Operations Research Letters, 2017, vol. 45, iss. 3, pp. 269—274.

7. Shapley L. S. Cores of convex games. Intern. Journal of Game Theory, 1971, vol. 1, pp. 11—26.

8. Shapley L. S. A value for ra-person games. Contributions to the theory of games, vol. II. Eds by H. W. Kuhn, A. W. Tucker. Princeton, Princeton University Press, 1953, vol. 28, pp. 307-317.

9. Petrosyan L., Bulgakova M., Sedakov A. Time-consistent solutions for two-stage network games with pairwise interactions. Game Theory for Networking Applications, 2018, pp. 15-23.

Received: May 8, 2018.

Accepted: December 18, 2018.

Author's information:

Maria A. Bulgakova — Postgraduate Student; maribulgakova@mail.ru