This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A338881 #6 Nov 14 2020 11:03:39 %S A338881 1,1,2,3,4,5,6,7,1,1,2,3,4,5,6,7,2,2,4,6,8,10,12,14,3,3,6,9,12,15,18, %T A338881 21,4,4,8,12,16,20,24,28,5,5,10,15,20,25,30,35,6,6,12,18,24,30,36,42, %U A338881 7,7,14,21,28,35,42,49,1,1,2,3,4,5,6,7,1,1,2,3,4,5,6,7,2 %N A338881 Product of the nonzero digits of (n written in base 8). %F A338881 G.f. A(x) satisfies: A(x) = (1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 5*x^5 + 6*x^6 + 7*x^7) * A(x^8). %t A338881 Table[Times @@ DeleteCases[IntegerDigits[n, 8], 0], {n, 0, 80}] %t A338881 nmax = 80; A[_] = 1; Do[A[x_] = (1 + x + 2 x^2 + 3 x^3 + 4 x^4 + 5 x^5 + 6 x^6 + 7 x^7) A[x^8] + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x] %o A338881 (PARI) a(n) = vecprod(select(x->x, digits(n, 8))); \\ _Michel Marcus_, Nov 14 2020 %Y A338881 Cf. A007094, A051801, A117592, A309959. %K A338881 nonn,base %O A338881 0,3 %A A338881 _Ilya Gutkovskiy_, Nov 13 2020