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.

Showing 1-10 of 79 results. Next

A280777 Denominators of coefficients in asymptotic expansion of C_n (number of connected chord diagrams, A000699).

Original entry on oeis.org

1, 2, 8, 16, 128, 1280, 3072, 215040, 3440640, 2293760, 137625600, 201850880, 72666316800, 1889324236800, 52901078630400, 176336928768000, 135426761293824, 191854578499584000, 1593868805996544000, 787371190162292736000, 31494847606491709440000, 29395191099392262144000
Offset: 0

Views

Author

N. J. A. Sloane, Jan 19 2017

Keywords

Examples

			Coefficients are 1, -5/2, -43/8, -579/16, -44477/128, -5326191/1280, -180306541/3072, ...
		

Crossrefs

Programs

  • PARI
    A000699_seq(N) = {
      my(a = vector(N)); a[1] = 1;
      for (n=2, N, a[n] = sum(k=1, n-1, (2*k-1)*a[k]*a[n-k])); a;
    };
    seq(N) = my(C = 'x*Ser(A000699_seq(N))); Vec(x*exp(1-(2*C+C^2)/(2*x))/C);
    apply(denominator, seq(22)) \\ Gheorghe Coserea, Jan 22 2017

Extensions

More terms from Gheorghe Coserea, Jan 22 2017

A178685 G.f.: A(x) = x/Series_Reversion(G(x)) where x*A(G(x)) = G(x) is the g.f. of A000699.

Original entry on oeis.org

1, 1, 3, 17, 143, 1569, 20911, 324673, 5720319, 112365569, 2430174335, 57327084033, 1464259637503, 40255921184257, 1185312426797823, 37219582236845057, 1241715255479353343, 43867539372405948417, 1636223622007966705663, 64261630682976366944257
Offset: 0

Views

Author

Paul D. Hanna, Dec 31 2010

Keywords

Comments

The limit of (a(n+1)/a(n))/n = 2. - Paul D. Hanna, Aug 08 2014
The g.f. of related sequence A000699, where A000699(n) is the number of irreducible diagrams with 2n nodes, satisfies: G(x) = x + x^2*[d/dx G(x)^2/x].

Examples

			G.f.: A(x) = 1 + x + 3*x^2 + 17*x^3 + 143*x^4 + 1569*x^5 + 20911*x^6 +...
A(x) satisfies: A(G(x)) = G(x)/x where G(x) is the g.f. of A000699:
G(x) = x + x^2 + 4*x^3 + 27*x^4 + 248*x^5 + 2830*x^6 + 38232*x^7 +...
which in turn satisfies: G(x) = x + x^2*[d/dx G(x)^2/x].
...
The table of coefficients in the initial powers of g.f. A(x) begins:
A^1: [1, 1,  3,  17,  143,  1569,  20911,  324673,  5720319, ...];
A^2: [1, 2,  7,  40,  329,  3526,  46107,  705444, 12289245, ...];
A^3: [1, 3, 12,  70,  567,  5949,  76350, 1150920, 19818801, ...];
A^4: [1, 4, 18, 108,  867,  8928, 112524, 1671000, 28436085, ...];
A^5: [1, 5, 25, 155, 1240, 12566, 155650, 2277050, 38285275, ...];
A^6: [1, 6, 33, 212, 1698, 16980, 206902, 2982096, 49529871, ...];
A^7: [1, 7, 42, 280, 2254, 22302, 267624, 3801036, 62355195, ...];
A^8: [1, 8, 52, 360, 2922, 28680, 339348, 4750872, 76971171, ...];
...
In the above table, notice that the main diagonal:
[1, 2, 12, 108, 1240, 16980, 267624, 4750872, 93615408, ...]
is related to the secondary diagonal in the following way:
[2*1, 4*3, 6*18, 8*155, 10*1698, 12*22302, 14*339348, 16*5850963, ...].
		

Crossrefs

Cf. A000699.

Programs

  • PARI
    {a(n)=local(A=[1,1]);for(i=2,n,A=concat(A,0);A[#A]=(2*(#A-1)*Vec(Ser(A)^(#A))[#A-1]-Vec(Ser(A)^(#A))[#A])/(#A));A[n+1]}
    for(n=0,25,print1(a(n),", "))
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A = 1 + x*A*(A + x*A')/(A - x*A' +x*O(x^n))); polcoeff(A,n)}
    for(n=0,25,print1(a(n),", "))

Formula

G.f. satisfies:
(1) [x^(n+1)] A(x)^(n+2) = 2(n+1)*[x^n] A(x)^(n+2).
(2) A(x) = 1 + x*A(x) * (A(x) + x*A'(x)) / (A(x) - x*A'(x)). - Paul D. Hanna, Aug 08 2014

A280776 Numerators of coefficients in asymptotic expansion of C_n (number of connected chord diagrams, A000699).

Original entry on oeis.org

1, -5, -43, -579, -44477, -5326191, -180306541, -203331297947, -58726239094693, -781618285277957, -1025587838964854273, -35763822710356866613, -330773478104531041960421, -237504847171108896327033959, -196526060612842999084524774697, -20633624138373135772483762873819
Offset: 0

Views

Author

N. J. A. Sloane, Jan 19 2017

Keywords

Examples

			Coefficients are 1, -5/2, -43/8, -579/16, -44477/128, -5326191/1280, -180306541/3072, ...
		

Crossrefs

Programs

  • PARI
    A000699_seq(N) = {
      my(a = vector(N)); a[1] = 1;
      for (n=2, N, a[n] = sum(k=1, n-1, (2*k-1)*a[k]*a[n-k])); a;
    };
    seq(N) = my(C = 'x*Ser(A000699_seq(N))); Vec(x*exp(1-(2*C+C^2)/(2*x))/C);
    apply(numerator, seq(16))  \\ Gheorghe Coserea, Jan 22 2017

Extensions

More terms from Gheorghe Coserea, Jan 22 2017

A212273 a(n) = n * A000699(n).

Original entry on oeis.org

1, 2, 12, 108, 1240, 16980, 267624, 4750872, 93615408, 2026018980, 47764893000, 1218621868200, 33459133580496, 983950866169544, 30862364045494800, 1028656519404191280, 36312294166359724128, 1353561646509338629572, 53133041579508472510728
Offset: 1

Views

Author

Michael Somos, May 12 2012

Keywords

Examples

			G.f. = x + 2*x^2 + 12*x^3 + 108*x^4 + 1240*x^5 + 16980*x^6 + 267624*x^7 + 4750872*x^8 + ...
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A); A = O(x); for(i=1, n, A = x + A * (2 * x * A' - A)); n * polcoeff(A, n)};

Formula

a(n) = [x^n] A178685(x)^n.
a(n) ~ n^(n+1) * 2^(n+1/2) / exp(n+1). - Vaclav Kotesovec, Feb 22 2014

A304685 a(n) = A000699(n) (mod 3).

Original entry on oeis.org

1, 1, 1, 0, 2, 1, 0, 0, 1, 0, 0, 0, 0, 1, 2, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0
Offset: 1

Views

Author

John M. Campbell, May 16 2018

Keywords

Comments

This is a 3-automatic sequence.

Examples

			We have that A000699(5) = 248, with 248 == 2 (mod 3), and the ternary expansion of 5 is given by the tuple (1, 2), so according to the above formula we have that a(5) = 2.
		

Crossrefs

Cf. A000699.

Programs

  • GAP
    a:=[1];;  for n in [2..90] do a[n]:=(n-1)*Sum([1..n-1],i->a[i]*a[n-i]) mod 3; od; a; # Muniru A Asiru, Aug 15 2018
  • Maple
    a:=proc(n) option remember; if n<=1 then 1; else
    add((2*i-1)*procname(i)*procname(n-i),i=1..n-1) mod 3; end if; end proc:
    seq(a(n), n=1..90); # Muniru A Asiru, Aug 15 2018
  • PARI
    a(n) = {my(A); A = O(x) ; for( i=1, n, A = x + A * (2 * x * A' - A)); polcoeff(A, n) % 3}; \\ Michel Marcus, Jul 04 2018; after A000699
    

Formula

For a natural number n, we have that a(n) = 1 if the ternary expansion of n is of the form 100...0 or is of the form 11...1200...0 for an even number of ones in this latter case, allowing runs of integers to be of length 0; a(n) = 2 if the ternary expansion of n is of the form 11...1200...0 if the length of the run of ones is odd; otherwise, a(n) = 0.

A051863 Duplicate of A000699.

Original entry on oeis.org

1, 1, 4, 27, 248, 2830, 38232, 593859, 10401712, 202601898, 4342263000
Offset: 1

Views

Author

Keywords

A174432 Partial sums of A000699.

Original entry on oeis.org

1, 2, 6, 33, 281, 3111, 41343, 635202, 11036914, 213638812, 4555901812, 106107724162, 2679887230354, 72962091956750, 2130453028323070, 66421485491085025, 2202438789394598209, 77400308039913410963
Offset: 1

Views

Author

Jonathan Vos Post, Mar 19 2010

Keywords

Comments

Partial sums of number of irreducible diagrams (in the sense of perturbation expansion in quantum field theory: spinor case in 4 spacetime dimensions) with 2n nodes. The subsequence of primes in this partial sum begins: 2, 281, and there never seem to be any more, as the underlying sequence is all even after A000699(8).

Examples

			a(18) = 1 + 1 + 4 + 27 + 248 + 2830 + 38232 + 593859 + 10401712 + 202601898 + 4342263000 + 101551822350 + 2573779506192 + 70282204726396 + 2057490936366320 + 64291032462761955 + 2136017303903513184 + 75197869250518812754.
		

Crossrefs

Formula

a(n) = SUM[i=1..n] A000699(i), where those A000699(n) = (n-1)*SUM[i=1..n-1] A000699(i)*A000699(n-i).

A036987 Fredholm-Rueppel sequence.

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Keywords

Comments

Binary representation of the Kempner-Mahler number Sum_{k>=0} 1/2^(2^k) = A007404.
a(n) = (product of digits of n; n in binary notation) mod 2. This sequence is a transformation of the Thue-Morse sequence (A010060), since there exists a function f such that f(sum of digits of n) = (product of digits of n). - Ctibor O. Zizka, Feb 12 2008
a(n-1), n >= 1, the characteristic sequence for powers of 2, A000079, is the unique solution of the following formal product and formal power series identity: Product_{j>=1} (1 + a(j-1)*x^j) = 1 + Sum_{k>=1} x^k = 1/(1-x). The product is therefore Product_{l>=1} (1 + x^(2^l)). Proof. Compare coefficients of x^n and use the binary representation of n. Uniqueness follows from the recurrence relation given for the general case under A147542. - Wolfdieter Lang, Mar 05 2009
a(n) is also the number of orbits of length n for the map x -> 1-cx^2 on [-1,1] at the Feigenbaum critical value c=1.401155... . - Thomas Ward, Apr 08 2009
A054525 (Mobius transform) * A001511 = A036987 = A047999^(-1) * A001511 = the inverse of Sierpiński's gasket * the ruler sequence. - Gary W. Adamson, Oct 26 2009 [Of course this is only vaguely correct depending on how the fuzzy indexing in these formulas is made concrete. - R. J. Mathar, Jun 20 2014]
Characteristic function of A000225. - Reinhard Zumkeller, Mar 06 2012
Also parity of the Catalan numbers A000108. - Omar E. Pol, Jan 17 2012
For n >= 2, also the largest exponent k >= 0 such that n^k in binary notation does not contain both 0 and 1. Unlike for the decimal version of this sequence, A062518, where the terms are only conjectural, for this sequence the values of a(n) can be proved to be the characteristic function of A000225, as follows: n^k will contain both 0 and 1 unless n^k = 2^r-1 for some r. But this is a special case of Catalan's equation x^p = y^q-1, which was proved by Preda Mihăilescu to have no nontrivial solution except 2^3 = 3^2 - 1. - Christopher J. Smyth, Aug 22 2014
Image, under the coding a,b -> 1; c -> 0, of the fixed point, starting with a, of the morphism a -> ab, b -> cb, c -> cc. - Jeffrey Shallit, May 14 2016
Number of nonisomorphic Boolean algebras of order n+1. - Jianing Song, Jan 23 2020

Examples

			G.f. = 1 + x + x^3 + x^7 + x^15 + x^31 + x^63 + x^127 + x^255 + x^511 + ...
a(7) = 1 since 7 = 2^3 - 1, while a(10) = 0 since 10 is not of the form 2^k - 1 for any integer k.
		

Crossrefs

The first row of A073346. Occurs for first time in A073202 as row 6 (and again as row 8).
Congruent to any of the sequences A000108, A007460, A007461, A007463, A007464, A061922, A068068 reduced modulo 2. Characteristic function of A000225.
If interpreted with offset=1 instead of 0 (i.e., a(1)=1, a(2)=1, a(3)=0, a(4)=1, ...) then this is the characteristic function of 2^n (A000079) and as such occurs as the first row of A073265. Also, in that case the INVERT transform will produce A023359.
This is Guy Steele's sequence GS(1, 3), also GS(3, 1) (see A135416).
Cf. A054525, A047999. - Gary W. Adamson, Oct 26 2009

Programs

  • Haskell
    a036987 n = ibp (n+1) where
       ibp 1 = 1
       ibp n = if r > 0 then 0 else ibp n' where (n',r) = divMod n 2
    a036987_list = 1 : f [0,1] where f (x:y:xs) = y : f (x:xs ++ [x,x+y])
    -- Same list generator function as for a091090_list, cf. A091090.
    -- Reinhard Zumkeller, May 19 2015, Apr 13 2013, Mar 13 2013
    
  • Maple
    A036987:= n-> `if`(2^ilog2(n+1) = n+1, 1, 0):
    seq(A036987(n), n=0..128);
  • Mathematica
    RealDigits[ N[ Sum[1/10^(2^n), {n, 0, Infinity}], 110]][[1]]
    (* Recurrence: *)
    t[n_, 1] = 1; t[1, k_] = 1;
    t[n_, k_] := t[n, k] =
      If[n < k, If[n > 1 && k > 1, -Sum[t[k - i, n], {i, 1, n - 1}], 0],
       If[n > 1 && k > 1, Sum[t[n - i, k], {i, 1, k - 1}], 0]];
    Table[t[n, k], {k, n, n}, {n, 104}]
    (* Mats Granvik, Jun 03 2011 *)
    mb2d[n_]:=1 - Module[{n2 = IntegerDigits[n, 2]}, Max[n2] - Min[n2]]; Array[mb2d, 120, 0] (* Vincenzo Librandi, Jul 19 2019 *)
    Table[PadRight[{1},2^k,0],{k,0,7}]//Flatten (* Harvey P. Dale, Apr 23 2022 *)
  • PARI
    {a(n) =( n++) == 2^valuation(n, 2)}; /* Michael Somos, Aug 25 2003 */
    
  • PARI
    a(n) = !bitand(n, n+1); \\ Ruud H.G. van Tol, Apr 05 2023
    
  • Python
    from sympy import catalan
    def a(n): return catalan(n)%2 # Indranil Ghosh, May 25 2017
    
  • Python
    def A036987(n): return int(not(n&(n+1))) # Chai Wah Wu, Jul 06 2022

Formula

1 followed by a string of 2^k - 1 0's. Also a(n)=1 iff n = 2^m - 1.
a(n) = a(floor(n/2)) * (n mod 2) for n>0 with a(0)=1. - Reinhard Zumkeller, Aug 02 2002 [Corrected by Mikhail Kurkov, Jul 16 2019]
Sum_{n>=0} 1/10^(2^n) = 0.110100010000000100000000000000010...
1 if n=0, floor(log_2(n+1)) - floor(log_2(n)) otherwise. G.f.: (1/x) * Sum_{k>=0} x^(2^k) = Sum_{k>=0} x^(2^k-1). - Ralf Stephan, Apr 28 2003
a(n) = 1 - A043545(n). - Michael Somos, Aug 25 2003
a(n) = -Sum_{d|n+1} mu(2*d). - Benoit Cloitre, Oct 24 2003
Dirichlet g.f. for right-shifted sequence: 2^(-s)/(1-2^(-s)).
a(n) = A000108(n) mod 2 = A001405(n) mod 2. - Paul Barry, Nov 22 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n,k)*Sum_{j=0..k} binomial(k, 2^j-1). - Paul Barry, Jun 01 2006
A000523(n+1) = Sum_{k=1..n} a(k). - Mitch Harris, Jul 22 2011
a(n) = A209229(n+1). - Reinhard Zumkeller, Mar 07 2012
a(n) = Sum_{k=1..n} A191898(n,k)*cos(Pi*(n-1)*(k-1))/n; (conjecture). - Mats Granvik, Mar 04 2013
a(n) = A000035(A000108(n)). - Omar E. Pol, Aug 06 2013
a(n) = 1 iff n=2^k-1 for some k, 0 otherwise. - M. F. Hasler, Jun 20 2014
a(n) = ceiling(log_2(n+2)) - ceiling(log_2(n+1)). - Gionata Neri, Sep 06 2015
From John M. Campbell, Jul 21 2016: (Start)
a(n) = (A000168(n-1) mod 2).
a(n) = (A000531(n+1) mod 2).
a(n) = (A000699(n+1) mod 2).
a(n) = (A000891(n) mod 2).
a(n) = (A000913(n-1) mod 2), for n>1.
a(n) = (A000917(n-1) mod 2), for n>0.
a(n) = (A001142(n) mod 2).
a(n) = (A001246(n) mod 2).
a(n) = (A001246(n) mod 4).
a(n) = (A002057(n-2) mod 2), for n>1.
a(n) = (A002430(n+1) mod 2). (End)
a(n) = 2 - A043529(n). - Antti Karttunen, Nov 19 2017
a(n) = floor(1+log(n+1)/log(2)) - floor(log(2n+1)/log(2)). - Adriano Caroli, Sep 22 2019
This is also the decimal expansion of -Sum_{k>=1} mu(2*k)/(10^k - 1), where mu is the Möbius function (A008683). - Amiram Eldar, Jul 12 2020

Extensions

Edited by M. F. Hasler, Jun 20 2014

A000698 A problem of configurations: a(0) = 1; for n>0, a(n) = (2n-1)!! - Sum_{k=1..n-1} (2k-1)!! a(n-k). Also the number of shellings of an n-cube, divided by 2^n n!.

Original entry on oeis.org

1, 1, 2, 10, 74, 706, 8162, 110410, 1708394, 29752066, 576037442, 12277827850, 285764591114, 7213364729026, 196316804255522, 5731249477826890, 178676789473121834, 5925085744543837186, 208256802758892355202, 7734158085942678174730
Offset: 0

Views

Author

Keywords

Comments

Also number of nonisomorphic unlabeled connected Feynman diagrams of order 2n-2 for the electron propagator of quantum electrodynamics (QED), including vanishing diagrams. [Corrected by Charles R Greathouse IV, Jan 24 2014][Clarified by Robert Coquereaux, Sep 14 2014]
a(n+1) is the moment of order 2*n for the probability density function rho(x) = (1/sqrt(2*Pi))*exp(x^2/2)/[(u(x))^2+Pi/2], with u(x) = Integral_{t=0..x} exp(t*t/2) dt, on the real interval -infinity..infinity. - Groux Roland, Jan 13 2009
Starting (1, 2, 10, 74, ...) = INVERTi transform of A001147: (1, 3, 15, 105, ...). - Gary W. Adamson, Oct 21 2009
The Cvitanovic et al. paper relates this sequence to A005411 and A005413. - Robert Munafo, Jan 24 2010
Hankel transform of a(n+1) is A168467. - Paul Barry, Nov 26 2009
a(n) = number of labeled Dyck (n-1)-paths (A000108) in which each vertex that terminates an upstep is labeled with an integer i in [0,h], where h is the height of the vertex . For example UDUD contributes 4 labeled paths--0D0D, 0D1D, 1D0D, 1D1D where upsteps are replaced by their labels--and UUDD contributes 6 labeled paths to a(3)=10. The Deléham (Mar 24 2007) formula below counts these labeled paths by number of "0" labels. - David Callan, Aug 23 2011
a(n) is the number of indecomposable perfect matchings on [2n]. A perfect matching on [2n] is decomposable if a nonempty subset of the edges forms a perfect matching on [2k] for some kDavid Callan, Nov 29 2012
From Robert Coquereaux, Sep 12 2014: (Start)
QED diagrams are graphs with two kinds of edges (lines): a (non-oriented), f (oriented), and only one kind of (internal) vertex: aff. They may have internal and external (i.e., pendant) lines. The order is the number of (internal) vertices. Vanishing diagrams: QED diagrams containing loops of type f with an odd number of vertices are set to 0 (Furry theorem). Proper diagrams: connected QED diagrams that remain connected when an arbitrary internal line is cut.
The number of Feynman diagrams of order 2n for the electron propagator (2-point function of QED), vanishing or not, proper or not, of order 2n, starting from n = 0, is given by 1, 2, 10, 74, 706, 8162, ..., i.e., this sequence A000698, with the first term (equal to 1) dropped. Call Sf the associated g.f.
The number of non-vanishing Feynman diagrams, for the same 2-point function, is given by 1, 1, 4, 25, 208, 2146, ..., i.e., by the sequence A005411, with a first term of order 0, equal to 1, added. Call S the associated g.f.
If one does not remove the vanishing diagram, but, at the same time, considers only those graphs that are proper, one obtains the Feynman diagrams (vanishing and non-vanishing) for the self-energy function of QED, 0, 1, 3, 21, 207, 2529, ..., i.e., the sequence A115974 with a first term of order 0, equal to 0, added. A115974 is twice A167872. Call Sigmaf the associated g.f.
If one removes the vanishing diagrams and, at the same time, considers only those graphs that are proper, one obtains the Feynman diagrams for the self-energy function of QED given by 0, 1, 3, 18, 153, 1638, ..., i.e., by the sequence A005412, with a first term of order 0, equal to 0, added. Call Sigma the associated g.f.
Then Sf = 1/(1-Sigmaf) and S = 1/(1-Sigma). (End)
For n>0 sum over all Dyck paths of semilength n-1 of products over all peaks p of (x_p+y_p)/y_p, where x_p and y_p are the coordinates of peak p. - Alois P. Heinz, May 22 2015
Also, counts certain isomorphism classes of closed normal linear lambda terms. [N. Zeilberger, 2015]. - N. J. A. Sloane, Sep 18 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
For n >= 2, a(n) is the number of coalescent histories for a pair consisting of a matching lodgepole gene tree and species tree with 2n-1 leaves. - Noah A Rosenberg, Jun 21 2022

Examples

			G.f. = 1 + x + 2*x^2 + 10*x^3 + 74*x^4 + 706*x^5 + 8162*x^6 + 110410*x^7 + ...
		

References

  • Dubois C., Giorgetti A., Genestier R. (2016) Tests and Proofs for Enumerative Combinatorics. In: Aichernig B., Furia C. (eds) Tests and Proofs. TAP 2016. Lecture Notes in Computer Science, vol 9762. Springer.
  • R. W. Robinson, Counting irreducible Feynman diagrams exactly and asymptotically, Abstracts Amer. Math. Soc., 2002, #975-05-270.
  • 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).

Crossrefs

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
Column k=1 of A258219, A258222.
Row sums of A322398.

Programs

  • Maple
    A006882 := proc(n) option remember; if n <= 1 then 1 else n*procname(n-2); fi; end;
    A000698:=proc(n) option remember; global df; local k; if n=0 then RETURN(1); fi; A006882(2*n-1) - add(A006882(2*k-1)*A000698(n-k),k=1..n-1); end;
    A000698 := proc(n::integer) local resul,fac,pows,c,c1,p,i ; if n = 0 then RETURN(1) ; else pows := combinat[partition](n) ; resul := 0 ; for p from 1 to nops(pows) do c := combinat[permute](op(p,pows)) ; c1 := op(1,c) ; fac := nops(c) ; for i from 1 to nops(c1) do fac := fac*doublefactorial(2*op(i,c1)-1) ; od ; resul := resul-(-1)^nops(c1)*fac ; od : fi ; RETURN(resul) ; end; # R. J. Mathar, Apr 24 2006
    # alternative Maple program:
    b:= proc(x, y, t) option remember; `if`(y>x or y<0, 0,
          `if`(x=0, 1, b(x-1, y-1, false)*`if`(t, (x+y)/y, 1) +
                       b(x-1, y+1, true)  ))
        end:
    a:= n-> `if`(n=0, 1, b(2*n-2, 0, false)):
    seq(a(n), n=0..25);  # Alois P. Heinz, May 23 2015
    a_list := proc(len) local n, A; if len=1 then return [1] fi: A := Array(-1..len-2); A[-1] := 1; A[0] := 1; for n to len-2 do A[n] := (2*n-1)*A[n-1]+add(A[j]*A[n-j-1], j=0..n-1) od: convert(A, list) end: a_list(20); # Peter Luschny, Jul 18 2017
  • Mathematica
    a[n_] := a[n] = (2n - 1)!! - Sum[ a[n - k](2k - 1)!!, {k, n-1}]; Array[a, 18, 0] (* Ignacio D. Peixoto, Jun 23 2006 *)
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ 2 - 1 / Sum[ (2 k - 1)!! x^k, {k, 0, n}], {x, 0, n}]]; (* Michael Somos, Nov 16 2011 *)
    a[n_]:= SeriesCoefficient[1+x(1/x+(E^((1/2)/x) Sqrt[2/\[Pi]] Sqrt[-(1/x)])/Erfc[Sqrt[-(1/x)]/Sqrt[2]]), {x,0,n}, Assumptions -> x >0](* Robert Coquereaux, Sep 14 2014 *)
    max = 20; g = t/Fold[1 - ((t + #2)*z)/#1 &, 1, Range[max, 1, -1]]; T[n_, k_] := SeriesCoefficient[g, {z, 0, n}, {t, 0, k}]; a[0] = 1; a[n_] := Sum[T[n-1, k], {k, 0, n}]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jan 31 2016, after Philippe Deléham *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 2 - 1 / sum( k=0, n, x^k * (2*k)! /(2^k * k!), x * O(x^n)), n))}; /* Michael Somos, Feb 08 2011 */
    
  • PARI
    {a(n) = my(A); if( n<1, n==0, A = vector(n); A[1] = 1; for( k=2, n, A[k] = (2*k - 3) * A[k-1] + sum( j=1, k-1, A[j] * A[k-j])); A[n])}; /* Michael Somos, Jul 24 2011 */
    
  • Python
    from sympy import factorial2, cacheit
    @cacheit
    def a(n): return 1 if n == 0 else factorial2(2*n - 1) - sum(factorial2(2*k - 1)*a(n - k) for k in range(1, n))
    [a(n) for n in range(51)]  # Indranil Ghosh, Jul 18 2017

Formula

G.f.: 2 - 1/(1 + Sum_{n>=1} (2*n-1)!! * x^n ).
a(n+1) = Sum_{k=0..n} A089949(n, k)*2^k. - Philippe Deléham, Aug 15 2005
a(n+1) = Sum_{k=0..n} A053979(n,k). - Philippe Deléham, Mar 24 2007
From Paul Barry, Nov 26 2009: (Start)
G.f.: 1+x/(1-2x/(1-3x/(1-4x/(1-5x/(1-6x/(1-... (continued fraction).
G.f.: 1+x/(1-2x-6x^2/(1-7x-20x^2/(1-11x-42x^2/(1-15x-72x^2/(1-19x-110x^2/(1-... (continued fraction). (End)
G.f.: 1 + x * B(x) * C(x) where B(x) is the g.f. for A001147 and C(x) is the g.f. for A005416. - Michael Somos, Feb 08 2011
G.f.: 1+x/W(0); where W(k)=1+x+x*2k-x*(2k+3)/W(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 17 2011
From Peter Bala, Dec 22 2011: (Start)
Recurrence relation: a(n+1) = (2*n-1)*a(n) + Sum_{k = 1..n} a(k)*a(n+1-k) for n >= 0 and a(1) = 1.
The o.g.f. B(x) = Sum_{n>=1} a(n)*x^(2*n-1) = x + 2*x^3 + 10*x^5 + 74*x^7 + ... satisfies the Riccati differential equation y'(x) = -1/x^2 + (1/x^3)*y(x) - (1/x^2)*y(x)^2 with initial condition y(0) = 0 (cf. A005412). The solution is B(x) = 1/z(x) + 1/x, where z(x) = -Sum_{n>=0} A001147(n) * x^(2*n+1) = -(x + x^3 + 3*x^5 + 15*x^7 + ...). The function b(x) = -B(1/x) satisfies b'(x) = -1 - (x + b(x))*b(x). Hence the differential operator (D^2 + x*D + 1), where D = d/dx, factorizes as (D - a(x))*(D - b(x)), where a(x) = -(x + b(x)), as conjectured by [Edgar, Problem 4.32]. For a refinement of this sequence see A053979. (End)
From Sergei N. Gladkovskii, Aug 19 2012, Oct 24 2012, Mar 19 2013, May 20 2013, May 29 2013, Aug 04 2013, Aug 05 2013: (Start)
Continued fractions:
G.f.: 2 - G(0) where G(k) = 1 - (k+1)*x/G(k+1).
G.f.: 2 - U(0) where U(k) = 1 - (2*k+1)*x/(1 - (2*k+2)*x/U(k+1)).
G.f.: 2 - U(0) where U(k) = 1 - (4*k+1)*x - (2*k+1)*(2*k+2)*x^2/U(k+1).
G.f.: 1/Q(0) where Q(k) = 1 - x*(2*k+2)/(1 - x*(2*k+3)/Q(k+1)).
G.f.: 1 + x/Q(0) where Q(k) = 1 - x*(k+2)/Q(k+1).
G.f.: 2 - G(0)/2 where G(k) = 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) - 1 + 2*x*(2*k+2)/ G(k+1))).
G.f.: 1 + x*G(0) where G(k) = 1 - x*(k+2)/(x*(k+2) - 1/G(k+1)).
G.f.: 2 - 1/B(x) where B(x) is the g.f. of A001147.
G.f.: 1 + x/(1-2*x*B(x)) where B(x) is the g.f. of A167872. (End)
a(n) ~ 2^(n+1/2) * n^n / exp(n). - Vaclav Kotesovec, Mar 10 2014
G.f.: 1 + x*(1/x + (sqrt(2/Pi) * exp(1/(2*x)) * sqrt(-1/x))/Erfc(sqrt(-1/x)/sqrt(2))) where Erfc(z) = 1 - Erf(z) is the complementary error function, and Erf(z) is the integral of the Gaussian distribution. This generating function is obtained from the generating functional of (4-dimensional) QED, evaluated in dimension 0 for the 2-point function, without the modification implementing Furry theorem. - Robert Coquereaux, Sep 14 2014
From Peter Bala, May 23 2017: (Start)
G.f. A(x) = 1 + x/(1 + x - 3*x/(1 + 3*x - 5*x/(1 + 5*x - 7*x/(1 + 7*x - ...)))).
A(x) = 1 + x/(1 + x - 3*x/(1 - 2*x/(1 - 5*x/(1 - 4*x/(1 - 7*x/(1 - 6*x/(1 - ...))))))). (End)

Extensions

Formula corrected by Ignacio D. Peixoto, Jun 23 2006
More terms from Sean A. Irvine, Feb 27 2011

A000260 Number of rooted simplicial 3-polytopes with n+3 nodes; or rooted 3-connected triangulations with 2n+2 faces; or rooted 3-connected trivalent maps with 2n+2 vertices.

Original entry on oeis.org

1, 1, 3, 13, 68, 399, 2530, 16965, 118668, 857956, 6369883, 48336171, 373537388, 2931682810, 23317105140, 187606350645, 1524813969276, 12504654858828, 103367824774012, 860593023907540, 7211115497448720, 60776550501588855
Offset: 0

Views

Author

Keywords

Comments

Number of rooted loopless planar maps with n edges. E.g., there are a(2)=3 loopless planar maps with 2 edges: two rooted paths (.-.-.) and one digon (.=.). - Valery A. Liskovets, Sep 25 2003
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Tamari lattice (rotation lattice of binary trees) of size n (see Pallo and Chapoton references). - Ralf Stephan, May 08 2007, Jean Pallo (Jean.Pallo(AT)u-bourgogne.fr), Sep 11 2007
Number of rooted triangulations of type [n, 0] (see Brown paper eq (4.8)). - Michel Marcus, Jun 23 2013
Equivalently, number of rooted bridgeless planar maps with n edges. - Noam Zeilberger, Oct 06 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
Number of uniquely sorted permutations of [2n+1] that avoid the pattern 231. Also the number of uniquely sorted permutations of [2n+1] that avoid 132. - Colin Defant, Jun 13 2019
The sequence 1,3,13,68,... appears naturally in integral geometry, namely in the algebra of unitarily invariant valuations on complex space forms. - Andreas Bernig, Feb 02 2020

Examples

			G.f. = 1 + x + 3*x^2 + 13*x^3 + 68*x^4 + 399*x^5 + 2530*x^6 + 16965*x^7 + ...
		

References

  • C. F. Earl and L. J. March, Architectural applications of graph theory, pp. 327-355 of R. J. Wilson and L. W. Beineke, editors, Applications of Graph Theory. Academic Press, NY, 1979.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 714.
  • Handbook of Combinatorics, North-Holland '95, p. 891.
  • 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).
  • W. T. Tutte, The enumerative theory of planar maps, in A Survey of Combinatorial Theory (J. N. Srivastava et al. eds.), pp. 437-448, North-Holland, Amsterdam, 1973.

Crossrefs

Row sums of A342981.
Column 0 of A146305 and A341856; Column 2 of A255918.
Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.

Programs

  • Magma
    [Binomial(4*n+1, n+1)-9*Binomial(4*n+1, n-1): n in [0..25]]; // Vincenzo Librandi, Nov 24 2016
  • Maple
    A000260 := proc(n)
        2*(4*n+1)!/((n+1)!*(3*n+2)!) ;
    end proc:
  • Mathematica
    Table[Binomial[4n+1,n+1]-9*Binomial[4n+1,n-1],{n,0,25}] (* Harvey P. Dale, Aug 23 2011 *)
    a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1/2, 3/4, 1, 5/4}, {4/3, 5/3, 2}, 256/27 x], {x, 0, n}]; (* Michael Somos, Dec 23 2014 *)
    terms = 22; G[] = 0; Do[G[x] = 1+x*G[x]^4 + O[x]^terms, terms];
    CoefficientList[(2-G[x])*G[x]^2, x] (* Jean-François Alcover, Jan 13 2018, after Mark van Hoeij *)
  • PARI
    {a(n) = if( n<0, 0, 2 * (4*n + 1)! / ((n + 1)! * (3*n + 2)!))}; /* Michael Somos, Sep 07 2005 */
    
  • PARI
    {a(n) = binomial( 4*n + 2, n + 1) / ((2*n + 1) * (3*n + 2))}; /* Michael Somos, Mar 28 2012 */
    
  • Sage
    def a(n):
        n = ZZ(n)
        return (4*n + 2).binomial(n + 1) // ((2*n + 1) * (3*n + 2))
    # F. Chapoton, Aug 06 2015
    

Formula

a(n) = 2*(4*n+1)! / ((n+1)!*(3*n+2)!) = binomial(4*n+1, n+1) - 9*binomial(4*n+1, n-1).
G.f.: (2-g)*g^2 where g = 1 + x*g^4 is the g.f. of A002293. - Mark van Hoeij, Nov 10 2011
G.f.: hypergeom([1,1/2,3/4,5/4],[2,4/3,5/3],256*x/27) = 1 + 120*x/(Q(0)-120*x); Q(k) = 8*x*(2*k+1)*(4*k+3)*(4*k+5) + 3*(k+2)*(3*k+4)*(3*k+5) - 24*x*(k+2)*(2*k+3)*(3*k+4)*(3*k+5)*(4*k+7)*(4*k+9)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 25 2011
a(n) = binomial(4*n + 2, n + 1) / ((2*n + 1) * (3*n + 2)). - Michael Somos, Mar 28 2012
a(n) * (n+1) = A069271(n). - Michael Somos, Mar 28 2012
0 = F(a(n), a(n+1), ..., a(n+8)) for all n in Z where a(-1) = 3/4 and F() is a polynomial of degree 2 with integer coefficients and 29 monomials. - Michael Somos, Dec 23 2014
D-finite with recurrence: 3*(3*n+2)*(3*n+1)*(n+1)*a(n) - 8*(4*n+1)*(2*n-1)*(4*n-1)*a(n-1) = 0. - R. J. Mathar, Oct 21 2015
a(n) = Sum_{k=1..A000108(n)} k * A263191(n,k). - Alois P. Heinz, Nov 16 2015
a(n) ~ 2^(8*n+7/2) / (sqrt(Pi) * n^(5/2) * 3^(3*n+5/2)). - Vaclav Kotesovec, Feb 26 2016
E.g.f.: 3F3(1/2,3/4,5/4; 4/3,5/3,2; 256*x/27). - Ilya Gutkovskiy, Feb 01 2017
From Gheorghe Coserea, Aug 17 2017: (Start)
G.f. y(x) satisfies:
0 = x^3*y^4 + 3*x^2*y^3 + x*(8*x+3)*y^2 - (20*x-1)*y + 16*x-1.
0 = x*(256*x - 27)*deriv(y,x) - 8*x^2*y^3 - 25*x*y^2 + 4*(24*x-11)*y + 44.
(End)
From Karol A. Penson, Apr 06 2022: (Start)
a(n) = Integral_{x=0...256/27} x^n*W(x), where
W(x) = (sqrt(2)/Pi)*(h1(x) - h2(x) + h3(x)) and
h1(x) = 3F2([-6/12,-2/12, 2/12], [ 3/12, 9/12], (27*x)/256)/((x/2)^(1/2));
h2(x) = 3F2([-3/12, 1/12, 5/12], [ 6/12, 15/12], (27*x)/256)/(x^(1/4));
h3(x) = 3F2([ 3/12, 7/12, 11/12], [18/12, 21/12], (27*x)/256)/(x^(-1/4)*32).
This integral representation is unique as the solution of n-th Hausdorff power moment of the function W. Using only the definition of a(n), W(x) can be proven to be positive. W(x) is singular at x = 0 and for x > 0 is monotonically decreasing to zero at x = 256/27. (End)
a(n) = (1/27^n) * Product_{1 <= i <= j <= 3*n} (3*i + j + 3)/(3*i + j - 1). Cf. A006013. - Peter Bala, Feb 21 2023

Extensions

Edited by F. Chapoton, Feb 03 2011
Showing 1-10 of 79 results. Next