|
%I #17 Mar 04 2021 16:22:15
%S 6,10,14,7,8,31,35,17,43,40,229,248,212,818,799,733,151,2191,1139,
%T 20894,877,6835,20528,34627,19687,91790,34502,367558,85336,46375,
%U 1342349,134683,109057,2758327,5921086,1655564,18147329,11934733,1315376
%N Least k >= 6 such that A087666(k) = n.
%C A087666: Consider the recurrence b(0) = n/3, b(n) = b(n-1)*floor(b(n-1)); sequence gives number of steps to reach an integer, or -1 if no integer is ever reached. - _Robert G. Wilson v_, Oct 10 2003 & Mar 10 2004
%H Chai Wah Wu, <a href="/A087710/b087710.txt">Table of n, a(n) for n = 0..60</a>
%p See the Maple program in A087666 for the best way to compute this sequence. - _N. J. A. Sloane_
%t f[n_] := Block[{c = 1, k = n/3}, If[ IntegerQ[k], 0, While[tn = Floor[k]; tf = k - tn; tn = Mod[tn, 3^100]; k = tn(tn + tf); ! IntegerQ[k], c++ ]; c++ ]]; a = Table[0, {50}]; Do[ b = f[n]; If[ a[[b + 1]] == 0, a[[b + 1]] = n; Print[b, " = ", n]], {n, 6, 10^7}]; a (* _Robert G. Wilson v_, Mar 10 2004, using the idea from _N. J. A. Sloane_'s Maple code in A087666)
%Y Cf. A087666, A086336, A087663.
%K nonn
%O 0,1
%A _Benoit Cloitre_, Sep 30 2003
%E More terms from _Robert G. Wilson v_, Oct 10 2003, Mar 10 2004
|
