login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A035227
Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 45.
1
1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 2, 1, 0, 0, 1, 1, 0, 0, 2, 1, 0, 0, 0, 0, 1, 0, 1, 0, 2, 0, 2, 0, 2, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 2, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 2, 0, 2, 0, 2, 1, 2, 0, 0, 1, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1, 2, 0, 0, 2, 1, 1
OFFSET
1,11
LINKS
FORMULA
From Amiram Eldar, Nov 20 2023: (Start)
a(n) = Sum_{d|n} Kronecker(45, d).
Multiplicative with a(p^e) = 1 if Kronecker(45, p) = 0 (p = 3 or 5), a(p^e) = (1+(-1)^e)/2 if Kronecker(45, p) = -1 (p is in A003631 \ {3}), and a(p^e) = e+1 if Kronecker(45, p) = 1 (p is in A045468).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 8*log(phi)/(3*sqrt(5)) = 0.573878587952..., where phi is the golden ratio (A001622) . (End)
MATHEMATICA
a[n_] := DivisorSum[n, KroneckerSymbol[45, #] &]; Array[a, 100] (* Amiram Eldar, Nov 20 2023 *)
PROG
(PARI) my(m = 45); direuler(p=2, 101, 1/(1-(kronecker(m, p)*(X-X^2))-X))
(PARI) a(n) = sumdiv(n, d, kronecker(45, d)); \\ Amiram Eldar, Nov 20 2023
CROSSREFS
KEYWORD
nonn,easy,mult
STATUS
approved