OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: Product_{i>=1} (1-x^(10*i))^10/(1-x^i).
Euler transform of period 10 sequence [1,1,1,1,1,1,1,1,1,-9, .. ].
MAPLE
with(numtheory): A:= proc(n, t) option remember; local d, j; `if`(n=0, 1, add(add(`if`(t=0 or irem(d, t)=0, d-d*t, d), d=divisors(j)) *A(n-j, t), j=1..n)/n) end: seq(A(n, 10), n=0..50);
MATHEMATICA
A[n_, t_] := A[n, t] = Module[{d, j}, If[n == 0, 1, Sum[Sum[If[t == 0 || Mod[d, t] == 0, d - d t, d], {d, Divisors[j]}] A[n - j, t], {j, 1, n}]/n]];
Table[A[n, 10], {n, 0, 50}] (* Jean-François Alcover, Dec 06 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Dec 03 2010
STATUS
approved