Let $a,b$ be two positive integers. Let the sequence $\{s_n\}_n$ be the set of all possible sums of squares $a^2+b^2$, such that they are in ascending order \begin{align*} & n=1 & s_1=1^2+1^2=2 \\ & n=2 & s_2=1^2+2^2=5 \\ & n=3 & s_3=2^2+1^2=5 \\ & n=4 & s_4=2^2+2^2=8 \end{align*} etc. Here $n$ serves as a counting index for the sequence and does not serve any other purpose. I have found that the number of identical values that each sum can have is equal to the number of lattice points lying on a quarter of circumference of a circle of radius $n$. For example, $s_2=s_3=5$ we have two values, but there are other cases, such as \begin{align*} & n=31 & s_{31}=1^2+7^2=50 \\ & n=32 & s_{32}=5^2+5^2=50 \\ & n=33 & s_{33}=7^2+1^2=50 \end{align*}
To the best of my knowledge, there are no references to the asymptotic form of $\{s_n\}_n$ as $n\to \infty$, if indeed it exists. A brute-force approach with $n\le10^6$ led to the expression $$ \{s_n\}_n\sim \beta n,\qquad\beta\approx1.276\ldots. $$ Therefore my question is: is there any known closed form of the asymptotics of this particular succession? A more profound and interesting question is, why does it seem to be linear?