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 A359938 #30 Aug 18 2024 20:48:39
%S A359938 1,2,1,2,2,2,1,2,1,4,1,2,1,2,2,2,2,2,1,4,1,2,1,2,2,3,1,2,1,4,1,2,1,3,
%T A359938 2,2,2,2,1,4,1,2,1,2,2,2,1,2,1,5,2,3,1,2,2,2,1,2,1,4,1,2,1,2,3,2,1,3,
%U A359938 1,4,1,2,1,3,2,2,1,3,1,4,1,3,1,2,3,2,1,2,1,4
%N A359938 Number of divisors d of n such that d-1 is square.
%C A359938 The Cartesian equation for the Witch of Agnesi is given as y = 8*k^3/(x^2+4*k^2). If we set x = n, then a(n)-1 is the number of positive integer solutions in k such that y is a positive integer. Let d = m^2+1 be a divisor of n, then k = m*n/2 is a solution. - _Thomas Scheuerle_, Aug 07 2024
%F A359938 G.f.: Sum_{k>=0} x^(k^2+1)/(1 - x^(k^2+1)).
%F A359938 From _Vaclav Kotesovec_, Jan 19 2023: (Start)
%F A359938 a(n) = Sum_{k=0..n} (floor(n/(k^2 + 1)) - floor((n-1)/(k^2 + 1))).
%F A359938 Sum_{k=1..n} a(k) = Sum_{k=0..n} floor(n/(k^2 + 1)).
%F A359938 Sum_{k=1..n} a(k) ~ c*n, where c = A113319 = (1 + Pi*coth(Pi))/2 = 2.0766... (End)
%t A359938 nmax = 100; Rest[CoefficientList[Series[Sum[x^(k^2 + 1)/(1 - x^(k^2 + 1)), {k, 0, Sqrt[nmax]}], {x, 0, nmax}], x]] (* _Vaclav Kotesovec_, Jan 19 2023 *)
%t A359938 Table[Sum[Floor[n/(k^2 + 1)] - Floor[(n-1)/(k^2 + 1)], {k, 0, Sqrt[n]}], {n, 1, 100}] (* _Vaclav Kotesovec_, Jan 19 2023 *)
%o A359938 (PARI) a(n) = sumdiv(n, d, issquare(d-1));
%o A359938 (PARI) my(N=100, x='x+O('x^N)); Vec(sum(k=0, sqrtint(N), x^(k^2+1)/(1-x^(k^2+1))))
%Y A359938 Cf. A002522, A046951, A113319, A359937, A359944.
%K A359938 nonn
%O A359938 1,2
%A A359938 _Seiichi Manyama_, Jan 19 2023