|
%I #25 Aug 01 2018 09:21:25
%S 1,1,1,5,1,5,1,6,1,7,1,12,8,10,1,17,8,1,13,13,13,5,10,11,5,9,8,19,10,
%T 11,7,11,5,9,27,9,13,5,23,5,9,17,9,11,11,7,21,9,7,5,17,27,11,7,9,17,5,
%U 13,9,21,11,7,13,9,9
%N Lengths of runs of identical terms in A251416.
%C It would be nice to have an alternative description of this sequence, one that is not based on A098550.
%C It appears (conjecture) that a(n)>1 for n>18. - _Alexander R. Povolotsky_, Dec 07 2014
%C Conjecture: a(n) = A247253(n-5) for n>12. - _Reinhard Zumkeller_, Dec 07 2014
%C The previous conjecture is equivalent to the statement that A251416(n) lists all primes and only primes after a(30)=18. - _M. F. Hasler_, Dec 08 2014
%H Reinhard Zumkeller, <a href="/A251417/b251417.txt">Table of n, a(n) for n = 1..10000</a>
%F Let f(n)=A098551(A251595(n)). Then one can prove that A251417(n) = f(n) - f(n-1), n>=2. - _Vladimir Shevelev_, Dec 09 2014
%e See A251595.
%t termsOfA251416 = 700;
%t f[lst_List] := Block[{k = 4}, While[GCD[lst[[-2]], k] == 1 || GCD[lst[[-1]], k] > 1 || MemberQ[lst, k], k++]; Append[lst, k]];
%t A098550 = Nest[f, {1, 2, 3}, termsOfA251416 - 3];
%t b[1] = 2;
%t b[n_] := b[n] = For[k = b[n-1], True, k++, If[FreeQ[A098550[[1 ;; n]], k], Return[k]]];
%t A251416 = Array[b, termsOfA251416];
%t Length /@ Split[A251416] (* _Jean-François Alcover_, Aug 01 2018, after _Robert G. Wilson v_ *)
%o (Haskell)
%o import Data.List (group)
%o a251417 n = a251417_list !! (n-1)
%o a251417_list = map length $ group a251416_list
%o -- _Reinhard Zumkeller_, Dec 05 2014
%Y Cf. A098550, A251416, A251595, A247253.
%K nonn
%O 1,4
%A _N. J. A. Sloane_, Dec 03 2014
|
