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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A153003 Toothpick sequence in the first three quadrants. 10
0, 1, 4, 7, 10, 16, 25, 31, 34, 40, 49, 58, 70, 91, 115, 127, 130, 136, 145, 154, 166, 187, 211, 226, 238, 259, 286, 316, 361, 427, 487, 511, 514, 520, 529, 538, 550, 571, 595, 610, 622, 643, 670, 700, 745, 811, 871, 898, 910, 931 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

From Omar E. Pol, Oct 01 2011: (Start)

On the infinite square grid, consider only the first three quadrants and count only the toothpicks of length 2.

At stage 0, we start from a vertical half toothpick at [(0,0),(0,1)]. This half toothpick represents one of the two components of the first toothpick placed in the toothpick structure of A139250, so a(0) = 0.

At stage 1, we place an orthogonal toothpick of length 2 centered at the end, so a(1) = 1. Also we place half toothpick at [(0,-1),(1,-1)]. This last half toothpick represents one of the two components of the third toothpick placed in the toothpick structure of A139250.

At stage 2, we place three toothpicks, so a(2) = 1+3 = 4.

In each subsequent stage, for every exposed toothpick end, place an orthogonal toothpick centered at that end.

The sequence gives the number of toothpicks after n stages. A153004 (the first differences) gives the number of toothpicks added to the structure at n-th stage.

Note that this sequence is different from the toothpick "corner" sequence A153006. For more information see A139250. (End)

LINKS

Table of n, a(n) for n=0..49.

David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata

N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS

Index entries for sequences related to toothpick sequences

FORMULA

a(n) = (A139250(n+1)-3)*3/4 + 1, if n >= 1.

From Omar E. Pol, Oct 01 2011: (Start)

a(n) = A139250(n+1) - A152998(n) + A153000(n-1) - 1, if n >= 1.

a(n) = A139250(n+1) - A153000(n-1) - 2, if n >= 1.

a(n) = A152998(n) + A153000(n-1), if n >= 1.

(End)

PROG

(Python)

def msb(n):

    t=0

    while n>>t>0: t+=1

    return 2**(t - 1)

def a139250(n):

    k=(2*msb(n)**2 + 1)/3

    return 0 if n==0 else k if n==msb(n) else k + 2*a139250(n - msb(n)) + a139250(n - msb(n) + 1) - 1

def a(n): return 0 if n==0 else (a139250(n + 1) - 3)*3/4 + 1

print [a(n) for n in xrange(101)] # Indranil Ghosh, Jul 01 2017

CROSSREFS

Cf. A139250, A139251, A152968, A152978, A152998, A153000, A153004.

Sequence in context: A299475 A310713 A180080 * A213484 A128429 A191154

Adjacent sequences:  A153000 A153001 A153002 * A153004 A153005 A153006

KEYWORD

nonn

AUTHOR

Omar E. Pol, Jan 02 2009

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 20 01:17 EDT 2018. Contains 316378 sequences. (Running on oeis4.)