A026419
a(n) = least positive integer > a(n-1) and not a(i)*a(j) for 1<=i
1, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 31, 33, 34, 37, 38, 39, 41, 43, 46, 47, 49, 51, 53, 57, 58, 59, 61, 62, 67, 69, 71, 73, 74, 79, 81, 82, 83, 86, 87, 89, 93, 94, 97, 101, 103, 106
Offset: 1
Keywords
Links
- Ivan Neretin, Table of n, a(n) for n = 1..1000
Crossrefs
Cf. A026416 and references therein.
Programs
-
Mathematica
a = {1, 4}; used = {a[[1]]*a[[2]]}; Do[k = a[[-1]] + 1; While[MemberQ[used, k], k++]; used = Union[used, k*a]; AppendTo[a, k], {n, 3, 60}]; a (* Ivan Neretin, Mar 07 2016 *)
Extensions
Name clarified by Robert C. Lyons, Feb 08 2025