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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A135106 Number of physical trees of alkane structures with n carbon vertices. 1
1, 1, 2, 6, 24, 118, 686, 4598, 34872, 295044, 2753958, 28103804, 311216626, 3716341042, 47597786154, 650812077852, 9461423560788, 145724617925326, 2370293673319292 (list; graph; refs; listen; history; internal format)
OFFSET

1,3

COMMENTS

Similar to A000602 (alkane trees with n carbon atoms) but keeping track of the history of attaching carbon atoms (methyls) to the backbone, as if these had been labeled.

LINKS

J. V. Knop, W. R. Muller, K. Szymanski et al., Computer enumeration and generation of physical trees, J. Comput. Chem. vol. 8 no. 4 (1987) pp 549-554.

EXAMPLE

Starting with a(1)=1, one C1 methane, we get a(2)=1, the C1-C2 backbone.

The third can be attached to either C1 ending up with C3-C1-C2, or to C2 ending up with C1-C2-C3, yielding a(3)=2 different propanes.

C4 may be attached to any of C1 to C3 in these two propanes, yielding a(4)=6 different butanes, four of which are linear and two of which are stars.

MAPLE

A135106 := proc(n) local numb, stack, istack, N, i ; numb := array(1..n) ; for i from 1 to n do numb[i] := 0 ; od: stack := array(1..7, 1..100) ; stack[1, 1]:=2 ; stack[2, 1]:=0 ; stack[3, 1]:=0 ; stack[4, 1]:=0 ; stack[5, 1]:=1 ; stack[6, 1]:=0 ; stack[7, 1]:=1 ; istack := 1 ; while istack <> 0 do for i from 1 to 7 do stack[i, istack+1] := stack[i, istack] ; od: if stack[6, istack] = 3 then istack := istack-1 ; else stack[6, istack] := stack[6, istack]+1 ; stack[1, istack+1] := stack[1, istack]+1 ; N := stack[6, istack] ; if stack[N, istack] <> 0 then stack[N, istack+1] := stack[N, istack+1]-1 ; stack[N+1, istack+1] := stack[N+1, istack+1]+1 ; stack[5, istack+1] := stack[N, istack]*stack[5, istack] ; stack[6, istack+1] := 0 ; stack[7, istack+1] := stack[7, istack]+1 ; numb[stack[7, istack+1]]:=numb[stack[7, istack+1]]+stack[5, istack+1] ; if stack[7, istack+1] <> n then istack := istack+1 ; fi ; fi ; fi ; od: numb[n] ; end: for n from 2 do print( A135106(n)) ; end: - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 18 2008

CROSSREFS

Sequence in context: A079106 A088713 A193938 * A005394 A095818 A052397

Adjacent sequences:  A135103 A135104 A135105 * A135107 A135108 A135109

KEYWORD

more,nonn

AUTHOR

R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 12 2008

EXTENSIONS

More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 18 2008

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 15 21:56 EST 2012. Contains 205860 sequences.