|
%I #8 Mar 07 2016 07:41:26
%S 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,
%T 33,34,37,38,39,41,43,46,47,49,51,53,57,58,59,61,62,67,69,71,73,74,79,
%U 81,82,83,86,87,89,93,94,97,101,103,106
%N a(n) = least positive integer > a(n-1) and not a(i)*a(j) for 1<=i<j<=n.
%H Ivan Neretin, <a href="/A026419/b026419.txt">Table of n, a(n) for n = 1..1000</a>
%t 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 *)
%Y Cf. A026416 and references therein.
%K nonn
%O 1,2
%A _Clark Kimberling_
|
