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 A112058 #10 May 06 2021 12:25:35 %S A112058 48,72,120,168,192,288,312,360,408,432,480,528,552,648,672,720,768, %T A112058 792,840,888,912,960,1008,1032,1128,1152,1200,1248,1272,1320,1368, %U A112058 1392,1488,1512,1560,1608,1632,1680,1728,1752,1800,1848,1872,1968 %N A112058 Mean of A112056 and A112057. %F A112058 a(n) = 4*A112054(n). %t A112058 a112046[n_]:=Block[{i=1}, While[JacobiSymbol[i, 2n + 1]==1, i++]; i]; 4*Select[Range[1000], a112046[2#] - a112046[2# - 1] != 0 &] (* _Indranil Ghosh_, May 24 2017 *) %o A112058 (Python) %o A112058 from sympy import jacobi_symbol as J %o A112058 def a112046(n): %o A112058 i=1 %o A112058 while True: %o A112058 if J(i, 2*n + 1)!=1: return i %o A112058 else: i+=1 %o A112058 def a(n): return a112046(2*n) - a112046(2*n - 1) %o A112058 print([4*n for n in range(1, 1001) if a(n)!=0]) # _Indranil Ghosh_, May 24 2017 %Y A112058 a(n) = A112056(n)+1 = A112057(n)-1 = (A112056(n)+A112057(n))/2. %K A112058 nonn %O A112058 1,1 %A A112058 _Antti Karttunen_, Aug 27 2005