A338854 Product of the nonzero digits of (n written in base 4).
1, 1, 2, 3, 1, 1, 2, 3, 2, 2, 4, 6, 3, 3, 6, 9, 1, 1, 2, 3, 1, 1, 2, 3, 2, 2, 4, 6, 3, 3, 6, 9, 2, 2, 4, 6, 2, 2, 4, 6, 4, 4, 8, 12, 6, 6, 12, 18, 3, 3, 6, 9, 3, 3, 6, 9, 6, 6, 12, 18, 9, 9, 18, 27, 1, 1, 2, 3, 1, 1, 2, 3, 2, 2, 4, 6, 3, 3, 6, 9, 1
Offset: 0
Programs
-
Mathematica
Table[Times @@ DeleteCases[IntegerDigits[n, 4], 0], {n, 0, 80}] nmax = 80; A[] = 1; Do[A[x] = (1 + x + 2 x^2 + 3 x^3) A[x^4] + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
-
PARI
a(n) = vecprod(select(x->x, digits(n, 4))); \\ Michel Marcus, Nov 12 2020