login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A007061 A maximally unpredictable sequence.
(Formerly M0075)
5
1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2, 1, 2, 1, 2, 1, 2, 2, 1, 1, 1, 2, 2, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2 (list; graph; refs; listen; history; internal format)
OFFSET

1,2

COMMENTS

Klaus Sutner remarks (Jun 26 2006) that this sequence is very similar to the Kimberling sequence A079101. Both sequences have every finite binary word as a factor; in fact, essentially the same proof works for both sequences.

Sutner continues: All words of length k seem to appear in the first 2^{k+2} bits. This is true for the first billion bits of the sequence, but no proof is known. The main open problem is whether the limiting density of 0's is 1/2. It seems to require a large amount of effort just to show that it is bounded away from 0, never mind some of the more exotic properties of the sequence (see the Sutner reference).

Start with a single bit 0. If the first n bits U(n) = a(1)a(2)...a(n) have already been chosen, let v be the longest suffix of U(n) that already appears in U(n-1). Find the last occurrence of v in U(n-1) and let b the bit that occurs immediately after. Then a(n+1) is the complement of b. (The entry gives the bits as 1's and 2s instead of 0's and 1's - compare A038219) - Joshua Zucker (joshua.zucker(AT)stanfordalumni.org), Aug 11 2006

REFERENCES

A. Ehrenfeucht and J. Mycielski, A pseudorandom sequence - how random is it?, Amer. Math. Monthly, 99 (1992), 373-375.

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

LINKS

Joshua Zucker, Table of n, a(n) for n = 1..1999

K. Sutner, The Ehrenfeucht-Mycielski sequence

PROG

(Haskell)

a007061 n = a007061_list !! (n-1)

a007061_list = 1 : f [1] where

   f us = a' : f (us ++ [a']) where

     a' = b $ reverse $ map (`splitAt` us) [0..length us - 1] where

        b ((xs, ys):xyss) | vs `isSuffixOf` xs = 3 - head ys

                         | otherwise          = b xyss

     vs = fromJust $ find (`isInfixOf` init us) $ tails us

-- Reinhard Zumkeller, Dec 05 2011

CROSSREFS

Cf. A038219 (0-1 version), A079101.

Cf. A201881 (run lengths).

Sequence in context: A100428 A206719 A093914 * A001817 A091954 A183953

Adjacent sequences:  A007058 A007059 A007060 * A007062 A007063 A007064

KEYWORD

nonn

AUTHOR

N. J. A. Sloane (njas(AT)research.att.com), Mira Bernstein (mira(AT)math.berkeley.edu)

EXTENSIONS

More terms from Joshua Zucker (joshua.zucker(AT)stanfordalumni.org), Aug 11 2006

Offset changed from 0 to 1, Aug 18 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 13 06:53 EST 2012. Contains 205451 sequences.