OFFSET
1,4
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
MAPLE
N:= 10^9: # for 3^i * 7^j <= N
R:= sort([seq(seq(3^i*7^j, j=0 .. ilog[7](N/3^i)), i=0..ilog[3](N))]):
map(padic:-ordp, R, 3); # Robert Israel, Jun 04 2025
MATHEMATICA
seq[lim_] := IntegerExponent[Flatten[Table[3^i*7^j, {i, 0, Log[3, lim]}, {j, 0, Log[7, lim/3^i]}]] // Sort, 3]; seq[10^9] (* Amiram Eldar, Nov 08 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
