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

A258234 Decimal expansion of a constant related to A107379.

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, May 24 2015

Keywords

Comments

Limit n->infinity (Integral_{x=0..1} Product_{k=1..n} x^k*(1-x^k) dx)^(1/n) = Limit n->infinity (A258191(n)/A258192(n))^(1/n) = 1/A258234 = 0.18515528932235959464731321119795428527382236445907508398560553036... .

Examples

			5.4008719041181541524660911191042700520294...
		

Crossrefs

Programs

  • Mathematica
    r^2/(r-1) /.FindRoot[-PolyLog[2, 1-r] == Log[r]^2, {r, E}, WorkingPrecision->117] (* Vaclav Kotesovec, Jun 11 2015 *)

Formula

Equals limit n->infinity A107379(n)^(1/n).
Equals limit n->infinity A173519(n)^(1/n).

Extensions

More terms from Vaclav Kotesovec, Jun 09 2015

A152146 Triangle read by rows: T(n,k) (n >= 0, 0 <= k <= n) = number of partitions of 2n into 2k odd parts.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 1, 1, 0, 2, 2, 1, 1, 0, 3, 3, 2, 1, 1, 0, 3, 5, 3, 2, 1, 1, 0, 4, 6, 5, 3, 2, 1, 1, 0, 4, 9, 7, 5, 3, 2, 1, 1, 0, 5, 11, 11, 7, 5, 3, 2, 1, 1, 0, 5, 15, 14, 11, 7, 5, 3, 2, 1, 1, 0, 6, 18, 20, 15, 11, 7, 5, 3, 2, 1, 1, 0, 6, 23, 26, 22, 15, 11, 7, 5, 3, 2, 1, 1
Offset: 0

Views

Author

R. J. Mathar, Sep 25 2009, indices corrected Jul 09 2012

Keywords

Comments

In both this and A152157, reading columns downwards "converges" to A000041.
Also the number of strict integer partitions of 2n with alternating sum 2k. Also the number of normal integer partitions of 2n of which 2k parts are odd, where a partition is normal if it covers an initial interval of positive integers. - Gus Wiseman, Jun 20 2021

Examples

			Triangle begins:
  1
  0  1
  0  1  1
  0  2  1   1
  0  2  2   1   1
  0  3  3   2   1   1
  0  3  5   3   2   1   1
  0  4  6   5   3   2   1  1
  0  4  9   7   5   3   2  1  1
  0  5 11  11   7   5   3  2  1  1
  0  5 15  14  11   7   5  3  2  1  1
  0  6 18  20  15  11   7  5  3  2  1  1
  0  6 23  26  22  15  11  7  5  3  2  1  1
  0  7 27  35  29  22  15 11  7  5  3  2  1  1
  0  7 34  44  40  30  22 15 11  7  5  3  2  1 1
  0  8 39  58  52  42  30 22 15 11  7  5  3  2 1 1
  0  8 47  71  70  55  42 30 22 15 11  7  5  3 2 1 1
  0  9 54  90  89  75  56 42 30 22 15 11  7  5 3 2 1 1
  0  9 64 110 116  97  77 56 42 30 22 15 11  7 5 3 2 1 1
  0 10 72 136 146 128 100 77 56 42 30 22 15 11 7 5 3 2 1 1
From _Gus Wiseman_, Jun 20 2021: (Start)
For example, row n = 6 counts the following partitions (B = 11):
  (75)  (3333)  (333111)  (33111111)  (3111111111)  (111111111111)
  (93)  (5331)  (531111)  (51111111)
  (B1)  (5511)  (711111)
        (7311)
        (9111)
The corresponding strict partitions are:
  (7,5)      (8,4)      (9,3)    (10,2)   (11,1)  (12)
  (6,5,1)    (5,4,3)    (7,3,2)  (9,2,1)
  (5,4,2,1)  (6,4,2)    (8,3,1)
             (7,4,1)
             (6,3,2,1)
The corresponding normal partitions are:
  43221    33321     3321111    321111111   21111111111  111111111111
  322221   332211    32211111   2211111111
  2222211  432111    222111111
           3222111
           22221111
(End)
		

Crossrefs

Cf. A035294 (row sums), A107379, A152140, A152157.
Column k = 1 is A004526.
Column k = 2-8 is A026810 - A026816.
The non-strict version is A239830.
The reverse non-strict version is A344610.
The reverse version is A344649
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A067659 counts strict partitions of odd length.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A124754 gives alternating sum of standard compositions (reverse: A344618).
A316524 is the alternating sum of the prime indices of n (reverse: A344616).
A344611 counts partitions of 2n with reverse-alternating sum >= 0.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-2)+`if`(i>n, 0, expand(sqrt(x)*b(n-i, i)))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(2*n, 2*n-1)):
    seq(T(n), n=0..12);  # Alois P. Heinz, Jun 21 2021
  • Mathematica
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&ats[#]==k&]],{n,0,30,2},{k,0,n,2}] (* Gus Wiseman, Jun 20 2021 *)

Formula

T(n,k) = A152140(2n,2k).

A206226 Number of partitions of n^2 into parts not greater than n.

Original entry on oeis.org

1, 1, 3, 12, 64, 377, 2432, 16475, 116263, 845105, 6292069, 47759392, 368379006, 2879998966, 22777018771, 181938716422, 1465972415692, 11902724768574, 97299665768397, 800212617435074, 6617003142869419, 54985826573015541, 458962108485797208, 3846526994743330075
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2012

Keywords

Comments

Also the number of partitions of n^2 using n or fewer numbers. Thus for n=3 one has: 9; 1,8; 2,7; 3,6; 4,5; 1,1,7; 1,2,6; 1,3,5; 1,4,4; 2,2,5; 2,3,4; 3,3,3. - J. M. Bergot, Mar 26 2014 [computations done by Charles R Greathouse IV]
The partitions in the comments above are the conjugates of the partitions in the definition. By conjugation we have: "partitions into parts <= m" are equinumerous with "partitions into at most m parts". - Joerg Arndt, Mar 31 2014
From Vaclav Kotesovec, May 25 2015: (Start)
In general, "number of partitions of j*n^2 into parts that are at most n" is (for j>0) asymptotic to c(j) * d(j)^n / n^2, where c(j) and d(j) are a constants.
-------
j c(j)
1 0.1582087202672504149766310999238...
2 0.0794245035465730707705885572860...
3 0.0530017980244665552354063060738...
4 0.0397666338404544208556554596295...
5 0.0318193213988281353709268311928...
...
17 0.0093617308583114626385718275875...
c(j) for big j asymptotically approaches 1 / (2*Pi*j).
---------
j d(j)
1 9.15337019245412246194853029240... = A258268
2 16.57962120993269533568313969522...
3 23.98280768122086592445663786762...
4 31.37931997386325137074644287711...
5 38.77298550971449870728474612568...
...
17 127.45526806942537991146993713837...
d(j) for big j asymptotically approaches j * exp(2).
(End)
d(j) = r^(2*j+1)/(r-1), where r is the root of the equation polylog(2, 1-r) + (j+1/2)*log(r)^2 = 0. - Vaclav Kotesovec, Jun 11 2015

Crossrefs

Column k=2 of A238016.
Cf. A258296 (j=2), A258293 (j=3), A258294 (j=4), A258295 (j=5).

Programs

  • 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(T(n^2, n), n=0..20); # Vaclav Kotesovec, May 25 2015 after Alois P. Heinz
  • Mathematica
    Table[SeriesCoefficient[Product[1/(1-x^k),{k,1,n}],{x,0,n^2}],{n,0,20}] (* Vaclav Kotesovec, May 25 2015 *)
    (* A program to compute the constants d(j) *) Table[r^(2*j+1)/(r-1) /.FindRoot[-PolyLog[2,1-r] == (j+1/2)*Log[r]^2, {r, E}, WorkingPrecision->60], {j, 1, 5}] (* Vaclav Kotesovec, Jun 11 2015 *)
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-x^k+x*O(x^(n^2)))),n^2)}
    for(n=0,25,print1(a(n),", "))

Formula

a(n) = [x^(n^2)] Product_{k=1..n} 1/(1 - x^k).
a(n) ~ c * d^n / n^2, where d = 9.1533701924541224619485302924013545... = A258268, c = 0.1582087202672504149766310999238742... . - Vaclav Kotesovec, Sep 07 2014

A206240 Number of partitions of n^2-n into parts not greater than n.

Original entry on oeis.org

1, 1, 2, 7, 34, 192, 1206, 8033, 55974, 403016, 2977866, 22464381, 172388026, 1341929845, 10573800028, 84192383755, 676491536028, 5479185281572, 44692412971566, 366844007355202, 3028143252035976, 25123376972033392, 209401287806758273, 1752674793617241002
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2012

Keywords

Comments

Also the number of partitions of n^2 into exactly n parts. - Seiichi Manyama, May 07 2018

Examples

			From _Seiichi Manyama_, May 07 2018: (Start)
n | Partitions of n^2 into exactly n parts
--+-------------------------------------------------------
1 | 1.
2 | 3+1 = 2+2.
3 | 7+1+1 = 6+2+1 = 5+3+1 = 5+2+2 = 4+4+1 = 4+3+2 = 3+3+3. (End)
		

Crossrefs

Programs

  • 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(T(n^2-n, n), n=0..20); # Vaclav Kotesovec, May 25 2015 after Alois P. Heinz
  • Mathematica
    Table[SeriesCoefficient[Product[1/(1-x^k),{k,1,n}],{x,0,n*(n-1)}],{n,0,20}] (* Vaclav Kotesovec, May 25 2015 *)
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-x^k+x*O(x^(n^2-n)))),n^2-n)}
    for(n=0,30,print1(a(n),", "))

Formula

a(n) = [x^(n^2-n)] Product_{k=1..n} 1/(1 - x^k).
a(n) ~ c * d^n / n^2, where d = 9.153370192454122461948530292401354540073... = A258268, c = 0.07005383646855329845970382163053268... . - Vaclav Kotesovec, Sep 07 2014

A152157 Triangle read by rows: T(n,k) (n>=0, 0<=k<=n) = number of partitions of 2n+1 into 2k+1 odd parts.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 1, 4, 3, 2, 1, 1, 1, 5, 5, 3, 2, 1, 1, 1, 7, 7, 5, 3, 2, 1, 1, 1, 8, 10, 7, 5, 3, 2, 1, 1, 1, 10, 13, 11, 7, 5, 3, 2, 1, 1, 1, 12, 18, 15, 11, 7, 5, 3, 2, 1, 1, 1, 14, 23, 21, 15, 11, 7, 5, 3, 2, 1, 1, 1, 16, 30, 28, 22, 15, 11, 7, 5, 3, 2, 1, 1
Offset: 0

Views

Author

R. J. Mathar, Sep 25 2009

Keywords

Comments

This number triangle satisfies T(n,k) = A008284(n+k+1,2*k+1), n,k >= 0. This means that T(n,k) is also the number of partitions of N:=n+k+1 into M:=2*k+1 parts. For the proof add 1 to every odd part of each partition of n':=2*n+1 into m:=2*k+1 parts which are all odd, and divide each part by a factor of 2, thus obtaining a partition of n+k+1 into m=2*k+1 parts. All partitions of N,for N>=1, into an odd number of parts M (M from {1,...,N}) are reached: just take k=(M-1)/2 and n=N-1-k. Each partition of N into an odd number of parts can only arise once from the given recipe (for given N and M the k and n values are unique). See also a comment by Franklin T. Adams-Watters on A152140. - Wolfdieter Lang, Jul 09 2012

Examples

			Triangle begins:
1
1  1
1  1   1
1  2   1   1
1  3   2   1   1
1  4   3   2   1   1
1  5   5   3   2   1   1
1  7   7   5   3   2   1  1
1  8  10   7   5   3   2  1  1
1 10  13  11   7   5   3  2  1  1
1 12  18  15  11   7   5  3  2  1  1
1 14  23  21  15  11   7  5  3  2  1  1
1 16  30  28  22  15  11  7  5  3  2  1  1
1 19  37  38  30  22  15 11  7  5  3  2  1  1
1 21  47  49  41  30  22 15 11  7  5  3  2  1 1
1 24  57  65  54  42  30 22 15 11  7  5  3  2 1 1
1 27  70  82  73  56  42 30 22 15 11  7  5  3 2 1  1
1 30  84 105  94  76  56 42 30 22 15 11  7  5 3 2  1 1
1 33 101 131 123  99  77 56 42 30 22 15 11  7 5 3  2 1 1
1 37 119 164 157 131 101 77 56 42 30 22 15 11 7 5  3 2 1 1
From _Wolfdieter Lang_, Jul 09 2012 (Start)
T(5,1) = 4 from the four partitions of 11 into 3 parts, all of which are odd: [1,1,9], [1,3,7], [1,5,5] and [3,3,5].
T(5,1) = 4 from the four partitions of 7 = 5+1+1 into 3 parts:
[1,1,5], [1,2,4], [1,3,3] and [2,2,3].
(End)
		

Crossrefs

Cf. A078408 (row sums), A107379, A152140, A152146, A008284.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1/sqrt(x), `if`(i<1, 0,
          b(n, i-2)+`if`(i>n, 0, expand(sqrt(x)*b(n-i, i)))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(2*n+1, 2*n+1)):
    seq(T(n), n=0..12);  # Alois P. Heinz, Jun 21 2021
  • Mathematica
    (* p = A008284 *) p[n_, 1] = 1; p[n_, k_] := p[n, k] = If[n >= k, Sum[p[n - i, k - 1], {i, 1, n - 1}] - Sum[p[n - i, k], {i, 1, k - 1}], 0];
    T[n_, k_] := p[n + k + 1, 2 k + 1];
    Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Oct 28 2019, after Wolfdieter Lang *)

Formula

T(n,k) = A152140(2n+1,2k+1).
T(n,k) = p(n+k+1,2*k+1), n >= 0, k >= 0, with p(N,M)= A008284(N,M), the number of partitions of N into M parts. See the sketch of the proof given above as a comment. - Wolfdieter Lang, Jul 09 2012
O.g.f. for column k: (x^k)/product(1-x^j,j=1..(2*k+1)), k>=0.
From the o.g.f.s of A008284. - Wolfdieter Lang, Jul 10 2012

Extensions

Indices corrected by R. J. Mathar, Jul 09 2012

A152140 Triangle read by rows: T(n,k) (n>=0, 0<=k<=n) = number of partitions of n into k odd parts.

Original entry on oeis.org

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

Views

Author

R. J. Mathar, Sep 25 2009, offset corrected Jul 09 2012

Keywords

Comments

The number of partitions of n into k odd parts is equal to the number of partitions of (n+k)/2 into k parts; or equivalently the number of partitions of (n-k)/2 into at most k parts. - Franklin T. Adams-Watters, Sep 25 2009

Examples

			n= 0, k= 0: [];
n= 1, k= 1: [1] ;
n= 2, k= 2: [1, 1] ;
n= 3, k= 1: [3] ;
n= 3, k= 3: [1, 1, 1] ;
n= 4, k= 2: [1, 3] ;
n= 4, k= 4: [1, 1, 1, 1];
n= 5, k= 1: [5];
n= 5, k= 3: [1, 1, 3];
n= 5, k= 5: [1, 1, 1, 1, 1];
n= 6, k= 2: [3, 3] or [1, 5];
n= 6, k= 4: [1, 1, 1, 3];
n= 6, k= 6: [1, 1, 1, 1, 1, 1];
Triangle begins:
1
0 1
0 0 1
0 1 0 1
0 0 1 0 1
0 1 0 1 0 1
0 0 2 0 1 0 1
0 1 0 2 0 1 0 1
0 0 2 0 2 0 1 0 1
0 1 0 3 0 2 0 1 0 1
0 0 3 0 3 0 2 0 1 0 1
0 1 0 4 0 3 0 2 0 1 0 1
0 0 3 0 5 0 3 0 2 0 1 0 1
0 1 0 5 0 5 0 3 0 2 0 1 0 1
0 0 4 0 6 0 5 0 3 0 2 0 1 0 1
0 1 0 7 0 7 0 5 0 3 0 2 0 1 0 1
0 0 4 0 9 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 8 0 10 0 7 0 5 0 3 0 2 0 1 0 1
0 0 5 0 11 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 10 0 13 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 5 0 15 0 14 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 12 0 18 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 6 0 18 0 20 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 14 0 23 0 21 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 6 0 23 0 26 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 16 0 30 0 28 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 7 0 27 0 35 0 29 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 19 0 37 0 38 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 7 0 34 0 44 0 40 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 21 0 47 0 49 0 41 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 8 0 39 0 58 0 52 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 24 0 57 0 65 0 54 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 8 0 47 0 71 0 70 0 55 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 27 0 70 0 82 0 73 0 56 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 0 9 0 54 0 90 0 89 0 75 0 56 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
0 1 0 30 0 84 0 105 0 94 0 76 0 56 0 42 0 30 0 22 0 15 0 11 0 7 0 5 0 3 0 2 0 1 0 1
		

Crossrefs

Cf. A000009 (row sums), A097304, A107379, A152146, A152157.

Programs

  • Maple
    b:= proc(n, i) option remember; local j; if n=0 then 1
          elif i<1 then 0 elif irem(i, 2)=0 then b(n, i-1)
          else []; for j from 0 to n/i do zip((x, y)->x+y, %,
          [0$j, b(n-i*j, i-2)], 0) od; %[] fi
        end:
    T:= n-> b(n$2):
    seq(T(n), n=0..13);  # Alois P. Heinz, May 31 2013
  • Mathematica
    nn = 10; CoefficientList[
    Series[Product[1/(1 - y x^i), {i, 1, nn, 2}], {x, 0, nn}], {x, y}] (* Geoffrey Critzer, May 31 2013 *)

A258268 Decimal expansion of a constant related to A206226.

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, May 25 2015

Keywords

Examples

			9.153370192454122461948530292401354540073...
		

Crossrefs

Programs

  • Mathematica
    r^3/(r-1) /.FindRoot[-PolyLog[2, 1-r] == 3*Log[r]^2/2, {r, E}, WorkingPrecision->120] (* Vaclav Kotesovec, Jun 11 2015 *)

Formula

Equals limit n->infinity A206226(n)^(1/n).
Equals limit n->infinity A206227(n)^(1/n).
Equals limit n->infinity A206240(n)^(1/n).

Extensions

More digits from Vaclav Kotesovec, Jun 10 2015

A281489 Number of partitions of n^2 into distinct odd parts.

Original entry on oeis.org

1, 1, 1, 2, 5, 12, 33, 93, 276, 833, 2574, 8057, 25565, 81889, 264703, 861889, 2824974, 9311875, 30851395, 102676439, 343112116, 1150785092, 3872588051, 13071583810, 44245023261, 150145281903, 510721124972, 1741020966255, 5947081503460, 20352707950277
Offset: 0

Views

Author

Alois P. Heinz, Jan 22 2017

Keywords

Examples

			a(0) = 1: [], the empty partition.
a(1) = 1: [1].
a(2) = 1: [1,3].
a(3) = 2: [1,3,5], [9].
a(4) = 5: [1,3,5,7], [7,9], [5,11], [3,13], [1,15].
a(5) = 12: [1,3,5,7,9], [5,9,11], [5,7,13], [3,9,13], [1,11,13], [3,7,15], [1,9,15], [3,5,17], [1,7,17], [1,5,19], [1,3,21], [25].
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(d*
          [0, 1, -1, 1][1+irem(d, 4)], d=divisors(j)), j=1..n)/n)
        end:
    a:= n-> b(n^2):
    seq(a(n), n=0..30);
  • Mathematica
    b[n_] := b[n] = If[n==0, 1, Sum[b[n-j]*Sum[d*{0, 1, -1, 1}[[1+Mod[d, 4]]], {d, Divisors[j]}], {j, 1, n}]/n];
    a[n_] := b[n^2];
    Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 23 2017, translated from Maple *)

Formula

a(n) = [x^(n^2)] Product_{j>=0} (1 + x^(2*j+1)).
a(n) = A000700(A000290(n)).
a(n) ~ exp(Pi*n/sqrt(6)) / (2^(7/4) * 3^(1/4) * n^(3/2)). - Vaclav Kotesovec, Apr 10 2017

A202261 Number of n-element subsets that can be chosen from {1,2,...,2*n} having element sum n^2.

Original entry on oeis.org

1, 1, 1, 3, 7, 18, 51, 155, 486, 1555, 5095, 17038, 57801, 198471, 689039, 2415043, 8534022, 30375188, 108815273, 392076629, 1420064031, 5167575997, 18885299641, 69287981666, 255121926519, 942474271999, 3492314839349, 12977225566680, 48349025154154
Offset: 0

Views

Author

Alois P. Heinz, Jan 20 2012

Keywords

Comments

a(n) is the number of partitions of n^2 into n distinct parts <= 2*n.
Taking the complement of each set, a(n) is also the number of partitions of n^2+n into n distinct parts <= 2*n. - Franklin T. Adams-Watters, Jan 20 2012
Also the number of partitions of n*(n+1)/2 into at most n parts of size at most n. a(4) = 7: 433, 442, 3322, 3331, 4222, 4321, 4411. - Alois P. Heinz, May 31 2020

Examples

			a(0) = 1: {}.
a(1) = 1: {1}.
a(2) = 1: {1,3}.
a(3) = 3: {1,2,6}, {1,3,5}, {2,3,4}.
a(4) = 7: {1,2,5,8}, {1,2,6,7}, {1,3,4,8}, {1,3,5,7}, {1,4,5,6}, {2,3,4,7},{2,3,5,6}.
		

Crossrefs

Column k=1 of A185282.

Programs

  • Maple
    b:= proc(n, i, t) option remember;
          `if`(it*(2*i-t+1)/2, 0,
          `if`(n=0, 1, b(n, i-1, t) +`if`(n b(n^2, 2*n, n):
    seq(a(n), n=0..30);
  • Mathematica
    b[n_, i_, t_] := b[n, i, t] = If[it*(2*i-t+1)/2, 0, If[n == 0, 1, b[n, i-1, t] + If[nJean-François Alcover, Feb 05 2015, after Alois P. Heinz *)

Formula

a(n) ~ sqrt(3) * 4^n / (Pi * n^2). - Vaclav Kotesovec, Sep 10 2014

A206227 Number of partitions of n^2+n into parts not greater than n.

Original entry on oeis.org

1, 1, 4, 19, 108, 674, 4494, 31275, 225132, 1662894, 12541802, 96225037, 748935563, 5900502806, 46976736513, 377425326138, 3056671009814, 24930725879856, 204623068332997, 1688980598900228, 14012122025369431, 116784468316023069, 977437078888272796, 8212186058546599006
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2012

Keywords

Crossrefs

Programs

  • 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(T(n^2+n, n), n=0..20); # Vaclav Kotesovec, May 25 2015 after Alois P. Heinz
  • Mathematica
    Table[SeriesCoefficient[Product[1/(1-x^k),{k,1,n}],{x,0,n*(n+1)}],{n,0,20}] (* Vaclav Kotesovec, May 25 2015 *)
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-x^k+x*O(x^(n^2+n)))),n^2+n)}
    for(n=0,30,print1(a(n),", "))

Formula

a(n) = [x^(n^2+n)] Product_{k=1..n} 1/(1 - x^k).
a(n) ~ c * d^n / n^2, where d = 9.1533701924541224619485302924013545... = A258268, c = 0.3572966225745094270279188015952797... . - Vaclav Kotesovec, Sep 07 2014
Showing 1-10 of 13 results. Next