|
|
A000621
|
|
Number of monosubstituted alkanes C(n-1)H(2n-1)-X with n-1 carbon atoms that are not stereoisomers.
(Formerly M0697 N0258)
|
|
28
|
|
|
1, 1, 1, 2, 3, 5, 8, 14, 23, 39, 65, 110, 184, 310, 520, 876, 1471, 2475, 4159, 6996, 11759, 19775, 33244, 55902, 93984, 158030, 265696, 446746, 751128, 1262940, 2123444, 3570318, 6002983, 10093259, 16970431, 28533590, 47975381, 80664329, 135626284, 228037752
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
Also number of monosubstituted alkanes C(n)H(2n+1)-X of the form R-CH2-X (primary) that are not stereoisomers.
Let the entries in the nine columns of Blair and Henze's Table I (JACS 54 (1932), p. 1098) be denoted by Ps(n), Pn(n), Ss(n), Sn(n), Ts(n), Tn(n), As(n), An(n), T(n) respectively (here P = Primary, S = Secondary, T = Tertiary, s = stereoisomers, n = non-stereoisomers and the last column T(n) gives total).
|
|
REFERENCES
|
S. R. Finch, Mathematical Constants, Cambridge, 2003, p. 300.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: A(x) satisfies A(x) = 1/(1-x*A(x^2)), with offset 0. - Paul D. Hanna, Aug 16 2002
Given g.f. A(x), then B(x) = A(x) / x satisfies 0 = f(B(x), B(x^2), B(x^4)) where f(u, v, w) = (1 - u)^2 * w - u^2 * v * (v - 1). - Michael Somos, Sep 03 2007
G.f.: x / (1 - x / (1 - x^2 / (1 - x^4 / (1 - ...)))). - Michael Somos, Sep 03 2007
For offset 0 (as considered in the 1937 Polya reference) we have
G.f.: A(x) = 1 / (1 - x / (1 - x^2 / (1 - x^4 / (1 - ...)))) and
A(x) satisfies A(x) = 1 + x*A(x)*A(x^2) (equivalent to Hanna's functional equation).
(End)
a(n) ~ c * beta^n, where beta = 1.681367524441880255591... (see A239804), c = 0.214536139134648555630... (see A239806). Asymptotic formula a(n) ~ K * beta^n from reference (Analytic Combinatorics, p. 283), where K = 0.3607140971, beta = 1.6813675244^n is for offset 0 (beta is same, but K = c * beta = 0.360714097160142828...). - Vaclav Kotesovec, Mar 27 2014
a(n) = T(2*n-1,1), where T(n,m) = Sum_{i=1..n-m} binomial(i+m-1,i)*((1+(-1)^(n-m))/2)*T((n-m)/2,i), n > m, T(n,n)=1. - Vladimir Kruchinin, Mar 18 2015
|
|
EXAMPLE
|
G.f. = x + x^2 + x^3 + 2*x^4 + 3*x^5 + 5*x^6 + 8*x^7 + 14*x^8 + 23*x^9 + ...
|
|
MATHEMATICA
|
nmax=40; a=1-x; Do[a=1/(1-x (a/.x->x^2)), {Log[2, nmax]+2}]; CoefficientList[Series[a, {x, 0, nmax-1}], x] (* Jean-François Alcover, Jun 16 2011, after Michael Somos, fixed by Vaclav Kotesovec, Mar 28 2014 *)
max = 40; cf = Fold[Function[1 - x^#2/#1], 1, 2^Reverse[Range[0, Floor[Log[2, max]]]]]; List @@ (1-Series[cf, {x, 0, 2*max}] // Normal) /. x -> 1 (* Jean-François Alcover, Sep 24 2014 *)
|
|
PROG
|
(PARI) {a(n) = my(A, m); if( n<1, 0, n--; m = 1; A = 1 + O(x); while( m<=n, m *= 2; A = 1 / (1 - x * subst(A, x, x^2)) ); polcoeff( A, n )) }; /* Michael Somos, Sep 03 2007 */
(Maxima)
T(n, m):=if m=n then 1 else sum(binomial(i+m-1, i)*((1+(-1)^(n-m))/2)*T((n-m)/2, i), i, 1, n-m);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,nice,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Additional comments from Bruce Corrigan, Nov 04 2002
|
|
STATUS
|
approved
|
|
|
|