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

A026823 Essentially same as A025148.

Original entry on oeis.org

0, 1, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 15, 17, 20, 24, 28, 32, 38
Offset: 1

Views

Author

Keywords

A025147 Number of partitions of n into distinct parts >= 2.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 2, 3, 3, 5, 5, 7, 8, 10, 12, 15, 17, 21, 25, 29, 35, 41, 48, 56, 66, 76, 89, 103, 119, 137, 159, 181, 209, 239, 273, 312, 356, 404, 460, 522, 591, 669, 757, 853, 963, 1085, 1219, 1371, 1539, 1725, 1933, 2164, 2418, 2702, 3016, 3362, 3746, 4171, 4637, 5155
Offset: 0

Views

Author

Keywords

Comments

From R. J. Mathar, Jul 31 2008: (Start)
These "partitions of n into distinct parts >= k" and "partitions of n into distinct parts, the least being k-1" come in pairs of similar, almost shifted but not identical, sequences:
The distinction in the definitions is that "distinct parts >= k" sets a lower bound to all parts, whereas "the least being ..." means that the lower limit must be attained by one of the parts. (End)
From N. J. A. Sloane, Sep 28 2008: (Start)
Generating functions and Maple programs for the sequences in the first and second columns of the above list are respectively:
For A025147, A025148, etc.:
f:=proc(k) product(1+x^j, j=k..100): series(%,x,100): seriestolist(%); end;
For A096765, A096749, etc.:
g:=proc(k) x^(k-1)*product(1+x^j, j=k..100): series(%,x,100): seriestolist(%); end; (End)
Also number of partitions of n+1 into distinct parts, the least being 1.
Number of different sums from 1+[1,3]+[1,4]+...+[1,n]. - Jon Perry, Jan 01 2004
Also number of partitions of n such that if k is the largest part, then all parts from 1 to k occur, k occurring at least twice. Example: a(7)=3 because we have [2,2,2,1],[2,2,1,1,1] and [1,1,1,1,1,1,1]. - Emeric Deutsch, Apr 09 2006
Also number of partitions of n+1 such that if k is the largest part, then all parts from 1 to k occur, k occurring exactly once. Example: a(7)=3 because we have [3,2,2,1],[3,2,1,1,1] and [2,1,1,1,1,1,1] (there is a simple bijection with the partitions defined before). - Emeric Deutsch, Apr 09 2006
Also number of partitions of n+1 into distinct parts where the number of parts is itself a part. - Reinhard Zumkeller, Nov 04 2007
Partial sums give A038348 (observed by Jonathan Vos Post, proved by several correspondents).
Trivially, number of partitions of n into distinct parts (as ascending lists) such that the first part is not 1, the second not 2, the third not 3, etc., see example. - Joerg Arndt, Jun 10 2013
Convolution with A033999 gives A270144 (apart from the offset). - R. J. Mathar, Jun 18 2016

Examples

			a(7) = 3, from {{3, 4}, {2, 5}, {7}}
From _Joerg Arndt_, Jun 10 2013: (Start)
There are a(17) = 21 partitions of 17 into distinct parts >=2:
01:  [ 2 3 4 8 ]
02:  [ 2 3 5 7 ]
03:  [ 2 3 12 ]
04:  [ 2 4 5 6 ]
05:  [ 2 4 11 ]
06:  [ 2 5 10 ]
07:  [ 2 6 9 ]
08:  [ 2 7 8 ]
09:  [ 2 15 ]
10:  [ 3 4 10 ]
11:  [ 3 5 9 ]
12:  [ 3 6 8 ]
13:  [ 3 14 ]
14:  [ 4 5 8 ]
15:  [ 4 6 7 ]
16:  [ 4 13 ]
17:  [ 5 12 ]
18:  [ 6 11 ]
19:  [ 7 10 ]
20:  [ 8 9 ]
21:  [ 17 ]
(End)
		

References

  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education, Vol. 31, No. 1, pp. 24-28, Winter 1997. MathEduc Database (Zentralblatt MATH, 1997c.01891).
  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem II, Missouri Journal of Mathematical Sciences, Vol. 16, No. 1, Winter 2004, pp. 12-17. Zentralblatt MATH, Zbl 1071.05501.

Crossrefs

Programs

  • Haskell
    a025147 = p 2 where
       p _ 0 = 1
       p k m = if m < k then 0 else p (k + 1) (m - k) + p (k + 1) m
    -- Reinhard Zumkeller, Dec 28 2011
    
  • Maple
    g:=product(1+x^j,j=2..65): gser:=series(g,x=0,62): seq(coeff(gser,x,n),n=0..57); # Emeric Deutsch, Apr 09 2006
    with(combstruct):ZL := {L = PowerSet(Sequence(Z,card>=2)) },unlabeled:seq(count([L,ZL],size=i),i=0..57); # Zerinvary Lajos, Mar 09 2007
  • Mathematica
    CoefficientList[Series[Product[1+q^n, {n, 2, 60}], {q, 0, 60}], q]
    FoldList[ PartitionsQ[ #2+1 ]-#1&, 0, Range[ 64 ] ]
    (* also *)
    d[n_] := Select[IntegerPartitions[n], Max[Length /@ Split@#] == 1 && Min[#] >= 2 &]; Table[d[n], {n, 12}] (* strict partitions, parts >= 2 *)
    Table[Length[d[n]], {n, 40}] (* A025147 for n >= 1 *)
    (* Clark Kimberling, Mar 07 2014 *)
    p[, 0] = 1; p[k, m_] := p[k, m] = If[m < k, 0, p[k+1, m-k] + p[k+1, m]]; Table[p[2, m], {m, 0, 59}] (* Jean-François Alcover, Apr 17 2014, after Reinhard Zumkeller *)
  • PARI
    a(n)=if(n,my(v=partitions(n));sum(i=1,#v,v[i][1]>1&&v[i]==vecsort(v[i],,8)),1) \\ Charles R Greathouse IV, Nov 20 2012

Formula

G.f.: Product_{k>=2} (1+x^k).
a(n) = A000009(n)-a(n-1) = Sum_{0<=k<=n} (-1)^k*A000009(n-k). - Henry Bottomley, May 09 2002
a(n)=t(n, 1), where t(n, k)=1+Sum_{i>j>k and i+j=n} t(i, j), 2<=k<=n. - Reinhard Zumkeller, Jan 01 2003
G.f.: 1 + Sum_{k=1..infinity} (x^(k*(k+3)/2) / Product_{j=1..k} (1-x^j)). - Emeric Deutsch, Apr 09 2006
The previous g.f. is a special case of the g.f. for partitions into distinct parts >= L, Sum_{n>=0} ( x^(n*(n+2*L-1)/2) / Product_{k=1..n} (1-x^k) ). - Joerg Arndt, Mar 24 2011
G.f.: Sum_{n>=1} ( x^(n*(n+1)/2-1) / Product_{k=1..n-1} (1-x^k) ), a special case of the g.f. for partitions into distinct parts >= L, Sum_{n>=L-1} ( x^(n*(n+1)/2-L*(L-1)/2) / Product_{k=1..n-(L-1)} (1-x^k) ). - Joerg Arndt, Mar 27 2011
a(n) = Sum_{1A060016(n-k+1,k-1), for n>0. - Reinhard Zumkeller, Nov 04 2007
a(n) = A096765(n+1). - R. J. Mathar, Jul 31 2008
From Vaclav Kotesovec, Aug 16 2015: (Start)
a(n) ~ 1/2 * A000009(n).
a(n) ~ exp(Pi*sqrt(n/3)) / (8*3^(1/4)*n^(3/4)).
(End)

Extensions

Corrected and extended by Dean Hickerson, Oct 10 2001

A087897 Number of partitions of n into odd parts greater than 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 5, 6, 8, 8, 10, 12, 13, 15, 18, 20, 23, 27, 30, 34, 40, 44, 50, 58, 64, 73, 83, 92, 104, 118, 131, 147, 166, 184, 206, 232, 256, 286, 320, 354, 394, 439, 485, 538, 598, 660, 730, 809, 891, 984, 1088, 1196, 1318, 1454, 1596, 1756
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2003

Keywords

Comments

Also number of partitions of n into distinct parts which are not powers of 2.
Also number of partitions of n into distinct parts such that the two largest parts differ by 1.
Also number of partitions of n such that the largest part occurs an odd number of times that is at least 3 and every other part occurs an even number of times. Example: a(10) = 2 because we have [2,2,2,1,1,1,1] and [2,2,2,2,2]. - Emeric Deutsch, Mar 30 2006
Also difference between number of partitions of 1+n into distinct parts and number of partitions of n into distinct parts. - Philippe LALLOUET, May 08 2007
In the Berndt reference replace {a -> -x, q -> x} in equation (3.1) to get f(x). G.f. is 1 - x * (1 - f(x)).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Also number of symmetric unimodal compositions of n+3 where the maximal part appears three times. - Joerg Arndt, Jun 11 2013
Let c(n) = number of palindromic partitions of n whose greatest part has multiplicity 3; then c(n) = a(n-3) for n>=3. - Clark Kimberling, Mar 05 2014
From Gus Wiseman, Aug 22 2021: (Start)
Also the number of integer partitions of n - 1 whose parts cover an interval of positive integers starting with 2. These partitions are ranked by A339886. For example, the a(6) = 1 through a(16) = 5 partitions are:
32 222 322 332 432 3322 3332 4332 4432 5432 43332
2222 3222 22222 4322 33222 33322 33332 44322
32222 222222 43222 43322 333222
322222 332222 432222
2222222 3222222
(End)

Examples

			1 + x^3 + x^5 + x^6 + x^7 + x^8 + 2*x^9 + 2*x^10 + 2*x^11 + 3*x^12 + 3*x^13 + ...
q + q^73 + q^121 + q^145 + q^169 + q^193 + 2*q^217 + 2*q^241 + 2*q^265 + ...
a(10)=2 because we have [7,3] and [5,5].
From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(22)=13 symmetric unimodal compositions of 22+3=25 where the maximal part appears three times:
01:  [ 1 1 1 1 1 1 1 1 3 3 3 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 2 3 3 3 2 1 1 1 1 1 1 ]
03:  [ 1 1 1 1 1 5 5 5 1 1 1 1 1 ]
04:  [ 1 1 1 1 2 2 3 3 3 2 2 1 1 1 1 ]
05:  [ 1 1 1 2 5 5 5 2 1 1 1 ]
06:  [ 1 1 2 2 2 3 3 3 2 2 2 1 1 ]
07:  [ 1 1 3 5 5 5 3 1 1 ]
08:  [ 1 1 7 7 7 1 1 ]
09:  [ 1 2 2 5 5 5 2 2 1 ]
10:  [ 1 4 5 5 5 4 1 ]
11:  [ 2 2 2 2 3 3 3 2 2 2 2 ]
12:  [ 2 3 5 5 5 3 2 ]
13:  [ 2 7 7 7 2 ]
(End)
From _Gus Wiseman_, Feb 16 2021: (Start)
The a(7) = 1 through a(19) = 8 partitions are the following (A..J = 10..19). The Heinz numbers of these partitions are given by A341449.
  7  53  9    55  B    75    D    77    F      97    H      99      J
         333  73  533  93    553  95    555    B5    755    B7      775
                       3333  733  B3    753    D3    773    D5      955
                                  5333  933    5533  953    F3      973
                                        33333  7333  B33    5553    B53
                                                     53333  7533    D33
                                                            9333    55333
                                                            333333  73333
(End)
		

References

  • J. W. L. Glaisher, Identities, Messenger of Mathematics, 5 (1876), pp. 111-112. see Eq. I

Crossrefs

The ordered version is A000931.
Partitions with no ones are counted by A002865, ranked by A005408.
The even version is A035363, ranked by A066207.
The version for factorizations is A340101.
Partitions whose only even part is the smallest are counted by A341447.
The Heinz numbers of these partitions are given by A341449.
A000009 counts partitions into odd parts, ranked by A066208.
A025147 counts strict partitions with no 1's.
A025148 counts strict partitions with no 1's or 2's.
A026804 counts partitions whose smallest part is odd, ranked by A340932.
A027187 counts partitions with even length/maximum, ranks A028260/A244990.
A027193 counts partitions with odd length/maximum, ranks A026424/A244991.
A058695 counts partitions of odd numbers, ranked by A300063.
A058696 counts partitions of even numbers, ranked by A300061.
A340385 counts partitions with odd length and maximum, ranked by A340386.

Programs

  • Haskell
    a087897 = p [3,5..] where
       p [] _ = 0
       p _  0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 12 2011
    
  • Maple
    To get 128 terms: t4 := mul((1+x^(2^n)),n=0..7); t5 := mul((1+x^k),k=1..128): t6 := series(t5/t4,x,100); t7 := seriestolist(t6);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1,
          `if`(i<3, 0, b(n, i-2)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(n, n-1+irem(n, 2)):
    seq(a(n), n=0..80);  # Alois P. Heinz, Jun 11 2013
  • Mathematica
    max = 65; f[x_] := Product[ 1/(1 - x^(2k+1)), {k, 1, max}]; CoefficientList[ Series[f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 16 2011, after Emeric Deutsch *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<3, 0, b[n, i-2]+If[i>n, 0, b[n-i, i]]] ]; a[n_] := b[n, n-1+Mod[n, 2]]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Apr 01 2015, after Alois P. Heinz *)
    Flatten[{1, Table[PartitionsQ[n+1] - PartitionsQ[n], {n, 0, 80}]}] (* Vaclav Kotesovec, Dec 01 2015 *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&OddQ[Times@@#]&]],{n,0,30}] (* Gus Wiseman, Feb 16 2021 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - x) * eta(x^2 + A) / eta(x + A), n))} /* Michael Somos, Nov 13 2011 */
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A087897_T(n,k):
        if n==0: return 1
        if k<3 or n<0: return 0
        return A087897_T(n,k-2)+A087897_T(n-k,k)
    def A087897(n): return A087897_T(n,n-(n&1^1)) # Chai Wah Wu, Sep 23 2023, after Alois P. Heinz

Formula

Expansion of q^(-1/24) * (1 - q) * eta(q^2) / eta(q) in powers of q.
Expansion of (1 - x) / chi(-x) in powers of x where chi() is a Ramanujan theta function.
G.f.: 1 + x^3 + x^5*(1 + x) + x^7*(1 + x)*(1 + x^2) + x^9*(1 + x)*(1 + x^2)*(1 + x^3) + ... [Glaisher 1876]. - Michael Somos, Jun 20 2012
G.f.: Product_{k >= 1} 1/(1-x^(2*k+1)).
G.f.: Product_{k >= 1, k not a power of 2} (1+x^k).
G.f.: Sum_{k >= 1} x^(3*k)/Product_{j = 1..k} (1 - x^(2*j)). - Emeric Deutsch, Mar 30 2006
a(n) ~ exp(Pi*sqrt(n/3)) * Pi / (8 * 3^(3/4) * n^(5/4)) * (1 - (15*sqrt(3)/(8*Pi) + 11*Pi/(48*sqrt(3)))/sqrt(n) + (169*Pi^2/13824 + 385/384 + 315/(128*Pi^2))/n). - Vaclav Kotesovec, Aug 30 2015, extended Nov 04 2016
G.f.: 1/(1 - x^3) * Sum_{n >= 0} x^(5*n)/Product_{k = 1..n} (1 - x^(2*k)) = 1/((1 - x^3)*(1 - x^5)) * Sum_{n >= 0} x^(7*n)/Product_{k = 1..n} (1 - x^(2*k)) = ..., extending Deutsch's result dated Mar 30 2006. - Peter Bala, Jan 15 2021
G.f.: Sum_{n >= 0} x^(n*(2*n+1))/Product_{k = 2..2*n+1} (1 - x^k). (Set z = x^3 and q = x^2 in Mc Laughlin et al., Section 1.3, Entry 7.) - Peter Bala, Feb 02 2021
a(2*n+1) = Sum{j>=1} A008284(n+1-j,2*j - 1) and a(2*n) = Sum{j>=1} A008284(n-j, 2*j). - Gregory L. Simay, Sep 22 2023

A367405 Triangle read by rows where T(n,k) is the number of strict integer partitions of n with two distinct parts summing to k.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 18 2023

Keywords

Examples

			Triangle begins:
  1
  0  1
  0  0  2
  1  1  1  2
  1  0  1  1  3
  1  1  1  1  2  3
  1  1  1  2  2  2  4
  2  2  3  2  3  2  3  4
  2  2  3  2  3  3  3  3  5
  3  2  4  3  4  4  5  3  4  5
  3  3  5  4  4  5  5  5  4  4  6
  4  3  6  5  6  5  7  5  7  4  5  6
  5  5  7  7  8  7  8  8  7  7  5  5  7
  6  5  9  8 10  7 10  9 10  7  9  5  6  7
  7  7 10 10 12 11 11 11 12 10  9  9  6  6  8
  9  7 13 11 15 12 13 13 15 13 13  9 11  6  7  8
Row n = 9 counts the following strict partitions:
  (6,2,1)  (5,3,1)  (4,3,2)  (5,3,1)  (6,2,1)  (6,2,1)  (8,1)
                             (4,3,2)  (4,3,2)  (5,3,1)  (7,2)
                                                        (6,3)
                                                        (5,4)
Row n = 13 counts the following strict partitions (A=10, B=11, C=12):
  A21   931   841   751   652   751   841   931   A21  A21  C1
  7321  7321  832   742   643   7321  742   832   832  931  B2
  6421  5431  7321  6421  6421  652   7321  7321  742  841  A3
              6421  5431  5431  6421  643   643   652  751  94
              5431              5431  5431  6421            85
                                                            76
		

Crossrefs

Column n = k is A004526.
Column k = 3 is A025148.
For subsets instead of partitions we have A365541, non-binary A365381.
The non-binary version is A365661, non-strict A365543.
The non-binary complement is A365663, non-strict A046663.
Row sums are A366741, non-strict A366738.
The non-strict version is A367404.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], UnsameQ@@#&&MemberQ[Total/@Subsets[#,{2}], k]&]], {n,3,10}, {k,3,n}]

A339101 Number of compositions (ordered partitions) of n into distinct parts >= 3.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 1, 3, 3, 5, 5, 7, 13, 15, 21, 29, 35, 43, 79, 87, 123, 161, 221, 259, 349, 531, 645, 857, 1115, 1471, 1903, 2403, 2979, 4493, 5357, 7135, 9013, 11919, 14925, 19685, 23939, 30667, 42679, 52215, 67035, 86009, 109541, 137923, 177493, 222027, 277749
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 23 2020

Keywords

Examples

			a(7) = 3 because we have [7], [4, 3] and [3, 4].
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, p) option remember;
          `if`(n=0, p!, `if`((i-2)*(i+3)/2 b(n$2, 0):
    seq(a(n), n=0..50);  # Alois P. Heinz, Nov 23 2020
  • Mathematica
    nmax = 50; CoefficientList[Series[Sum[k! x^(k (k + 5)/2)/Product[1 - x^j, {j, 1, k}], {k, 0, nmax}], {x, 0, nmax}], x]

Formula

G.f.: Sum_{k>=0} k! * x^(k*(k + 5)/2) / Product_{j=1..k} (1 - x^j).

A026835 Triangular array read by rows: T(n,k) = number of partitions of n into distinct parts in which every part is >=k, for k=1,2,...,n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 1, 4, 2, 1, 1, 1, 1, 5, 3, 2, 1, 1, 1, 1, 6, 3, 2, 1, 1, 1, 1, 1, 8, 5, 3, 2, 1, 1, 1, 1, 1, 10, 5, 3, 2, 1, 1, 1, 1, 1, 1, 12, 7, 4, 3, 2, 1, 1, 1, 1, 1, 1, 15, 8, 5, 3, 2, 1, 1, 1, 1, 1, 1, 1, 18, 10, 6, 4, 3, 2, 1, 1, 1, 1, 1, 1, 1, 22, 12, 7, 4, 3, 2, 1, 1
Offset: 1

Views

Author

Keywords

Comments

T(n,1)=A000009(n), T(n,2)=A025147(n) for n>1, T(n,3)=A025148(n) for n>2, T(n,4)=A025149(n) for n>3.
A219922(n) = smallest number of row containing n. - Reinhard Zumkeller, Dec 01 2012

Examples

			From _Michael De Vlieger_, Aug 03 2020: (Start)
Table begins:
   1
   1   1
   2   1   1
   2   1   1   1
   3   2   1   1   1
   4   2   1   1   1   1
   5   3   2   1   1   1   1
   6   3   2   1   1   1   1   1
   8   5   3   2   1   1   1   1   1
  10   5   3   2   1   1   1   1   1   1
  12   7   4   3   2   1   1   1   1   1   1
  15   8   5   3   2   1   1   1   1   1   1   1
  ... (End)
		

Crossrefs

Programs

  • Haskell
    import Data.List (tails)
    a026835 n k = a026835_tabl !! (n-1) !! (k-1)
    a026835_row n = a026835_tabl !! (n-1)
    a026835_tabl = map
       (\row -> map (p $ last row) $ init $ tails row) a002260_tabl
       where p 0      _ = 1
             p _     [] = 0
             p m (k:ks) = if m < k then 0 else p (m - k) ks + p m ks
    -- Reinhard Zumkeller, Dec 01 2012
  • Mathematica
    Nest[Function[{T, n, r}, Append[T, Table[1 + Total[T[[##]] & @@@ Select[r, #[[-1]] > k + 1 &]], {k, 0, n}]]] @@ {#1, #2, Transpose[1 + {#2 - #3, #3}]} & @@ {#1, #2, Range[Ceiling[#2/2] - 1]} & @@ {#, Length@ #} &, {{1}}, 12] // Flatten (* Michael De Vlieger, Aug 03 2020 *)

Formula

G.f.: Sum_{k>=1} (y^k*(-1+Product_{i>=k} (1+x^i))). - Vladeta Jovovic, Aug 25 2003
T(n, k) = 1 + Sum(T(i, j): i>=j>k and i+j=n+1). - Reinhard Zumkeller, Jan 01 2003
T(n, k) > 1 iff 2*k < n. - Reinhard Zumkeller, Jan 01 2003

A096749 Number of partitions of n into distinct parts, the least being 2.

Original entry on oeis.org

0, 0, 1, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 15, 17, 20, 24, 28, 32, 38, 44, 51, 59, 68, 78, 91, 103, 118, 136, 155, 176, 201, 228, 259, 294, 332, 375, 425, 478, 538, 607, 681, 764, 858, 961, 1075, 1203, 1343, 1499, 1673, 1863, 2073, 2308, 2564, 2847, 3161, 3504
Offset: 0

Views

Author

N. J. A. Sloane, Sep 28 2008

Keywords

Comments

The old entry with this sequence number was a duplicate of A071569.
a(n), n>2 is the Euler transform of [0,0,1,1,1] joined with period [0,1]. - Georg Fischer, Aug 15 2020

Crossrefs

Cf. A096765 (least=1), A022824 (3), A022825 (4), A022826 (5), A022827 (6), A022828 (7), A022829 (8), A022830 (9), A022831 (10).

Programs

  • Maple
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`((i-2)*(i+3)/2 `if`(n<2, 0, b(n-2$2)):
    seq(a(n), n=0..60);  # Alois P. Heinz, Feb 07 2014
    # Using the function EULER from Transforms (see link at the bottom of the page).
    [0,0,1,op(EULER([0,0,1,1,seq(irem(n,2),n=1..57)]))]; # Peter Luschny, Aug 19 2020
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[(i-2)*(i+3)/2Jean-François Alcover, Oct 13 2014, after Alois P. Heinz *)
    Join[{0}, Table[Count[Last /@ Select[IntegerPartitions@n, DeleteDuplicates[#] == # &], 2], {n, 66}]] (* Robert Price, Jun 13 2020 *)

Formula

G.f.: x^2*Product_{j>=3} (1+x^j). - R. J. Mathar, Jul 31 2008
a(n) = A025148(n-2), n>1. - R. J. Mathar, Sep 30 2008
G.f.: Sum_{k>=1} x^(k*(k + 3)/2) / Product_{j=1..k-1} (1 - x^j). - Ilya Gutkovskiy, Nov 24 2020

A341283 Number of ways to partition n labeled elements into sets of different sizes of at least 3.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 1, 36, 57, 211, 331, 958, 29228, 64065, 294659, 1232479, 3549717, 11296603, 557617987, 1512758550, 8514685860, 41183585167, 251022906729, 838303110637, 4183056225010, 263978773601641, 887708421995331, 5813843897797861, 32212405278588967, 216518890998518716
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 28 2021

Keywords

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1,
         `if`(i>n, 0, b(n, i+1)+b(n-i, i+1)*binomial(n, i)))
        end:
    a:= n-> b(n, 3):
    seq(a(n), n=0..30);  # Alois P. Heinz, Apr 28 2021
  • Mathematica
    nmax = 29; CoefficientList[Series[Product[(1 + x^k/k!), {k, 3, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
    a[0] = 1; a[n_] := a[n] = -(n - 1)! Sum[DivisorSum[k, # (-#!)^(-k/#) &, # > 2 &] a[n - k]/(n - k)!, {k, 1, n}]; Table[a[n], {n, 0, 29}]

Formula

E.g.f.: Product_{k>=3} (1 + x^k/k!).

A343348 Irregular triangle read by rows where T(n,k) is the number of strict integer partitions of n with least gap k.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 2, 1, 0, 2, 1, 0, 1, 3, 1, 1, 0, 3, 2, 1, 0, 5, 2, 1, 0, 5, 3, 1, 0, 1, 7, 3, 1, 1, 0, 8, 4, 2, 1, 0, 10, 5, 2, 1, 0, 12, 6, 3, 1, 0, 15, 7, 3, 1, 0, 1, 17, 9, 4, 1, 1, 0, 21, 10, 4, 2, 1, 0, 25, 12, 6, 2, 1, 0, 29, 15, 6, 3, 1, 0, 35, 17, 8, 3, 1, 0
Offset: 0

Views

Author

Gus Wiseman, Apr 18 2021

Keywords

Comments

The least gap (or mex) of a partition is the least positive integer that is not a part.
Row lengths are chosen to be consistent with the non-strict case A264401.

Examples

			Triangle begins:
   1
   0   1
   1   0
   1   0   1
   1   1   0
   2   1   0
   2   1   0   1
   3   1   1   0
   3   2   1   0
   5   2   1   0
   5   3   1   0   1
   7   3   1   1   0
   8   4   2   1   0
  10   5   2   1   0
  12   6   3   1   0
  15   7   3   1   0   1
		

Crossrefs

Row sums are A000009.
Row lengths are A002024.
Column k = 1 is A025147.
Column k = 2 is A025148.
The non-strict version is A264401.
A000009 counts strict partitions.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A257993 gives the least gap of the partition with Heinz number n.
A339564 counts factorizations with a selected factor.
A342050 ranks partitions with even least gap.
A342051 ranks partitions with odd least gap.

Programs

  • Mathematica
    mingap[q_]:=Min@@Complement[Range[If[q=={},0,Max[q]]+1],q];
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&mingap[#]==k&]],{n,0,15},{k,Round[Sqrt[2*(n+1)]]}]
Showing 1-9 of 9 results.