OFFSET
1,2
COMMENTS
Equal to nonzero terms of A099308 when sorted into ascending order. In this order, which is dictated by the primorial base expansion of n (A049345) and mapped to products of prime powers by A276086, all terms of A099308 that are prime(k)-smooth appear before the terms that are not prime(k)-smooth.
Number of terms whose greatest prime factor (A006530) is prime(n) [in other words, that are prime(n)-smooth but not prime(n-1)-smooth] is given by A351071(n): 1, 4, 8, 44, 216, 1474, 11130, ...
For all n > 1, A003415(a(n)) is also a term of the sequence.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..12878 (all 17-smooth terms of this sequence)
Antti Karttunen, 105368 initial terms, without indices (all 19-smooth terms of this sequence, and also A276086(9699690) = 23, the first 23-smooth term)
PROG
(PARI)
A003415checked(n) = if(n<=1, 0, my(f=factor(n), s=0); for(i=1, #f~, if(f[i, 2]>=f[i, 1], return(0), s += f[i, 2]/f[i, 1])); (n*s));
A099307(n) = { my(s=1); while(n>1, n = A003415checked(n); s++); if(n, s, 0); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Feb 10 2022
STATUS
approved