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 23 results. Next

A003114 Number of partitions of n into parts 5k+1 or 5k+4.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 14, 17, 19, 23, 26, 31, 35, 41, 46, 54, 61, 70, 79, 91, 102, 117, 131, 149, 167, 189, 211, 239, 266, 299, 333, 374, 415, 465, 515, 575, 637, 709, 783, 871, 961, 1065, 1174, 1299, 1429, 1579, 1735, 1913, 2100, 2311, 2533, 2785
Offset: 0

Views

Author

Keywords

Comments

Expansion of Rogers-Ramanujan function G(x) in powers of x.
Same as number of partitions into distinct parts where the difference between successive parts is >= 2.
As a formal power series, the limit of polynomials S(n,x): S(n,x)=sum(T(i,x),0<=i<=n); T(i,x)=S(i-2,x).x^i; T(0,x)=1,T(1,x)=x; S(n,1)=A000045(n+1), the Fibonacci sequence. - Claude Lenormand (claude.lenormand(AT)free.fr), Feb 04 2001
The Rogers-Ramanujan identity is 1 + Sum_{n >= 1} t^(n^2)/((1-t)*(1-t^2)*...*(1-t^n)) = Product_{n >= 1} 1/((1-t^(5*n-1))*(1-t^(5*n-4))).
Coefficients in expansion of permanent of infinite tridiagonal matrix:
1 1 0 0 0 0 0 0 ...
x 1 1 0 0 0 0 0 ...
0 x^2 1 1 0 0 0 ...
0 0 x^3 1 1 0 0 ...
0 0 0 x^4 1 1 0 ...
................... - Vladeta Jovovic, Jul 17 2004
Also number of partitions of n such that the smallest part is greater than or equal to number of parts. - Vladeta Jovovic, Jul 17 2004
Also number of partitions of n such that if k is the largest part, then each of {1, 2, ..., k-1} occur at least twice. Example: a(9)=5 because we have [3, 2, 2, 1, 1], [2, 2, 2, 1, 1, 1], [2, 2, 1, 1, 1, 1, 1], [2, 1, 1, 1, 1, 1, 1, 1] and [1, 1, 1, 1, 1, 1, 1, 1, 1]. - Emeric Deutsch, Feb 27 2006
Also number of partitions of n such that if k is the largest part, then k occurs at least k times. Example: a(9)=5 because we have [3, 3, 3], [2, 2, 2, 2, 1], [2, 2, 2, 1, 1, 1], [2, 2, 1, 1, 1, 1, 1] and [1, 1, 1, 1, 1, 1, 1, 1, 1]. - Emeric Deutsch, Apr 16 2006
a(n) = number of NW partitions of n, for n >= 1; see A237981.
For more about the generalized Rogers-Ramanujan series G[i](x) see the Andrews-Baxter and Lepowsky-Zhu papers. The present series is G[1](x). - N. J. A. Sloane, Nov 22 2015
Convolution of A109700 and A109697. - Vaclav Kotesovec, Jan 21 2017

Examples

			G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 3*x^6 + 3*x^7 + 4*x^8 + 5*x^9 + ...
G.f. = 1/q + q^59 + q^119 + q^179 + 2*q^239 + 2*q^299 + 3*q^359 + 3*q^419 + ...
From _Joerg Arndt_, Dec 27 2012: (Start)
The a(16)=17 partitions of 16 where all parts are 1 or 4 (mod 5) are
  [ 1]  [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]
  [ 2]  [ 4 1 1 1 1 1 1 1 1 1 1 1 1 ]
  [ 3]  [ 4 4 1 1 1 1 1 1 1 1 ]
  [ 4]  [ 4 4 4 1 1 1 1 ]
  [ 5]  [ 4 4 4 4 ]
  [ 6]  [ 6 1 1 1 1 1 1 1 1 1 1 ]
  [ 7]  [ 6 4 1 1 1 1 1 1 ]
  [ 8]  [ 6 4 4 1 1 ]
  [ 9]  [ 6 6 1 1 1 1 ]
  [10]  [ 6 6 4 ]
  [11]  [ 9 1 1 1 1 1 1 1 ]
  [12]  [ 9 4 1 1 1 ]
  [13]  [ 9 6 1 ]
  [14]  [ 11 1 1 1 1 1 ]
  [15]  [ 11 4 1 ]
  [16]  [ 14 1 1 ]
  [17]  [ 16 ]
The a(16)=17 partitions of 16 where successive parts differ by at least 2 are
  [ 1]  [ 7 5 3 1 ]
  [ 2]  [ 8 5 3 ]
  [ 3]  [ 8 6 2 ]
  [ 4]  [ 9 5 2 ]
  [ 5]  [ 9 6 1 ]
  [ 6]  [ 9 7 ]
  [ 7]  [ 10 4 2 ]
  [ 8]  [ 10 5 1 ]
  [ 9]  [ 10 6 ]
  [10]  [ 11 4 1 ]
  [11]  [ 11 5 ]
  [12]  [ 12 3 1 ]
  [13]  [ 12 4 ]
  [14]  [ 13 3 ]
  [15]  [ 14 2 ]
  [16]  [ 15 1 ]
  [17]  [ 16 ]
(End)
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 109, 238.
  • G. E. Andrews, R. Askey and R. Roy, Special Functions, Cambridge University Press, 1999; Exercise 6(e), p. 591.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 669.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 107.
  • G. H. Hardy, Ramanujan, AMS Chelsea Publ., Providence, RI, 2002, pp. 90-92.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth ed., Clarendon Press, Oxford, 2003, pp. 290-291.
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A188216 (least part k occurs at least k times).
For the generalized Rogers-Ramanujan series G[1], G[2], G[3], G[4], G[5], G[6], G[7], G[8] see A003114, A003106, A006141, A264591, A264592, A264593, A264594, A264595. G[0] = G[1]+G[2] is given by A003113.
Row sums of A268187.

Programs

  • Haskell
    a003114 = p a047209_list where
       p _      0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Jan 05 2011
    
  • Haskell
    a003114 = p 1 where
       p _ 0 = 1
       p k m = if k > m then 0 else p (k + 2) (m - k) + p (k + 1) m
    -- Reinhard Zumkeller, Feb 19 2013
  • Maple
    g:=sum(x^(k^2)/product(1-x^j,j=1..k),k=0..10): gser:=series(g,x=0,65): seq(coeff(gser,x,n),n=0..60); # Emeric Deutsch, Feb 27 2006
  • Mathematica
    CoefficientList[ Series[Sum[x^k^2/Product[1 - x^j, {j, 1, k}], {k, 0, 10}], {x, 0, 65}], x][[1 ;; 61]] (* Jean-François Alcover, Apr 08 2011, after Emeric Deutsch *)
    Table[Count[IntegerPartitions[n], p_ /; Min[p] >= Length[p]], {n, 0, 24}] (* Clark Kimberling, Feb 13 2014 *)
    a[ n_] := SeriesCoefficient[ 1 / (QPochhammer[ x^1, x^5] QPochhammer[ x^4, x^5]), {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ Product[ (1 - x^k)^{-1, 0, 0, -1, 0}[[Mod[k, 5, 1]]], {k, n}], {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    nmax = 60; kmax = nmax/5;
    s = Flatten[{Range[0, kmax]*5 + 1}~Join~{Range[0, kmax]*5 + 4}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Aug 02 2020 *)
  • PARI
    {a(n) = my(t); if( n<0, 0, t = 1 + x * O(x^n); polcoeff( sum(k=1, sqrtint(n), t *= x^(2*k - 1) / (1 - x^k) * (1 + x * O(x^(n - k^2))), 1), n))}; /* Michael Somos, Oct 15 2008 */
    

Formula

G.f.: Sum_{k>=0} x^(k^2)/(Product_{i=1..k} 1-x^i).
The g.f. above is the special case D=2 of sum(n>=0, x^(D*n*(n+1)/2 - (D-1)*n) / prod(k=1..n, 1-x^k) ), the g.f. for partitions into distinct part where the difference between successive parts is >= D. - Joerg Arndt, Mar 31 2014
G.f.: 1 + sum(i=1, oo, x^(5i+1)/prod(j=1 or 4 mod 5 and j<=5i+1, 1-x^j) + x^(5i+4)/prod(j=1 or 4 mod 5 and j<=5i+4, 1-x^j)). - Jon Perry, Jul 06 2004
G.f.: (Product_{k>0} 1 + x^(2*k)) * (Sum_{k>=0} x^(k^2) / (Product_{i=1..k} 1 - x^(4*i))). - Michael Somos, Oct 19 2006
Euler transform of period 5 sequence [ 1, 0, 0, 1, 0, ...]. - Michael Somos, Oct 15 2008
Expansion of f(-x^5) / f(-x^1, -x^4) in powers of x where f(,) is the Ramanujan general theta function. - Michael Somos, May 17 2015
Expansion of f(-x^2, -x^3) / f(-x) in powers of x where f(,) is the Ramanujan general theta function. - Michael Somos, Jun 13 2015
a(n) ~ phi^(1/2) * exp(2*Pi*sqrt(n/15)) / (2 * 3^(1/4) * 5^(1/2) * n^(3/4)) * (1 - (3*sqrt(15)/(16*Pi) + Pi/(60*sqrt(15))) / sqrt(n)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Aug 23 2015, extended Jan 24 2017
a(n) = (1/n)*Sum_{k=1..n} A284150(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 21 2017

A003319 Number of connected permutations of [1..n] (those not fixing [1..j] for 0 < j < n). Also called indecomposable permutations, or irreducible permutations.

Original entry on oeis.org

1, 1, 1, 3, 13, 71, 461, 3447, 29093, 273343, 2829325, 31998903, 392743957, 5201061455, 73943424413, 1123596277863, 18176728317413, 311951144828863, 5661698774848621, 108355864447215063, 2181096921557783605, 46066653228356851631, 1018705098450570562877
Offset: 0

Views

Author

Keywords

Comments

Also the number of permutations with no global descents, as introduced by Aguiar and Sottile [Corollaries 6.3, 6.4 and Remark 6.5].
Also the dimensions of the homogeneous components of the space of primitive elements of the Malvenuto-Reutenauer Hopf algebra of permutations. This result, due to Poirier and Reutenauer [Theoreme 2.1] is stated in this form in the work of Aguiar and Sottile [Corollary 6.3] and also in the work of Duchamp, Hivert and Thibon [Section 3.3].
Related to number of subgroups of index n-1 in free group of rank 2 (i.e., maximal number of subgroups of index n-1 in any 2-generator group). See Problem 5.13(b) in Stanley's Enumerative Combinatorics, Vol. 2.
Also the left border of triangle A144107, with row sums = n!. - Gary W. Adamson, Sep 11 2008
Hankel transform is A059332. Hankel transform of aerated sequence is A137704(n+1). - Paul Barry, Oct 07 2008
For every n, a(n+1) is also the moment of order n for the probability density function rho(x) = exp(x)/(Ei(1,-x)*(Ei(1,-x) + 2*I*Pi)) on the interval 0..infinity, with Ei the exponential-integral function. - Groux Roland, Jan 16 2009
Also (apparently), a(n+1) is the number of rooted hypermaps with n darts on a surface of any genus (see Walsh 2012). - N. J. A. Sloane, Aug 01 2012
Also recurrent sequence A233824 (for n > 0) in Panaitopol's formula for pi(x), the number of primes <= x. - Jonathan Sondow, Dec 19 2013
Also the number of mobiles (cyclic rooted trees) with an arrow from each internal vertex to a descendant of that vertex. - Brad R. Jones, Sep 12 2014
Up to sign, Möbius numbers of the shard intersection orders of type A, see Theorem 1.3 in Reading reference. - F. Chapoton, Apr 29 2015
Also, a(n) is the number of distinct leaf matrices of complete non-ambiguous trees of size n. - Daniel Chen, Oct 23 2022

Examples

			G.f. = 1 + x + x^2 + 3*x^3 + 13*x^4 + 71*x^5 + 461*x^6 + 3447*x^7 + 29093*x^8 + ...
From _Peter Luschny_, Aug 03 2022: (Start)
A permutation p in [n] (where n >= 0) is reducible if there exists an i in 1..n-1 such that for all j in the range 1..i and all k in the range i+1..n it is true that p(j) < p(k). (Note that a range a..b includes a and b.) If such an i exists we say that i splits the permutation at i.
Examples:
* () is not reducible since there is no index i which splits (). (=> a(0) = 1)
* (1) is not reducible since there is no index i which splits (1). (=> a(1) = 1)
* (1, 2) is reducible since index 1 splits (1, 2) as p(1) < p(2).
* (2, 1) is not reducible since at the only potential splitting point i = 1 we have p(1) > p(2). (=> a(2) = 1)
* For n = 3 we have (1, 2, 3), (1, 3, 2), and (2, 1, 3) are reducible and (2, 3, 1), (3, 1, 2), and (3, 2, 1) are irreducible. (End)
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 25, Example 20.
  • E. W. Bowen, Letter to N. J. A. Sloane, Aug 27 1976.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 84 (#25), 262 (#14) and 295 (#16).
  • P. de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 23, N_{n,2}.
  • I. M. Gessel and R. P. Stanley, Algebraic Enumeration, chapter 21 in Handbook of Combinatorics, Vol. 2, edited by R. L. Graham et al., The MIT Press, Mass, 1995.
  • M. Kauers and P. Paule, The Concrete Tetrahedron, Springer 2011, p. 22.
  • H. P. Robinson, Letter to N. J. A. Sloane, Nov 19 1973.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, Chap. 1, Ex. 128; Vol. 2, 1999, see Problem 5.13(b).

Crossrefs

See A167894 for another version.
Bisections give A272656, A272657.
Row sums of A111184 and A089949.
Leading diagonal of A059438. A diagonal of A263484.
Cf. A090238, A000698, A356291 (reducible permutations).
Column k=0 of A370380 and A370381 (without pair of initial terms and with different offset).

Programs

  • Maple
    INVERTi([seq(n!,n=1..20)]);
    A003319 := proc(n) option remember; n! - add((n-j)!*A003319(j), j=1..n-1) end;
    [seq(A003319(n), n=0..50)]; # N. J. A. Sloane, Dec 28 2011
    series(2 - 1/hypergeom([1,1], [], x), x=0,50); # Mark van Hoeij, Apr 18 2013
  • Mathematica
    a[n_] := a[n] = n! - Sum[k!*a[n-k], {k, 1, n-1}]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Oct 11 2011, after given formula *)
    CoefficientList[Assuming[Element[x,Reals],Series[2-E^(1/x)* x/ExpIntegralEi[1/x],{x,0,20}]],x] (* Vaclav Kotesovec, Mar 07 2014 *)
    a[ n_] := If[ n < 2, 1, a[n] = (n - 2) a[n - 1] + Sum[ a[k] a[n - k], {k, n - 1}]]; (* Michael Somos, Feb 23 2015 *)
    Table[SeriesCoefficient[1 + x/(1 + ContinuedFractionK[-Floor[(k + 2)/2]*x, 1, {k, 1, n}]), {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Sep 29 2017 *)
  • PARI
    {a(n) = my(A); if( n<1, 1, A = vector(n); A[1] = 1; for( k=2, n, A[k] = (k - 2) * A[k-1] + sum( j=1, k-1, A[j] * A[k-j])); A[n])}; /* Michael Somos, Jul 24 2011 */
    
  • PARI
    {if(n<1,1,a(n)=local(A=x);for(i=1,n,A=x-x*A+A^2+x^2*A' +x*O(x^n));polcoeff(A,n))} /* Paul D. Hanna, Jul 30 2011 */
    
  • Sage
    def A003319_list(len):
        R, C = [1], [1] + [0] * (len - 1)
        for n in range(1, len):
            for k in range(n, 0, -1):
                C[k] = C[k - 1] * k
            C[0] = -sum(C[k] for k in range(1, n + 1))
            R.append(-C[0])
        return R
    print(A003319_list(21))  # Peter Luschny, Feb 19 2016

Formula

G.f.: 2 - 1/Sum_{k>=0} k!*x^k.
Also a(n) = n! - Sum_{k=1..n-1} k!*a(n-k) [Bowen, 1976].
Also coefficients in the divergent series expansion log Sum_{n>=0} n!*x^n = Sum_{n>=1} a(n+1)*x^n/n [Bowen, 1976].
a(n) = (-1)^(n-1) * det {| 1! 2! ... n! | 1 1! ... (n-1)! | 0 1 1! ... (n-2)! | ... | 0 ... 0 1 1! |}.
INVERTi transform of factorial numbers, A000142 starting from n=1. - Antti Karttunen, May 30 2003
Gives the row sums of the triangle [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, ...] where DELTA is the operator defined in A084938; this triangle A089949. - Philippe Deléham, Dec 30 2003
a(n+1) = Sum_{k=0..n} A089949(n,k). - Philippe Deléham, Oct 16 2006
L.g.f.: Sum_{n>=1} a(n)*x^n/n = log( Sum_{n>=0} n!*x^n ). - Paul D. Hanna, Sep 19 2007
G.f.: 1+x/(1-x/(1-2*x/(1-2*x/(1-3*x/(1-3*x/(1-4*x/(1-4*x/(1-...)))))))) (continued fraction). - Paul Barry, Oct 07 2008
a(n) = -Sum_{i=0..n} (-1)^i*A090238(n, i) for n > 0. - Peter Luschny, Mar 13 2009
From Gary W. Adamson, Jul 14 2011: (Start)
a(n) = upper left term in M^(n-1), M = triangle A128175 as an infinite square production matrix (deleting the first "1"); as follows:
1, 1, 0, 0, 0, 0, ...
2, 2, 1, 0, 0, 0, ...
4, 4, 3, 1, 0, 0, ...
8, 8, 7, 4, 1, 0, ...
16, 16, 15, 11, 5, 1, ...
... (End)
O.g.f. satisfies: A(x) = x - x*A(x) + A(x)^2 + x^2*A'(x). - Paul D. Hanna, Jul 30 2011
From Sergei N. Gladkovskii, Jun 24 2012: (Start)
Let A(x) be the g.f.; then
A(x) = 1/Q(0), where Q(k) = x + 1 + x*k - (k+2)*x/Q(k+1).
A(x) = (1-1/U(0))/x, when U(k) = 1 + x*(2*k+1)/(1 - 2*x*(k+1)/(2*x*(k+1) + 1/U(k+1))). (End)
From Sergei N. Gladkovskii, Aug 03 2013: (Start)
Continued fractions:
G.f.: 1 - G(0)/2, where G(k) = 1 + 1/(1 - x*(2*k+2)/(x*(2*k+2) - 1 + x*(2*k+2)/G(k+1))).
G.f.: (x/2)*G(0), where G(k) = 1 + 1/(1 - x*(k+1)/(x*(k+1/2) + 1/G(k+1))).
G.f.: x*G(0), where G(k) = 1 - x*(k+1)/(x - 1/G(k+1)).
G.f.: 1 - 1/G(0), where G(k) = 1 - x*(k+1)/(x*(k+1) - 1/(1 - x*(k+1)/(x*(k+1) - 1/G(k+1)))).
G.f.: x*W(0), where W(k) = 1 - x*(k+1)/(x*(k+1) - 1/(1 - x*(k+2)/(x*(k+2) - 1/W(k+1)))).
(End)
a(n) = A233824(n-1) if n > 0. (Proof. Set b(n) = A233824(n), so that b(n) = n*n! - Sum_{k=1..n-1} k!*b(n-k). To get a(n+1) = b(n) for n >= 0, induct on n, use (n+1)! = n*n! + n!, and replace k with k+1 in the sum.) - Jonathan Sondow, Dec 19 2013
a(n) ~ n! * (1 - 2/n - 1/n^2 - 5/n^3 - 32/n^4 - 253/n^5 - 2381/n^6 - 25912/n^7 - 319339/n^8 - 4388949/n^9 - 66495386/n^10), for coefficients see A260503. - Vaclav Kotesovec, Jul 27 2015
For n>0, a(n) = (A059439(n) - A259472(n))/2. - Vaclav Kotesovec, Aug 03 2015
From Peter Bala, May 23 2017: (Start)
G.f.: 1 + x/(1 + x - 2*x/(1 + 2*x - 3*x/(1 + 3*x - 4*x/(1 + 4*x - ...)))). Cf. A000698.
G.f.: 1/(1 - x/(1 + x - x/(1 - 2*x/(1 - 2*x/(1 - 3*x/(1 - 3*x/(1 - 4*x/(1 - 4*x/(1 - ...))))))))). (End)
Conjecture: a(n) = A370380(n-2, 0) = A370381(n-2, 0) for n > 1 with a(0) = a(1) = 1. - Mikhail Kurkov, Apr 26 2024

Extensions

More terms from Michael Somos, Jan 26 2000
Additional comments from Marcelo Aguiar (maguiar(AT)math.tamu.edu), Mar 28 2002
Added a(0)=0 (some of the formulas may now need adjusting). - N. J. A. Sloane, Sep 12 2012
Edited and set a(0) = 1 by Peter Luschny, Aug 03 2022

A003105 Schur's 1926 partition theorem: number of partitions of n into parts 6n+1 or 6n-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 3, 3, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20, 23, 26, 30, 34, 38, 42, 47, 53, 60, 67, 74, 82, 91, 102, 114, 126, 139, 153, 169, 187, 207, 228, 250, 274, 301, 331, 364, 399, 436, 476, 520, 569, 622, 679, 739, 804, 875, 953, 1038, 1128, 1224, 1327
Offset: 0

Views

Author

Keywords

Comments

There are many (at least 8) equivalent definitions of this sequence (besides the comments below, see also Schur, Alladi, Andrews). - N. J. A. Sloane, Jun 17 2011
Coefficients of replicable function number 72e. - N. J. A. Sloane, Jun 10 2015
Also number of partitions of n into odd parts in which no part appears more than twice, cf. A070048 and A096938. - Vladeta Jovovic, Jan 18 2005
Also number of partitions of n into distinct parts congruent to 1 or 2 modulo 3. (Follows from second g.f.) - N. Sato, Jul 20 2005
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A262928 and A261612. - Vaclav Kotesovec, Jan 13 2017
Convolution of A109702 and A109701. - Vaclav Kotesovec, Jan 21 2017

Examples

			G.f: A(x) = 1 + x + x^2 + x^3 + x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 3*x^8 + 3*x^9 + 4*x^10 + ...
T72e = 1/q + q^11 + q^23 + q^35 + q^47 + 2*q^59 + 2*q^71 + 3*q^83 + ...
The logarithm of the g.f. begins:
log(A(x)) = x + x^2/2 + x^3/3 + x^4/4 + 6*x^5/5 + x^6/6 + 8*x^7/7 + x^8/8 + x^9/9 + 6*x^10/10 + 12*x^11/11 + x^12/12 + ... + A186099(n)*x^n/n + ... . - _Paul D. Hanna_, Feb 17 2013
		

References

  • K. Alladi, Refinements of Rogers-Ramanujan type identities. In Special Functions, q-Series and Related Topics (Toronto, ON, 1995), 1-35, Fields Inst. Commun., 14, Amer. Math. Soc., Providence, RI, 1997.
  • G. E. Andrews, Schur's theorem, partitions with odd parts and the Al-Salam-Carlitz polynomials. In q-Series From a Contemporary Perspective (South Hadley, MA, 1998), 45-56, Contemp. Math., 254, Amer. Math. Soc., Providence, RI, 2000.
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • I. Schur, Zur Additiven Zahlentheorie, Ges. Abh., Vol. 2, Springer, pp. 43-50.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a003105 n = p 1 n where
       p k m | m == 0 = 1 | m < k = 0 | otherwise = q k (m-k) + p (k+2) m
       q k m | m == 0 = 1 | m < k = 0 | otherwise = p (k+2) (m-k) + p (k+2) m
    -- Reinhard Zumkeller, Nov 12 2011
  • Maple
    with(combinat);
    A:=proc(n) local i, j, t3, t2, t1;
        t2:=0;
        t1:=firstpart(n);
        for j from 1 to numbpart(n)+2 do
            t3:=1;
            for i from 1 to nops(t1) do
                if (t1[i] mod 6) <> 1 and (t1[i] mod 6) <> 5 then t3:=0; fi;
            od;
            if t3=1 then t2:=t2+1; fi;
            if nops(t1) = 1 then RETURN(t2); fi;
            t1:=nextpart(t1);
        od;
    end;
    # brute-force Maple program from N. J. A. Sloane, Jun 17 2011
  • Mathematica
    max = 63; f[x_] := 1/Product[1 - x^k + x^(2k), {k, 0, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 01 2011, after Vladeta Jovovic *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] / QPochhammer[ -x^3, x^3], {x, 0, n}]; (* Michael Somos, Jul 05 2014 *)
    nmax = 100; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[If[Mod[k, 3] != 0, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}];], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 13 2017 *)
    nmax = 63; kmax = nmax/6;
    s = Flatten[{Range[0, kmax]*6 + 1}~Join~{Range[kmax]*6 - 1}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Jul 31 2020 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A) / (eta(x + A) * eta(x^6 + A)), n))}; /* Michael Somos, Jan 09 2005 */
    
  • PARI
    {S(n,x)=sumdiv(n,d,d*(1-x^d)^(n/d))}
    {a(n)=polcoeff(exp(sum(k=1,n,S(k,x)*x^k/k)+x*O(x^n)),n)}
    for(n=0,60,print1(a(n),", "))
    /* Paul D. Hanna, Feb 17 2013 */
    

Formula

G.f.: 1/Product_{k>=0} (1-x^(6*k+1))*(1-x^(6*k+5)) = Product_{k>=0} (1+x^(3*k+1))*(1+x^(3*k+2)) = 1/Product_{k>=0} (1-x^k+x^(2*k)). - Vladeta Jovovic, Jun 08 2003
Expansion of chi(-x^3) / chi(-x) in powers of x where chi() is a Ramanujan theta function. - Michael Somos, Mar 04 2012
Expansion of f(x, x^2) / f(-x^3) = f(-x^6) / f(-x, -x^5) in powers of x where f() is Ramanujan theta function. - Michael Somos, Jul 05 2014
Expansion of q^(1/12) * eta(q^2) * eta(q^3) / (eta(q) * eta(q^6)) in powers of q. - Michael Somos, Jan 09 2005
Euler transform of period 6 sequence [1, 0, 0, 0, 1, 0, ...]. - Michael Somos, Jan 09 2005
Given g.f. A(x), then B(q) = (A(q^12) / q)^4 satisfies 0 = f(B(q), B(q^2)) where f(u, v) = u*v^4 + (1 - u^3) * v^3 + 6*u^2*v^2 + (u^4 - u)*v + u^3. - Michael Somos, Jan 09 2005
The logarithmic derivative equals A186099. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} A186099(n) * x^n/n ) where A186099(n) = sum of divisors of n congruent to 1 or 5 mod 6. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} S(n,x) * x^n/n ) where S(n,x) = Sum_{d|n} d*(1-x^d)^(n/d). - Paul D. Hanna, Feb 17 2013
a(n) ~ Pi*sqrt(2) / sqrt(3*(12*n-1)) * BesselI(1, Pi*sqrt(12*n-1) / (3*sqrt(6))) ~ exp(Pi*sqrt(2*n)/3) / (2^(5/4) * sqrt(3) * n^(3/4)) * (1 - (9/(8*Pi) + Pi/36)/sqrt(2*n) + (5 - 135/(4*Pi^2) + Pi^2/81)/(64*n)). - Vaclav Kotesovec, Aug 23 2015, extended Jan 09 2017
a(n) = (1/n)*Sum_{k=1..n} A186099(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 21 2017

Extensions

More terms from Vladeta Jovovic, Jun 08 2003

A002945 Continued fraction for cube root of 2.

Original entry on oeis.org

1, 3, 1, 5, 1, 1, 4, 1, 1, 8, 1, 14, 1, 10, 2, 1, 4, 12, 2, 3, 2, 1, 3, 4, 1, 1, 2, 14, 3, 12, 1, 15, 3, 1, 4, 534, 1, 1, 5, 1, 1, 121, 1, 2, 2, 4, 10, 3, 2, 2, 41, 1, 1, 1, 3, 7, 2, 2, 9, 4, 1, 3, 7, 6, 1, 1, 2, 2, 9, 3, 1, 1, 69, 4, 4, 5, 12, 1, 1, 5, 15, 1, 4
Offset: 0

Views

Author

Keywords

Examples

			2^(1/3) = 1.25992104989487316... = 1 + 1/(3 + 1/(1 + 1/(5 + 1/(1 + ...)))).
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002946, A002947, A002948, A002949, A002580 (decimal expansion).
Cf. A002351, A002352 (convergents).

Programs

  • Magma
    ContinuedFraction(2^(1/3)); // Vincenzo Librandi, Oct 08 2017
  • Maple
    N:= 100: # to get a(1) to a(N)
    a[1] := 1: p[1] := 1: q[1] := 0: p[2] := 1: q[2] := 1:
    for n from 2 to N do
      a[n] := floor((-1)^(n+1)*3*p[n]^2/(q[n]*(p[n]^3-2*q[n]^3)) - q[n-1]/q[n]);
      p[n+1] := a[n]*p[n] + p[n-1];
      q[n+1] := a[n]*q[n] + q[n-1];
    od:
    seq(a[i],i=1..N); # Robert Israel, Jul 30 2014
  • Mathematica
    ContinuedFraction[Power[2, (3)^-1],70] (* Harvey P. Dale, Sep 29 2011 *)
  • PARI
    allocatemem(932245000); default(realprecision, 21000); x=contfrac(2^(1/3)); for (n=1, 20000, write("b002945.txt", n-1, " ", x[n])); \\ Harry J. Smith, May 08 2009
    

Formula

From Robert Israel, Jul 30 2014: (Start)
Bombieri/van der Poorten give a complicated formula:
a(n) = floor((-1)^(n+1)*3*p(n)^2/(q(n)*(p(n)^3-2*q(n)^3)) - q(n-1)/q(n)),
p(n+1) = a(n)*p(n) + p(n-1),
q(n+1) = a(n)*q(n) + q(n-1),
with a(1) = 1, p(1) = 1, q(1) = 0, p(2) = 1, q(2) = 1. (End)

Extensions

BCMATH link from Keith R Matthews (keithmatt(AT)gmail.com), Jun 04 2006
Offset changed by Andrew Howroyd, Jul 04 2024

A034297 Number of ordered positive integer solutions (m_1, m_2, ..., m_k) (for some k) to Sum_{i=1..k} m_i=n with |m_i-m_{i-1}| <= 1 for i = 2 ... k.

Original entry on oeis.org

1, 1, 2, 4, 6, 11, 17, 29, 47, 78, 130, 215, 357, 595, 990, 1651, 2748, 4584, 7643, 12744, 21256, 35451, 59133, 98636, 164531, 274463, 457837, 763746, 1274060, 2125356, 3545491, 5914545, 9866602, 16459421, 27457549, 45804648, 76411272, 127469285, 212644336
Offset: 0

Views

Author

Keywords

Comments

Compositions of n where successive parts differ by at most 1, see example. [Joerg Arndt, Dec 10 2012]

Examples

			From _Joerg Arndt_, Dec 10 2012: (Start)
The a(6) = 17 such compositions of 6 are
[ #]     composition
[ 1]    [ 1 1 1 1 1 1 ]
[ 2]    [ 1 1 1 1 2 ]
[ 3]    [ 1 1 1 2 1 ]
[ 4]    [ 1 1 2 1 1 ]
[ 5]    [ 1 1 2 2 ]
[ 6]    [ 1 2 1 1 1 ]
[ 7]    [ 1 2 1 2 ]
[ 8]    [ 1 2 2 1 ]
[ 9]    [ 1 2 3 ]
[10]    [ 2 1 1 1 1 ]
[11]    [ 2 1 1 2 ]
[12]    [ 2 1 2 1 ]
[13]    [ 2 2 1 1 ]
[14]    [ 2 2 2 ]
[15]    [ 3 2 1 ]
[16]    [ 3 3 ]
[17]    [ 6 ]
(End)
		

Crossrefs

Column k=1 of A214246, A214248.
Row sums of A309939.

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=i, 1, `if`(n<0 or i<1, 0, add(b(n-i, i+j), j=-1..1)))
        end:
    a:= n-> add(b(n, k), k=0..n):
    seq(a(n), n=0..50);  # Alois P. Heinz, Jul 06 2012
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == i, 1, If[n<0 || i<1, 0, Sum[b[n-i, i+j], {j, -1, 1}] ]]; a[n_] := Sum[b[n, k], {k, 1, n}]; Array[a, 50] (* Jean-François Alcover, Mar 13 2015, after Alois P. Heinz *)
  • PARI
    N=70;  nil=-1;
    T = matrix(N, N, i, j, nil);
    doIt(last, left) = my(c); c = T[last, left]; if (c == nil, c = 0; for (i = max(1, last - 1), last + 1, c += b(i, left - i)); T[last, left] = c); c;
    b(last, left) = if (left == 0, return(1)); if (left < 0, return(0)); doIt(last, left);
    a(n) = sum (i = 1, n, b(i, n - i));
    vector(N,n,a(n))  \\ David Wasserman, Feb 02 2006
    
  • Python
    from sympy.core.cache import cacheit
    @cacheit
    def b(n, i): return 1 if n==i else 0 if n<0 or i<1 else sum(b(n - i, i + j) for j in range(-1, 2))
    def a(n): return sum(b(n, k) for k in range(n + 1))
    print([a(n) for n in range(51)]) # Indranil Ghosh, Aug 14 2017, after Maple code

Formula

a(n) ~ c * d^n, where d = 1.668202067018461116361070469945501401879811945303435230637248..., c = 0.762436680050402638439806786781869262562176911054246754543346... . - Vaclav Kotesovec, Sep 02 2014

Extensions

More terms from David Wasserman, Feb 02 2006
a(0)=1 prepended by Alois P. Heinz, Aug 14 2017

A002949 Continued fraction for cube root of 6.

Original entry on oeis.org

1, 1, 4, 2, 7, 3, 508, 1, 5, 5, 1, 1, 1, 2, 1, 1, 24, 1, 1, 1, 3, 3, 30, 4, 10, 158, 6, 1, 1, 2, 12, 1, 10, 1, 1, 3, 2, 1, 1, 89, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 7, 1, 2, 18, 1, 17, 2, 2, 10, 14, 3, 1, 2, 1, 2, 1, 5, 1, 1, 2, 26, 1, 4, 65, 1, 1, 1, 27, 1, 2, 1, 4
Offset: 0

Views

Author

Keywords

Examples

			6^(1/3) = 1.81712059283213965... = 1 + 1/(1 + 1/(4 + 1/(2 + 1/(7 + ...)))). - _Harry J. Smith_, May 08 2009
		

References

  • H. P. Robinson, Letter to N. J. A. Sloane, Nov 13 1973.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A005486 (decimal expansion).
Cf. A002359, A002360 (convergents).

Programs

  • Magma
    SetDefaultRealField(RealField(100)); ContinuedFraction(6^(1/3)); // G. C. Greubel, Nov 02 2018
  • Maple
    with(numtheory):
    cfrac(6^(1/3),100,'quotients'); # Muniru A Asiru, Nov 02 2018
  • Mathematica
    ContinuedFraction[6^(1/3), 100] (* G. C. Greubel, Nov 02 2018 *)
  • PARI
    { allocatemem(932245000); default(realprecision, 21000); x=contfrac(6^(1/3)); for (n=1, 20000, write("b002949.txt", n-1, " ", x[n])); } \\ Harry J. Smith, May 08 2009
    

Extensions

Offset changed by Andrew Howroyd, Jul 05 2024

A039924 G.f.: Sum_{k>=0} x^(k^2)*(-1)^k/(Product_{i=1..k} 1-x^i).

Original entry on oeis.org

1, -1, -1, -1, 0, 0, 1, 1, 2, 1, 2, 1, 1, 0, 0, -2, -1, -3, -3, -4, -3, -5, -3, -4, -2, -3, 0, -1, 3, 2, 5, 5, 9, 7, 11, 9, 13, 10, 13, 9, 12, 7, 9, 3, 5, -3, -1, -9, -7, -17, -15, -24, -21, -31, -27, -37, -31, -40, -33, -41, -31, -39, -27, -33, -18, -24, -6, -11, 9, 5, 26, 23
Offset: 0

Views

Author

Keywords

Comments

Ramanujan used the form Sum_{k>=0} x^(k^2) / (Product_{i=1..k} 1-(-x)^i), which is obtained by changing the sign of x. - Michael Somos, Jul 20 2003
Coefficients in expansion of determinant of infinite tridiagonal matrix shown below in powers of x^2 (Lehmer 1973):
1 x 0 0 0 0 ...
x 1 x^2 0 0 0 ...
0 x^2 1 x^3 0 0 ...
0 0 x^3 1 x^4 0 ...
... ... ... ... ... ... ...
Convolution inverse of A003116.

Examples

			G.f. = 1 - x - x^2 - x^3 + x^6 + x^7 + 2*x^8 + x^9 + 2*x^10 + x^11 + x^12 + ...
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 55, Eq. (26.11).
  • G. H. Hardy, P. V. Seshu Aiyar and B. M. Wilson, editors, Collected Papers of Srinivasa Ramanujan, Cambridge, 1923; p. 354.
  • D. H. Lehmer, Combinatorial and cyclotomic properties of certain tridiagonal matrices. Proceedings of the Fifth Southeastern Conference on Combinatorics, Graph Theory, and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1974), pp. 53-74. Congressus Numerantium, No. X, Utilitas Math., Winnipeg, Man., 1974. MR0441852.
  • Herman P. Robinson, personal communication to N. J. A. Sloane.

Crossrefs

Programs

  • Maple
    qq:=n->mul( 1-(-q)^i, i=1..n); add (q^(n^2)/qq(n),n=0..100): series(t1,q,99);
  • Mathematica
    CoefficientList[ Series[ Sum[x^k^2*(-1)^k / Product[1-x^i,{i,1,k}], {k,0,100}], {x,0,100}],x][[1 ;; 72]] (* Jean-François Alcover, Apr 08 2011 *)
    a[ n_] := If[n < 0, 0, SeriesCoefficient[ Sum[ (-1)^k x^k^2 / QPochhammer[ x, x, k], {k, 0, Sqrt[n]}], {x, 0, n}]] (* Michael Somos, Jan 04 2014 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, sqrtint(n), x^k^2 / prod(i=1, k, x^i - 1, 1 + x * O(x^n))), n))} /* Michael Somos, Jul 20 2003 */

Formula

a(n) = A286041(n) - A286316(n) (conjectured). - George Beck, May 05 2017
Proof from Doron Zeilberger, Aug 20 2018: (Start)
The generating function for partitions whose parts differ by at least 2 with exactly k parts is (famously) q^(k^2)/((1-q)*...*(1-q^k)).
Indeed, if you take any such partition and remove 1 from the smallest part, 3 from the second-smallest part, etc., you remove 1+3+...+(2k-1) = k^2 and are left with an ordinary partition whose number of parts is <= k whose generating function is 1/((1-q)*...*(1-q)^k).
Summing these up famously gives the generating function for partitions whose differences is >= 2. Sticking a (-1)^k in front gives the generating function for the difference between such partitions with an even number of parts and an odd number of parts, since (-1)^even=1 and (-1)^odd=-1. (End)

Extensions

More terms from Vladeta Jovovic, Mar 05 2001

A111537 Column 1 of triangle A111536.

Original entry on oeis.org

1, 2, 8, 44, 296, 2312, 20384, 199376, 2138336, 24936416, 314142848, 4252773824, 61594847360, 950757812864, 15586971531776, 270569513970944, 4959071121374720, 95721139472072192, 1941212789888952320, 41271304403571227648, 918030912312297752576, 21325054720042613565440
Offset: 0

Views

Author

Paul D. Hanna, Aug 06 2005

Keywords

Comments

Row sums of triangle in A200659. - Philippe Deléham, Nov 21 2011

References

  • A. N. Khovanskii. The Application of Continued Fractions and Their Generalizations to Problem in Approximation Theory. Groningen: Noordhoff, Netherlands, 1963. See p.141 (10.19).
  • H. P. Robinson, Letter to N. J. A. Sloane, Nov 19 1973.

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 1,
          n*(n+1)! -add((n-k+1)!*a(k), k=1..n-1))
        end:
    seq(a(n), n=0..30);  # Alois P. Heinz, May 06 2013
  • Mathematica
    a[n_] := a[n] = If[n==0, 1, n*(n+1)! - Sum[(n-k+1)!*a[k], {k, 1, n-1}]]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 13 2017, after Alois P. Heinz *)
  • PARI
    {a(n)=if(n<0,0,(matrix(n+2,n+2,m,j,if(m==j,1,if(m==j+1,-m+1, -(m-j-1)*polcoeff(log(sum(i=0,m,(i+1)!/1!*x^i)),m-j-1))))^-1)[n+2,2])}

Formula

a(n) = A111536(n+1, 1) = 2*A111536(n, 0) = 2*A111529(n) for n >= 1.
G.f.: log(Sum_{n>=0} (n+1)!*x^n) = Sum_{n>=1} a(n)*x^n/n.
a(n+1) = (n+3)! - 2*(n+2)! - Sum_{k=0..n-1} (n-k+1)!*a(k+1).
a(n+1) is the moment of order n for the measure of density x*exp(-x)/((x*exp(-x)*Ei(x)-1)^2+(Pi*x*exp(-x))^2) on the interval 0..infinity.
G.f.: 1/(1-2*x/(1-2*x/(1-3*x/(1-3*x/(1-4*x/(1-4*x/(1-5*x/(1-...(continued fraction). - Philippe Deléham, Nov 21 2011
G.f. (1-U(0))/x; where U(k) = 1-x*(k+1)/(1-x*(k+2)/U(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Jun 29 2012
G.f. -1 + 1/x + U(0)/x where U(k) = 2*x - 1 + 2*x*k - x^2*(k+1)*(k+2)/U(k+1), U(0)=x - W(1,1;-x)/W(1,2;-x), W(a,b,x)= 1 - a*b*x/1! + a*(a+1)*b*(b+1)*x^2/2! - ... + a*(a+1)*...*(a+n-1)*b*(b+1)*...*(b+n-1)*x^n/n! + ...; see [A. N. Khovanskii, p. 141 (10.19)]; (continued fraction, 1-step). - Sergei N. Gladkovskii, Aug 15 2012
G.f.: 1/Q(0), where Q(k) = 1 + k*x - x*(k+2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 03 2013
G.f.: 1/x - 1/( x*G(0)), where G(k) = 1 - x*(k+1)/(x - 1/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Aug 03 2013
a(0) = 1; a(n) = n * a(n-1) + Sum_{k=0..n-1} a(k) * a(n-k-1). - Ilya Gutkovskiy, Jul 05 2020

A002950 Continued fraction for fifth root of 2.

Original entry on oeis.org

1, 6, 1, 2, 1, 1, 1, 3, 25, 1, 4, 3, 3, 7, 52, 1, 2, 3, 2, 15, 2, 2, 4, 16, 2, 7, 1, 1, 1, 10, 21, 1, 1, 1, 141, 2, 4, 1, 4, 2, 1, 1, 17, 1, 3, 3, 4, 1, 3, 1, 3, 2, 1, 1, 2, 33, 1, 6, 6, 1, 2, 4, 1, 21, 1, 3, 3, 8, 10, 1, 46, 6, 1, 10, 1, 1, 1, 1, 2, 11, 1, 3, 1
Offset: 0

Views

Author

Keywords

Examples

			2^(1/5) = 1.148698354997035006798626946... = 1 + 1/(6 + 1/(1 + 1/(2 + 1/(1 + ...)))). - _Harry J. Smith_, May 12 2009
		

References

  • H. P. Robinson, Letter to N. J. A. Sloane, Nov 13 1973.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A005531 (decimal expansion).
Cf. A002361, A002362 (convergents).

Programs

  • Magma
    SetDefaultRealField(RealField(100)); ContinuedFraction(2^(1/5)); // G. C. Greubel, Nov 02 2018
  • Maple
    with(numtheory):
    cfrac(2^(1/5),100,'quotients'); # Muniru A Asiru, Nov 02 2018
  • Mathematica
    ContinuedFraction[2^(1/5), 100] (* G. C. Greubel, Nov 02 2018 *)
  • PARI
    { allocatemem(932245000); default(realprecision, 21000); x=contfrac(2^(1/5)); for (n=1, 20000, write("b002950.txt", n-1, " ", x[n])); } \\ Harry J. Smith, May 12 2009
    

Extensions

Offset changed by Andrew Howroyd, Jul 05 2024

A225084 Triangle read by rows: T(n,k) is the number of compositions of n with maximal up-step k; n>=1, 0<=k

Original entry on oeis.org

1, 2, 0, 3, 1, 0, 5, 2, 1, 0, 7, 6, 2, 1, 0, 11, 12, 6, 2, 1, 0, 15, 26, 14, 6, 2, 1, 0, 22, 50, 33, 14, 6, 2, 1, 0, 30, 97, 72, 34, 14, 6, 2, 1, 0, 42, 180, 156, 77, 34, 14, 6, 2, 1, 0, 56, 332, 328, 173, 78, 34, 14, 6, 2, 1, 0, 77, 600, 681, 378, 177, 78, 34, 14, 6, 2, 1, 0, 101, 1078, 1393, 818, 393, 178, 78, 34, 14, 6, 2, 1, 0
Offset: 1

Views

Author

Joerg Arndt, Apr 27 2013

Keywords

Comments

T(n,k) is the number of compositions [p(1), p(2), ..., p(k)] of n such that max(p(j) - p(j-1)) == k.
The first column is A000041 (partition numbers).
Sum of first and second column is A003116.
Sum of the first three columns is A224959.
The second columns deviates from A054454 after the term 600.
Row sums are A011782.

Examples

			Triangle starts:
01: 1,
02: 2, 0,
03: 3, 1, 0,
04: 5, 2, 1, 0,
05: 7, 6, 2, 1, 0,
06: 11, 12, 6, 2, 1, 0,
07: 15, 26, 14, 6, 2, 1, 0,
08: 22, 50, 33, 14, 6, 2, 1, 0,
09: 30, 97, 72, 34, 14, 6, 2, 1, 0,
10: 42, 180, 156, 77, 34, 14, 6, 2, 1, 0,
11: 56, 332, 328, 173, 78, 34, 14, 6, 2, 1, 0,
12: 77, 600, 681, 378, 177, 78, 34, 14, 6, 2, 1, 0,
13: 101, 1078, 1393, 818, 393, 178, 78, 34, 14, 6, 2, 1, 0,
14: 135, 1917, 2821, 1746, 863, 397, 178, 78, 34, 14, 6, 2, 1, 0,
15: 176, 3393, 5660, 3695, 1872, 877, 398, 178, 78, 34, 14, 6, 2, 1, 0,
...
The fifth row corresponds to the following statistics:
#:  M   composition
01:  0  [ 1 1 1 1 1 ]
02:  1  [ 1 1 1 2 ]
03:  1  [ 1 1 2 1 ]
04:  2  [ 1 1 3 ]
05:  1  [ 1 2 1 1 ]
06:  1  [ 1 2 2 ]
07:  2  [ 1 3 1 ]
08:  3  [ 1 4 ]
09:  0  [ 2 1 1 1 ]
10:  1  [ 2 1 2 ]
11:  0  [ 2 2 1 ]
12:  1  [ 2 3 ]
13:  0  [ 3 1 1 ]
14:  0  [ 3 2 ]
15:  0  [ 4 1 ]
16:  0  [ 5 ]
There are 7 compositions with no up-step (M=0), 6 with M=1, 2 with M=2, and 1 with M=3.
		

Crossrefs

Cf. A225085 (partial sums of rows).
T(2n,n) gives A059570.

Programs

  • Maple
    b:= proc(n, v) option remember; `if`(n=0, 1, add((p->
          `if`(i seq(coeff(b(n, 0), x, i), i=0..n-1):
    seq(T(n), n=1..14);  # Alois P. Heinz, Feb 22 2014
  • Mathematica
    b[n_, v_] := b[n, v] = If[n == 0, 1, Sum[Function[{p}, If[iJean-François Alcover, Feb 18 2015, after Alois P. Heinz *)
Showing 1-10 of 23 results. Next