A103929 -id:A103929 - cp's OEIS frontend

A008951 Array read by columns: number of partitions of n into parts of 2 kinds.

1, 1, 1, 2, 2, 3, 4, 1, 5, 7, 2, 7, 12, 5, 11, 19, 9, 1, 15, 30, 17, 2, 22, 45, 28, 5, 30, 67, 47, 10, 42, 97, 73, 19, 1, 56, 139, 114, 33, 2, 77, 195, 170, 57, 5, 101, 272, 253, 92, 10, 135, 373, 365, 147, 20, 176, 508, 525, 227, 35, 1, 231, 684, 738, 345, 62, 2, 297

Offset: 0

N. J. A. Sloane

Fine-Riordan array S_n(m) = a(n,m) with extra row for n=0 added.
Row n of this triangle has length floor(1/2 + sqrt(2*(n+1))), n>=0. This is sequence {A002024(n+1)} = [1,2,2,3,3,3,4,4,4,4,5,5,5,5,5,6,6,6,6,6,6,...].
Written as a triangle this becomes A103923.
a(n,m) also gives the number of partitions of n-t(m), where t(m):=A000217(m) (triangular numbers), with two kinds of parts 1,2,..m. See the column o.g.f.'s in table A103923.
In general, column m is asymptotic to exp(Pi*sqrt(2*n/3)) * 6^(m/2) * n^((m-2)/2) / (4*sqrt(3) * m! * Pi^m), equivalently to 6^(m/2) * n^(m/2) / (m! * Pi^m) * p(n), where p(n) is the partition function A000041. - Vaclav Kotesovec, Aug 28 2015

			Array begins:
m\n 0 1 2 3 4 .5 .6 .7 .8 ...
0 | 1 1 2 3 5 .7 11 15 22 ... (A000041)
1 | . 1 2 4 7 12 19 ... (A000070)
2 | . . . 1 2 .5 .9 ... (A000097)
3 | . . . . . .. .1 ... (A000098)
[1]; [1,1]; [2,2]; [3,4,1]; [5,7,2]; [7,12,5]; [11,19,9,1]...
a(3,1) = 4 because the partitions (3), (1,2) and (1^3) have q values 1,2 and 1 which sum to 4.
a(3,1) = 4 because the exponents of part 1 in the above given partitions of 3 are 0,1,3 and they sum to 4.
a(3,1) = 4 because the partitions of 3-t(1)=2 with two kinds of part 1, say 1 and 1' and one kind of part 2 are (2),(1^2), (1'^2) and (11').

H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 90.
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 199.

Alois P. Heinz, Columns n = 0..500, flattened
William K. Keith, Restricted k-color partitions, arXiv preprint arXiv:1408.4089 [math.CO], 2014.
Wolfdieter Lang, First 20 rows and comments.

The first column (m=0) gives A000041(n). Columns m=1..10 are A000070 (partial sums of partition numbers), A000097, A000098, A000710, A103924-A103929.

a:= proc(n, m) option remember; `if`(n<0, 0,
      `if`(m=0, combinat[numbpart](n), a(n-m, m-1) +a(n-m, m)))
    end:
seq(seq(a(n,m), m=0..round(sqrt(2*n+2))-1), n=0..20);  # Alois P. Heinz, Nov 16 2012

a[n_, 0] := PartitionsP[n]; a[n_, m_] /; (n >= m*(m+1)/2) := a[n, m] = a[n-m, m-1] + a[n-m, m]; a[n_, m_] = 0; Flatten[ Table[ a[n, m], {n, 0, 18}, {m, 0, Floor[1/2 + Sqrt[2*(n+1)]] - 1}]](* Jean-François Alcover, May 02 2012, after recurrence formula *)
DeleteCases[Flatten@Transpose@Table[ConstantArray[0, m (m + 1)/2]~Join~Table[Length@IntegerPartitions[n, All, Range@n~Join~Range@m], {n, 0, 21 - m (m + 1)/2}] , {m, 0, 6}], 0](* Robert Price, Jul 28 2020 *)

Riordan gives formula.
a(n, m) is the sum over partitions of n of Product_{j=1..m} k(j), where k(j) is the number of parts of size j (exponent of j in a given partition of n), if m>=1. If m=0 then a(n, 0)=p(n):=A000041(n) (number of partitions of n). O is counted as a part for n=0 and only for this n.
a(n, m) is the sum over partitions of n of binomial(q(partition), m), with q the number of distinct parts of a given partition. m>=0.
a(n, m) = a(n-m, m-1) + a(n-m, m), n >= t(m):=m*(m+1)/2 = A000217(m) (triangular numbers), otherwise 0, with input a(n, 0) = p(n):=A000041(n).

More terms from Robert G Bearden (nem636(AT)myrealbox.com), Apr 27 2004

Correction, comments and Riordan formulas from Wolfdieter Lang, Apr 28 2005

A103923 Triangle of partitions of n with parts of sizes 1,2,...,m, each of two different kinds, m>=1.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 3, 4, 2, 1, 5, 7, 5, 2, 1, 7, 12, 9, 5, 2, 1, 11, 19, 17, 10, 5, 2, 1, 15, 30, 28, 19, 10, 5, 2, 1, 22, 45, 47, 33, 20, 10, 5, 2, 1, 30, 67, 73, 57, 35, 20, 10, 5, 2, 1, 42, 97, 114, 92, 62, 36, 20, 10, 5, 2, 1, 56, 139, 170, 147, 102, 64, 36, 20, 10, 5, 2, 1, 77, 195

Offset: 0

Wolfdieter Lang, Mar 24 2005

The corresponding Fine-Riordan triangle is A008951.
This is the array p_2(n,m) of Gupta et al. written as a triangle. p_2(n,m) is defined on p. x of this reference as the number of partitions of n into parts consisting of two varieties of each of the integers 1 to m and one variety of each larger integer. Therefore a(n,m) gives these numbers for the partitions of n-m.
a(n,m)= sum over partitions of n+t(m)-m of binomial(q(partition),m), with t(m):=A000217(m) and q the number of distinct parts of a given partition. m>=0.
a(n,m)= number of partitions of 2*n-m with exactly m odd parts.
a(n,m)= sum over partitions of n+t(m)-m of product(k[j],j=1..m), with t(m):=A000217(m) and k[j]=number of parts of size j (exponent of j in a given partition of n), if m>=1. If m=0 then a(n,0)=p(n):=A000041(n) (number of partitions of n). 0 is counted as a part for n=0 and only for this n.

			Triangle starts:
[1];
[1,1];
[2,2,1];
[3,4,2,1];
[5,7,5,2,1];
...
a(4,2)=5 from the partitions of 4-2=2 with two varieties of parts 1 and of 2, namely (2),(2'),(1^2),(1'^2) and (1,1').
a(4,2)=5 from the partitions of 4+t(2)-2=5 which have products of the exponents of parts 1 and 2: 0*0,1*0,0*1,2*1,1*2,5*0 and sum to 4.
a(4,2)=5 from the partitions of 4+t(2)-2=5 which have number of distinct parts (q values) 1,2,2,2,2,2,1. The corresponding binomial(q,2) values are 0,1,1,1,1,0 and sum to 4.
a(4,2)=5 from the partitions of 2*4-2=6 with exactly two odd parts, namely (1,5), (3^2), (1^2,4), (1,2,3) and (1^2,2^2), which are 5 in number.

H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958 (reprinted 1962), pp. 90-121.
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 199.

Alois P. Heinz, Rows n = 0..140, flattened
J. Huang, A. Senger, P. Wear, T. Wu, Partition statistics equidistributed with the number of hook difference one cells, 2013. See Remark 5.7. - _N. J. A. Sloane_, May 20 2014
W. Lang: First 16 rows.

The column sequences (without leading 0's) are, for m=0..10: A000041, A000070, A000097, A000098, A000710, A103924-A103929.

Cf. A000712, A124577.

with(numtheory):
b:= proc(n, k) option remember; `if`(n=0, 1, add(add(d*
      `if`(d<=k, 2, 1), d=divisors(j)) *b(n-j, k), j=1..n)/n)
    end:
A:= (n, k)-> b(n, k) -`if`(k=0, 0, b(n, k-1)):
seq(seq(A(n, k), k=0..n), n=0..14);  # Alois P. Heinz, Sep 14 2014

a[n_, 0] := a[n, 0] = PartitionsP[n]; a[n_, m_] /; n= m >= 0 := a[n, m] = a[n-1, m-1] + a[n-m, m]; Table[a[n, m], {n, 0, 14}, {m, 0, n}] // Flatten (* Jean-François Alcover, Dec 09 2014 *)
Flatten@Table[Length@IntegerPartitions[n-m, All, Range@n~Join~Range@m],  {n, 0, 12}, {m, 0, n}] (* Robert Price, Jul 29 2020 *)

a(n, m) = a(n-1, m-1) + a(n-m, m), n>=m>=0, with a(n, 0)= A000041(n) (partition numbers), a(n, m)=0 if n

a(n, m) = sum(a(n-1-j*m, m-1), j=0..floor((n-m)/m)), m>=1, input a(n, 0)= A000041(n).

G.f. column m: product(1/(1-x^j), j=1..m)*P(x), with P(x)= product(1/(1-x^j), j=1..infty), the o.g.f. for the partition numbers A000041.

G.f. column m>=1: (product(1/(1-x^k), k=1..m)^2)*product(1/(1-x^j), j=(m+1)..infty). For m=0 put the first product equal to 1.

A062876 Numbers of lattice points corresponding to incrementally largest circle radii in A062875.

Original entry on oeis.org

4, 12, 20, 28, 44, 52, 68, 76, 92, 116, 124, 148, 164, 172, 188, 212, 236, 244, 268, 284, 292, 316, 332, 356, 388, 404, 412, 428, 436, 452, 508, 524, 548, 556, 596, 604, 628, 652, 668, 692, 716, 724, 764, 772, 788, 796, 844, 892, 908, 916, 932, 956, 964

Offset: 1

Eric W. Weisstein

For n = 1 and n >= 3, a(n) is the smallest nonsquarefree number divisible by prime(n). - David James Sycamore, Jun 15 2024

Ray Chandler, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Circle Lattice Points

Cf. A046112, A062875.
Cf. A017113, A111333.
Cf. A000040, A103929

[4] cat [4*NthPrime(n): n in [2..60]]; // Vincenzo Librandi, May 08 2015

Join[{4}, Table[4 Prime[n], {n, 2, 50}]] (* Vincenzo Librandi, May 08 2015 *)

a(n)=if(n>1,4*prime(n),4) \\ Charles R Greathouse IV, May 08 2015

from sympy import prime
def A062876(n): return prime(n)<<2 if n>1 else 4 # Chai Wah Wu, Aug 02 2024

a(n) = A017113(A111333(n)-1) = 8*A111333(n) - 4.

For n >= 2 a(n) = 4*A000040(n) (a term in A013929). - David James Sycamore, Jun 15 2024

Edited and extended by Ray Chandler, Jan 05 2012

Showing 1-3 of 3 results.

cp's OEIS Frontend

A008951 Array read by columns: number of partitions of n into parts of 2 kinds.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A103923 Triangle of partitions of n with parts of sizes 1,2,...,m, each of two different kinds, m>=1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A062876 Numbers of lattice points corresponding to incrementally largest circle radii in A062875.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Formula

Extensions