This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A363972 #13 Jun 30 2023 15:36:39 %S A363972 1,1,1,1,5,1,1,1,10,5,1,1,17,1,5,1,26,10,1,5,37,1,1,1,54,17,10,1,65,5, %T A363972 1,1,82,26,5,10,101,1,17,5,122,37,1,1,158,1,1,1,170,54,26,17,197,10,5, %U A363972 1,226,65,1,5,257,1,46,1,310,82,1,26,325,5,1,10,362,101,54,1,401,17,1,5,451,122,1 %N A363972 Expansion of Sum_{k>0} k^2 * x^(4*k-3) / (1 - x^(4*k-3)). %F A363972 a(n) = Sum_{d|n, d==1 mod 4} ((d+3)/4)^2. %t A363972 a[n_] := DivisorSum[n, ((#+3)/4)^2 &, Mod[#, 4] == 1 &]; Array[a, 100] (* _Amiram Eldar_, Jun 30 2023 *) %o A363972 (PARI) a(n) = sumdiv(n, d, (d%4==1)*((d+3)/4)^2); %Y A363972 Cf. A001826, A363903. %Y A363972 Cf. A050450, A363977. %K A363972 nonn,easy %O A363972 1,5 %A A363972 _Seiichi Manyama_, Jun 30 2023