cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A001677 Number of series-parallel networks with n edges.

Original entry on oeis.org

1, 2, 3, 6, 12, 26, 59, 146, 368, 976, 2667, 7482, 21440, 62622, 185637, 557680, 1694256, 5198142, 16086486, 50165218, 157510504, 497607008, 1580800091, 5047337994, 16190223624, 52153429218, 168657986843, 547389492416
Offset: 2

Views

Author

Keywords

Examples

			a(5) = 24 - (1/2)*(1*10+2*4+4*2+10*1) = 6.
		

References

  • 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).
  • B. D. H. Tellegen, Geometrical configurations and duality of electrical networks, Philips Technical Review, 5 (1940), 324-330.

Crossrefs

Programs

  • Mathematica
    m = 29; ClearAll[a, b, s]; a[1] = 1; a[2] = 2; a[3] = 4; b[1] = 1; b[n_ /; n >= 2] = a[n]/2; ex = Product[ 1/(1-x^k)^b[k], {k, 1, m}] - 1 - Sum[ a[k]*x^k, {k, 1, m}]; coes = CoefficientList[ Series[ ex, {x, 0, m}], x]; sol = Solve[ Thread[ coes == 0]][[1]]; Do[ s[k] = a[k] /. sol, {k, 1, m}]; a[2] = 1; a[3] = 2; a[n_] := s[n] - (1/2)*Sum[ s[i]*s[n-i], {i, 1, n-1}] - If[ OddQ[n], 0, s[n/2]/2]; Table[ a[n], {n, 2, m}] (* Jean-François Alcover, Feb 24 2012 *)

Formula

a(n) = s(n) - (1/2)*Sum_{i=1..n-1} s(i)*s(n-i) - (1/2)*s(n/2), where s() = A000084 and the last term is omitted if n is odd.

Extensions

More terms from David W. Wilson, Sep 20 2000