OFFSET
0,12
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = [x^n y^9] 1/Product_{j>=2} (1-y*x^A002113(j)).
MAPLE
p:= proc(n) option remember; local i, s; s:= ""||n;
for i to iquo(length(s), 2) do if
s[i]<>s[-i] then return false fi od; true
end:
h:= proc(n) option remember; `if`(n<1, 0,
`if`(p(n), n, h(n-1)))
end:
b:= proc(n, i, t) option remember; `if`(n=0, 1, `if`(t*i<n,
0, b(n, h(i-1), t)+b(n-i, h(min(n-i, i)), t-1)))
end:
a:= n-> (k-> b(n, h(n), k)-b(n, h(n), k-1))(9):
seq(a(n), n=0..100);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Alois P. Heinz, Sep 19 2018
STATUS
approved