|
|
A026425
|
|
a(n) = least positive integer > a(n-1) and not a(i)*a(j) for 1<=i<=j<=n, n >= 2.
|
|
2
|
|
|
1, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 21, 22, 23, 26, 27, 29, 31, 33, 34, 37, 38, 39, 41, 43, 46, 47, 51, 53, 57, 58, 59, 61, 62, 67, 69, 71, 73, 74, 79, 82, 83, 86, 87, 89, 93, 94, 97, 101, 103, 106, 107, 109, 111, 113
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
MATHEMATICA
|
a = {1, 4}; used = {a[[1]]^2, a[[1]]*a[[2]], a[[2]]^2}; Do[k = a[[-1]] + 1; While[MemberQ[used, k], k++]; AppendTo[a, k]; used = Union[used, k*a], {n, 3, 60}]; a (* Ivan Neretin, Mar 07 2016 *)
|
|
CROSSREFS
|
Cf. A026416 and references therein.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|