A097901 Duplicate of A026427.
3, 4, 5, 6, 7, 8, 10, 11, 13, 14, 17, 19, 22, 23, 26, 27, 29, 31, 34, 37, 38, 41, 43, 45, 46
Offset: 1
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
a(8) is not 10 because we already have 10 = 2*5. Of course all primes appear. 16 appears because 16 is not a product of earlier terms.
a[1]=1; a[2]=2; a[n_] := a[n] = For[k = a[n-1] + 1, True, k++, If[ FreeQ[ Table[ a[i]*a[j], {i, 1, n-2}, {j, i+1, n-1}], k], Return[k]]]; Table[a[n], {n, 1, 101}] (* Jean-François Alcover, May 16 2013 *)
from itertools import count, islice def agen(): # generator of terms a, products = [1, 2], {2} yield from a for k in count(3): if k not in products: yield k products.update(k*a[i] for i in range(len(a))) a.append(k) products.discard(k) print(list(islice(agen(), 62))) # Michael S. Branicky, Jun 09 2025
Comments