Сумма обратных квадратов гипотенузы пифагоровых треугольников

Воспользуемся обозначениями из

http://mathworld.wolfram.com/DoubleSeries.html

Зафиксируем положительно определенную бинарную квадратичную форму$q$данный$q(m,n)=am^2+bmn+cn^2$,$a,b,c$целые числа. Мы используем суммирование по набору индексов

$$J=\Bbb Z\times\Bbb Z-\{(0,0)\}\ .$$ Мы определяем $$\begin{aligned} S(q;s) = S(a,b,c;s) &=\sum_{(m,n)\in J} q(m,n)^{-s}=\sum_{(m,n)\in J} (am^2+bmn+cn^2)^{-s}\ ,\\ S_1(q;s) = S_1(a,b,c;s) &=\sum_{(m,n)\in J} \color{blue}{(-1)^m}\; q(m,n)^{-s}\ ,\\ S_2(q;s) = S_2(a,b,c;s) &=\sum_{(m,n)\in J} \color{blue}{(-1)^n}\; q(m,n)^{-s}\ ,\\ S_{12}(q;s) = S_{12}(a,b,c;s) &=\sum_{(m,n)\in J}\color{blue}{(-1)^{m+n}}\; q(m,n)^{-s}\ . \end{aligned}$$ Последние три суммы являются «подкрученными версиями» первой суммы, «подкрутка» происходит за счет использования символа для первого параметра, для второго, для обоих. В нашем случае$q(m,n)=m^2 +n^2$, и$(a,b,c)=(1,0,1)$, имеем симметричный случай (относительно замены$a\leftrightarrow c$).

мы бросим$q$ниже из обозначений в$S_?(q,s)$, так как мы используем только приведенную выше квадратичную форму$q$. Я решил во время операции редактирования, которая должна быстро привести нас к числам, которые мы можем вычислить, что для проверки лучше ввести версии.$S^+$для всех сумм, где плюсовой индекс указывает на дальнейшее ограничение на$(m,n)\in J$с

$$(+)\qquad m,n>0\ .$$

Из лок. цит. извлекаем следующие соотношения:

$$\begin{aligned} S(s) &= \sum_{(m,n)\in J}(m^2+n^2)^{-s} \\ &= 4\beta(s)\;\zeta(s)\ ,\\ %S_1(s) =S_2(s) &= \sum_{(m,n)\in J}(-1)^m(m^2+n^2)^{-s} =\dots %= 2^{-s}b_2(2s) = -2^{-s}\cdot 4\beta(2s)\; \eta(2s) %\ , %\\ S_{12}(s) &= \sum_{(m,n)\in J}(-1)^{m+n}(m^2+n^2)^{-s} \\ &= -4 \beta(s) \;\eta(s)=-4\beta(s) \;(1-2^{1-s})\; \zeta(s)\ . \\[2mm] &\qquad\text{ Then the plus versions are:} \\[3mm] S^+(s) &= \beta(s)\;\zeta(s) - \zeta(2s)\ , \\ -S_{12}^+(s) &= \beta(s)\;\eta(s) - \eta(2s) \\ &= \beta(s)\;(1-2^{1-s})\zeta(s) - (1-2^{1-2s})\zeta(2s)\ , \\ &\qquad \text{ which gives} \\ S^+(s)-S_{12}^+(s) &=2\beta(s)\;(1-2^{-s})\zeta(s) - 2(1-2^{-2s})\zeta(2s)\ . \end{aligned}$$

Теперь найдем линейную комбинацию вышеуказанных сумм, соответствующую суммированию$q(m,n)^{-s}$по множеству$K$из всех$(m,n)$с положительной (составляющей) разной четности. Это

$$\frac 12(\ S^+(s)-S^+_{12}(s)\ )\ .$$ Пока мы можем написать: $$\begin{aligned} &\beta(s)\;(1-2^{-s})\zeta(s) - (1-2^{-2s})\zeta(2s) \\ &\qquad=\frac 12(\ S^+(s)-S^+_{12}(s)\ ) \\ &\qquad=\sum_{\substack{(m,n)\in K\\m,n> 0}} q(m,n)^{-s}\\ &\qquad=2\sum_{\substack{(m,n)\in K\\m>n> 0}} q(m,n)^{-s}\\ &\qquad=2\sum_{\substack{(m,n)\in K\\m>n> 0\\ d=(m,n)\text{ odd}}} q(m,n)^{-s}\qquad\text{ and with }M=m/d,\ N=n/d\\ &\qquad=2\sum_{d>0\text{ odd}}d^{-2s} \sum_{\substack{(M,N)\in K\\M>N> 0\\ (M,N)=1}} q(M,N)^{-s}\\ &\qquad= 2(1-2^{-2s})\; \zeta(2s)\; \color{blue}{ \sum_{\substack{(M,N)\in K\\M>N> 0\\ (M,N)=1}} q(M,N)^{-s} } \ . \end{aligned}$$ Выделенная сумма в последнем выражении и есть нужная нам сумма, возьмем ее за$s=2$.

Значение, которое мы получаем:

$$\color{brown}{ \frac{\beta(2)\;\zeta(2)}{2(1+2^{-2})\zeta(4)} -\frac 12\ = \frac{6C}{\pi^2} - \frac 12.}$$

$$\color{brown}{ \frac{\beta(2)\zeta(2)}{2(1+2^{-2})\zeta(4)} - \frac{1}{2} = \frac{6C}{\pi^2} - \frac 1 2.}$$

где$C$каталонская постоянная. Численно:

sage: E = catalan * zeta(2) / 2 / (1+2^-2) / zeta(4) - 1/2 
sage: E.n()
0.0568403090661582