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.

Previous Showing 51-60 of 579 results. Next

A351293 Number of non-Look-and-Say partitions of n. Number of integer partitions of n such that there is no way to choose a disjoint strict integer partition of each multiplicity.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 4, 5, 9, 14, 21, 28, 44, 56, 80, 111, 148, 192, 264, 335, 447, 575, 743, 937, 1213, 1513, 1924, 2396, 3011, 3715, 4646, 5687, 7040, 8600, 10556, 12804, 15650, 18897, 22930, 27593, 33296, 39884, 47921, 57168, 68360, 81295, 96807, 114685
Offset: 0

Views

Author

Gus Wiseman, Feb 16 2022

Keywords

Comments

First differs from A336866 (non-Wilf partitions) at a(9) = 14, A336866(9) = 15, the difference being the partition (2,2,2,1,1,1).
See A239455 for the definition of Look-and-Say partitions.

Examples

			The a(3) = 1 through a(9) = 14 partitions:
  (21)  (31)  (32)  (42)    (43)    (53)     (54)
              (41)  (51)    (52)    (62)     (63)
                    (321)   (61)    (71)     (72)
                    (2211)  (421)   (431)    (81)
                            (3211)  (521)    (432)
                                    (3221)   (531)
                                    (3311)   (621)
                                    (4211)   (3321)
                                    (32111)  (4221)
                                             (4311)
                                             (5211)
                                             (32211)
                                             (42111)
                                             (321111)
		

Crossrefs

The complement is counted by A239455, ranked by A351294.
These are all non-Wilf partitions (counted by A336866, ranked by A130092).
A variant for runs is A351203, complement A351204, ranked by A351201.
These partitions appear to be ranked by A351295.
Non-Wilf partitions in the complement are counted by A351592.
A000569 = graphical partitions, complement A339617.
A032020 = number of binary expansions with all distinct run-lengths.
A044813 = numbers whose binary expansion has all distinct run-lengths.
A098859 = Wilf partitions (distinct multiplicities), ranked by A130091.
A181819 = Heinz number of the prime signature of n (prime shadow).
A329738 = compositions with all equal run-lengths.
A329739 = compositions with all distinct run-lengths, for all runs A351013.
A351017 = binary words with all distinct run-lengths, for all runs A351016.
A351292 = patterns with all distinct run-lengths, for all runs A351200.

Programs

  • Mathematica
    disjointFamilies[y_]:=Select[Tuples[IntegerPartitions/@Length/@Split[y]],UnsameQ@@Join@@#&];
    Table[Length[Select[IntegerPartitions[n],Length[disjointFamilies[#]]==0&]],{n,0,15}] (* Gus Wiseman, Aug 13 2025 *)

Formula

a(n) = A000041(n) - A239455(n).

Extensions

Edited by Gus Wiseman, Aug 12 2025

A359901 Triangle read by rows where T(n,k) is the number of integer partitions of n with median k = 1..n.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 2, 2, 0, 1, 3, 1, 0, 0, 1, 4, 2, 3, 0, 0, 1, 6, 3, 1, 0, 0, 0, 1, 8, 6, 2, 4, 0, 0, 0, 1, 11, 7, 3, 1, 0, 0, 0, 0, 1, 15, 10, 4, 2, 5, 0, 0, 0, 0, 1, 20, 13, 7, 3, 1, 0, 0, 0, 0, 0, 1, 26, 19, 11, 4, 2, 6, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 21 2023

Keywords

Comments

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Examples

			Triangle begins:
   1
   1  1
   1  0  1
   2  2  0  1
   3  1  0  0  1
   4  2  3  0  0  1
   6  3  1  0  0  0  1
   8  6  2  4  0  0  0  1
  11  7  3  1  0  0  0  0  1
  15 10  4  2  5  0  0  0  0  1
  20 13  7  3  1  0  0  0  0  0  1
  26 19 11  4  2  6  0  0  0  0  0  1
  35 24 14  5  3  1  0  0  0  0  0  0  1
  45 34 17  8  4  2  7  0  0  0  0  0  0  1
  58 42 23 12  5  3  1  0  0  0  0  0  0  0  1
For example, row n = 9 counts the following partitions:
  (7,1,1)              (5,2,2)      (3,3,3)  (4,4,1)  .  .  .  .  (9)
  (6,1,1,1)            (6,2,1)      (4,3,2)
  (3,3,1,1,1)          (3,2,2,2)    (5,3,1)
  (4,2,1,1,1)          (4,2,2,1)
  (5,1,1,1,1)          (4,3,1,1)
  (3,2,1,1,1,1)        (2,2,2,2,1)
  (4,1,1,1,1,1)        (3,2,2,1,1)
  (2,2,1,1,1,1,1)
  (3,1,1,1,1,1,1)
  (2,1,1,1,1,1,1,1)
  (1,1,1,1,1,1,1,1,1)
		

Crossrefs

Column k=1 is A027336(n+1).
For mean instead of median we have A058398, see also A008284, A327482.
Row sums are A325347.
The mean statistic is ranked by A326567/A326568.
Including half-steps gives A359893.
The odd-length case is A359902.
The median statistic is ranked by A360005(n)/2.
First appearances of medians are ranked by A360006, A360007.
A000041 counts partitions, strict A000009.
A027193 counts odd-length partitions, strict A067659, ranked by A026424.
A067538 counts partitions w/ integer mean, strict A102627, ranks A316413.
A240219 counts partitions w/ the same mean as median, complement A359894.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Median[#]==k&]],{n,15},{k,n}]

A069905 Number of partitions of n into 3 positive parts.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120, 127, 133, 140, 147, 154, 161, 169, 176, 184, 192, 200, 208, 217, 225, 234, 243, 252, 261, 271, 280, 290, 300, 310, 320, 331, 341
Offset: 0

Views

Author

N. J. A. Sloane, May 04 2002

Keywords

Comments

Number of binary bracelets of n beads, 3 of them 0. For n >= 3, a(n-3) is the number of binary bracelets of n beads, 3 of them 0, with 00 prohibited. - Washington Bomfim, Aug 27 2008
Also number of partitions of n-3 into parts 1, 2, and 3. - Joerg Arndt, Sep 05 2013
Number of incongruent triangles with integer sides that have perimeter 2n-3 (see the Jordan et al. link). - Freddy Barrera, Aug 18 2018
Number of ordered triples (x,y,z) of nonnegative integers such that x+y+z=n and xDennis P. Walsh, Apr 19 2019
Number of incongruent triangles formed from any 3 vertices of a regular n-gon. - Frank M Jackson, Sep 11 2022
Also a(n-3) for n > 2, otherwise 0 is the number of incongruent scalene triangles formed from the vertices of a regular n-gon. - Frank M Jackson, Nov 27 2022

Examples

			G.f. = x^3 + x^4 + 2*x^5 + 3*x^6 + 4*x^7 + 5*x^8 + 7*x^9 + 8*x^10 + 10*x^11 + ...
		

References

  • Ross Honsberger, Mathematical Gems III, Math. Assoc. Amer., 1985, p. 39.
  • Donald E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, Section 7.2.1.4, p. 410.
  • Donald E. Knuth, The Art of Computer Programming, vol. 4,fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.

Crossrefs

Another version of A001399, which is the main entry for this sequence.
Cf. A005044, A008284, A008615, A026810 (4 positive parts).

Programs

  • GAP
    List([0..70],n->NrPartitions(n,3)); # Muniru A Asiru, May 17 2018
    
  • Haskell
    a069905 n = a069905_list !! n
    a069905_list = scanl (+) 0 a008615_list
    -- Reinhard Zumkeller, Apr 28 2014
    
  • Magma
    [(n^2+6) div 12: n in [0..70]]; // Vincenzo Librandi, Oct 14 2015
    
  • Maple
    A069905 := n->round(n^2/12): seq(A069905(n), n=0..70);
  • Mathematica
    a[ n_]:= Round[ n^2 / 12] (* Michael Somos, Sep 04 2013 *)
    CoefficientList[Series[x^3/((1-x)(1-x^2)(1-x^3)), {x, 0, 70}], x] (* Vincenzo Librandi, Oct 14 2015 *)
    Drop[LinearRecurrence[{1,1,0,-1,-1,1}, Append[Table[0,{5}],1],70],2] (* Robert A. Russell, May 17 2018 *)
  • PARI
    a(n) = floor((n^2+6)/12);  \\ Washington Bomfim, Jul 03 2012
    
  • PARI
    my(x='x+O('x^70)); concat([0, 0, 0], Vec(x^3/((1-x)*(1-x^2)*(1-x^3)))) \\ Altug Alkan, Oct 14 2015
    
  • SageMath
    [round(n^2/12) for n in range(70)] # G. C. Greubel, Apr 03 2019

Formula

G.f.: x^3/((1-x)*(1-x^2)*(1-x^3)) = x^3/((1-x)^3*(1+x+x^2)*(1+x)).
a(n) = round(n^2/12).
a(n) = floor((n^2+6)/12). - Washington Bomfim, Jul 03 2012
a(-n) = a(n). - Michael Somos, Sep 04 2013
a(n) = a(n-1) + A008615(n-1) for n > 0. - Reinhard Zumkeller, Apr 28 2014
Let n = 6k + m. Then a(n) = n^2/12 + a(m) - m^2/12. Also, a(n) = 3*k^2 + m*k + a(m). Example: a(35) = a(6*5 + 5) = 35^2/12 + a(5) - 5^2/12 = 102 = 3*5^2 + 5*5 + a(5). - Gregory L. Simay, Oct 13 2015
a(n) = a(n-1) +a(n-2) -a(n-4) -a(n-5) +a(n-6), n>5. - Wesley Ivan Hurt, Oct 16 2015
a(n) = A008284(n,3). - Robert A. Russell, May 13 2018
a(n) = A005044(2*n) = A005044(2*n - 3). - Freddy Barrera, Aug 18 2018
a(n) = floor((n^2+k)/12) for all integers k such that 3 <= k <= 7. - Giacomo Guglieri, Apr 03 2019
From Wesley Ivan Hurt, Apr 19 2019: (Start)
a(n) = Sum_{k=1..floor(n/3)} Sum_{i=k..floor((n-k)/2)} 1.
a(n) = Sum_{i=1..floor(n/3)} floor((n-i)/2) - i + 1. (End)
Sum_{n>=3} 1/a(n) = 15/4 + Pi^2/18 - Pi/(2*sqrt(3)) + tanh(Pi/(2*sqrt(3))) * Pi/sqrt(3). - Amiram Eldar, Sep 27 2022
E.g.f.: (8*exp(-x/2)*cos(sqrt(3)*x/2) + (3*x^2 + 3*x - 8)*cosh(x) + (3*x^2 + 3*x + 1)*sinh(x))/36. - Stefano Spezia, Apr 05 2023
From Ridouane Oudra, Dec 12 2024: (Start)
a(n) = (n^2 + 2*gcd(n,3) - 3*gcd(n,2))/12.
a(n) = (A198442(n) + A079978(n))/3.
a(n) = A000212(n) - A002620(n).
a(n) = A008133(n+1) - A307018(n+1). (End)
a(n) = (A309511(n) + A309513(n))/3. - Ray Chandler, Mar 13 2025

A115994 Triangle read by rows: T(n,k) is number of partitions of n with Durfee square of size k (n>=1; 1<=k<=floor(sqrt(n))).

Original entry on oeis.org

1, 2, 3, 4, 1, 5, 2, 6, 5, 7, 8, 8, 14, 9, 20, 1, 10, 30, 2, 11, 40, 5, 12, 55, 10, 13, 70, 18, 14, 91, 30, 15, 112, 49, 16, 140, 74, 1, 17, 168, 110, 2, 18, 204, 158, 5, 19, 240, 221, 10, 20, 285, 302, 20, 21, 330, 407, 34, 22, 385, 536, 59, 23, 440, 698, 94, 24, 506, 896, 149, 25
Offset: 1

Views

Author

Emeric Deutsch, Feb 11 2006

Keywords

Comments

Row n has floor(sqrt(n)) terms. Row sums yield A000041. Column 2 yields A006918. sum(k*T(n,k),k=1..floor(sqrt(n)))=A115995.
T(n,k) is number of partitions of n-k^2 into parts of 2 kinds with at most k of each kind.
The limit of the diagonals is A000712 (partitions into parts of two kinds). In particular, if 0<=m<=n, T(n(n+1)/2 + m, n) = A000712(m). These partitions in this range can be viewed as an equilateral right triangle of side n, with one partition appended on the top (at the left) and another appended on the right. - Franklin T. Adams-Watters, Jan 11 2006
Successive columns approach closer and closer to A000712. - N. J. A. Sloane, Mar 10 2007

Examples

			T(5,2) = 2 because the only partitions of 5 having Durfee square of size 2 are [3,2] and [2,2,1]; the other five partitions ([5], [4,1], [3,1,1], [2,1,1,1] and [1,1,1,1,1]) have Durfee square of size 1.
Triangle starts:
  1;
  2;
  3;
  4,  1;
  5,  2;
  6,  5;
  7,  8;
  8, 14;
  9, 20,  1;
  ...
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (pp. 27-28).
  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004 (pp. 75-78).

Crossrefs

For another version see A115720. Row lengths A000196.

Programs

  • Maple
    g:=sum(t^k*q^(k^2)/product((1-q^j)^2,j=1..k),k=1..40): gser:=series(g,q=0,32): for n from 1 to 27 do P[n]:=coeff(gser,q^n) od: for n from 1 to 27 do seq(coeff(P[n],t^j),j=1..floor(sqrt(n))) od; # yields sequence in triangular form
    # second Maple program:
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    T:= (n, k)-> add(b(m, k)*b(n-k^2-m, k), m=0..n-k^2):
    seq(seq(T(n, k), k=1..floor(sqrt(n))), n=1..30); # Alois P. Heinz, Apr 09 2012
  • Mathematica
    Map[Select[#,#>0&]&,Drop[Transpose[Table[CoefficientList[ Series[x^(n^2)/Product[1-x^i,{i,1,n}]^2,{x,0,nn}],x],{n,1,10}]],1]] //Grid (* Geoffrey Critzer, Sep 27 2013 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; T[n_, k_] := Sum[b[m, k]*b[n-k^2-m, k], {m, 0, n-k^2}]; Table[T[n, k], {n, 1, 30}, {k, 1, Sqrt[n]}] // Flatten (* Jean-François Alcover, Dec 25 2015, after Alois P. Heinz *)

Formula

G.f.: sum(k>=1, t^k*q^(k^2)/product(j=1..k, (1-q^j)^2 ) ).
T(n,k) = Sum_{i=0}^{n-k^2} P*(i,k)*P*(n-k^2-i), where P*(n,k) = P(n+k,k) is the number of partitions of n objects into at most k parts.

Extensions

Edited and verified by Franklin T. Adams-Watters, Mar 11 2006

A364350 Number of strict integer partitions of n such that no part can be written as a nonnegative linear combination of the others.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 3, 2, 3, 3, 5, 3, 6, 5, 7, 6, 9, 7, 11, 10, 14, 12, 16, 15, 20, 17, 24, 22, 27, 29, 32, 30, 41, 36, 49, 45, 50, 52, 65, 63, 70, 77, 80, 83, 104, 98, 107, 116, 126, 134, 152, 148, 162, 180, 196, 195, 227, 227, 238, 272, 271, 293, 333, 325
Offset: 0

Views

Author

Gus Wiseman, Aug 15 2023

Keywords

Comments

A way of writing n as a (presumed nonnegative) linear combination of a finite sequence y is any sequence of pairs (k_i,y_i) such that k_i >= 0 and Sum k_i*y_i = n. For example, the pairs ((3,1),(1,1),(1,1),(0,2)) are a way of writing 5 as a linear combination of (1,1,1,2), namely 5 = 3*1 + 1*1 + 1*1 + 0*2. Of course, there are A000041(n) ways to write n as a linear combination of (1..n).

Examples

			The a(16) = 6 through a(22) = 12 strict partitions:
  (16)     (17)     (18)     (19)     (20)      (21)      (22)
  (9,7)    (9,8)    (10,8)   (10,9)   (11,9)    (12,9)    (13,9)
  (10,6)   (10,7)   (11,7)   (11,8)   (12,8)    (13,8)    (14,8)
  (11,5)   (11,6)   (13,5)   (12,7)   (13,7)    (15,6)    (15,7)
  (13,3)   (12,5)   (14,4)   (13,6)   (14,6)    (16,5)    (16,6)
  (7,5,4)  (13,4)   (7,6,5)  (14,5)   (17,3)    (17,4)    (17,5)
           (14,3)   (8,7,3)  (15,4)   (8,7,5)   (19,2)    (18,4)
           (15,2)            (16,3)   (9,6,5)   (11,10)   (19,3)
           (7,6,4)           (17,2)   (9,7,4)   (8,7,6)   (12,10)
                             (8,6,5)  (11,5,4)  (9,7,5)   (9,7,6)
                             (9,6,4)            (10,7,4)  (9,8,5)
                                                (10,8,3)  (7,6,5,4)
                                                (11,6,4)
                                                (11,7,3)
		

Crossrefs

For sums of subsets instead of combinations of partitions we have A151897.
For sums instead of combinations we have A237667, binary A236912.
For subsets instead of partitions we have A326083, complement A364914.
The complement in strict partitions is A364839, non-strict A364913.
A more strict variation is A364915.
The case of all positive coefficients is A365006.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A108917 counts knapsack partitions, ranks A299702.
A116861 and A364916 count linear combinations of strict partitions.
A323092 (ranks A320340) and A120641 count double-free partitions.
A364912 counts linear combinations of partitions of k.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&And@@Table[combs[#[[k]],Delete[#,k]]=={},{k,Length[#]}]&]],{n,0,15}]
  • Python
    from sympy.utilities.iterables import partitions
    def A364350(n):
        if n <= 1: return 1
        alist, c = [set(tuple(sorted(set(p))) for p in partitions(i)) for i in range(n)], 1
        for p in partitions(n,k=n-1):
            if max(p.values(),default=0)==1:
                s = set(p)
                if not any(set(t).issubset(s-{q}) for q in s for t in alist[q]):
                    c += 1
        return c # Chai Wah Wu, Sep 23 2023

Extensions

More terms and offset corrected by Martin Fuller, Sep 11 2023

A131689 Triangle of numbers T(n,k) = k!*Stirling2(n,k) = A000142(k)*A048993(n,k) read by rows, T(n, k) for 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 6, 6, 0, 1, 14, 36, 24, 0, 1, 30, 150, 240, 120, 0, 1, 62, 540, 1560, 1800, 720, 0, 1, 126, 1806, 8400, 16800, 15120, 5040, 0, 1, 254, 5796, 40824, 126000, 191520, 141120, 40320, 0, 1, 510, 18150, 186480, 834120, 1905120, 2328480, 1451520, 362880
Offset: 0

Views

Author

Philippe Deléham, Sep 14 2007

Keywords

Comments

Triangle T(n,k), 0 <= k <= n, read by rows given by [0,1,0,2,0,3,0,4,0,5,0,6,0,7,0,...] DELTA [1,1,2,2,3,3,4,4,5,5,6,6,...] where DELTA is the operator defined in A084938; another version of A019538.
See also A019538: version with n > 0 and k > 0. - Philippe Deléham, Nov 03 2008
From Peter Bala, Jul 21 2014: (Start)
T(n,k) gives the number of (k-1)-dimensional faces in the interior of the first barycentric subdivision of the standard (n-1)-dimensional simplex. For example, the barycentric subdivision of the 1-simplex is o--o--o, with 1 interior vertex and 2 interior edges, giving T(2,1) = 1 and T(2,2) = 2.
This triangle is used when calculating the face vectors of the barycentric subdivision of a simplicial complex. Let S be an n-dimensional simplicial complex and write f_k for the number of k-dimensional faces of S, with the usual convention that f_(-1) = 1, so that F := (f_(-1), f_0, f_1,...,f_n) is the f-vector of S. If M(n) denotes the square matrix formed from the first n+1 rows and n+1 columns of the present triangle, then the vector F*M(n) is the f-vector of the first barycentric subdivision of the simplicial complex S (Brenti and Welker, Lemma 2.1). For example, the rows of Pascal's triangle A007318 (but with row and column indexing starting at -1) are the f-vectors for the standard n-simplexes. It follows that A007318*A131689, which equals A028246, is the array of f-vectors of the first barycentric subdivision of standard n-simplexes. (End)
This triangle T(n, k) appears in the o.g.f. G(n, x) = Sum_{m>=0} S(n, m)*x^m with S(n, m) = Sum_{j=0..m} j^n for n >= 1 as G(n, x) = Sum_{k=1..n} (x^k/(1 - x)^(k+2))*T(n, k). See also the Eulerian triangle A008292 with a Mar 31 2017 comment for a rewritten form. For the e.g.f. see A028246 with a Mar 13 2017 comment. - Wolfdieter Lang, Mar 31 2017
T(n,k) = the number of alignments of length k of n strings each of length 1. See Slowinski. An example is given below. Cf. A122193 (alignments of strings of length 2) and A299041 (alignments of strings of length 3). - Peter Bala, Feb 04 2018
The row polynomials R(n,x) are the Fubini polynomials. - Emanuele Munarini, Dec 05 2020
From Gus Wiseman, Feb 18 2022: (Start)
Also the number of patterns of length n with k distinct parts (or with maximum part k), where we define a pattern to be a finite sequence covering an initial interval of positive integers. For example, row n = 3 counts the following patterns:
(1,1,1) (1,2,2) (1,2,3)
(2,1,2) (1,3,2)
(2,2,1) (2,1,3)
(1,1,2) (2,3,1)
(1,2,1) (3,1,2)
(2,1,1) (3,2,1)
(End)
Regard A048994 as a lower-triangular matrix and divide each term A048994(n,k) by n!, then this is the matrix inverse. Because Sum_{k=0..n} (A048994(n,k) * x^n / n!) = A007318(x,n), Sum_{k=0..n} (A131689(n,k) * A007318(x,k)) = x^n. - Natalia L. Skirrow, Mar 23 2023
T(n,k) is the number of ordered partitions of [n] into k blocks. - Alois P. Heinz, Feb 21 2025

Examples

			The triangle T(n,k) begins:
  n\k 0 1    2     3      4       5        6        7        8        9      10 ...
  0:  1
  1:  0 1
  2:  0 1    2
  3:  0 1    6     6
  4:  0 1   14    36     24
  5:  0 1   30   150    240     120
  6:  0 1   62   540   1560    1800      720
  7:  0 1  126  1806   8400   16800    15120     5040
  8:  0 1  254  5796  40824  126000   191520   141120    40320
  9:  0 1  510 18150 186480  834120  1905120  2328480  1451520   362880
  10: 0 1 1022 55980 818520 5103000 16435440 29635200 30240000 16329600 3628800
  ... reformatted and extended. - _Wolfdieter Lang_, Mar 31 2017
From _Peter Bala_, Feb 04 2018: (Start)
T(4,2) = 14 alignments of length 2 of 4 strings of length 1. Examples include
  (i) A -    (ii) A -    (iii) A -
      B -         B -          - B
      C -         - C          - C
      - D         - D          - D
There are C(4,1) = 4 alignments of type (i) with a single gap character - in column 1, C(4,2) = 6 alignments of type (ii) with two gap characters in column 1 and C(4,3) = 4 alignments of type (iii) with three gap characters in column 1, giving a total of 4 + 6 + 4 = 14 alignments. (End)
		

Crossrefs

Case m=1 of the polynomials defined in A278073.
Cf. A000142 (diagonal), A000670 (row sums), A000012 (alternating row sums), A210029 (central terms).
Cf. A008292, A028246 (o.g.f. and e.g.f. of sums of powers).
A version for partitions is A116608, or by maximum A008284.
A version for compositions is A235998, or by maximum A048004.
Classes of patterns:
- A000142 = strict
- A005649 = anti-run, complement A069321
- A019536 = necklace
- A032011 = distinct multiplicities
- A060223 = Lyndon
- A226316 = (1,2,3)-avoiding, weakly A052709, complement A335515
- A296975 = aperiodic
- A345194 = alternating, up/down A350354, complement A350252
- A349058 = weakly alternating
- A351200 = distinct runs
- A351292 = distinct run-lengths

Programs

  • Julia
    function T(n, k)
        if k < 0 || k > n return 0 end
        if n == 0 && k == 0 return 1 end
        k*(T(n-1, k-1) + T(n-1, k))
    end
    for n in 0:7
        println([T(n, k) for k in 0:n])
    end
    # Peter Luschny, Mar 26 2020
    
  • Maple
    A131689 := (n,k) -> Stirling2(n,k)*k!: # Peter Luschny, Sep 17 2011
    # Alternatively:
    A131689_row := proc(n) 1/(1-t*(exp(x)-1)); expand(series(%,x,n+1)); n!*coeff(%,x,n); PolynomialTools:-CoefficientList(%,t) end:
    for n from 0 to 9 do A131689_row(n) od; # Peter Luschny, Jan 23 2017
  • Mathematica
    t[n_, k_] := k!*StirlingS2[n, k]; Table[t[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 25 2014 *)
    T[n_, k_] := If[n <= 0 || k <= 0, Boole[n == 0 && k == 0], Sum[(-1)^(i + k) Binomial[k, i] i^(n + k), {i, 0, k}]]; (* Michael Somos, Jul 08 2018 *)
  • PARI
    {T(n, k) = if( n<0, 0, sum(i=0, k, (-1)^(k + i) * binomial(k, i) * i^n))};
    /* Michael Somos, Jul 08 2018 */
    
  • SageMath
    @cached_function
    def F(n): # Fubini polynomial
        R. = PolynomialRing(ZZ)
        if n == 0: return R(1)
        return R(sum(binomial(n, k)*F(n - k)*x for k in (1..n)))
    for n in (0..9): print(F(n).list()) # Peter Luschny, May 21 2021

Formula

T(n,k) = k*(T(n-1,k-1) + T(n-1,k)) with T(0,0)=1. Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A000629(n), A033999(n), A000007(n), A000670(n), A004123(n+1), A032033(n), A094417(n), A094418(n), A094419(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6 respectively. [corrected by Philippe Deléham, Feb 11 2013]
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A000142(n), A000670(n), A122704(n) for x=-1, 0, 1, 2 respectively. - Philippe Deléham, Oct 09 2007
Sum_{k=0..n} (-1)^k*T(n,k)/(k+1) = Bernoulli numbers A027641(n)/A027642(n). - Peter Luschny, Sep 17 2011
G.f.: F(x,t) = 1 + x*t + (x+x^2)*t^2/2! + (x+6*x^2+6*x^3)*t^3/3! + ... = Sum_{n>=0} R(n,x)*t^n/n!.
The row polynomials R(n,x) satisfy the recursion R(n+1,x) = (x+x^2)*R'(n,x) + x*R(n,x) where ' indicates differentiation with respect to x. - Philippe Deléham, Feb 11 2013
T(n,k) = [t^k] (n! [x^n] (1/(1-t*(exp(x)-1)))). - Peter Luschny, Jan 23 2017
The n-th row polynomial has the form x o x o ... o x (n factors), where o denotes the black diamond multiplication operator of Dukes and White. See also Bala, Example E8. - Peter Bala, Jan 08 2018

A364272 Number of strict integer partitions of n containing the sum of some subset of the parts. A variation of sum-full strict partitions.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 3, 1, 4, 3, 8, 6, 11, 10, 17, 16, 26, 25, 39, 39, 54, 60, 82, 84, 116, 126, 160, 177, 222, 242, 302, 337, 402, 453, 542, 601, 722, 803, 936, 1057, 1234, 1373, 1601, 1793, 2056, 2312, 2658, 2950, 3395, 3789, 4281, 4814, 5452, 6048
Offset: 0

Views

Author

Gus Wiseman, Aug 01 2023

Keywords

Comments

First differs from A316402 at a(16) = 11 due to (7,5,3,1).

Examples

			The a(6) = 1 through a(16) = 11 partitions (A=10):
  (321) . (431) . (532)  (5321) (642)  (5431) (743)  (6432)  (853)
                  (541)         (651)  (6421) (752)  (6531)  (862)
                  (4321)        (5421) (7321) (761)  (7431)  (871)
                                (6321)        (5432) (7521)  (6532)
                                              (6431) (9321)  (6541)
                                              (6521) (54321) (7432)
                                              (7421)         (7621)
                                              (8321)         (8431)
                                                             (8521)
                                                             (A321)
                                                             (64321)
		

Crossrefs

The non-strict complement is A237667, ranks A364531.
The non-strict version is A237668, ranks A364532.
The complement in strict partitions is A364349, binary A364533.
The linear combination-free version is A364350.
For subsets of {1..n} we have A364534, complement A151897.
The binary version is A364670, allowing re-used parts A363226.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A108917 counts knapsack partitions, strict A275972, ranks A299702.
A236912 counts binary sum-free partitions, complement A237113.
A323092 counts double-free partitions, ranks A320340.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], UnsameQ@@#&&Intersection[#, Total/@Subsets[#,{2,Length[#]}]]!={}&]],{n,0,30}]

A008483 Number of partitions of n into parts >= 3.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 2, 2, 3, 4, 5, 6, 9, 10, 13, 17, 21, 25, 33, 39, 49, 60, 73, 88, 110, 130, 158, 191, 230, 273, 331, 391, 468, 556, 660, 779, 927, 1087, 1284, 1510, 1775, 2075, 2438, 2842, 3323, 3872, 4510
Offset: 0

Views

Author

T. Forbes (anthony.d.forbes(AT)googlemail.com)

Keywords

Comments

a(0) = 1 because the empty partition vacuously has each part >= 3. - Jason Kimberley, Jan 11 2011
Number of partitions where the largest part occurs at least three times. - Joerg Arndt, Apr 17 2011
By removing a single part of size 3, an A026796 partition of n becomes an A008483 partition of n - 3.
For n >= 3 the sequence counts the isomorphism classes of authentication codes AC(2,n,n) with perfect secrecy and with largest probability 0.5 that an interceptor could deceive with a substituted message. - E. Keith Lloyd (ekl(AT)soton.ac.uk).
For n >= 1, also the number of regular graphs of degree 2. - Mitch Harris, Jun 22 2005
(1 + 0*x + 0*x^2 + x^3 + x^4 + x^5 + 2*x^6 + ...) = (1 + x + 2*x^2 + 3*x^3 + 5*x^4 + ...) * 1 / (1 + x + 2*x^2 + 2*x^3 + 3*x^4 + 3*x^5 + 4*x^6 + 4*x^7 + ...). - Gary W. Adamson, Jun 30 2009
Because the triangle A051031 is symmetric, a(n) is also the number of (n-3)-regular graphs on n vertices. Since the disconnected (n-3)-regular graph with minimum order is 2K_{n-2}, then for n > 4 there are no disconnected (n-3)-regular graphs on n vertices. Therefore for n > 4, a(n) is also the number of connected (n-3)-regular graphs on n vertices. - Jason Kimberley, Oct 05 2009
Number of partitions of n+2 such that 2*(number of parts) is a part. - Clark Kimberling, Feb 27 2014
For n >= 1, a(n) is the number of (1,1)-separable partitions of n, as defined at A239482. For example, the (1,1)-separable partitions of 11 are [10,1], [7,1,2,1], [6,1,3,1], [5,1,4,1], [4,1,2,1,2,1], [3,1,3,1,2,1], so that a(11) = 6. - Clark Kimberling, Mar 21 2014
From Peter Bala, Dec 01 2024: (Start)
Let P(3, n) denote the set of partitions of n into parts k >= 3. Then A000041(n) = (1/2) * Sum_{parts k in all partitions in P(3, n+3)} phi(k), where phi(k) is the Euler totient function (see A000010). For example, with n = 5, there are 3 partitions of n + 3 = 8 into parts greater then 3, namely, 8, 5 + 3 and 4 + 4, and (1/2)*(phi(8) + phi(5) + phi(3) + 2*phi(4)) = 7 = A000041(5). (End)

Crossrefs

Essentially the same sequence as A026796 and A281356.
From Jason Kimberley, Nov 07 2009 and Jan 05 2011 and Feb 03 2011: (Start)
Not necessarily connected simple regular graphs: A005176 (any degree), A051031 (triangular array), specified degree k: A000012 (k=0), A059841 (k=1), this sequence (k=2), A005638 (k=3), A033301 (k=4), A165626 (k=5), A165627 (k=6), A165628 (k=7).
2-regular simple graphs: A179184 (connected), A165652 (disconnected), this sequence (not necessarily connected).
2-regular not necessarily connected graphs without multiple edges [partitions without 2 as a part]: this sequence (no loops allowed [without 1 as a part]), A027336 (loops allowed [parts may be 1]).
Not necessarily connected 2-regular graphs with girth at least g [partitions into parts >= g]: A026807 (triangle); chosen g: A000041 (g=1 -- multigraphs with loops allowed), A002865 (g=2 -- multigraphs with loops forbidden), this sequence (g=3), A008484 (g=4), A185325 (g=5), A185326 (g=6), A185327 (g=7), A185328 (g=8), A185329 (g=9).
Not necessarily connected 2-regular graphs with girth exactly g [partitions with smallest part g]: A026794 (triangle); chosen g: A002865 (g=2), A026796 (g=3), A026797 (g=4), A026798 (g=5), A026799 (g=6), A026800(g=7), A026801 (g=8), A026802 (g=9), A026803 (g=10), ... (End)
Cf. A008284.

Programs

  • Magma
    p := NumberOfPartitions; A008483 :=  func< n | n eq 0 select 1 else n le 2 select 0 else p(n) - p(n-1) - p(n-2) + p(n-3)>; // Jason Kimberley, Jan 11 2011
    
  • Maple
    series(1/product((1-x^i),i=3..50),x,51);
    ZL := [ B,{B=Set(Set(Z, card>=3))}, unlabeled ]: seq(combstruct[count](ZL, size=n), n=0..46); # Zerinvary Lajos, Mar 13 2007
    with(combstruct):ZL2:=[S,{S=Set(Cycle(Z,card>2))}, unlabeled]:seq(count(ZL2,size=n),n=0..46); # Zerinvary Lajos, Sep 24 2007
    with(combstruct):a:=proc(m) [A,{A=Set(Cycle(Z,card>m))},unlabeled]; end: A008483:=a(2):seq(count(A008483,size=n),n=0..46); # Zerinvary Lajos, Oct 02 2007
  • Mathematica
    f[1, 1] = 1; f[n_, k_] := f[n, k] = If[n < 0, 0, If[k > n, 0, If[k == n, 1, f[n, k + 1] + f[n - k, k]]]]; Table[ f[n, 3], {n, 49}] (* Robert G. Wilson v, Jan 31 2011 *)
    Rest[Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, 2*Length[p]]], {n, 50}]]  (* Clark Kimberling, Feb 27 2014 *)
  • PARI
    a(n) = numbpart(n)-numbpart(n-1)-numbpart(n-2)+numbpart(n-3) \\ Charles R Greathouse IV, Jul 19 2011
    
  • Python
    from sympy import partition
    def A008483(n): return partition(n)-partition(n-1)-partition(n-2)+partition(n-3) # Chai Wah Wu, Jun 10 2025

Formula

a(n) = p(n) - p(n - 1) - p(n - 2) + p(n - 3) where p(n) is the number of unrestricted partitions of n into positive parts (A000041).
G.f.: Product_{m>=3} 1/(1-x^m).
G.f.: (Sum_{n>=0} x^(3*n)) / (Product_{k=1..n} (1 - x^k)). - Joerg Arndt, Apr 17 2011
a(n) = A121081(n+3) - A121659(n+3). - Reinhard Zumkeller, Aug 14 2006
Euler transformation of A179184. a(n) = A179184(n) + A165652(n). - Jason Kimberley, Jan 05 2011
a(n) ~ Pi^2 * exp(Pi*sqrt(2*n/3)) / (12*sqrt(3)*n^2). - Vaclav Kotesovec, Feb 26 2015
G.f.: exp(Sum_{k>=1} x^(3*k)/(k*(1 - x^k))). - Ilya Gutkovskiy, Aug 21 2018
a(n) = Sum_{j=0..floor(n/2)} A008284(n-2*j,j). - Gregory L. Simay, Apr 27 2023
G.f.: 1 + Sum_{n >= 1} x^(n+2)/Product_{k = 0..n-1} (1 - x^(k+3)). - Peter Bala, Dec 01 2024

A237113 Number of partitions of n such that some part is a sum of two other parts.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 3, 3, 8, 10, 17, 22, 37, 47, 71, 91, 133, 170, 236, 301, 408, 515, 686, 860, 1119, 1401, 1798, 2232, 2829, 3495, 4378, 5381, 6682, 8165, 10060, 12238, 14958, 18116, 22018, 26533, 32071, 38490, 46265, 55318, 66193, 78843, 93949, 111503, 132326
Offset: 0

Views

Author

Clark Kimberling, Feb 04 2014

Keywords

Comments

These are partitions containing the sum of some 2-element submultiset of the parts, a variation of binary sum-full partitions where parts cannot be re-used, ranked by A364462. The complement is counted by A236912. The non-binary version is A237668. For re-usable parts we have A363225. - Gus Wiseman, Aug 10 2023

Examples

			Of the 11 partitions of 6, only these 3 include a part that is a sum of two other parts: [3,2,1], [2,2,1,1], [2,1,1,1,1].  Thus, a(6) = 3.
From _Gus Wiseman_, Aug 09 2023: (Start)
The a(0) = 0 through a(9) = 10 partitions:
  .  .  .  .  (211)  (2111)  (321)    (3211)    (422)      (3321)
                             (2211)   (22111)   (431)      (4221)
                             (21111)  (211111)  (3221)     (4311)
                                                (4211)     (5211)
                                                (22211)    (32211)
                                                (32111)    (42111)
                                                (221111)   (222111)
                                                (2111111)  (321111)
                                                           (2211111)
                                                           (21111111)
(End)
		

Crossrefs

The complement for subsets is A085489, with re-usable parts A007865.
For subsets of {1..n} we have A088809, with re-usable parts A093971.
The complement is counted by A236912, ranks A364461.
The non-binary complement is A237667, ranks A364531.
The non-binary version is A237668, ranks A364532.
With re-usable parts we have A363225, ranks A364348.
The complement with re-usable parts is A364345, ranks A364347.
These partitions have ranks A364462.
The strict case is A364670, with re-usable parts A363226.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A108917 counts knapsack partitions, ranks A299702.
A323092 counts double-free partitions, ranks A320340.

Programs

  • Mathematica
    z = 20; t = Map[Count[Map[Length[Cases[Map[Total[#] &, Subsets[#, {2}]],  Apply[Alternatives, #]]] &, IntegerPartitions[#]], 0] &, Range[z]] (* A236912 *)
    u = PartitionsP[Range[z]] - t  (* A237113, Peter J. C. Moses, Feb 03 2014 *)
    Table[Length[Select[IntegerPartitions[n],Intersection[#,Total/@Subsets[#,{2}]]!={}&]],{n,0,30}] (* Gus Wiseman, Aug 09 2023 *)

Formula

a(n) = A000041(n) - A236912(n).

Extensions

a(0)=0 prepended by Alois P. Heinz, Sep 17 2023

A359902 Triangle read by rows where T(n,k) is the number of odd-length integer partitions of n with median k.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 2, 1, 0, 0, 1, 2, 2, 0, 0, 0, 1, 4, 2, 1, 0, 0, 0, 1, 4, 3, 2, 0, 0, 0, 0, 1, 7, 4, 3, 1, 0, 0, 0, 0, 1, 8, 6, 3, 2, 0, 0, 0, 0, 0, 1, 12, 8, 4, 3, 1, 0, 0, 0, 0, 0, 1, 14, 11, 5, 4, 2, 0, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 21 2023

Keywords

Comments

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Examples

			Triangle begins:
  1
  0  1
  1  0  1
  1  0  0  1
  2  1  0  0  1
  2  2  0  0  0  1
  4  2  1  0  0  0  1
  4  3  2  0  0  0  0  1
  7  4  3  1  0  0  0  0  1
  8  6  3  2  0  0  0  0  0  1
 12  8  4  3  1  0  0  0  0  0  1
 14 11  5  4  2  0  0  0  0  0  0  1
 21 14  8  4  3  1  0  0  0  0  0  0  1
 24 20 10  5  4  2  0  0  0  0  0  0  0  1
 34 25 15  6  5  3  1  0  0  0  0  0  0  0  1
For example, row n = 9 counts the following partitions:
  (7,1,1)              (5,2,2)      (3,3,3)  (4,4,1)  .  .  .  .  (9)
  (3,3,1,1,1)          (6,2,1)      (4,3,2)
  (4,2,1,1,1)          (2,2,2,2,1)  (5,3,1)
  (5,1,1,1,1)          (3,2,2,1,1)
  (2,2,1,1,1,1,1)
  (3,1,1,1,1,1,1)
  (1,1,1,1,1,1,1,1,1)
		

Crossrefs

Column k=1 is A002865(n-1).
Row sums are A027193 (odd-length ptns), strict A067659.
This is the odd-length case of A359901, with half-steps A359893.
The median statistic is ranked by A360005(n)/2.
First appearances of medians are ranked by A360006, A360007.
A000041 counts partitions, strict A000009.
A058398 counts partitions by mean, see also A008284, A327482.
A067538 counts partitions w/ integer mean, strict A102627, ranked by A316413.
A240219 counts partitions w/ the same mean as median, complement A359894.
A325347 counts partitions w/ integer median, complement A307683.
A326567/A326568 gives mean of prime indices.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OddQ[Length[#]]&&Median[#]==k&]],{n,15},{k,n}]
Previous Showing 51-60 of 579 results. Next