OFFSET
1,2
COMMENTS
Numbers k such that k | A002895(k).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..132
EXAMPLE
MATHEMATICA
seq[kmax_] := Module[{d0 = 1, d1 = 4, d2, s = {1}}, Do[d2 = ((20*k^3 - 30*k^2 + 18*k - 4)*d1 - 64*(k-1)^3*d0)/k^3; If[Divisible[d2, k], AppendTo[s, k]]; d0 = d1; d1 = d2, {k, 2, kmax}]; s]; seq[5000]
PROG
(PARI) lista(kmax) = {my(d0 = 1, d1 = 4, d2); print1("1, "); for(k = 2, kmax, d2 = ((20*k^3 - 30*k^2 + 18*k - 4)*d1 - 64*(k-1)^3*d0)/k^3; if(!(d2 % k), print1(k, ", ")); d0 = d1; d1 = d2); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 17 2024
STATUS
approved