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.

A285401 Positions of 0 in A285177; complement of A285402.

This page as a plain text file.
%I A285401 #21 Jul 27 2024 03:55:11
%S A285401 1,2,4,5,7,8,10,11,17,18,20,21,23,24,26,27,29,30,36,37,39,40,42,43,45,
%T A285401 46,48,49,55,56,58,59,61,62,64,65,67,68,74,75,81,82,88,89,95,96,102,
%U A285401 103,105,106,108,109,111,112,114,115,121,122,124,125,127,128
%N A285401 Positions of 0 in A285177; complement of A285402.
%C A285401 Conjecture: a(n)/n -> (61-sqrt(3))/26 = 2.279...
%C A285401 From _Michel Dekking_, Feb 10 2021: (Start)
%C A285401 This conjecture is false. In fact,
%C A285401        a(n)/n --> (5+sqrt(17))/4 = 2.28077...
%C A285401 Let mu be the defining morphism for A285177, i.e,
%C A285401       mu(0) = 11,   mu(1) = 001.
%C A285401 The sequence A285177 is the fixed point x = 0010010010011111... starting with 0 of mu^2:
%C A285401       mu^2(0) = 001001,   mu^2(1) = 1111001.
%C A285401 The 0's in x are at positions a(1)=1, a(2)=2, a(3)=4, etc.
%C A285401 Now suppose that N_0(K) = n is the number of 0's in a prefix x[1,K] of length K of x. Then obviously a(n) = K +/- 6.
%C A285401 Also N_0(K) + N_1(K) = K, where N_1(K) is the number of 1's in x[1,K].
%C A285401 So
%C A285401       K/N_0(K) = a(n)/n +/- 6/n.
%C A285401 Letting n tend to infinity, we find that
%C A285401       a(n)/n --> 1/f0,
%C A285401 where f0 is the frequency of 0's in x.
%C A285401 It is well known that these exist and are equal to the normalized eigenvector of the Perron-Frobenius eigenvalue of the incidence matrix of the morphism mu.
%C A285401 A simple computation yields that f0 = 4/(5+sqrt(17)).
%C A285401 It follows that a(n)/n --> (5+sqrt(17))/4.
%C A285401 (End)
%H A285401 Clark Kimberling, <a href="/A285401/b285401.txt">Table of n, a(n) for n = 1..10000</a>
%e A285401 As a word, A285177 = 001001..., in which 0 is in positions 1,2,4,5,7,...
%t A285401 s = Nest[Flatten[# /. {0 -> {1, 1}, 1 -> {0, 0, 1}}] &, {0}, 10] (* A285177 *)
%t A285401 Flatten[Position[s, 0]]  (* A285401 *)
%t A285401 Flatten[Position[s, 1]]  (* A285402 *)
%Y A285401 Cf. A285177, A285402, A285403.
%K A285401 nonn,easy
%O A285401 1,2
%A A285401 _Clark Kimberling_, Apr 26 2017