Could be solved easily with https://en.wikipedia.org/wiki/Fermat%27s_theorem_on_sums_of_two_squares and where = number of decompositions of in two squares. The French version of the above page gives this formula and the meaning of and .