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-2 of 2 results.

A026820 Euler's table: triangular array T read by rows, where T(n,k) = number of partitions in which every part is <= k for 1 <= k <= n. Also number of partitions of n into at most k parts.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 3, 4, 5, 1, 3, 5, 6, 7, 1, 4, 7, 9, 10, 11, 1, 4, 8, 11, 13, 14, 15, 1, 5, 10, 15, 18, 20, 21, 22, 1, 5, 12, 18, 23, 26, 28, 29, 30, 1, 6, 14, 23, 30, 35, 38, 40, 41, 42, 1, 6, 16, 27, 37, 44, 49, 52, 54, 55, 56, 1, 7, 19, 34, 47, 58, 65, 70, 73, 75, 76, 77
Offset: 1

Views

Author

Keywords

Examples

			Triangle starts:
  1;
  1, 2;
  1, 2,  3;
  1, 3,  4,  5;
  1, 3,  5,  6,  7;
  1, 4,  7,  9, 10, 11;
  1, 4,  8, 11, 13, 14, 15;
  1, 5, 10, 15, 18, 20, 21, 22;
  1, 5, 12, 18, 23, 26, 28, 29, 30;
  1, 6, 14, 23, 30, 35, 38, 40, 41, 42;
  1, 6, 16, 27, 37, 44, 49, 52, 54, 55, 56;
  ...
		

References

  • G. Chrystal, Algebra, Vol. II, p. 558.
  • D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section XIV.2, p. 493.

Crossrefs

Partial sums of rows of A008284, row sums give A058397, central terms give A171985, mirror is A058400.
T(n,n) = A000041(n), T(n,1) = A000012(n), T(n,2) = A008619(n) for n>1, T(n,3) = A001399(n) for n>2, T(n,4) = A001400(n) for n>3, T(n,5) = A001401(n) for n>4, T(n,6) = A001402(n) for n>5, T(n,7) = A008636(n) for n>6, T(n,8) = A008637(n) for n>7, T(n,9) = A008638(n) for n>8, T(n,10) = A008639(n) for n>9, T(n,11) = A008640(n) for n>10, T(n,12) = A008641(n) for n>11, T(n,n-2) = A007042(n-1) for n>2, T(n,n-1) = A000065(n) for n>1.

Programs

  • Haskell
    import Data.List (inits)
    a026820 n k = a026820_tabl !! (n-1) !! (k-1)
    a026820_row n = a026820_tabl !! (n-1)
    a026820_tabl = zipWith
       (\x -> map (p x) . tail . inits) [1..] $ tail $ inits [1..] where
       p 0 _ = 1
       p _ [] = 0
       p m ks'@(k:ks) = if m < k then 0 else p (m - k) ks' + p m ks
    -- Reinhard Zumkeller, Dec 18 2013
    
  • Maple
    T:= proc(n, k) option remember;
          `if`(n=0 or k=1, 1, T(n, k-1) + `if`(k>n, 0, T(n-k, k)))
        end:
    seq(seq(T(n, k), k=1..n), n=1..12); # Alois P. Heinz, Apr 21 2012
  • Mathematica
    t[n_, k_] := Length@ IntegerPartitions[n, k]; Table[ t[n, k], {n, 12}, {k, n}] // Flatten
    (* Second program: *)
    T[n_, k_] := T[n, k] = If[n==0 || k==1, 1, T[n, k-1] + If[k>n, 0, T[n-k, k]]]; Table[T[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Jean-François Alcover, Sep 22 2015, after Alois P. Heinz *)
  • PARI
    T(n,k)=my(s); forpart(v=n,s++,,k); s \\ Charles R Greathouse IV, Feb 27 2018
    
  • SageMath
    from sage.combinat.partition import number_of_partitions_length
    from itertools import accumulate
    for n in (1..11):
        print(list(accumulate([number_of_partitions_length(n, k) for k in (1..n)])))
    # Peter Luschny, Jul 28 2022

Formula

T(T(n,n),n) = A134737(n). - Reinhard Zumkeller, Nov 07 2007
T(A000217(n),n) = A173519(n). - Reinhard Zumkeller, Feb 20 2010
T(n,k) = T(n,k-1) + T(n-k,k). - Thomas Dybdahl Ahle, Jun 13 2011
T(n,k) = Sum_{i=1..min(k,floor(n/2))} T(n-i,i) + Sum_{j=1+floor(n/2)..k} A000041(n-j). - Bob Selcoe, Aug 22 2014 [corrected by Álvar Ibeas, Mar 15 2018]
O.g.f.: Product_{i>=0} 1/(1-y*x^i). - Geoffrey Critzer, Mar 11 2012
T(n,k) = A008284(n+k,k). - Álvar Ibeas, Jan 06 2015

A119620 Number of partitions of floor(3n/2) into n parts each from {1,2,...,n}.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 5, 5, 7, 7, 11, 11, 15, 15, 22, 22, 30, 30, 42, 42, 56, 56, 77, 77, 101, 101, 135, 135, 176, 176, 231, 231, 297, 297, 385, 385, 490, 490, 627, 627, 792, 792, 1002, 1002, 1255, 1255, 1575, 1575, 1958, 1958, 2436, 2436, 3010, 3010, 3718, 3718
Offset: 0

Views

Author

John W. Layman, Jun 07 2006

Keywords

Comments

The bisection {1,1,2,3,5,7,11,15,22,...} agrees with the initial terms of A008641, Number of partitions of n into at most 12 parts and also A008635, Molien series for A_12.
a(2n+1)=a(2n) for all n>0. If the partition {...,1} is a member of a(2n) then the partition {...,1,1} is a member of a(2n+1). - Robert G. Wilson v, Jun 09 2006
Number of partitions of n where all parts (except for possibly the first part) are even; see example. - Joerg Arndt, Apr 22 2013
For n >= 2, a(n) = number of partitions p of n such that floor(n/2) is a part of p. For n >= 1, a(n) = number of partitions p of n such that ceiling(n/2) is a part of p. - Clark Kimberling, Feb 28 2014
From Gus Wiseman, Oct 28 2021: (Start)
If we insert zeros every three terms, this counts partitions of n such that n = floor(3*k/2), where k is the number of parts. This counts by sum rather than length. These partitions are ranked by A347452.
Also the number of integer partitions of n with alternating product 1, where the alternating product of a sequence (y_1,...,y_k) is Product_i y_i^((-1)^(i-1)). These are the conjugates of the partitions (ranked by A336119) described in Arndt's comment above. For example, the a(2) = 1 through a(10) = 7 partitions are:
11 111 22 221 33 331 44 441 55
1111 11111 2211 22111 2222 22221 3322
111111 1111111 3311 33111 4411
221111 2211111 222211
11111111 111111111 331111
22111111
1111111111
These partitions are ranked by A028982. The odd-length case is A035363 (shifted), which is also the version for sum instead of product. The multiplicative version (factorizations) is A347438.
(End)

Examples

			For n=8, floor(3*n/2) is 12 and there are five partitions of 12 into 8 parts each in the range 1-8 inclusive, namely: {5,1,1,1,1,1,1,1}, {4,2,1,1,1,1,1,1}, {3,3,1,1,1,1,1,1}, {3,2,2,1,1,1,1,1} and {2,2,2,2,1,1,1,1}. Thus a(8)=5.
From _Joerg Arndt_, Apr 22 2013: (Start)
a(8) = a(9) = 5, counting the following partitions where all parts (except for possibly the first part) are even:
01:  [ 2 2 2 2 ]
02:  [ 4 2 2 ]
03:  [ 4 4 ]
04:  [ 6 2 ]
05:  [ 8 ]
and
01:  [ 3 2 2 2 ]
02:  [ 5 2 2 ]
03:  [ 5 4 ]
04:  [ 7 2 ]
05:  [ 9 ]
(End)
G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 3*x^6 + 3*x^7 + 5*x^8 + 5*x^9 + 7*x^10 + ...
		

Crossrefs

Both bisections are A000041.
An adjoint version is A108711.
A027187 counts partitions of even length.
A027193 counts partitions of odd length.
A325534 counts separable partitions.
A325535 counts inseparable partitions.

Programs

  • Maple
    # Using the function EULER from Transforms (see link at the bottom of the page).
    [1, op(EULER([1,0,seq(irem(n,2),n=2..55)]))]; # Peter Luschny, Aug 19 2020
  • Mathematica
    (* first do *) Needs["DiscreteMath`Combinatorica`"] (* then *) f[n_] := f[n] = Length@ Select[ Partitions[ Floor[3n/2], n], Length@# == n &]; Table[ If[n > 1, f[2Floor[n/2]], f[n]], {n, 57}] (* Robert G. Wilson v, Jun 09 2006 *)
    Table[ PartitionsP[ Floor[n/2]], {n, 57}] (* Robert G. Wilson v, Jun 09 2006 *)
    Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, Ceiling[n/2]]], {n, 50}] (* Clark Kimberling, Feb 28 2014 *)
    a[ n_] := SeriesCoefficient[ (1 + x) / QPochhammer[x^2], {x, 0, n}]; (* Michael Somos, Mar 01 2014 *)
  • PARI
    a(n)=numbpart(n\2); \\ Joerg Arndt, Apr 22 2013

Formula

a(n) = A000041(floor(n/2)). - Vladeta Jovovic, Jun 10 2006
G.f.: (Sum_{n>=0} x^(4*n) / Product_{k=1..n} (1-x^(2*k))) / (1 - x). - Michael Somos, Mar 01 2014 [corrected by Jason Yuen, Jan 24 2025]

Extensions

More terms from Robert G. Wilson v, Jun 09 2006
Added a(0)=1. - Michael Somos, Mar 01 2014
Showing 1-2 of 2 results.