В самом начале Книги I «Начал» Евклид Александрийский дал два основополагающих определения: «Точка есть то, что не имеет частей; линия же — длина без ширины (\(\mu\tilde{\eta}\kappa\omicron\varsigma\ \dot{\alpha}\pi\lambda\alpha\tau\acute{\varepsilon}\varsigma\)).» Согласно нашей повседневной геометрической интуиции, идеальная математическая кривая имеет строго нулевую толщину (\(\delta = 0\)). На протяжении более двух тысячелетий казалось самоочевидным, что одномерная нить нулевой ширины никогда не сможет без просветов заполнить внутренность двумерного квадрата \([0,1]^2\).
Однако в 1890 году итальянский математик Джузеппе Пеано (Giuseppe Peano, 1858–1932) разрушил эту двухтысячелетнюю уверенность в пятистраничной статье в Mathematische Annalen, построив непрерывное сюръективное отображение \(f: [0,1] \to [0,1]^2\) единичного отрезка на единичный квадрат. Год спустя (1891) Давид Гильберт (David Hilbert, 1862–1943) предложил прозрачную геометрическую рекурсию — кривую Гильберта. Ещё более удивительно то, что объект, который в конце XIX века называли «патологическим монстром» чистой математики, столетие спустя стал одним из важнейших алгоритмов современной информатики, лежащим в основе Google Maps (библиотека Google S2), пространственных баз данных (Hilbert R-Tree), кэш-независимых (cache-oblivious) алгоритмов и балансировки нагрузки суперкомпьютеров.
1. Шок Кантора (1877) и теорема Нетто (1879)
В 1877 году создатель теории множеств Георг Кантор обнаружил, что между точками одномерного отрезка \([0,1]\) и точками двумерного квадрата \([0,1]^2\) существует биекция (взаимно однозначное соответствие). Потрясённый собственным доказательством, он написал Рихарду Дедекинду знаменитую фразу: «Je le vois, mais je ne le crois pas!» (Я вижу это, но не верю этому!)
Идея Кантора состояла в чередовании (interleaving) десятичных цифр (или неделимых блоков цифр, оканчивающихся ненулевой цифрой, чтобы строго устранить неоднозначность конечных дробей вида \(0{,}5 = 0{,}4999\dots\)):
\[x = 0.a_1 a_2 a_3 \dots, \quad y = 0.b_1 b_2 b_3 \dots \quad \Longrightarrow \quad t = \phi(x, y) = 0.a_1 b_1 a_2 b_2 a_3 b_3 \dots\]\[|[0,1]| = |[0,1]^2| = |[0,1]^n| = |\mathbb{R}| = \mathfrak{c} = 2^{\aleph_0}\]Тем самым было доказано, что отрезок \([0,1]\), квадрат \([0,1]^2\) и \(n\)-мерный гиперкуб \([0,1]^n\) имеют одну и ту же мощность континуума \(\mathfrak{c} = 2^{\aleph_0}\). Однако Дедекинд сразу заметил, что биекция Кантора всюду разрывна. В 1879 году Ойген Нетто доказал, что между \([0,1]\) и \([0,1]^2\) не существует непрерывной биекции (гомеоморфизма): так как \([0,1]\) компактно, а \([0,1]^2\) хаусдорфово, любая непрерывная биекция была бы гомеоморфизмом; но при удалении внутренней точки \(t_0 = 1/2\) множество \([0,1]\setminus\{1/2\} = [0,1/2) \cup (1/2,1]\) становится несвязным, тогда как проколотый квадрат \([0,1]^2 \setminus \{f(1/2)\}\) остаётся линейно связным — противоречие!
2. Кривая Пеано (1890): отказ от инъективности открывает путь к непрерывной сюръекции!
Пеано увидел логическую лазейку в теореме Нетто: невозможно лишь одновременное выполнение непрерывности, сюръективности и инъективности. Если отказаться от инъективности (разрешив самопересечения, т. е. кратные точки), то можно построить непрерывную сюръекцию \(f: [0,1] \to [0,1]^2\)!
Пеано использовал троичное разложение \(t = 0{,}t_1 t_2 t_3 \dots_3\) (\(t_m \in \{0,1,2\}\)) и оператор дополнения \(k(d) = 2 - d\):
\[t = \sum_{m=1}^{\infty} \frac{t_m}{3^m} \;\; (t_m \in \{0,1,2\}), \qquad k(d) = 2 - d, \qquad k^n(d) = \begin{cases} d & (n \text{ even}) \\ 2 - d & (n \text{ odd}) \end{cases}\]\[x_n = k^{t_2 + t_4 + \cdots + t_{2n-2}}(t_{2n-1}), \qquad y_n = k^{t_1 + t_3 + \cdots + t_{2n-1}}(t_{2n})\]Поскольку \(k^n(d)\) зависит только от чётности показателя \(n\), в точках неоднозначности троичной записи (например, \(t = 1/3 = 0{,}0222\dots_3 = 0{,}1000\dots_3\)) смена \(t_1 = 0\) (чётное) на \(t_1 = 1\) (нечётное) включает оператор \(k^1\), который обращает все последующие нули в двойки и идеально склеивает оба предела в точке \((1/3, 1)\)!
3. Геометрическая рекурсия Гильберта (1891) и строгие аналитические доказательства
В 1891 году Давид Гильберт заменил разбиение Пеано \(3 \times 3\) двоичной рекурсией на \(2 \times 2 = 4\) квадрантах: на шаге \(n\) отрезок \([0,1]\) делится на \(4^n\) замкнутых подотрезков \(I_{n,k}\) длины \(4^{-n}\), а квадрат \([0,1]^2\) — на \(4^n\) малых квадратов \(Q_{n,k}\) со стороной \(2^{-n}\), удовлетворяющих Аксиоме A (Смежность: \(Q_{n,k}\) и \(Q_{n,k+1}\) всегда имеют общую сторону) и Аксиоме B (Вложенность: \(Q_{n,k} = \bigcup_{j=0}^3 Q_{n+1, 4k-j}\)).
flowchart TD
A["Единичный отрезок I = [0,1] и единичный квадрат Q = [0,1]²"] --> B["Разбиение шага n: 4ⁿ отрезков I_{n,k} (длина 4⁻ⁿ)<br/>и 4ⁿ подквадратов Q_{n,k} (сторона 2⁻ⁿ, диагональ √2·2⁻ⁿ)"]
B --> C["Аксиома A (Смежность): Q_{n,k} и Q_{n,k+1} всегда имеют общее ребро"]
B --> D["Аксиома B (Вложенность): При m > n и t ∈ I_{n,k} точки H_n(t), H_m(t) ∈ Q_{n,k}"]
C --> E["1/2-Гёльдерова непрерывность: ||H(t₁) - H(t₂)||₂ ≤ √6 · |t₁ - t₂|^{1/2}"]
D --> F["Равномерная оценка Коши: sup_t ||H_m(t) - H_n(t)||₂ ≤ √2 · 2⁻ⁿ → 0"]
F --> G["Равномерная сходимость в полном пространстве C([0,1], ℝ²): H(t) = lim_{n→∞} H_n(t)"]
G --> H["Теорема 1 (Непрерывность): По теореме Вейерштрасса H(t) всюду непрерывна"]
G --> I["Теорема 2 (Сюръективность): Вложенные отрезки и компактность ⇒ H([0,1]) = [0,1]²"]
H --> J["Теорема 3 (Неинъективность): Кратные точки, напр. H(1/6)=H(1/2)=H(5/6)=(1/2,1/2)"]
E --> K["Теорема 4 (Нигде не дифференцируемость): Разностное отношение ≥ 2^{n-1} → ∞"]
- Теорема 1 (Равномерная сходимость и непрерывность): Для любых \(m > n \ge 1\) и \(t \in I_{n,k}\) по Аксиоме B точки \(H_n(t)\) и \(H_m(t)\) лежат в одном подквадрате \(Q_{n,k}\) диаметра \(\sqrt{2}\cdot 2^{-n}\):
Следовательно, \(\{H_n\}\) — фундаментальная последовательность в банаховом пространстве \(\mathcal{C}([0,1], \mathbb{R}^2)\), и её равномерный предел \(H(t)\) непрерывен на \([0,1]\). 2. Теорема 2 (Сюръективность): Для любой точки \(p \in [0,1]^2\) система вложенных подквадратов \(Q_{n,k_n} \ni p\) задаёт систему вложенных отрезков \(I_{n,k_n}\) длины \(4^{-n} \to 0\) с единственной общей точкой \(t^* \in [0,1]\), где \(H(t^*) = p\):
\[\bigcap_{n=1}^{\infty} I_{n, k_n} = \{t^*\}, \qquad \|H(t^*) - p\|_2 \le \frac{\sqrt{2}}{2^n} \xrightarrow{n \to \infty} 0 \quad \Longrightarrow \quad H([0,1]) = [0,1]^2\]- Теорема 3 (Кратные точки): Центр \((1/2, 1/2)\) проходится трижды (\(H(1/6) = H(1/2) = H(5/6) = (1/2, 1/2)\)), а в двоично-рациональных узлах сетки возникают четырёхкратные точки.
- Теорема 4 (Бесконечная длина, \(\frac{1}{2}\)-гёльдеровость и нигде не дифференцируемость): Длина ломаной \(L(H_n) = 2^n - 2^{-n} \to \infty\). Для любого \(t \in I_{n,k}\) хотя бы одна вершина \(v_n = H(s_n)\) квадрата \(Q_{n,k}\) удалена от \(H(t)\) минимум на \(2^{-(n+1)}\) при \(|s_n - t| \le 4^{-n}\), поэтому разностное отношение расходится как \(2^{n-1} \to \infty\) — \(H(t)\) не дифференцируема ни в одной точке \(t \in [0,1]\):
4. Разрешение кризиса размерности: топологическая размерность \(\dim_T\) против размерности Хаусдорфа \(D_H\)
В 1911 году Л. Э. Я. Брауэр доказал теорему об инвариантности размерности (\(\mathbb{R}^n \cong \mathbb{R}^m \iff n = m\)) через локальные группы гомологий. Большая индуктивная / лебегова размерность покрытия \(\dim_T\) определяется минимальной кратностью пересечения \(d + 1\) открытых покрытий (\(\dim_T([0,1]) = 1\), но \(\dim_T([0,1]^2) = 2\) подобно тройным стыкам кирпичной кладки, что топологически объясняет неизбежность минимум тройных точек самопересечения!). Напротив, размерность Хаусдорфа (фрактальная размерность) \(D_H\) измеряет масштабное заполнение:
\[H_k(U, U \setminus \{p\}; \mathbb{Z}) \cong \tilde{H}_{k-1}(S^{n-1}; \mathbb{Z}) \cong \begin{cases} \mathbb{Z} & (k = n) \\ 0 & (k \ne n) \end{cases}, \qquad \dim_T([0,1]) = 1 < \dim_T([0,1]^2) = 2\]\[D_H = \frac{\log N}{\log(1/r)} \quad \Longrightarrow \quad D_H(\text{Peano}) = \frac{\log 9}{\log 3} = 2, \qquad D_H(\text{Hilbert}) = \frac{\log 4}{\log 2} = 2\]Кривая Гильберта за счёт бесконечного складывания повышает свою фрактальную размерность до \(D_H = 2\), а за счёт самопересечений поднимает и топологическую размерность образа до \(\dim_T = 2\).
| Тип кривой / отображения | Непрерывность | Инъективность (без самопересечений) | Сюръективность на \([0,1]^2\) | Дифференцируемость | Топологическая разм. \(\dim_T(\text{Образ})\) | Разм. Хаусдорфа \(D_H\) | 2D мера Лебега \(\mu_2\) |
|---|---|---|---|---|---|---|---|
| Гладкая простая кривая (отрезок/дуга) | Да (\(C^\infty\)) | Да | Нет | Всюду дифференцируема | \(1\) | \(1\) | \(0\) |
| Кривая Коха (снежинка) | Да | Да | Нет | Нигде не дифференцируема | \(1\) | \(\frac{\log 4}{\log 3} \approx 1{,}262\) | \(0\) |
| Кривая Осгуда (1903) | Да | Да (жорданова дуга) | Нет (пустая внутренность) | Почти всюду не дифф. | \(1\) | \(2\) | \(> 0\) (положительная площадь!) |
| Чередование цифр Кантора | Нет (разрывна) | Да (биекция) | Да | Нигде не непр./дифф. | \(2\) | \(2\) | \(1\) |
| Кривая Пеано (база 3, 1890) | Да | Нет (кратные точки) | Да (заполняет квадрат) | Нигде не дифференцируема | \(2\) | \(\frac{\log 9}{\log 3} = 2\) | \(1\) |
| Кривая Гильберта (база 2, 1891) | Да (\(\frac{1}{2}\)-Гёльдер) | Нет (до 4-кратных точек) | Да (заполняет квадрат) | Нигде не дифференцируема | \(2\) | \(\frac{\log 4}{\log 2} = 2\) | \(1\) |
5. Приложения в компьютерных науках: индексация многомерного мира в 1D-памяти
Поскольку компьютерная память (адреса DRAM, блоки SSD, ключи B+-деревьев) одномерна, а геопространственные и физические данные многомерны, дискретная кривая Гильберта обеспечивает оптимальное сохранение локальности (Locality Preservation) благодаря \(\frac{1}{2}\)-гёльдеровости и квадродеревянной вложенности:
\[z = \text{Morton}(x, y) = (y_{n-1} x_{n-1} \dots y_1 x_1 y_0 x_0)_2, \qquad C_{\text{dilation}} = \sup_{t_1 \ne t_2} \frac{\|\gamma(t_1) - \gamma(t_2)\|_2^2}{|t_1 - t_2|} = 6 \;\; (\text{Hilbert})\]В отличие от кривой Мортона (Z-кривой), которая просто чередует биты координат и совершает длинные диагональные скачки на границах квадрантов (\(C_{\text{dilation}} = \infty\)), кривая Гильберта на каждом шаге \(t \to t+1\) смещается ровно на манхэттенское расстояние \(1\) (\(C_{\text{dilation}} = 6\)), сокращая число кластеров поиска в B+-дереве на \(25\%\text{–}40\%\).
- Google S2 Geometry Library (Google Maps, Uber, MongoDB): Проецирует земную сферу \(S^2\) на \(6\) граней куба, выравнивает площади ячеек нелинейным преобразованием и покрывает каждую грань 30-уровневой кривой Гильберта \(H_{30}\). Каждый участок Земли размером \(0{,}7\,\text{см} \times 0{,}7\,\text{см}\) кодируется одним 64-битным целым числом (
S2CellId), превращая двумерные геозапросы в молниеносное 1D-сканирование диапазонов в B+-деревьях. - Hilbert R-Tree и суперкомпьютерные расчёты (HPC): Сортировка полигонов по индексу Гильберта минимизирует перекрытие узлов R-дерева (PostGIS, Apache Sedona), а 3D-кривая Гильберта разрезает октодеревья в космологических симуляциях (Gadget-4, AMReX) для идеальной динамической балансировки нагрузки.