Аннотация:
В работе исследуется распределение суммарной длины хорд единичной окружности, соединяющих вершины правильного $q$-угольника с номерами $k$ и $ak^{2} \ (\operatorname{mod} q)$, $k=1,2,\dots, q$, в случае, когда $q$ – простое число, а величина $a$ пробегает полную систему вычетов по модулю $q$.
Библиография: 13 названий.
Ключевые слова:
хорда окружности, вычеты по модулю, суммы Гаусса, наименьший невычет, символ Лежандра, плотность распределения.
Зададимся целым числом $q\geqslant 5$ и рассмотрим правильный $q$-угольник, образованный корнями степени $q$ из единицы. Зафиксируем целое число $a$ и соединим хордой вершину $e(k/q)$ с вершиной $e(m/q)$, где $m\equiv ak\ (\operatorname{mod} q)$, а $e(u) = e^{2\pi iu}$. Проделав эту операцию для каждого $k=1,2,3,\dots, q$, получим некоторое множество хорд, которое будем называть конструкцией на единичной окружности, образованной вычетами по модулю $q$.
Такие конструкции обладают несомненной визуальной красотой (см. рис. 1) и в ряде случаев допускают красивые физические интерпретации. Их геометрические и арифметические свойства хорошо известны и нашли отражение в ряде статей, по большей части научно-популярных (см., например [1]–[7])1[x]1К сожалению, мне не удалось ознакомиться с работой [3].. Однако лишь недавно (2023 г., см. [8]) был отмечен тот факт, что суммарная длина $\mathscr{L}(q,a)$ хорд, отвечающих заданным $q$ и $a$, в случае, когда числа $a-1$ и $q$ взаимно просты, не зависит от $a$ и равна
Подобные конструкции можно рассматривать и в случае, когда линейное отображение $k\mapsto ak\pmod{q}$, сопоставляющее началу хорды с номером $k$ ее конец, заменяется другим: $k\mapsto af(k)\pmod{q}$. Визуальная красота таких конструкций, как правило, исчезает (см. рис. 2), но они становятся источником содержательных теоретико-числовых задач. Одна из них рассматривается в настоящей работе.
Именно, обозначим через $\mathscr{L}(f;q,a)$ суммарную длину хорд такой конструкции, отвечающей заданным числам $q$ и $a$. Вопрос, которому посвящена статья, состоит в следующем: как распределены величины $\mathscr{L}(f;q,a)$ при изменении $a$ в пределах $1\leqslant a\leqslant q$, в случае квадратичной функции $f(k)=k^{2}$? При этом мы ограничиваемся случаем, когда $q$ принимает значения простых чисел.
Несложно показать (теорема 1), что в первом приближении величины $\mathscr{L}(f;q,a)$ не зависят от $a$: главный член соответствующей формулы совпадает с (1.1). Напротив, остаточный член такой формулы имеет порядок $\sqrt{q}$ и существенным образом зависит от $a$. При этом оказывается, что величина
хорошо приближается некоторым тригонометрическим многочленом $\mathcal{R}_{q}(\varphi)$ и ее распределение зависит от того, какие значения принимает символ Лежандра $\bigl(\frac{p}{q}\bigr)$ на первых простых числах $p$, не превосходящих некоторой границы. Однако формализовать это утверждение удается лишь в двух специальных случаях (теоремы 2, 3), когда простое число $q$ лежит в особых арифметических прогрессиях с большой разностью.
В общем же случае удается лишь получить нижние оценки на число тех вычетов $a$, для которых остаточный член (1.2) принимает значения, близкие к заданному (теорема 4). Некоторые свойства полинома $\mathcal{R}_{q}(\varphi)$, связанного с (1.2), исследуются в теореме 5.
Обозначения
Всюду в статье $q$ – нечетное простое число, $\bigl(\frac{\cdot}{q}\bigr)$ – символ Лежандра по модулю $q$; $p_{n}$, $n= 1,2,3,\dots$, – все подряд идущие простые числа, занумерованные по возрастанию; $\psi(z)=\Gamma'(z)/\Gamma(z)$ – дигамма-функция; $\theta, \theta_{1}, \theta_{2},\dots$ – комплексные числа, не превосходящие по модулю единицы, в разных соотношениях, вообще говоря, разные. Поскольку функция $f(k)=k^{2}$ зафиксирована, вместо $\mathscr{L}(f;q,a)$ условимся писать $\mathscr{L}(q,a)$. Для натурального $n\geqslant 2$ через $P^{+}(n)$ и $P^{-}(n)$ обозначаются его наибольший и, соответственно, наименьший простой делитель. Число $n$ называется $y$-гладким (для $y\geqslant 2$), если $n=1$ либо $P^{+}(n)\leqslant y$, и называется $y$-грубым, если $P^{-}(n)>y$. Наконец, символом $\Omega(n)$ обозначается число простых делителей $n$ с учетом кратности.
Пусть $\alpha$ – точка максимума $Q(\varphi)$ на отрезке $0\leqslant\varphi\leqslant 0.5$, и $\beta=Q(\alpha)=(1/2)\operatorname{tg}{\pi\alpha}$; приближенное вычисление дает: $\alpha=0.1863009\dots$, $\beta=0.331372\dots$ . Несложно проверить, что производная $Q'(\varphi)$ имеет на отрезке $0\leqslant\varphi\leqslant 1$ нули в точках $\alpha$ и $1-\alpha$, монотонно убывает на интервале $(0,0.5)$ от $+\infty$ до значения $(-\pi/2)= Q'(0.5)$, и монотонно возрастает на интервале $(0.5,1)$ от $(-\pi/2)$ до $+\infty$. Далее, вторая производная $Q''(\varphi)$ монотонно возрастает на интервале $(0,1)$ от $-\infty$ до $+\infty$, обращаясь в нуль в точке $\varphi=0.5$. Наконец, третья производная $Q^{(3)}(\varphi)$ монотонно убывает на $(0,0.5)$ от $+\infty$ до значения $\pi^{3}=Q^{(3)}(0.5)$ и монотонно возрастает на $(0.5,1)$ до $+\infty$.
Для всякого $x$ c условиями $0\leqslant x\leqslant \beta$ определим величины $\xi(x)$ и $\eta(x)$ как наименьший и наибольший корни уравнения $Q(\varphi)=x$, лежащие на отрезке $0\leqslant \varphi\leqslant 0.5$. Тогда $\xi(x)$ монотонно возрастает на отрезке $0\leqslant x\leqslant \beta$ от значения $\xi(0)=0$ до значения $\xi(\beta)=\alpha$, а $\eta(x)$ монотонно убывает от $\eta(0)=0.5$ до $\eta(\beta)=\alpha$. Очевидно, если $0\leqslant x\leqslant\beta$, то уравнение $Q(\varphi)=-x$ имеет на отрезке $0.5\leqslant \varphi\leqslant 1$ два корня: $\xi_{1}$ и $\eta_{1}$, причем $\xi_{1}=1-\xi(x)$, $\eta_{1}=1-\eta(x)$.
Обозначим для каждой из функций $F=P,Q$ и любого вещественного $x$ через $\mathscr{E}(F;x)$ меру множества тех $\varphi$, $0\leqslant\varphi\leqslant 1$, для которых $F(\varphi)\leqslant x$. Определим, далее, для вещественного $x$ величину
В случае $F=Q$ при $|x|>\beta+|h|$ имеем, очевидно: $\omega(F;x,h)\equiv 0$. Следующая лемма дает оценку величины $\omega(F;x,h)$ при малых $h$.
Лемма 3. Существует абсолютная постоянная $c_{0}$ такая, что при любом $h$, $|h|\leqslant c_{0}$, и любом $x$ для каждой из функций $F=P,Q$ справедлива оценка
Доказательство. В случае $F=P$ несложно заключить, что наибольшее значение величины $\omega(F;x,h)$ равно длине отрезка $AD$ горизонтальной касательной к графику $P(\varphi)$ в точке $\varphi=1/2$, заключенного между двумя синусоидами (см. рис. 5). Обозначая длину $AB$ через $\delta$ (так что $AD=2\delta$), будем иметь
Докажем теперь утверждение в случае $F=Q$. Определим (см. рис. 6) для разных положений горизонтальной прямой $x=x_{i}$ отрезки $A_{i}B_{i}$, $C_{i}D_{i}$, предполагая, что в случаях, когда $x_{i}$ лежит между $\beta$ и $\beta+h$ либо между $(-\beta)$ и $(-\beta+h)$, точки $B_{i}$, $C_{i}$ совпадают и имеют одинаковую абсциссу (равную $\alpha$ либо $1-\alpha$). Тогда величина $\omega=\omega(F;x,h)$ равна суммарной длине указанных отрезков.
Будем считать, что $16h\leqslant \min\bigl\{Q(\alpha/2),Q(3\alpha/2)\bigr\}=0.276\dots$ . В силу симметрии рис. 6 относительно точки $(0.5,0.5h)$ достаточно доказать утверждение в случае $h>0$, $0.5h<x\leqslant \beta+h$.
Пусть сначала $0.5h<x\leqslant 16h$. Положим тогда $\xi_{0}=\xi(16h)$, $\xi_{1}=\xi(h)$, так что $0<\xi_{1}<\xi_{0}$. Введя обозначения как на рис. 7, (а), будем, очевидно, иметь $BE\leqslant B'E'= \xi_{0}$. Далее, если $A'$ – точка с координатами $(-\xi_{1},0)$, то $AE\leqslant A'O= \xi_{1}\leqslant A''E''=\xi_{0}$. Следовательно,
Таким образом, $\omega\leqslant 2(1.61\sqrt{h}+2.59\sqrt{h})=8.4\sqrt{h}$.
Пусть теперь $16h<x\leqslant \beta$, и пусть $u,v$ – абсциссы точек $C$ и $D$ (см. рис. 7, (b)). Очевидно, $\alpha\leqslant u\leqslant \eta_{0}$, $\alpha<v<0.5$. Поскольку $Q(u)=x$, $Q(v)+h=x=Q(u)$, то $h=Q(u)-Q(v)$. Полагая $\Delta=v-u$, будем иметь
Допустим теперь, что $\alpha \leqslant u < 3\alpha/2$. Покажем, что $v<2\alpha$. Действительно, положив $\eta_{1}=\eta(15h)$, будем иметь: $Q(\eta_{1})+h=16h<x=Q(v)+h$, т.е. $Q(\eta_{1})<Q(v)$ и, следовательно, $v<\eta_{1}$.
Таким образом, $v=u+\Delta<3\alpha/2+\Delta < 3\alpha/2+\alpha/2=2\alpha$, что и требовалось. Поэтому $\zeta_{3}\leqslant v < 2\alpha$, $|Q''(\zeta_{3})|> |Q''(2\alpha)|$, так что окончательно находим
Итак, $\Delta=CD\leqslant 0.72\sqrt{h}$ при $16h<x\leqslant \beta$. Но эта же оценка остается справедливой и в случае $\beta < x\leqslant \beta+h$, так как длина отрезка $C'D'$, отвечающего указанному $x$, не превосходит длины отрезка $CD$, отвечающего $x=\beta$.
Оценим, наконец, длину $\delta$ отрезка $AB$. Как и выше, достаточно ограничиться случаем $16h<x\leqslant\beta$. Пусть $Q(u)+h=Q(v)$, $v= u+\delta$. Тогда
Доказательство. Для упрощения выкладок вместо функции $f(k)=k^{2}$ будем рассматривать функцию $f(k)=-k^{2}$. Пусть $A, B$ – точки на единичной окружности, и пусть $\pi\varphi$ – длина любой из дуг, соединяющих эти точки. Тогда $|AB|=2|\sin{\pi\varphi}|$. Полагая $\varphi=(ak^{2}+k)/q$, и пользуясь разложением
Выражение в фигурных скобках равно $2\cos(2\pi\,\overline{4}\,\overline{a}n/q)$ либо $\bigl(-2i\sin(2\pi\,\overline{4}\,\overline{a}n/q)\bigl)$ в зависимости от того, сравнимо ли $q$ с $1$ или $3$ по модулю $4$. Следовательно, второе слагаемое в (3.1) представится в виде
Если $a$ пробегает полную систему вычетов по модулю $q$, то множитель при $(-8/\pi)\sqrt{q}$ в правой части последнего равенства пробегает то же множество значений, что и величина
Несложно проверить, что произведение $\bigl(\frac{\ell}{q}\bigr)\mathcal{H}\bigl(\ell a/q\bigr)$ не меняется от замены $\ell$ на $q-\ell$. Следовательно,
4. Распределение величин $r(q;a)$ в случае специальных модулей $q$
Результаты вычисления значений величин $r(q;a)$, $a=1,2,\dots,q$, для заданных $q$ можно представить в виде гистограмм, примеры которых приведены на рис. 8. В свете теоремы 1 естественно предположить, что распределение величин $r(q;a)$ существенно зависит от свойств полинома $\mathcal{R}_{q}(\varphi)$. В свою очередь, поведение этого полинома во многом определяется тем, какие значения принимает символ Лежандра $\bigl(\frac{n}{q}\bigr)$ на начальном отрезке натурального ряда. По этой причине отыскать предельную форму гистограмм удается лишь в очень частном случае, когда наименьший квадратичный невычет $n_{q}$ по модулю $q$ достаточно велик, т.е. когда полином $\mathcal{R}_{q}(\varphi)$ с хорошей точностью приближается функцией
Последняя же совпадает с функцией $P(\varphi)$ в случае $q\equiv 1\pmod{4}$ и с функцией $Q(\varphi)$ при $q\equiv 3\pmod{4}$. Справедливы следующие утверждения.
Тогда при любом $x$, $-0.5\leqslant x\leqslant 0.5$, для количества $\mathcal{N}_{1}(q;x)$ вычетов $a$, $0<a\leqslant q$, удовлетворяющих неравенству $r(q;a)\leqslant x$, справедлива формула
Тогда при любом $x$, $-\beta\leqslant x\leqslant \beta$, для количества $\mathcal{N}_{2}(q;x)$ тех вычетов $a$, $0<a\leqslant q$, что удовлетворяют условию $r(q;a)\leqslant x$, справедлива формула
Доказательство. Докажем теорему 2. Из условий (4.2) следует, что $p_{N}<n_{q}$, где $n_{q}$ – наименьший квадратичный невычет по модулю $q$. Очевидно, $n_{q}<q$, так что $(p_{n},q)=1$ для всех $n\leqslant N$. Поэтому $\bigl(\frac{\ell}{q}\bigr)=+1$ для всякого $p_{N}$-гладкого числа $\ell$. Обозначим через $\mathcal{A}$ множество всех $p_{N}$-гладких чисел отрезка $1\leqslant \ell\leqslant t$, $t=(q-1)/2$, а через $\mathcal{B}$ – дополнение к нему. Тогда
Далее, всякое $\ell\in \mathcal{B}$ имеет вид $n=km$, где $m\in \mathcal{A}$, а $k$ является $p_{N}$-грубым. Следовательно,
$$
\begin{equation*}
\begin{aligned} \, |S_{2}(\varphi)| &\leqslant 2\sum_{\substack{\ell\leqslant t \\ \ell=km\in \mathcal{B}}}\frac{1}{4m^{2}k^2-1} \leqslant 2\sum_{\substack{\ell\leqslant t \\ \ell=km\in \mathcal{B}}}\frac{1}{(4m^{2}-1)k^2} \\ &< 2\sum_{m=1}^{+\infty}\frac{1}{4m^{2}-1}\sum_{\substack{p_{N}<k\leqslant t \\ P^{-}(k)>p_{N}}}\frac{1}{k^{2}} < \sum_{\substack{k>p_{N}\\ P^{-}(k)>p_{N}}}\frac{1}{k^{2}}. \end{aligned}
\end{equation*}
\notag
$$
Обозначим через $\Phi(x,y)$ количество $y$-грубых чисел, не превосходящих $x$. Применяя к сумме по $k$ преобразование Абеля наряду с оценкой $\Phi(x,y)\leqslant x/\ln{y}$, справедливой при всех $2\leqslant y\leqslant x$ (см. [10]), будем иметь:
Зафиксируем теперь произвольное $x$, $-0.5\leqslant x\leqslant 0.5$. Тогда количество $\nu(q;x)$ чисел $a$ ряда $1,2,\dots, q-1$ с условием $s(q;a)\leqslant x$ складывается из количества $\nu_{1}(q;a)$ чисел того же ряда с условиями
имеет на отрезке $0\leqslant \varphi\leqslant 1$ два решения: $\varphi(x)$ и $1-\varphi(x)$, где $\varphi(x)$ определена выше. Условия (4.6) в этом случае принимают вид
Иными словами, $\nu_{1}(q;x)$ есть число квадратичных вычетов на отрезке $q\varphi(x)\leqslant a\leqslant q(1-\varphi(x))$. В силу классической оценки суммы символов Лежандра, принадлежащей Виноградову (см. [13], а также [12; гл. V, вопрос 11b ($\varkappa$)]), имеют место равенства
Аналогично, $\nu_{2}(q;x)$ совпадает с числом квадратичных невычетов, которые лежат в объединении промежутков $0<a\leqslant q\varphi(-x)$, $q(1-\varphi(-x))\leqslant a <q$ или, что то же, в промежутке $q(1-\varphi(-x))\leqslant a\leqslant q(1+\varphi(-x))$. Так находим
Эта же формула остается справедливой и в случае $-(\pi/4-1/2)\leqslant x<0$. Случай $|x|>\pi/4-1/2$ рассматривается аналогично. В итоге получаем равенство
По доказанному выше, график полинома $\mathcal{R}_{q}(\varphi)$ заключен в полосе ширины $2h$, где $h$ определено в (4.5), окружающей график $P(\varphi)$ (ее границы получаются сдвигом графика $P(\varphi)$ по вертикали на $\pm h$).
Согласно лемме 3 суммарная длина промежутков, на которых выполняется неравенство $\mathcal{R}_{q}(\varphi)\leqslant x$, отличается от величины $\mathscr{E}(P;x)$ самое большее на $8.4\sqrt{h}$. Количество дробей вида $a/q$, содержащихся в таких промежутках, не превосходит числа дробей $a/q$, попавших в проекции промежутков $AB$, $CD$ (см. рис. 6) на горизонтальную ось. Последнее же число не превосходит, очевидно
Из условий, которым удовлетворяет $p_{N}$, заключаем: $|\Delta(q;x)|<10\varepsilon$. Теорема 2 доказана.
Доказательство теоремы 3 проводится с помощью тех же рассуждений, с заменой (4.5) аналогичным неравенством для разности $Q(\varphi)-\mathcal{R}_{q}(\varphi)$.
На рис. 9, 10 представлены графики функций распределения $F_{i}(x)$ и отвечающих им плотностей $f_{i}(x)=F_{i}'(x)$, $i=1,2$.
В свою очередь, графики плотности интересно сопоставить с гистограммами, отвечающими значениям $q=9601$ и $q=18191$ (см. рис. 11). В первом случае $q\equiv 1\pmod{4}$, $n_{q}=13$, однако при этом $\bigl(\frac{17}{q}\bigr)=\bigl(\frac{19}{q}\bigr)=+1$. Во втором случае $q\equiv 3 \pmod{4}$ и при этом $n_{q}=23$.
Тогда $\alpha-h=\xi(\beta-\delta)$, так что для разложения в ряд по степеням $\delta$ функции $\xi(\beta-\delta)$ достаточно найти разложение для $h$. В свою очередь, согласно (4.13) имеем: $\delta=\beta-Q(\alpha-h)=Q(\alpha)-Q(\alpha-h)$, откуда
В частности, беря $\varepsilon$ и $\delta$ достаточно малыми, причем так, чтобы $\varepsilon=o(\sqrt{\delta})$, для количества чисел $a$ с условиями $0<a<q$, $\beta-\delta<r(q;a)\leqslant \beta$ получим выражение вида
можно получить нижнюю оценку на количество тех квадратичных вычетов (невычетов) $a$, для которых величина $r(q;a)$ близка к $f_{q}(\varphi)$ (соответственно, к $-f_{q}(\varphi)$), где $\varphi$ – произвольное число отрезка $0\leqslant \varphi\leqslant 1$.
Теорема 4. Пусть $q$ – достаточно большое простое число, $2q^{-1/4}\sqrt{\ln{q}}<\varepsilon<10^{-2}$, $N$ – произвольное целое с условиями $3\varepsilon^{-1}<N\leqslant \varepsilon\sqrt{q}/(\ln{q})$. Пусть, далее, $\varphi$ – произвольное фиксированное число из промежутка $-0.5\leqslant\varphi\leqslant 0.5$. Тогда по крайней мере для
Доказательство. Ограничимся случаем $q\equiv 1\,(\operatorname{mod}{4})$, когда $\mathcal{H}(x)=\cos{(2\pi x)}$; для $q\equiv 3\pmod{4}$ рассуждения аналогичны.
Пусть $\varphi$ фиксировано, и пусть $a$ принимает значения $[q\varphi]+\ell$, где $|\ell|\leqslant \varepsilon q/N$. Полагая $\vartheta= \{q\varphi\}$, будем иметь
Замечание 2. Пусть $\varepsilon>0$ фиксировано. Беря $N=\bigl[\varepsilon\sqrt{q}/\ln{q}\bigr]$, получим не менее $(0.5-o(1))\sqrt{q}\ln{q}$ квадратичных вычетов, для которых
Далее, если $\psi(t)$ – произвольная положительная монотонно возрастающая и неограниченная сверху при $t\to+\infty$ функция с условием $\psi(t)<t$, то, полагая $N=[\varepsilon^{-1}\psi(\varepsilon^{-1})]$, в случае $q\geqslant q_{0}(\varepsilon)$ получим не менее
Аналогичные утверждения справедливы и для квадратичных невычетов (с заменой $f_{q}(\varphi)$ на $-f_{q}(\varphi)$).
6. Некоторые свойства полиномов $\mathcal{R}_{q}(\varphi)$
Как отмечалось выше, основным препятствием, которое не позволяет получить аналоги теорем 2,3 в случае произвольного простого $q$, является сложное поведение тригонометрических полиномов $\mathcal{R}_{q}(\varphi)$ и, в частности, отсутствие в общем случае аналога леммы 3, т.е. оценки типа $\omega(\mathcal{R}_{q};h,x)\leqslant g(h)$, где $g(h)\to 0$ при $h\to 0$. В этом разделе помещена теорема, в некоторой степени указывающая на различие в поведении $\mathcal{R}_{q}(\varphi)$ при различных $q$. Именно, оказывается, что для положительной доли простых чисел $q$ значение $\mathcal{R}_{q}(0.5)$ может быть сколь угодно близко к минимуму $\mathcal{R}_{q}(\varphi)$ на отрезке $0 \leqslant \varphi\leqslant 1$, и в то же время для положительной доли простых $q$ полином $\mathcal{R}_{q}(\varphi)$ в точке $\varphi_{0}\approx 0.4$ принимает значение меньшее, чем в точке $\varphi=0.5$ (см. рис. 12). Более точно, справедливо следующее утверждение.
Теорема 5. Пусть $0<\varepsilon<10^{-2}$ – сколь угодно малое фиксированное число. Тогда существует постоянная $c=c(\varepsilon)$ такая, что при $x\to+\infty$ для не менее чем $c\pi(x)$ простых чисел $q$ с условием $q\equiv 1\pmod{4}$ неравенство $\mathcal{R}_{q}(\varphi)>\mathcal{R}_{q}(0.5)-\varepsilon$ выполнено на всем промежутке $0\leqslant\varphi\leqslant 1$. Вместе с тем, имеется не менее $(2^{-8}+o(1))\pi(x)$ простых чисел $q$ с условием $q\equiv 1\pmod{4}$, для которых можно указать число $\varphi_{0}$ с условиями
Обозначая, как и выше, через $\mathcal{A}$ множество всех $p_{N}$-гладких чисел, не превосходящих $(q-1)/2$, а через $\mathcal{B}$ – дополнение к нему, будем иметь
Если $n=2^{2s+1}m$, где $P^{+}(m)\leqslant p_{N}$, то слагаемое, отвечающее значению $2n=2^{2(s+1)}m$, встретится во второй сумме; верно и обратное. Следовательно,
Остается заметить, что всякое проcтое $q$, удовлетворяющее (6.1), принадлежит прогрессии $q\equiv 5\pmod{8}$ и, сверх того, одной из $(p_{k}-1)/2=\varphi(p_{k})/2$ прогрессий вида $q\equiv a\pmod{p_{k}}$, где $a$ – произвольный квадратичный вычет по модулю $p_{k}$. Ясно, что при фиксированном $p_{N}$ и $x\to+\infty$ общее число таких $q$ равно
Для завершения доказательства остается заметить, что условию (6.2) удовлетворяют $(2^{-8}+o(1))\pi(x)$ простых чисел $q$, $q\leqslant x$. Теорема 5 доказана.
Автор благодарит М. Р. Габдуллина за помощь в подготовке статьи.
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1.
A. J. Picard, “Graphing in modular systems”, Math. Teacher, 64:5 (1971), 459–466
2.
P. Locke, “Residue designs”, Math. Teacher, 65:3 (1972), 260–263
3.
T. E. Moore, “Aspects of group theory in residue designs”, Two-year College Mathematics Readings, eds. W. Page, Math. Association Amer., Washington, DC, 1981, 254–257
4.
I. D. Johnson, “Paving the way to algebraic thought using residue designs”, Math. Teacher, 91:4 (1998), 326–332
5.
T. Koshy, Elementary Number Theory with Applications, Elsevier, 2007
6.
D. Richeson, “Residue designs, string art, and number theory”, Proceedings of Bridges 2023: Mathematics, Art, Music, Architecture, Culture, eds. J. Holdener, E. Torrence, C. Fong, K. Seaton, Tessellations Publ., Phoenix, 2023, 365–368
7.
К. П. Кохась, “Волшебная арифметика хорд”, Квант, 2024, № 3, 13–23