A088178 Sequence of distinct products b(n)*b(n+1), n=1,2,3,..., of the terms b(n) of A088177.
1, 2, 4, 6, 3, 5, 10, 8, 12, 9, 15, 20, 16, 24, 18, 21, 7, 11, 22, 14, 28, 32, 40, 25, 30, 36, 42, 35, 45, 27, 33, 44, 48, 60, 50, 70, 49, 56, 64, 72, 54, 66, 55, 65, 13, 17, 34, 26, 39, 51, 68, 52, 78, 84, 98, 63, 81, 90, 80, 88, 77, 91, 104, 96, 108, 99, 110, 100, 120, 132
Offset: 1
Links
- Ivan Neretin, Table of n, a(n) for n = 1..10000 (first 1000 terms from Michael De Vlieger)
Programs
-
Mathematica
a088177[n_Integer] := Module[{t = {1, 1}}, Do[AppendTo[t, 1]; While[Length[Union[Most[t]*Rest[t]]] < i - 1, t[[-1]]++], {i, 3, n}]; t]; a088178[n_Integer] := Last[a088177[n]]*Last[a088177[n + 1]]; a088178 /@ Range[120] (* Michael De Vlieger, Aug 30 2014, based on T. D. Noe's script at A088177 *)
-
Python
from itertools import islice def A088178(): # generator of terms yield 1 p, a = {1}, 1 while True: n, na = 1, a while na in p: n += 1 na += a p.add(na) a = n yield na A088178_list = list(islice(A088178(),20)) # Chai Wah Wu, Oct 21 2021
Formula
a(m) < a(n)^2 for m < n. - Thomas Ordowski, Sep 02 2014
Extensions
Edited by N. J. A. Sloane, Oct 18 2021
Comments