login
a(1) = 1; for n >= 2, a(n) = noz(a(n-1) + 3*n - 2), where noz(n) = A004719(n).
4

%I #14 Jun 07 2024 16:36:35

%S 1,5,12,22,35,51,7,29,54,82,113,147,184,224,267,313,362,414,469,527,

%T 588,652,719,789,862,938,117,199,284,372,463,557,654,754,857,963,172,

%U 284,399,517,638,762,889,119,252,388,527,669,814,962,1113,1267,1424,1584

%N a(1) = 1; for n >= 2, a(n) = noz(a(n-1) + 3*n - 2), where noz(n) = A004719(n).

%C Zeroless analog of the positive pentagonal numbers.

%H Paolo Xausa, <a href="/A373171/b373171.txt">Table of n, a(n) for n = 1..10000</a>

%t noz[n_] := FromDigits[DeleteCases[IntegerDigits[n], 0]];

%t Block[{n = 1}, NestList[noz[++n*3 - 2 + #] &, 1, 100]]

%o (PARI) noz(n) = fromdigits(select(sign, digits(n))); \\ A004719

%o lista(nn) = my(va=vector(nn)); for (n=1, nn, va[n] = if (n==1, 1, noz(va[n-1] + 3*n - 2))); va; \\ _Michel Marcus_, Jun 03 2024

%Y Row n = 5 of A373169.

%Y Cf. A000326, A004719, A243658, A370812, A373172.

%K nonn,base,easy

%O 1,2

%A _Paolo Xausa_, May 28 2024