OFFSET
1,4
LINKS
FORMULA
a(1) = 1; a(n) = Sum_{d|n, d<n} A136522(n/d) * a(d).
MAPLE
ispali:= proc(n) local L;
L:= convert(n, base, 10);
evalb(L = ListTools:-Reverse(L))
end proc:
N:= 100: # for a(1)..a(N)
Palis:= select(ispali, {$2..N}):
A:= Vector(N):
A[1]:= 1:
for n from 2 to N do
A[n]:= add(A[n/d], d= numtheory:-divisors(n) intersect Palis)
od:
convert(A, list); # Robert Israel, Nov 15 2020
MATHEMATICA
a[n_] := If[n == 1, n, Sum[If[(d < n && PalindromeQ[n/d]), a[d], 0], {d, Divisors[n]}]]; Table[a[n], {n, 1, 90}]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Ilya Gutkovskiy, Nov 02 2020
STATUS
approved