Математик научила таксистов объезжать Манхэттен по ровному кругу и измерять число π

Математик Мишель Рудольф-Лилит из Национального центра научных исследований Франции описала особенности окружностей, начерченных в дискретном пространстве, в качестве примера которого ученый рассмотрела пересечения улиц и проспектов Манхэттена — центрального района Нью-Йорка. Оказалось, что можно аналитически описать несколько алгоритмов, следуя которым, гипотетический таксист проедет вдоль линии, максимально приближенной к идеальной окружности, а при достаточно большом ее радиусе можно с хорошей точностью измерить число π. Исследование в виде препринта выложено на arXiv.org.

В математике есть понятие «норма», простейший смысл которого — это расстояние (или, точнее, длина вектора) в некотором пространстве. В жизни мы чаще всего пользуемся евклидовой нормой, которая для плоскости очень похожа на теорему Пифагора. Чтобы найти расстояние от некоторого центра с координатами (0,0) до точки с координатами А (x,y), проведем в эту точку вектор OA. Квадрат его нормы ||OA||² = x² + y² и будет равен квадрату расстояния до точки А.

Но как быть, если мы находимся в дискретном пространстве? Например, как рассказывает Рудольф-Лилит, все в том же Манхэттене расстояние от одной точки до другой по прямой может быть 4,3 мили, но на самом деле кратчайший маршрут составит около 6 миль, ведь ехать можно только по улицам и проспектам, которые пересекаются под прямым углом. В этом случае пространство можно задать как множество точек, соответствующее перекресткам. Тогда любой маршрут будет задаваться линией, которую мы проведем по этим точкам, причем использовать диагонали нельзя.

Для такого пространства уже не имеет смысла использовать евклидову (или L2) норму, так как она не дает представления о реальной длине маршрута. Вместо этого подойдет норма L1 (ее, кстати, называют «манхэттенской»), которая для нашего же вектора OA будет определена как |OA| = |x| + |y|. У этой нормы есть ряд отличительных особенностей. Например, в евклидовом пространстве (с L2 нормой) для любой пары точек есть только один кратчайший маршрут — отрезок, соединяющий эти две точки. В дискретном пространстве вы такой отрезок провести не можете, и кратчайших маршрутов в этом случае будет множество. Например, можно сначала переместиться на x по проспекту, потом повернуть, и сместиться на y по улице, а можно двигаться зигзагом и сворачивать на каждом перекрестке влево-вправо. Легко убедиться, что пройденное расстояние будет одинаковым.

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

Прежде, чем перейти к выводу необходимых соотношений, Рудольф-Лилит напомнила, что результат может зависеть от того, как именно мы определим понятие окружности. Так, если следовать «школьному» определению (на самом деле, авторства Евклида), согласно которому окружность — это множество точек, лежащих на равном расстоянии от выбранного центра, то для дискретного пространства эта линия будет выглядеть очень непривычно: это будет квадрат, повернутый на 45 градусов. Тогда L1 норма векторов, проведенных во все его точки, действительно будет одинакова и если для такой «окружности» определить число π как отношение длины окружности к диаметру, оно будет равно 4, а не 3,14... .

Чтобы избежать непонимания, автор подчеркнула, что изучаться будет именно привычная нам окружность, которую для центра в (0,0) и радиуса R можно задать параметрически как множество точек, удовлетворяющих x² + y² = R². «Дискретная» окружность тогда будет задаваться как множество точек в дискретном пространстве, лежащих как можно ближе к «настоящей» окружности. Автор приводит широко используемые алгоритмы, в частности mid-point и  signum, которые давно и успешно используются для черчения таких дискретных окружностей. Рудольф-Лилит отмечает, что «таксисту будет несложно прокатить пассажира „по кругу“, в уме вычисляя, куда повернуть на том или ином перекрестке».

С точки зрения математика гораздо интереснее понять, можно ли с помощью построения такой окружности точно вычислить знаменитое число π. Кажется, что при увеличении радиуса окружности, мы должны все точнее приближать π, однако в математике известны похожие задачи, где в любом бесконечном пределе все равно не удается приблизить одно расстояние другим: сколько бы перекрестков не было по пути, вам никогда не превратить 6 миль в 4,3, так уж устроено дискретное пространство.

Тем не менее, автор показывает, что с увеличением радиуса воображаемой окружности точность определения числа π через отношение длины окружности к диаметру растет. Для окружности радиусом в десять перекрестков можно определить π с точностью до третьего знака после запятой, а для радиуса в 1000000 перекрестков — до девятого знака. Рудольф-Лилит не только сделала численное моделирование, но получила необходимые аналитические выражения, которые показывают справедливость ее утверждений в предельном переходе к бесконечному радиусу.

Математик в заключении признается, что удивлена, что несмотря на особенности дискретного пространства улиц Нью-Йорка, удается воспроизвести π — фундаментальную константу для евклидовой окружности. «Теперь, сидя в такси, мы будем не только слушать размышления водителя вслух, но и гадать, какие еще сюрпризы нам приподнесет город, который никогда не спит». 

Число π — один из любимых объектов математиков. Эту иррациональную константу можно измерить многими, часто неординарными способами: например, бросая иголки на стол или стреляя 200 раз из дробовика Mossberg 500. А четырнадцатого марта (03/14 в формате дат, принятом в США) по всему миру традиционно отмечается день числа π. 

Тарас Молотилин