A121162 -id:A121162 - cp's OEIS frontend

A121331 Number of bridged bicyclic skeletons with n carbon atoms (see Parks et al. for precise definition).

1, 2, 6, 15, 39, 99, 258, 671, 1762, 4657, 12372, 33036, 88590, 238483, 644045, 1744542, 4737341, 12894158, 35165994, 96083192, 262951511, 720685274, 1977846334, 5434588909, 14949284828, 41163690109, 113451949753, 312955174089, 863965424349, 2386874582238

Offset: 5

N. J. A. Sloane, Aug 27 2006

nonn

Equivalently, the number of connected graphs on n unlabeled nodes with exactly 2 cycles of the same even length joined along half their length and all nodes having degree at most 4. The resulting graph will have three equal length cycles. - Andrew Howroyd, May 25 2018

			From _Andrew Howroyd_, May 25 2018: (Start)
Illustration of graphs for n=5 and n=6:
    o          o--o       o
   /|\        /|\        /|\
  o o o      o o o      o o o--o
   \|/        \|/        \|/
    o          o          o
.
Illustration of graphs for n=7:
    o--o      o--o--o   o--o        o        o          o   o
   /|\       /|\       /|\         /|\      /|\        /|\ /
  o o o     o o o     o o o--o    o o o    o o o--o   o o o
   \|/       \|/       \|/       / \|/ \    \|/   |    \|/ \
    o--o      o         o       o   o   o    o    o     o   o
(End)

Andrew Howroyd, Table of n, a(n) for n = 5..200
Camden A. Parks and James B. Hendrickson, Enumeration of monocyclic and bicyclic carbon skeletons, J. Chem. Inf. Comput. Sci., vol. 31, 334-339 (1991).

Cf. A121158, A121162, A121165, A305132.

G[n_] := Module[{g}, g[] = 0; Do[g[x] = 1 + x*(g[x]^3/6 + g[x^2]*g[x]/2 + g[x^3]/3) + O[x]^n // Normal, {n}]; g[x]];
C1[n_] := Sum[(d1^(3*k)+3*d1^k*d2^k + 2*d3^k), {k, 1, Quotient[n, 3]}]/12;
C2[n_] := Sum[(d1^Mod[k, 2]*d2^Quotient[k, 2])^3 + 3*d1^Mod[k, 2]* d2^(Quotient[k, 2] + k) + 2*d3^Mod[k, 2]*d6^Quotient[k, 2], {k, 1, Quotient[n, 3]}]/12;
seq[n_] := Module[{s, d, g}, s = G[n]; d = x*(s^2 + (s /. x -> x^2))/2; g[p_, e_] := Normal[(p+O[x]^(Quotient[n, e]+1))] /. x :> x^e; g[s, 1]^2* (C1[n-2] /. Thread[{d1, d2, d3} :> {g[d, 1], g[d, 2], g[d, 3]}]) + g[s, 2]*(C2[n-2] /. Thread[{d1, d2, d3, d6} :> {g[d, 1], g[d, 2], g[d, 3], g[d, 6]}]) + O[x]^n] // CoefficientList[#, x]& // Drop[#, 3]&;
seq[33] (* Jean-François Alcover, Sep 08 2019, after Andrew Howroyd *)

\\ here G is A000598 as series
G(n)={my(g=O(x)); for(n=1, n, g = 1 + x*(g^3/6 + subst(g, x, x^2)*g/2 + subst(g, x, x^3)/3) + O(x^n)); g}
C1(n)={sum(k=1, n\3, (d1^(3*k) + 3*d1^k*d2^k + 2*d3^k))/12}
C2(n)={sum(k=1, n\3, (d1^(k%2)*d2^(k\2))^3 + 3*d1^(k%2)*d2^(k\2+k) + 2*d3^(k%2)*d6^(k\2))/12}
seq(n)={my(s=G(n)); my(d=x*(s^2+subst(s, x, x^2))/2); my(g(p,e)=subst(p + O(x*x^(n\e)), x, x^e)); Vec(O(x^n/x) + g(s,1)^2*substvec(C1(n-2),[d1,d2,d3],[g(d,1), g(d,2), g(d,3)]) + g(s,2)*substvec(C2(n-2), [d1,d2,d3,d6], [g(d,1), g(d,2), g(d,3), g(d,6)]))} \\ Andrew Howroyd, May 25 2018

a(n) ~ c * d^n / sqrt(n), where d = 1/A261340 = 2.815460033176150746526616778..., c = 0.0064202170754... . - Vaclav Kotesovec, Sep 08 2019

Corrected by Franklin T. Adams-Watters and T. D. Noe, Oct 25 2006

a(24) corrected and terms a(26) and beyond from Andrew Howroyd, May 25 2018

A125669 Number of bicyclic skeletons with n carbon atoms and the parameter 'alpha' having the value of 0. See the paper by Hendrickson and Parks for details.

Original entry on oeis.org

1, 4, 20, 76, 288, 1005, 3433, 11324, 36712, 116809, 367076, 1140226, 3510491, 10722708, 32539364, 98178211, 294767639, 881147521, 2623934079, 7787024985, 23039064263, 67977412951, 200072442611, 587532484513, 1721812143140, 5036454320102, 14706743476128

Offset: 6

Parthasarathy Nambi, Jan 29 2007

nonn

Here 'alpha' is the number of atoms the two rings have in common.

Equivalently, the number of connected graphs on n unlabeled nodes with exactly 2 cycles without any shared node and all nodes having degree at most 4. See A121162 for the special case of the two cycles having the same length. - Andrew Howroyd, May 25 2018

			If n=6 then the number of bicyclics when 'alpha' = zero is 1.
If n=7 then the number of bicyclics when 'alpha' = zero is 4.
If n=8 then the number of bicyclics when 'alpha' = zero is 20.
If n=9 then the number of bicyclics when 'alpha' = zero is 76.
From _Andrew Howroyd_, May 25 2018: (Start)
Case n=7: illustrations of the 4 graphs:
     o   o   o       o   o   o       o   o---o       o   o---o
    / \ / \ / \     / \ /   / \     / \     / \     / \   \   \
   o---o   o---o   o---o---o---o   o---o---o---o   o---o---o---o
(End)

Andrew Howroyd, Table of n, a(n) for n = 6..200
J. B. Hendrickson and C. A. Parks, Generation and Enumeration of Carbon skeletons, J. Chem. Inf. Comput. Sci., 31 (1991), pp. 101-107. See Table VII column 2 on page 104.

Cf. A121162, A121330, A121331, A121332, A121333, A125064, A125670.

\\ here G is A000598 as series
G(n)={my(g=O(x)); for(n=1, n, g = 1 + x*(g^3/6 + subst(g, x, x^2)*g/2 + subst(g, x, x^3)/3) + O(x^n)); g}
C1(n)={sum(i=1, n\2-1, sum(j=1, n\2-i, (d1^(2*(i+j)) + 2*d1^(2*i)*d2^j + d2^(i+j))*(1 + d1)^2))/(8*(1-d1))}
C2(n)={sum(k=1, n\4,  2*(d2^(2*k) + d4^k)*(1 + d2))*(1+d1)/(8*(1-d2))}
seq(n)={my(s=G(n)); my(d=x*(s^2+subst(s, x, x^2))/2); my(g(p,e)=subst(p + O(x*x^(n\e)), x, x^e)); Vec(O(x^n/x) + g(s,1)^2*substvec(C1(n-2),[d1,d2],[g(d,1),g(d,2)]) + g(s,2)*substvec(C2(n-2), [d1,d2,d4], [g(d,1),g(d,2),g(d,4)]))} \\ Andrew Howroyd, May 25 2018

Terms a(16) and beyond from Andrew Howroyd, May 25 2018

cp's OEIS Frontend

A121331 Number of bridged bicyclic skeletons with n carbon atoms (see Parks et al. for precise definition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions