cp's OEIS Frontend

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.

A105600 Assume the conjectured terms of A105594 are the correct beginnings of the trajectories described in A003508. a(n) is a record length of b(n) iterations to arrive at the collected trajectories. This sequence cites the a(n)'s.

This page as a plain text file.
%I A105600 #6 May 22 2013 15:12:15
%S A105600 1,5,9,16,25,43,91,105,427,463,484,4085,4306,4413,5583,6273,10172,
%T A105600 18105,24946,31686,31886
%N A105600 Assume the conjectured terms of A105594 are the correct beginnings of the trajectories described in A003508. a(n) is a record length of b(n) iterations to arrive at the collected trajectories. This sequence cites the a(n)'s.
%C A105600 The trajectory in A003508, etc., is defined as a(1)=n, for n>1, a(n) = a(n-1) + 1 + sum of distinct prime factors of a(n-1) that are < a(n-1).
%t A105600 a[1] = 1; a[n_] := a[n] = a[n - 1] + 1 + Plus @@ Select[ Flatten[ Table[ #[[1]], {1}] & /@ FactorInteger[ a[n - 1]]], # < a[n - 1] &]; t = Table[ a[n], {n, 1500}]; f[n_] := Module[{b, k = 1}, b[1] = n; b[m_] := b[m] = b[m - 1] + 1 + Plus @@ Select[ Flatten[ Table[ #[[1]], {1}] & /@ FactorInteger[ b[m - 1]]], # < b[m - 1] &]; While[ Position[t, b[k]] == {} && k < 1000, k++ ]; If[ k == 1000, t = Select[ Union[ Join[t, Table[ b[i], {i, 2, k}]]], # > n &]; -1, k - 1]]; lst = {{1, 0}}; Do[d = f[n]; If[d > lst[[ -1, 2]], AppendTo[lst, {n, d}]], {n, 60000}]; Transpose[ lst][[1]]
%Y A105600 Cf. A105593, the b(n)'s are in A105600.
%K A105600 nonn
%O A105600 0,2
%A A105600 _R. K. Guy_ and _Robert G. Wilson v_, Apr 15 2005