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

A098859 Number of partitions of n into parts each of which is used a different number of times.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 7, 10, 13, 15, 21, 28, 31, 45, 55, 62, 82, 105, 116, 153, 172, 208, 251, 312, 341, 431, 492, 588, 676, 826, 905, 1120, 1249, 1475, 1676, 2003, 2187, 2625, 2922, 3409, 3810, 4481, 4910, 5792, 6382, 7407, 8186, 9527, 10434
Offset: 0

Views

Author

David S. Newman, Oct 11 2004

Keywords

Comments

Fill, Janson and Ward refer to these partitions as Wilf partitions. - Peter Luschny, Jun 04 2012

Examples

			a(6)=7 because 6= 4+1+1= 3+3= 3+1+1+1= 2+2+2= 2+1+1+1+1= 1+1+1+1+1+1. Four unrestricted partitions of 6 are not counted by a(6): 5+1, 4+2, 3+2+1 because at least two different summands are each used once; 2+2+1+1 because each summand is used twice.
From _Gus Wiseman_, Apr 19 2019: (Start)
The a(1) = 1 through a(9) = 15 partitions are the following. The Heinz numbers of these partitions are given by A130091.
  1   2    3     4      5       6        7         8          9
      11   111   22     221     33       322       44         333
                 211    311     222      331       332        441
                 1111   2111    411      511       422        522
                        11111   3111     2221      611        711
                                21111    4111      2222       3222
                                111111   22111     5111       6111
                                         31111     22211      22221
                                         211111    41111      33111
                                         1111111   221111     51111
                                                   311111     411111
                                                   2111111    2211111
                                                   11111111   3111111
                                                              21111111
                                                              111111111
(End)
		

Crossrefs

Row sums of A182485.
Cf. A047966 (each part appears the same number of times), A090858, A116608, A130091, A325242.

Programs

Formula

log(a(n)) ~ N*log(N) where N = (6*n)^(1/3) (see Fill, Janson and Ward). - Peter Luschny, Jun 04 2012

Extensions

Corrected and extended by Vladeta Jovovic, Oct 22 2004

A325280 Triangle read by rows where T(n,k) is the number of integer partitions of n with adjusted frequency depth k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 1, 2, 3, 0, 0, 1, 3, 4, 3, 0, 0, 0, 1, 1, 4, 8, 1, 0, 0, 0, 1, 3, 6, 9, 3, 0, 0, 0, 0, 1, 2, 8, 12, 7, 0, 0, 0, 0, 0, 1, 3, 11, 17, 10, 0, 0, 0, 0, 0, 0, 1, 1, 11, 26, 17, 0, 0, 0, 0, 0, 0, 0, 1, 5, 19, 25, 27
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

The adjusted frequency depth of an integer partition is 0 if the partition is empty, and otherwise it is one plus the number of times one must take the multiset of multiplicities to reach a singleton. For example, the partition (32211) has adjusted frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2).
The term "frequency depth" appears to have been coined by Clark Kimberling in A225485 and A225486, and can be applied to both integers (A323014) and integer partitions (this sequence).

Examples

			Triangle begins:
  1
  0  1
  0  1  1
  0  1  1  1
  0  1  2  1  1
  0  1  1  2  3  0
  0  1  3  4  3  0  0
  0  1  1  4  8  1  0  0
  0  1  3  6  9  3  0  0  0
  0  1  2  8 12  7  0  0  0  0
  0  1  3 11 17 10  0  0  0  0  0
  0  1  1 11 26 17  0  0  0  0  0  0
  0  1  5 19 25 27  0  0  0  0  0  0  0
  0  1  1 17 44 38  0  0  0  0  0  0  0  0
  0  1  3 25 53 52  1  0  0  0  0  0  0  0  0
  0  1  3 29 63 76  4  0  0  0  0  0  0  0  0  0
  0  1  4 37 83 98  8  0  0  0  0  0  0  0  0  0  0
Row n = 9 counts the following partitions:
  (9)  (333)        (54)      (441)       (3321)
       (111111111)  (63)      (522)       (4221)
                    (72)      (711)       (4311)
                    (81)      (3222)      (5211)
                    (432)     (6111)      (32211)
                    (531)     (22221)     (42111)
                    (621)     (33111)     (321111)
                    (222111)  (51111)
                              (411111)
                              (2211111)
                              (3111111)
                              (21111111)
		

Crossrefs

Row sums are A000041. Column k = 2 is A032741. Column k = 3 is A325245.
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or this sequence (length/frequency depth).

Programs

  • Mathematica
    fdadj[ptn_List]:=If[ptn=={},0,Length[NestWhileList[Sort[Length/@Split[#]]&,ptn,Length[#]>1&]]];
    Table[Length[Select[IntegerPartitions[n],fdadj[#]==k&]],{n,0,16},{k,0,n}]
  • PARI
    \\ depth(p) gives adjusted frequency depth of partition.
    depth(p)={if(!#p, 0, my(r=1); while(#p > 1, my(L=List(), k=0); for(i=1, #p, if(i==#p||p[i]<>p[i+1], listput(L,i-k); k=i)); listsort(L); p=L; r++); r)}
    row(n)={my(v=vector(1+n)); forpart(p=n, v[1+depth(Vec(p))]++); v}
    { for(n=0, 10, print(row(n))) } \\ Andrew Howroyd, Jan 18 2023

A336866 Number of integer partitions of n without all distinct multiplicities.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 4, 5, 9, 15, 21, 28, 46, 56, 80, 114, 149, 192, 269, 337, 455, 584, 751, 943, 1234, 1527, 1944, 2422, 3042, 3739, 4699, 5722, 7100, 8668, 10634, 12880, 15790, 19012, 23093, 27776, 33528, 40102, 48264, 57469, 68793, 81727, 97372, 115227
Offset: 0

Views

Author

Gus Wiseman, Aug 09 2020

Keywords

Examples

			The a(0) = 0 through a(9) = 15 partitions (empty columns shown as dots):
  .  .  .  (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)
                                                      (222111)
                                                      (321111)
		

Crossrefs

A098859 counts the complement.
A130092 gives the Heinz numbers of these partitions.
A001222 counts prime factors with multiplicity.
A013929 lists nonsquarefree numbers.
A047966 counts uniform partitions.
A047967 counts non-strict partitions.
A071625 counts distinct prime multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A327498 gives the maximum divisor with distinct prime multiplicities.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!UnsameQ@@Length/@Split[#]&]],{n,0,30}]

Formula

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

A225485 Number of partitions of n that have frequency depth k, an array read by rows.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 3, 4, 3, 1, 1, 4, 8, 1, 1, 3, 6, 9, 3, 1, 2, 8, 12, 7, 1, 3, 11, 17, 10, 1, 1, 11, 26, 17, 1, 5, 19, 25, 27, 1, 1, 17, 44, 38, 1, 3, 25, 53, 52, 1, 1, 3, 29, 63, 76, 4
Offset: 1

Views

Author

Clark Kimberling, May 08 2013

Keywords

Comments

Let S = {x(1),...,x(k)} be a multiset whose distinct elements are y(1),...,y(h). Let f(i) be the frequency of y(i) in S. Define F(S) = {f(1),..,f(h)}, F(1,S) = F(S), and F(m,S) = F(F(m-1),S) for m>1. Then lim(F(m,S)) = {1} for every S, so that there is a least positive integer i for which F(i,S) = {1}, which we call the frequency depth of S.
Equivalently, the frequency depth of an integer partition is the number of times one must take the multiset of multiplicities to reach (1). For example, the partition (32211) has frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2) -> (1). - Gus Wiseman, Apr 19 2019
From Clark Kimberling, Sep 26 2023: (Start)
Below, m^n abbreviates the sum m+...+m of n terms. In the following list, the numbers p_1,...,p_k are distinct, m >= 1, and k >= 1. The forms of the partitions being counted are as follows:
column 1: [n],
column 2: [m^k],
column 3: [p_1^m,...,p_k^m],
column 4: [(p_1^m_1)^m,..., (p_k^m_k)^m], distinct numbers m_i.
Column 3 is of special interest. Assume first that m = 1, so that the form of partition being counted is p = [p_1,...,p_k], with conjugate given by [q_1,...,q_m] where q_i is the number of parts of p that are >= i. Since the p_i are distinct, the distinct parts of q are the integers 1,2,...,k. For the general case that m >= 1, the distinct parts of q are the integers m,...,km. Let S(n) denote the set of partitions of n counted by column 3. Then if a and b are in the set S*(n) of conjugates of partitions in S(n), and if a > b, then a - b is also in S*(n). Call this the subtraction property. Conversely, if a partition q has the subtraction property, then q must consist of a set of numbers m,..,km for some m. Thus, column 3 counts the partitions of n that have the subtraction property. (End)

Examples

			The first 9 rows:
  n = 1 .... 0
  n = 2 .... 1..1
  n = 3 .... 1..1..1
  n = 4 .... 1..2..1..1
  n = 5 .... 1..1..2..3
  n = 6 .... 1..3..4..3
  n = 7 .... 1..1..4..8..1
  n = 8 .... 1..3..6..9..3
  n = 9 .... 1..2..8.12..7
For the 7 partitions of 5, successive frequencies are shown here:
  5 -> 1 (depth 1)
  41 -> 11 -> 2 -> 1 (depth 3)
  32 -> 11 -> 2 -> 1 (depth 3)
  311 -> 12 -> 11 -> 2 -> 1 (depth 4)
  221 -> 12 -> 11 -> 2 -> 1 (depth 4)
  2111 -> 13 -> 11 -> 2 -> 1 (depth 4)
  11111 -> 5 -> 1 (depth 2)
Summary: 1 partition has depth 1; 1 has depth 2; 2 have 3; and 3 have 4, so that the row for n = 5 is 1..1..2..3 .
		

Crossrefs

Row sums are A000041.
Column k = 2 is A032741.
Column k = 3 is A325245.
a(n!) = A325272(n).
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or A325280 (length/frequency depth).

Programs

  • Mathematica
    c[s_] := c[s] = Select[Table[Count[s, i], {i, 1, Max[s]}], # > 0 &]
    f[s_] := f[s] = Drop[FixedPointList[c, s], -2]
    t[s_] := t[s] = Length[f[s]]
    u[n_] := u[n] = Table[t[Part[IntegerPartitions[n], i]],
      {i, 1, Length[IntegerPartitions[n]]}];
    Flatten[Table[Count[u[n], k], {n, 2, 25}, {k, 1, Max[u[n]]}]]

A183558 Number of partitions of n containing a clique of size 1.

Original entry on oeis.org

0, 1, 1, 2, 3, 6, 7, 13, 16, 25, 33, 49, 61, 90, 113, 156, 198, 269, 334, 448, 556, 726, 902, 1163, 1428, 1827, 2237, 2817, 3443, 4302, 5219, 6478, 7833, 9632, 11616, 14197, 17031, 20712, 24769, 29925, 35688, 42920, 50980, 61059, 72318, 86206, 101837, 120941
Offset: 0

Views

Author

Alois P. Heinz, Jan 05 2011

Keywords

Comments

All parts of a number partition with the same value form a clique. The size of a clique is the number of elements in the clique.

Examples

			a(5) = 6, because 6 partitions of 5 contain (at least) one clique of size 1: [1,1,1,2], [1,2,2], [1,1,3], [2,3], [1,4], [5].
From _Gus Wiseman_, Apr 19 2019: (Start)
The a(1) = 1 through a(8) = 16 partitions are the following. The Heinz numbers of these partitions are given by A052485 (weak numbers).
  (1)  (2)  (3)   (4)    (5)     (6)      (7)       (8)
            (21)  (31)   (32)    (42)     (43)      (53)
                  (211)  (41)    (51)     (52)      (62)
                         (221)   (321)    (61)      (71)
                         (311)   (411)    (322)     (332)
                         (2111)  (3111)   (331)     (422)
                                 (21111)  (421)     (431)
                                          (511)     (521)
                                          (2221)    (611)
                                          (3211)    (3221)
                                          (4111)    (4211)
                                          (31111)   (5111)
                                          (211111)  (32111)
                                                    (41111)
                                                    (311111)
                                                    (2111111)
(End)
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0], `if`(i<1, [0, 0],
          add((l->`if`(j=1, [l[1]$2], l))(b(n-i*j, i-1)), j=0..n/i)))
        end:
    a:= n-> b(n$2)[2]:
    seq(a(n), n=0..50);
  • Mathematica
    max = 50; f = (1 - Product[1 - x^j + x^(2*j), {j, 1, max}])/Product[1 - x^j, {j, 1, max}]; s = Series[f, {x, 0, max}]; CoefficientList[s, x] (* Jean-François Alcover, Oct 01 2014. Edited by Gus Wiseman, Apr 19 2019 *)

Formula

G.f.: (1-Product_{j>0} (1-x^(j)+x^(2*j))) / (Product_{j>0} (1-x^j)).
From Vaclav Kotesovec, Nov 15 2016: (Start)
a(n) = A000041(n) - A007690(n).
a(n) ~ exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*n). (End)

Extensions

a(0)=0 prepended by Gus Wiseman, Apr 19 2019

A175804 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the n-th term of the k-th differences of partition numbers A000041.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, -1, 0, 1, 3, 2, 1, 1, 2, 5, -4, -2, -1, 0, 2, 7, 9, 5, 3, 2, 2, 4, 11, -21, -12, -7, -4, -2, 0, 4, 15, 49, 28, 16, 9, 5, 3, 3, 7, 22, -112, -63, -35, -19, -10, -5, -2, 1, 8, 30, 249, 137, 74, 39, 20, 10, 5, 3, 4, 12, 42, -539, -290, -153, -79, -40, -20, -10, -5, -2, 2, 14, 56
Offset: 0

Views

Author

Alois P. Heinz, Dec 04 2010

Keywords

Comments

Odlyzko showed that the k-th differences of A000041(n) alternate in sign with increasing n up to a certain index n_0(k) and then stay positive.
Are there any zeros after the first four, which all lie in columns k = 1, 2? - Gus Wiseman, Dec 15 2024

Examples

			Square array A(n,k) begins:
   1,  0,  1, -1,  2,  -4,   9,  ...
   1,  1,  0,  1, -2,   5, -12,  ...
   2,  1,  1, -1,  3,  -7,  16,  ...
   3,  2,  0,  2, -4,   9, -19,  ...
   5,  2,  2, -2,  5, -10,  20,  ...
   7,  4,  0,  3, -5,  10, -20,  ...
  11,  4,  3, -2,  5, -10,  22,  ...
		

Crossrefs

Columns k=0-5 give: A000041, A002865, A053445, A072380, A081094, A081095.
Main diagonal gives A379378.
For primes we have A095195 or A376682.
Row n = 0 is A281425.
Row n = 1 is A320590 except first term.
For composites we have A377033.
For squarefree numbers we have A377038.
For nonsquarefree numbers we have A377046.
For prime powers we have A377051.
Antidiagonal sums are A377056, absolute value version A378621.
The version for strict partitions is A378622, first column A293467.
A000009 counts strict integer partitions, differences A087897, A378972.

Programs

  • Maple
    A41:= combinat[numbpart]:
    DD:= proc(p) proc(n) option remember; p(n+1) -p(n) end end:
    A:= (n,k)-> (DD@@k)(A41)(n):
    seq(seq(A(n, d-n), n=0..d), d=0..11);
  • Mathematica
    max = 11; a41 = Array[PartitionsP, max+1, 0]; a[n_, k_] := Differences[a41, k][[n+1]]; Table[a[n, k-n], {k, 0, max}, {n, 0, k}] // Flatten (* Jean-François Alcover, Aug 29 2014 *)
    nn=5;Table[Table[Sum[(-1)^(k-i)*Binomial[k,i]*PartitionsP[n+i],{i,0,k}],{k,0,nn}],{n,0,nn}] (* Gus Wiseman, Dec 15 2024 *)

Formula

A(n,k) = (Delta^(k) A000041)(n).
A(n,k) = Sum_{i=0..k} (-1)^(k-i) * binomial(k,i) * A000041(n+i). In words, row x is the inverse zero-based binomial transform of A000041 shifted left x times. - Gus Wiseman, Dec 15 2024

A325268 Triangle read by rows where T(n,k) is the number of integer partitions of n with omicron k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 3, 0, 1, 0, 1, 5, 0, 0, 1, 0, 1, 7, 2, 0, 0, 1, 0, 1, 12, 1, 0, 0, 0, 1, 0, 1, 17, 2, 1, 0, 0, 0, 1, 0, 1, 24, 4, 0, 0, 0, 0, 0, 1, 0, 1, 33, 5, 1, 1, 0, 0, 0, 0, 1, 0, 1, 44, 9, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 57, 14, 3, 0, 1
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

The omega-sequence of an integer partition is the sequence of lengths of the multisets obtained by repeatedly taking the multiset of multiplicities until a singleton is reached. The omicron of the partition is 0 if the omega-sequence is empty, 1 if it is a singleton, and otherwise the second-to-last part. For example, the partition (32211) has chain of multisets of multiplicities {1,1,2,2,3} -> {1,2,2} -> {1,2} -> {1,1} -> {2}, so its omega-sequence is (5,3,2,2,1), and its omicron is 2.

Examples

			Triangle begins:
  1
  0  1
  0  1  1
  0  1  1  1
  0  1  3  0  1
  0  1  5  0  0  1
  0  1  7  2  0  0  1
  0  1 12  1  0  0  0  1
  0  1 17  2  1  0  0  0  1
  0  1 24  4  0  0  0  0  0  1
  0  1 33  5  1  1  0  0  0  0  1
  0  1 44  9  1  0  0  0  0  0  0  1
  0  1 57 14  3  0  1  0  0  0  0  0  1
  0  1 76 20  3  0  0  0  0  0  0  0  0  1
Row n = 8 counts the following partitions.
  (8)  (44)       (431)  (2222)  (11111111)
       (53)       (521)
       (62)
       (71)
       (332)
       (422)
       (611)
       (3221)
       (3311)
       (4211)
       (5111)
       (22211)
       (32111)
       (41111)
       (221111)
       (311111)
       (2111111)
		

Crossrefs

Row sums are A000041. Column k = 2 is A325267.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or A325280 (length/frequency depth).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Switch[#,{},0,{},1,,NestWhile[Sort[Length/@Split[#]]&,#,Length[#]>1&]//First]==k&]],{n,0,10},{k,0,n}]
  • PARI
    omicron(p)={if(!#p, 0, my(r=1); while(#p > 1, my(L=List(), k=0); r=#p; for(i=1, #p, if(i==#p||p[i]<>p[i+1], listput(L,i-k); k=i)); listsort(L); p=L); r)}
    row(n)={my(v=vector(1+n)); forpart(p=n, v[1 + omicron(Vec(p))]++); v}
    { for(n=0, 10, print(row(n))) } \\ Andrew Howroyd, Jan 18 2023

A225486 Maximal frequency depth for the partitions of n.

Original entry on oeis.org

0, 2, 3, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
Offset: 1

Views

Author

Clark Kimberling, May 08 2013

Keywords

Comments

See A225485 for the definition of frequency depth.
The frequency depth of an integer partition is the number of times one must take the multiset of multiplicities to reach (1). For example, the partition (32211) has frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2) -> (1). Differs from A325282 at a(0) and a(1). - Gus Wiseman, Apr 19 2019

Examples

			(See A225485.)
		

Crossrefs

Run lengths are A325258, i.e., first differences of Levine's sequence A011784.

Programs

  • Mathematica
    c[s_] := c[s] = Select[Table[Count[s, i], {i, 1, Max[s]}], # > 0 &]
    f[s_] := f[s] = Drop[FixedPointList[c, s], -2]
    t[s_] := t[s] = Length[f[s]]
    u[n_] := u[n] = Table[t[Part[IntegerPartitions[n], k]],
        {k, 1, Length[IntegerPartitions[n]]}];
    Prepend[Table[Max[u[n]], {n, 2, 10}], 0]
    (* second program *)
    grw[q_]:=Join@@Table[ConstantArray[i,q[[Length[q]-i+1]]],{i,Length[q]}];
    Join@@MapIndexed[ConstantArray[#2[[1]]-1,#1]&,Length[#]-Last[#]&/@NestList[grw,{1,1},6]] (* Gus Wiseman, Apr 19 2019 *)

Formula

a(n) = number of terms in row n of the array in A225485, for n > 0.

Extensions

More terms from Gus Wiseman, Apr 19 2019

A056556 First tetrahedral coordinate; repeat m (m+1)*(m+2)/2 times.

Original entry on oeis.org

0, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
Offset: 0

Views

Author

Henry Bottomley, Jun 26 2000

Keywords

Comments

If {(X,Y,Z)} are triples of nonnegative integers with X >= Y >= Z ordered by X, Y and Z, then X=A056556(n), Y=A056557(n) and Z=A056558(n).
From Gus Wiseman, Jul 03 2019: (Start)
Also the maximum number of distinct multiplicities among integer partitions of n. For example, random partitions of 56 realizing each number of distinct multiplicities are:
1: (24,17,6,5,3,1)
2: (10,9,9,5,5,4,4,3,3,2,1,1)
3: (6,5,5,5,4,4,4,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
4: (28,5,5,3,3,3,2,2,1,1,1,1,1)
5: (13,4,4,4,4,4,3,3,3,2,2,2,2,2,2,1,1)
6: (6,5,5,4,4,4,3,3,3,3,2,2,2,2,2,1,1,1,1,1,1)
The maximum number of distinct multiplicities is 6, so a(56) = 6.
(End)

Examples

			3 is (3+1) * (3+2)/2 = 10 times in the sequence all these occurrences are in consecutive places. The first 3 is at position binomial(3 + 2, 3) = 10, the last one at binomial((3 + 1) + 2, 3) - 1. - _David A. Corneth_, Oct 14 2022
		

Crossrefs

Programs

  • Mathematica
    Table[Table[m, {(m+1)(m+2)/2}], {m, 0, 7}] // Flatten (* Jean-François Alcover, Feb 28 2019 *)
  • PARI
    a(n)=my(t=polrootsreal(x^3+3*x^2+2*x-6*n)); t[#t]\1 \\ Charles R Greathouse IV, Feb 22 2017
    
  • Python
    from math import comb
    from sympy import integer_nthroot
    def A056556(n): return (m:=integer_nthroot(6*(n+1),3)[0])-(nChai Wah Wu, Nov 04 2024

Formula

a(n) = floor(x) where x is the (largest real) solution to x^3 + 3x^2 + 2x - 6n = 0; a(A000292(n)) = n+1.
a(n+1) = a(n)+1 if a(n) = A056558(n), otherwise a(n). - Graeme McRae, Jan 09 2007
a(n) = floor(t/3 + 1/t - 1), where t = (81*n + 3*sqrt(729*n^2 - 3))^(1/3). - Ridouane Oudra, Mar 21 2021
a(n) = floor(t + 1/(3*t) - 1), where t = (6*n)^(1/3), for n>=1. - Ridouane Oudra, Nov 04 2022
a(n) = m if n>=binomial(m+2,3) and a(n) = m-1 otherwise where m = floor((6n+6)^(1/3)). - Chai Wah Wu, Nov 04 2024

Extensions

Incorrect formula deleted by Ridouane Oudra, Nov 04 2022

A117571 Expansion of (1+2*x^2)/((1-x)*(1-x^3)).

Original entry on oeis.org

1, 1, 3, 4, 4, 6, 7, 7, 9, 10, 10, 12, 13, 13, 15, 16, 16, 18, 19, 19, 21, 22, 22, 24, 25, 25, 27, 28, 28, 30, 31, 31, 33, 34, 34, 36, 37, 37, 39, 40, 40, 42, 43, 43, 45, 46, 46, 48, 49, 49, 51, 52, 52, 54, 55, 55, 57, 58, 58, 60, 61, 61, 63, 64, 64, 66, 67, 67, 69, 70, 70, 72
Offset: 0

Views

Author

Paul Barry, Mar 29 2006

Keywords

Comments

Row sums of A116948.
Place n+2 equally-spaced points around a circle, labeled 0,1,2,...,n+1. For each i = 0..n+1 such that 2i != i mod n+2, draw an (undirected) chord from i to 2i mod n+2. Then a(n) is the number of distinct chords. - Kival Ngaokrajang, May 13 2016 (Edited by N. J. A. Sloane, Jun 23 2016)
From Gus Wiseman, Apr 19 2019: (Start)
Also the number of integer partitions of n + 2 with 1 fewer distinct multiplicities than (not necessarily distinct) parts. These are partitions of the form (x,x), (x,y), (x,x,y), or (x,y,y). For example, the a(0) = 1 through a(8) = 9 partitions are the following. The Heinz numbers of these partitions are given by A325270.
(11) (21) (22) (32) (33) (43) (44) (54) (55)
(31) (41) (42) (52) (53) (63) (64)
(211) (221) (51) (61) (62) (72) (73)
(311) (411) (322) (71) (81) (82)
(331) (332) (441) (91)
(511) (422) (522) (433)
(611) (711) (442)
(622)
(811)
(End)

Crossrefs

Programs

  • Magma
    [1 + Floor(2*n/3) + Floor((n+1)/3) : n in [0..100]]; // Wesley Ivan Hurt, Jul 25 2016
  • Maple
    A117571:=n->1 + floor(2*n/3) + floor((n+1)/3): seq(A117571(n), n=0..100); # Wesley Ivan Hurt, Jul 25 2016
  • Mathematica
    CoefficientList[Series[(1 + 2 x^2)/((1 - x) (1 - x^3)), {x, 0, 71}], x] (* Michael De Vlieger, May 13 2016 *)

Formula

G.f.: (1+2*x^2)/((1-x)*(1-x^3)).
a(n) = a(n-1) + a(n-3) - a(n-4) for n>3.
a(n) = cos(2*Pi*n/3+Pi/6)/sqrt(3)-sin(2*Pi*n/3+Pi/6)/3+(3n+2)/3.
a(n) = Sum_{k=0..n} 2*A001045(L((n-k+2)/3)) where L(j/p) is the Legendre symbol of j and p.
a(n) = 1 + floor((n+1)/3) + floor(2*n/3). - Wesley Ivan Hurt, Jul 25 2016
a(n) = n+sign((n-1) mod 3). - Wesley Ivan Hurt, Sep 25 2017
Showing 1-10 of 39 results. Next