A359227 Number of divisors of 4*n-3 of form 4*k+1.
1, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 4, 2, 2, 2, 2, 4, 2, 2, 2, 3, 4, 2, 2, 2, 2, 4, 2, 2, 4, 2, 4, 2, 2, 2, 2, 4, 2, 4, 2, 2, 4, 3, 2, 2, 2, 4, 4, 2, 2, 2, 4, 2, 2, 2, 4, 6, 2, 2, 2, 2, 4, 2, 2, 2, 4, 4, 2, 4, 2, 2, 4, 3, 2, 4, 2, 4, 2, 2, 2, 2, 6, 2, 4, 2, 2, 4, 2, 2, 4
Offset: 1
Programs
-
Mathematica
a[n_] := DivisorSum[4*n-3, 1 &, Mod[#, 4] == 1 &]; Array[a, 100] (* Amiram Eldar, Aug 23 2023 *)
-
PARI
a(n) = sumdiv(4*n-3, d, d%4==1);
-
PARI
my(N=100, x='x+O('x^N)); Vec(sum(k=1, N, x^k/(1-x^(4*k-3))))
Formula
a(n) = A001826(4*n-3).
G.f.: Sum_{k>0} x^k/(1 - x^(4*k-3)).