О вычислении классических сумм Якобсталя

В статье рассказывается о быстрых алгоритмах вычисления классических сумм Якобсталя. Эти алгоритмы основаны на нетривиальной связи сумм Якобсталя с представлением простых чисел $p\equiv 1 \pmod 4$ и $p\equiv 1 \pmod 3$ бинарными квадратичными формами $A^2+B^2$ и $A^2+3B^2$ соответственно.