A026422 a(n) = least positive integer > a(n-1) and not a(i)*a(j) for 1 <= i <= j < n.

1, 2, 3, 5, 7, 8, 11, 12, 13, 17, 18, 19, 20, 23, 27, 28, 29, 30, 31, 32, 37, 41, 42, 43, 44, 45, 47, 48, 50, 52, 53, 59, 61, 63, 66, 67, 68, 70, 71, 72, 73, 75, 76, 78, 79, 80, 83, 89, 92, 97, 98, 99, 101, 102, 103, 105, 107, 108, 109, 110

Offset

1,2

Links

Vincenzo Librandi, Table of n, a(n) for n = 1..5000

Mathematica

a[1]=1; a[n_] := a[n] = (t= Union[Flatten[Table[a[i] a[j], {i, 1, n-1}, {j, i, n-1}]]]; Do[If[FreeQ[t, k], an = k; Break[]], {k, a[n-1]+1, Last[t]+1}]; an); Array[a, 60] (* Jean-François Alcover, May 06 2011 *)
Select[Range[110], OddQ[Total[FactorInteger[#]][[2]]] &] (* T. D. Noe, May 07 2011 *)
g = 110; t = Array[1 &, g];
Table[If[t[[j]] == 1, t[[j*i]] = 0, t[[i*j]] = 1], {j, 2, g/2}, {i, 2, g/j}]; Flatten[Position[t, 1]] (* Horst H. Manninger, Mar 15 2023 *)

Prog

(PARI) is(n)=bigomega(n)%2 || n==1 \\ Charles R Greathouse IV, Sep 16 2015

Crossrefs

Apart from initial term, same as A026424.

Cf. A026416 and references therein.

Sequence in context: A245303 A359766 A166982 * A026424 A298207 A347453

Adjacent sequences: A026419 A026420 A026421 * A026423 A026424 A026425

Keyword

nonn

Author

Clark Kimberling

Extensions

Name corrected by Charles R Greathouse IV, Sep 16 2015

Status

approved