2024 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 11 (69). Вып. 1
МАТЕМАТИКА. МЕХАНИКА. АСТРОНОМИЯ
МАТЕМАТИКА
УДК 519.2 МБС 60F99
Непрерывный вариант задачи об эгоистичной парковке
С. М. Ананьевский1, А. П. Чен2
1 Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
2 Университет Западной Онтарио,
Канада, Онтарио, Лондон, ул. Ричмонд, 1151
Для цитирования: Ананьевский С. М., Чен А. П. Непрерывный вариант задачи об эгоистичной парковке // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1. С. 84-95. https://doi.org/10.21638/spbu01.2024.104
Работа посвящена исследованию новой модели случайного заполнения отрезка большой длины интервалами меньшей длины. Рассмотрена новая постановка задачи. Изучается модель, в которой единичные интервалы размещаются на отрезке только в том случае, если заполняемый отрезок имеет длину не меньше 2. При этом положение размещаемого интервала подчинено равномерному закону распределения. В работе исследуется поведение среднего числа размещенных интервалов в зависимости от длины заполняемого отрезка. Получено точное выражение для аналога константы Реньи. Ключевые слова: случайное заполнение отрезка, задача о парковке, асимптотическое поведение математического ожидания.
Задача случайного заполнения отрезка интервалами, также известная как задача о парковке, имеет давнюю историю, которая берет свое начало в 1958 г. Именно тогда вышла работа Реньи [1], в которой была исследована асимптотика среднего числа случайно размещенных интервалов единичной длины на отрезке, длина которого бесконечно увеличивается.
В работе Реньи постановка задачи была следующей. На отрезке [0, х] для некоторого фиксированного х > 1 случайным образом размещается интервал + 1), который разбивает первоначальный отрезок на два отрезка [0, ¿] и [£ + 1, х], которые в дальнейшем продолжают заполняться по тому же правилу. Размещение интервала
© Санкт-Петербургский государственный университет, 2024 84 https://doi.org/10.21638/spbu01.2024.104
(£,£ + 1) случайным образом означает, что £ является случайной величиной с равномерным законом распределения на отрезке [0, х — 1]. Если какой-то из отрезков [0,£] или + 1, х] будет иметь длину меньше 1, то он в дальнейшем не заполняется. Указанный способ заполнения продолжается до тех пор, пока длина каждого из всех незаполненных отрезков станет меньше 1. После этого подсчитывается количество размещенных интервалов. Его мы обозначим Ых. Если в начале процесса заполнения отрезка х < 1, то Ых = 0.
Реньи предложил для этой задачи следующую интерпретацию. На улице длины х случайным образом паркуются автомобили единичной длины. Тогда Ых означает число запаркованных на этой улице автомобилей. Эта интерпретация и дала название задаче как задача о парковке.
В работе Реньи [1] было показано, что при любом п > 1
Е[Мх] = Ах + Л — 1 + 0(х-п), (х ^ (1)
Для константы А также было получено следующее выражение:
' -и
Г -2 [ йи
А = / е о ¿г. (2)
В последние годы появился ряд работ по теме задачи о парковке (см. [2-9]), в которых рассматриваются различные модели этой задачи, отличные от классической постановки. В нескольких работах изучаются дискретные аналоги этой задачи (см., например, работы [2, 7-9]). В работе [7] изучается дискретная модель, в которой длина заполняемого отрезка х — целое число, а £ — случайная величина, принимающая с равными вероятностями значения 0,1,...х — 1. При этом процесс размещения интервалов заканчивается, когда длина всех незаполненных отрезков становится меньше 2. Эта модель была названа задачей об эгоистичной парковке. Так же, как и в классической задаче, Ых означает число размещенных интервалов на отрезке длины х. В этой модели выполняются равенства
2х — 1
ЕНХ = —-— при х > 2 и Ых = 0 при х < 2. (3)
В этой работе мы рассматриваем новую модель задачи о парковке, которую можно считать аналогом задачи об эгоистичной парковке в непрерывном случае.
На отрезке [0, х] для некоторого фиксированного х > 1 случайным образом размещается интервал (£,£ + 1), который разбивает первоначальный отрезок на два отрезка [0,£] и [£ + 1,х], которые в дальнейшем продолжают заполняться по тому же правилу. Размещение интервала (£,£ + 1) случайным образом означает, что £ является случайной величиной с равномерным законом распределения на отрезке [0, х — 1]. Если какой-то из отрезков [0, £] или [£ + 1,х] будет иметь длину меньше 2, то он в дальнейшем не заполняется. Указанный способ заполнения продолжается до тех пор, пока длина всех незаполненных отрезков окажется меньше 2. После этого подсчитывается количество размещенных интервалов. Его мы будем обозначать за Мх. Если в начале процесса заполнения длина первоначального отрезка х < 2, то Мх = 0.
Теорема. Для определенной выше случайной величины Мх справедливо соотношение з
ЕМх = кх + к-1 + 0^^ (4)
где
ж / ь
1 — ехр(— и)
СО ft
к = J (t + 1) exp I -t - 2 J ■
0 \ 0
du\dt& 0.531417.
При доказательстве теоремы мы будем придерживаться схемы, предложенной в работах Дворецкого и Роббинса [10] и Реньи [1]. Докажем вначале следующую лемму.
Лемма. Пусть /(х) — функция, определенная для х > 0 и удовлетворяющая условию
х
2
/(х + 1) = - [ №<й+р(х+1), (5)
х
о
где р(х) — непрерывная функция, определенная для х > 1. И пусть выполнено
22 т<~.
i=2
где через px обозначено
Px = sup |p(i)|.
x<t<x+1
Тогда существует константа X такая, что
2n 2n n j!
sup \f(x) - Xx - X\ < — sup \f(x) - Xx - X\ +— У—Rj,
x
<n+2 n! к x< 2 n! . 23
- - - - 3=1
где
2j +1 , 2(j + l)(j + 3) ^ Pi
= —p+---
J J i=j+2
Доказательство леммы. Из равенства (5) для положительных x, y и y > x имеем
x y y
2 i' 2 f 1 2 f f(y+1) = - / f(t)dt+- / /(¿)<Й+р(г,+ 1) = -(х/(х+1)-а^(х+1))+- / /(*)<Й+р(гИ) У J У J y y J
0
x 2 { s
f(y + 1) = ~f{x + 1) + - / f(t)dt +P(y + 1) - -p(x + 1). (6)
y y y
x
Введем следующие обозначения:
г ■ f fit) я fit)
Ix= Jrß ГТТ' S*= SUP гтт-
x<i<x+l t + 1 x<t<x+1 t +1
u
или
Можно легко увидеть, что функция ](х) = х + 1 удовлетворяет условию (5) с р(х) = 0, и, следовательно,
у
х2
у + 2=-(х + 2) + - (7)
У У .1
х
Вычтем из равенства (5) равенство (7), умноженное на 1х, и получим
у
х 2 х
/(у+1)-1х-(у+2) = ~(/(х+1)-1х-(х+2))+- / (т-1х-(г+1))(М+р(у+1)—р(х+1).
У У У
х
Следовательно, в силу определения 1х и свойств функции р(х) для х < у < х + 1 /(у + 1) — 1х(у + 2) > 0 + 0 — рх+1 — Рх+1 = —2рх+1.
Из чего следует
т > т 2рх+1
Вводя обозначение
Ах = 2
ОО
, Px+i
f-f x + i + 1
i=1
и последовательно заменяя x в неравенстве выше на x + 1,x + 2,..., получаем, что
Iy > Ix Ax.
Следовательно, справедливо неравенство
lim inf Iy > Ix — Ax.
y—ж
Из определения Ax и из условия леммы имеем Ax = о(1). Поэтому справедливо неравенство
lim inf Iy > lim sup Ix.
у — ж x—ж
Вспоминая неравенство Iy > Ix — Ax и полагая x =1, делаем вывод:
3 lim Ix = , > —ю.
x—ж
Аналогичным образом получаем соответствующие утверждения для Sy:
Sy < Sx + Ax,
3 lim Sx = ЬЖ, Sx < ю.
x—^^O
Так как Ix < Sx, то выполнены неравенства
ю < ^ < ^^ж < ю. (8)
Вспоминая равенство (6), получаем для x > 0 и y > 0
y
f(y+l)-f(x+l) = ^-f(x+l) + - [ f(t)dt+p(y + l)--p(x + l).
y y y
x
В силу неравенства (8) имеем f (x) = O(x).
Поэтому из последнего равенства следует соотношение
sup f (y +1) - f (x +1)| = O(1) + 2px,
x<y<x+l
w
из которого следует (учитывая ^ < оо ), что
i=2 1
Sx — Ix = 0(1) при x
Поэтому
Iw = Sw =
Теперь определим константу:
к = lim Ix = lim Sx = lim ^ ^ .
x^w x^w x^w x +1
Применяя уже полученные неравенства, можем записать
Ix — Ax < к < Sx + Ax. Заметим, что для любого x > 1 найдется x' такое, что
f (x')
x < x < x + 1,
— к
< Ax
x' + 1
Это верно, поскольку в противном случае выполнялось бы одно из неравенств
Ix >к + ДХ или Sx <к — Ах,
что противоречит (8).
Для каждого x = n такое x' обозначим через xn. Для n > 2 имеем
\f (xn) - фп + 1)\<(n + 2)Дп, (n < xn < n +1). (9)
Введем функцию f *(x) = f (x) — k(x + 1). Заметим, функция f * удовлетворяет условию леммы. Применяя равенство (6) для n < y < n+1 и x = xn+i — 1 и равенство (9), получаем
n+1 2 n+1 \f*(y + 1)1 < —^—(n + 3)Д„+1 + - sup \f*(t)\+pn+1 + -!—pn+1.
n n n<t<n+1 n
Введем обозначение Tx = sup \f *(i)\ и преобразуем последнее неравенство:
x<t<x+1
т - . 2n+1 , (n+l)(n + 3)A 2
j-n+1 S —J-n H--Pn+1 H--^n+1 — —J-n + rt-n,
n n n n
где Rn удовлетворяет условию леммы.
Последовательно применяя данное неравенство для натуральных п, получаем требуемое утверждение:
2п 2п /1! 2! п! Тпц--\ < —-Т-] —- [ —Д-1 + —ттИ1? +... + —ап + п\ п\ V 2 22 2™
□
Доказательство теоремы. Рассмотрим отрезок [0, х+1] и предположим, что первый случайно расположенный интервал занял место + 1). Заметим, что число интервалов по окончании процесса заполнения отрезка, расположенных слева от первого интервала — М4, и число интервалов по окончании процесса заполнения отрезка, расположенных справа от первого интервала — Мх-^ являются независимыми случайными величинами при фиксированном значении Ь. Общее число размещенных интервалов будет равно М.1 + Мх^ + 1. Таким образом,
Е (Мх\г)= ЕМг + ЕМ— + 1, (0 < Ь < х).
Введем следущее обозначение: ¡л(х) = ЕМх; учитывая, что случайная величина Ь равномерно распределена на [0,х], имеем равенство
х
ц{х + 1) = - [ /лШЛ + 1. х
Определим функцию I(х) = ¡л(х) + 1. Тогда последнее равенство примет вид
х
/(х + 1) = - [ №<и. х
Начальные данные для функции I(х) будут следующими:
I (х) = 1 (0 < х < 2), и I(2) = 2.
Мы видим, что функция I (х) последовательно определена на промежутках 2 < х < 3, 3 < х < 4,...:
2 < х < 3 :
3 < х < 4 :
I (х)
с-1
№ = / №№ + / №№ = 4 "
х1
4 < х < 5 :
I (х)
с-1
х1
I(1)31 + 11(1)31 + у I(г)зг
8 -
16
х 1 х 1
8 , (х - 2 ■ 1п
и т. д.
2
4
3
2
Воспользуемся приведенной леммой для получения оценки ¡j,(x). Как мы уже отметили, функция f (x) = ¡j,(x) + 1 удовлетворяет равенству (5) с p(x) = 0. Поэтому применима лемма и
M(x)
Зк = lim -.
х^ж x
И в силу равенства (8) получаем для любых x > 0
. ju(t) + 1 ju(t) + 1 mf -= lx < к < Ьх = sup -.
x<i<x+l t +1 x<t<x+1 t + 1
Вспомним, что при 3 < x < 4
2 x — 1
2
= / №№ + J №№ = 4 "
x — 1 W J I x — 1
\0 2
Поэтому в случае x = 3 имеют место неравенства
0.5 < к < 0.5359
Отметим, что можно получить и более точную оценку, если положить x = 4:
16 8 , (x — 2
4 < ж < 5 : f(x) = 8----In
x— 1 x— 1
тогда
0.5306 < к < 0.5333.
Из равенства ц(х) = 0 при 1 < х < 2 и из даже более грубой оценки ^ < к < | следует
sup l^-(x) + 1 — kx — к| = max |1 — kx — к| = |1 — 3к| < 1.
1<x<2 1<ж<2
Итак, вновь воспользовавшись леммой для p(x) = 0, получаем, что
1 2" Зк, - < к < 1 : sup |/1,(х) + 1 — кх — к| < —-.
2 n+1<x<n+2 n!
По формуле Стирлинга о приближении факториала справедливо соотношение
(/2е\ж^\ /1,(х) — кх — к + 1 = О I ( — J I,
что частично доказывает теорему.
Далее укажем явный вид константы к. У нас имеется интегральное уравнение
x
ц(х+ 1) = - f n(t)dt+ 1 (10)
x
0
с начальными данными
при 0 < x < 2 yU,(x) = 0, при 2 < x < 3 ¡j,(x) = 1,
при 3 < х < 4 /л(х) = 3 —
x — 1
16 8 , /x - 2\
при 4 < ж < 5 /х(ж) = 7----In - . (11)
x- 1 x-1 2
Найдем
lim
х^ж x
Умножим обе части равенства (10) на x, продифференцируем по x и получим дифференциальное уравнение:
x^'(x + 1) + n(x + 1) = 2^(x) + 1 при x > 1. (12)
Рассмотрим преобразование Лапласа функции ^(x):
ж
ф(а)=!Kx)e-sxdx. 0
Преобразование Лапласа существует, поскольку для функции ¡j,(x) выполнены неравенства: 0 < ¡j,(x) < x. Умножая обе части равенства (12) на e-sx и интегрируя от 1 до то, приходим к следующему равенству:
ж ж
J XI)!{х + 1 )e-sxdx + J ц{х + l)e-sx = 2ф{а) + —. (13)
1 1
Учитывая начальные данные (11), перепишем каждое слагаемое в уравнении (13):
ж ж
j*x +1)e-sxdx = esj = ✓ (14)
1
/^" = ' l/" + 1У" + 1KS<X+"dx^"" S<X+1,dxj =
= es ( / tf)'(t)e-stdt - I f)'(t)e-stdt | = ~ ( es / /л'(t)e-stdt | ,
то есть
ж
A. i rs 1
ds
I хц'{х + l)e-sxdx = ~ ^е8 ^ . (15)
Интегрируя по частям правую часть равенства (15), получаем
ж
J ц'(г)в-згзг = -е-2я + вф(в). (16)
2
и
Из равенств (13)—(16) следует, что функция ф(в) удовлетворяет следующему дифференциальному уравнению:
а
-(зе>ф(з)) = ф(з)(е>- 2)-е~а--.
ав в
Введем новое обозначение ш(в) = еяф(в). Тогда последнее уравнение перепишется в следующем виде:
е~а
8ш'{8) = -2ш(в)е-8 -е~а--. (17)
в
Это дифференциальное уравнение относительно ш можно решить, например, методом Лагранжа. Прежде чем выписывать решение уравнения (17), заметим, что в силу неравенства 0 < ц(х) < х и условия ¡л(х) = 0 при 0 < х < 1 выполнено соотношение
сю
1 1
11
О < Cü(s) < es / xe-sxdx = — + J s2 s
из которого следует, что
lim w(s) = 0.
в^ж
Тогда решение уравнения (17) имеет вид
cj(s) = j(t + 1) ехр ( -t - 2 [ ---—du | dt.
Следовательно,
ф{й) = ^-у- I (t + 1) ехр I -t - 2 [ ---—du I dt,
из чего следует
lim s2ф(s) = C0,
s^+0
где
сю / Ь
С0 = J(t+l)exp{-t-2 J —-—du I dt « 0.531417.
о V о )
Далее воспользуемся тауберовой теоремой: если функция а(х) монотонно возрастает при х > 0 и если при некотором в > 0
ж
lim sß i e—sxda(x) = C,
в^+0 J
тогда
lim ^ = er.
Функция
a(x) = j j(t)dt
0
монотонно возрастает, и
ж
aß e-SXda(X) = А
'i e--Xda(X)=
0
Поэтому согласно тауберовой теореме при ß = 2
a(x) к lim —г— = —.
х^ж x2 2
Разделив обе части уравнения (10) на x +1 и устремив x к бесконечности, получаем, что
х^ж x +1
или
к = lim ^ = С0 « 0.531417,
х^ж x
что полностью доказывает теорему.
Приведем сравнение поведения среднего числа разместившихся единичных интервалов на отрезке большой длины для модели эгоистичной парковки, рассматриваемой в работе [7], и для модели, рассматриваемой в настоящей работе. Напомним, в работе [7] было установлено, что
22
ENX = ¡j,(x) = —х + — — 1 при натуральных х > 2,
и значение константы А = |. В настоящей работе
EMx=fj,(x) = Kx + K-l + 0^j^j j,
где к « 0.531417.
Таким образом, мы можем сделать вывод, что в непрерывном случае случайное размещение единичных интервалов будет менее плотным или, что то же самое, более разреженным.
Литература
1. Renyi A. On a one-dimensional problem concerning space-filling. Publ. of the Math. Inst. of Hungarian Acad. of Sciences 3, 109—127 (1958).
2. Pinsky R. G., Problems from the Discrete to the Continuous. Springer International Publishing Switzerland. Chapter 3, 21-34 (2014).
3. 1льенко А. Б., Фатенко В. В. Узагальнення задачi Рень! про паркування. Науковi вiстi НТУУ КПЬ- мiжнародний науково-техтчний журнал 4 (114), 54-60 (2017).
4. Clay M.P., Simanyi N.J. Renyi's parking problem revisited. 29 Dec. 2014. ArXiv:1406.1781v1 [math. PR].
5. Ананьевский С. М. Задача парковки для отрезков различной длины. Записки научных семинаров ПОМИ 228, 16-23 (1996).
6. Ананьевский С.М., Крюков Н. А. Об асимптотической нормальности в одном обобщении задачи Реньи. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 6 (64), вып. 3, 353-362 (2019). https://doi.org/10.21638/11701/spbu01.2019.301
7. Ананьевский С.М., Крюков Н.А. Задача об эгоистичной парковке. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 5 (63), вып. 4, 549-555 (2018). https://doi.org/10.21638/11701/spbu01.2018.402
8. Ананьевский С.М. Некоторые обобщения задачи о «парковке». Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 3 (61), вып. 4, 525-532 (2016). https://doi.org/10.21638/11701/spbu01.2016.401
9. Ананьевский С. М., Чен А. П. Обобщение задачи об эгоистичной парковке. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 9 (67), вып.3, 464-472 (2022). https://doi.org/10.21638/spbu01.2022.307
10. Dvoretzky A., Robbins H. On the "parking' problem. Publ. of the Math. Inst. of Hungarian Acad. of Sciences 9, 209-226 (1964).
Статья поступила в редакцию 31 марта 2023 г.;
доработана 28 апреля 2023 г.; рекомендована к печати 31 августа 2023 г.
Контактная информация:
Ананьевский Сергей Михайлович — канд. физ.-мат. наук, доц.; ananjevskii@mail.ru Чен Александр Петрович — ассистент преподавателя; sasha.24chen@mail.ru
A continuous version of the selfish parking problem
S. M. Ananjevskii1, A. P. Chen2
1 St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
2 Western University,
1151, ul. Richmond, London, Ontario, Canada
For citation: Ananjevskii S. M., Chen A. P. A continuous version of the selfish parking problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 1, pp. 84-95. https://doi.org/10.21638/spbu01.2024.104 (In Russian)
The work is devoted to the study of a new model of random filling of a segment of large length with intervals of smaller length. A new formulation of the problem is considered. We study a model in which unit intervals are placed on a segment only if the segment being filled has a length of at least 2. In this case, the position of the placed interval is subject to a uniform distribution law. The paper investigates the behavior of the average number of placed intervals depending on the length of the filled segment. An exact expression is obtained for the analog of the Renyi constant.
Keywords: random filling of a segment, parking problem, asymptotic behavior of the expectation.
References
1. Renyi A. On a one-dimensional problem concerning space-filling. Publ. of the Math. Inst. of Hungarian Acad. of Sciences 3, 109-127 (1958).
2. Pinsky R. G., Problems from the Discrete to the Continuous. Springer International Publishing Switzerland. Chapter 3, 21-34 (2014).
3. Ilienko A.B., Fatenko V. V. A Generalization of Renyi's Parking Problem. Research Bulletin of the National Technical University of Ukraine. Kyiv Politechnic Institute. Physics and Mathematics 4 (114), 54-60 (2017). (In Ukrainian)
4. Clay M.P., Simanyi N.J. Renyi's parking problem revisited. 29 Dec. 2014. ArXiv:1406.1781v1 [math. PR].
5. Ananjevskii S. M. The parking problem for segments of different length. Journal of Mathematical Sciences. Zapiski nauchnyh semitnarov POMI 228, 16—23 (1996). (In Russian). [Eng. transl.: J. Soviet Math. 93, 125-130 (1999)].
6. Ananjevskii S.M., Kryukov N.A. On Asymptotic Normality in One Generalization of the Renyi Problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 6 (64), iss.3, 353-362 (2019). https://doi.org/10.21638/11701/spbu01.2019.301 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 52, iss.3, 227-233 (2019). https://doi.org/10.21638/11701/spbu01.2019.301].
7. Ananjevskii S. M., Kryukov N.A. The problem of selfish parking. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 5 (63), iss. 4, 549-555 (2018). https://doi.org/10.21638/11701/spbu01.2018.402 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 51, iss.4, 322-326 (2018). https://doi.org/10.3103/S1063454118040039].
8. Ananjevskii S.M. Some generalizations of "parking" problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 3 (61), iss.4, 525-532 (2016). https://doi.org/10.21638 /11701/spbu01.2016.401
9. Ananjevskii S. M., Chen A. P. Generalization of the selfish parking problem. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 9 (67), iss.3, 464-473 (2022). https://doi.org/10.21638/spbu01.2022.307 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 55, iss.3, 290-296 (2022). https://doi.org/10.1134/S1063454122030025].
10. Dvoretzky A., Robbins H. On the "parking' problem. Publ. of the Math. Inst. of Hungarian Acad. of Sciences 9, 209-226 (1964).
Received: March 31, 2023 Revised: April 28, 2023 Accepted: August 31, 2023
Authors' information:
Sergey M. Ananjevskii — ananjevskii@mail.ru Alexander P. Chen — sasha.24chen@mail.ru