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 31-40 of 3902 results. Next

A048574 Self-convolution of 1 2 3 5 7 11 15 22 30 42 56 77 ... (A000041).

Original entry on oeis.org

1, 4, 10, 22, 43, 80, 141, 240, 397, 640, 1011, 1568, 2395, 3604, 5360, 7876, 11460, 16510, 23588, 33418, 47006, 65640, 91085, 125596, 172215, 234820, 318579, 430060, 577920, 773130, 1030007, 1366644, 1806445, 2378892, 3121835, 4082796
Offset: 2

Views

Author

Keywords

Comments

Number of proper partitions of n into parts of two kinds (i.e. both kinds must be present). - Franklin T. Adams-Watters, Feb 08 2006

Examples

			a(4) = 22 because (1,2,3,5)*(5,3,2,1) = 5 + 6 + 6 + 5 = 22
		

Crossrefs

Essentially the same as A052837.
Cf. A122768.
Column k=2 of A060642.

Programs

  • Haskell
    a048574 n = a048574_list !! (n-2)
    a048574_list = f (drop 2 a000041_list) [1] where
    f (p:ps) rs = (sum $ zipWith (*) rs $ tail a000041_list) : f ps (p : rs)
    -- Reinhard Zumkeller, Nov 09 2015
    
  • Maple
    spec := [S,{C=Sequence(Z,1 <= card),B=Set(C,1 <= card),S=Prod(B,B)},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20); # Franklin T. Adams-Watters, Feb 08 2006
    # second Maple program:
    a:= n-> (p-> add(p(j)*p(n-j), j=1..n-1))(combinat[numbpart]):
    seq(a(n), n=2..40);  # Alois P. Heinz, May 26 2018
  • Mathematica
    a[n_] := First[ ListConvolve[ pp = Array[ PartitionsP, n], pp]]; Table[ a[n], {n, 1, 36}] (* Jean-François Alcover, Oct 21 2011 *)
    Table[ListConvolve[PartitionsP[Range[n]],PartitionsP[Range[n]]],{n,40}]// Flatten (* Harvey P. Dale, Oct 29 2020 *)
  • PARI
    a(n) = sum(k=1, n-1, numbpart(k)*numbpart(n-k)); \\ Michel Marcus, Dec 11 2016

Formula

From Franklin T. Adams-Watters, Feb 08 2006: (Start)
a(0) = 0, a(n) = A000712(n)-2*A000041(n) for n>0.
a(n) = Sum_{k=1..n-1} A000041(k)*A000041(n-k).
G.f.: ((Product_{k>0} 1/(1-x^k))-1)^2 = (exp(Sum_{k>0} (x^k/(1-x^k)/k))-1)^2. (End)
a(n) ~ exp(2*Pi*sqrt(n/3)) / (4*3^(3/4)*n^(5/4)). - Vaclav Kotesovec, Mar 10 2018

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Sep 29 2000

A221530 Triangle read by rows: T(n,k) = A000005(k)*A000041(n-k).

Original entry on oeis.org

1, 1, 2, 2, 2, 2, 3, 4, 2, 3, 5, 6, 4, 3, 2, 7, 10, 6, 6, 2, 4, 11, 14, 10, 9, 4, 4, 2, 15, 22, 14, 15, 6, 8, 2, 4, 22, 30, 22, 21, 10, 12, 4, 4, 3, 30, 44, 30, 33, 14, 20, 6, 8, 3, 4, 42, 60, 44, 45, 22, 28, 10, 12, 6, 4, 2, 56, 84, 60, 66, 30, 44, 14, 20, 9, 8, 2, 6
Offset: 1

Views

Author

Omar E. Pol, Jan 19 2013

Keywords

Comments

T(n,k) is the number of partitions of n that contain k as a part multiplied by the number of divisors of k.
It appears that T(n,k) is also the total number of appearances of k in the last k sections of the set of partitions of n multiplied by the number of divisors of k.
T(n,k) is also the number of partitions of k into equal parts multiplied by the number of ones in the j-th section of the set of partitions of n, where j = (n - k + 1).
For another version see A245095. - Omar E. Pol, Jul 15 2014

Examples

			For n = 6:
  -------------------------
  k   A000005        T(6,k)
  1      1  *  7   =    7
  2      2  *  5   =   10
  3      2  *  3   =    6
  4      3  *  2   =    6
  5      2  *  1   =    2
  6      4  *  1   =    4
  .         A000041
  -------------------------
So row 6 is [7, 10, 6, 6, 4, 2]. Note that the sum of row 6 is 7+10+6+6+2+4 = 35 equals A006128(6).
.
Triangle begins:
  1;
  1,   2;
  2,   2,  2;
  3,   4,  2,  3;
  5,   6,  4,  3,  2;
  7,  10,  6,  6,  2,  4;
  11, 14, 10,  9,  4,  4,  2;
  15, 22, 14, 15,  6,  8,  2,  4;
  22, 30, 22, 21, 10, 12,  4,  4,  3;
  30, 44, 30, 33, 14, 20,  6,  8,  3,  4;
  42, 60, 44, 45, 22, 28, 10, 12,  6,  4,  2;
  56, 84, 60, 66, 30, 44, 14, 20,  9,  8,  2,  6;
  ...
		

Crossrefs

Similar to A221529.
Columns 1-2: A000041, A139582. Leading diagonals 1-3: A000005, A000005, A062011. Row sums give A006128.

Programs

  • Mathematica
    A221530row[n_]:=DivisorSigma[0,Range[n]]PartitionsP[n-Range[n]];Array[A221530row,10] (* Paolo Xausa, Sep 04 2023 *)
  • PARI
    row(n) = vector(n, i, numdiv(i)*numbpart(n-i)); \\ Michel Marcus, Jul 18 2014

Formula

T(n,k) = d(k)*p(n-k) = A000005(k)*A027293(n,k).

A233346 Primes of the form p(k)^2 + q(m)^2 with k > 0 and m > 0, where p(.) is the partition function (A000041), and q(.) is the strict partition function (A000009).

Original entry on oeis.org

2, 5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 101, 109, 113, 137, 149, 157, 193, 229, 241, 349, 373, 509, 709, 733, 1033, 1049, 1213, 1249, 1453, 1493, 1669, 1789, 2141, 2237, 2341, 2917, 3037, 3137, 3361, 4217, 5801, 5897, 6029, 6073, 8821, 10301, 10937, 11057, 18229, 18289, 19249, 20173, 20341, 20389, 21017, 24001, 30977, 36913, 42793
Offset: 1

Views

Author

Zhi-Wei Sun, Dec 07 2013

Keywords

Comments

Conjecture: The sequence contains infinitely many terms.
This follows from part (i) of the conjecture in A233307. Similarly, the conjecture in A232504 implies that there are infinitely many primes of the form p(k) + q(m) with k and m positive integers.

Examples

			a(1) = 2 since p(1)^2 + q(1)^2 = 1^2 + 1^2 = 2.
a(2) = 5 since p(1)^2 + q(3)^2 = 1^2 + 2^2 = 5.
		

Crossrefs

Programs

  • Mathematica
    SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]]
    n=0
    Do[If[Mod[Prime[m]+1,4]>0,Do[If[PartitionsP[j]>=Sqrt[Prime[m]],Goto[aa],
    If[SQ[Prime[m]-PartitionsP[j]^2]==False,Goto[bb],Do[If[PartitionsQ[k]^2==Prime[m]-PartitionsP[j]^2,
    n=n+1;Print[n," ",Prime[m]];Goto[aa]];If[PartitionsQ[k]^2>Prime[m]-PartitionsP[j]^2,Goto[bb]];Continue,{k,1,2*Sqrt[Prime[m]]}]]];
    Label[bb];Continue,{j,1,Sqrt[Prime[m]]}]];
    Label[aa];Continue,{m,1,4475}]

A234470 Number of ways to write n = k + m with k > 0 and m > 2 such that p(k + phi(m)/2) is prime, where p(.) is the partition function (A000041) and phi(.) is Euler's totient function.

Original entry on oeis.org

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

Views

Author

Zhi-Wei Sun, Dec 26 2013

Keywords

Comments

Conjecture: a(n) > 0 if n > 3 is not among 27, 34, 50, 61, 74, 78, 115, 120, 123, 127.
This implies that there are infinitely many primes in the range of the partition function p(n).

Examples

			a(26) = 1 since 26 = 2 + 24 with p(2 + phi(24)/2) = p(6) = 11 prime.
a(54) = 1 since 54 = 27 + 27 with p(27 + phi(27)/2) = p(36) = 17977 prime.
a(73) = 1 since 73 = 1 + 72 with p(1 + phi(72)/2) = p(36) = 17977 prime.
a(110) = 1 since 110 = 65 + 45 with p(65 + phi(45)/2) = p(77) = 10619863 prime.
a(150) = 1 since 150 = 123 + 27 with p(123 + phi(27)/2) = p(132) = 6620830889 prime.
a(170) = 1 since 170 = 167 + 3 with p(167 + phi(3)/2) = p(168) = 228204732751 prime.
		

Crossrefs

Programs

  • Mathematica
    f[n_,k_]:=PartitionsP[k+EulerPhi[n-k]/2]
    a[n_]:=Sum[If[PrimeQ[f[n,k]],1,0],{k,1,n-3}]
    Table[a[n],{n,1,100}]

A339106 Triangle read by rows: T(n,k) = A000203(n-k+1)*A000041(k-1), n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 3, 1, 4, 3, 2, 7, 4, 6, 3, 6, 7, 8, 9, 5, 12, 6, 14, 12, 15, 7, 8, 12, 12, 21, 20, 21, 11, 15, 8, 24, 18, 35, 28, 33, 15, 13, 15, 16, 36, 30, 49, 44, 45, 22, 18, 13, 30, 24, 60, 42, 77, 60, 66, 30, 12, 18, 26, 45, 40, 84, 66, 105, 88, 90, 42, 28, 12, 36, 39, 75, 56, 132, 90, 154, 120, 126, 56
Offset: 1

Views

Author

Omar E. Pol, Nov 23 2020

Keywords

Comments

Conjecture 1: T(n,k) is the sum of all divisors of all (n - k + 1)'s in the n-th row of triangle A176206, assuming that A176206 has offset 1. The same for the triangle A340061.
Conjecture 2: the sum of row n equals A066186(n), the sum of all parts of all partitions of n.

Examples

			Triangle begins:
   1;
   3,  1;
   4,  3,  2;
   7,  4,  6,  3;
   6,  7,  8,  9,  5;
  12,  6, 14, 12, 15,  7;
   8, 12, 12, 21, 20, 21,  11;
  15,  8, 24, 18, 35, 28,  33,  15;
  13, 15, 16, 36, 30, 49,  44,  45,  22;
  18, 13, 30, 24, 60, 42,  77,  60,  66,  30;
  12, 18, 26, 45, 40, 84,  66, 105,  88,  90,  42;
  28, 12, 36, 39, 75, 56, 132,  90, 154, 120, 126, 56;
...
For n = 6 the calculation of every term of row 6 is as follows:
-------------------------
k   A000041        T(6,k)
1      1  *  12  =   12
2      1  *  6   =    6
3      2  *  7   =   14
4      3  *  4   =   12
5      5  *  3   =   15
6      7  *  1   =    7
.         A000203
-------------------------
The sum of row 6 is 12 + 6 + 14 + 12 + 15 + 7 = 66, equaling A066186(6).
		

Crossrefs

Mirror of A221529.
Row sums give A066186 (conjectured).
Main diagonal gives A000041.
Columns 1 and 2 give A000203.
Column 3 gives A074400.
Column 4 gives A272027.
Column 5 gives A274535.
Column 6 gives A319527.

Programs

  • Mathematica
    T[n_, k_] := DivisorSigma[1, n - k + 1] * PartitionsP[k - 1]; Table[T[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Amiram Eldar, Jan 08 2021 *)
  • PARI
    T(n, k) = sigma(n-k+1)*numbpart(k-1); \\ Michel Marcus, Jan 08 2021

Formula

T(n,k) = sigma(n-k+1)*p(k-1), n >= 1, 1 <= k <= n.

A087787 a(n) = Sum_{k=0..n} (-1)^(n-k)*A000041(k).

Original entry on oeis.org

1, 0, 2, 1, 4, 3, 8, 7, 15, 15, 27, 29, 48, 53, 82, 94, 137, 160, 225, 265, 362, 430, 572, 683, 892, 1066, 1370, 1640, 2078, 2487, 3117, 3725, 4624, 5519, 6791, 8092, 9885, 11752, 14263, 16922, 20416, 24167, 29007, 34254, 40921, 48213, 57345, 67409
Offset: 0

Views

Author

Vladeta Jovovic, Oct 07 2003

Keywords

Comments

Essentially first differences of A024786 (see the formula). Also, a(n) is the number of 2's in the last section of the set of partitions of n+2 (see A135010). - Omar E. Pol, Sep 10 2008
From Gus Wiseman, May 20 2024: (Start)
Also the number of integer partitions of n containing an even number of ones, ranked by A003159. The a(0) = 1 through a(8) = 15 partitions are:
() . (2) (3) (4) (5) (6) (7) (8)
(11) (22) (32) (33) (43) (44)
(211) (311) (42) (52) (53)
(1111) (222) (322) (62)
(411) (511) (332)
(2211) (3211) (422)
(21111) (31111) (611)
(111111) (2222)
(3311)
(4211)
(22211)
(41111)
(221111)
(2111111)
(11111111)
Also the number of integer partitions of n + 1 containing an odd number of ones, ranked by A036554.
(End)

Crossrefs

The unsigned version is A000070, strict A036469.
For powers of 2 instead number of partitions we have A001045.
The strict or odd version is A025147 or A096765.
The ordered version (compositions instead of partitions) is A078008.
For powers of 2 instead of -1 we have A259401, cf. A259400.
A002865 counts partitions with no ones, column k=0 of A116598.
A072233 counts partitions by sum and length.

Programs

  • Mathematica
    Table[Sum[(-1)^(n-k)*PartitionsP[k], {k,0,n}], {n,0,50}] (* Vaclav Kotesovec, Aug 16 2015 *)
    (* more efficient program *) sig = 1; su = 1; Flatten[{1, Table[sig = -sig; su = su + sig*PartitionsP[n]; Abs[su], {n, 1, 50}]}] (* Vaclav Kotesovec, Nov 06 2016 *)
    Table[Length[Select[IntegerPartitions[n], EvenQ[Count[#,1]]&]],{n,0,30}] (* Gus Wiseman, May 20 2024 *)
  • Python
    from sympy import npartitions
    def A087787(n): return sum(-npartitions(k) if n-k&1 else npartitions(k) for k in range(n+1)) # Chai Wah Wu, Oct 25 2023

Formula

G.f.: 1/(1+x)*1/Product_{k>0} (1-x^k).
a(n) = 1/n*Sum_{k=1..n} (sigma(k)+(-1)^k)*a(n-k).
a(n) = A024786(n+2)-A024786(n+1). - Omar E. Pol, Sep 10 2008
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*sqrt(3)*n) * (1 + (11*Pi/(24*sqrt(6)) - sqrt(3/2)/Pi)/sqrt(n) - (11/16 + (23*Pi^2)/6912)/n). - Vaclav Kotesovec, Nov 05 2016
a(n) = A000041(n) - a(n-1). - Jon Maiga, Aug 29 2019
Alternating partial sums of A000041. - Gus Wiseman, May 20 2024

A115131 Waring numbers for power sums functions in terms of elementary symmetric functions; irregular triangle T(n,k), read by rows, for n >= 1 and 1 <= k <= A000041(n).

Original entry on oeis.org

1, -2, 1, 3, -3, 1, -4, 4, 2, -4, 1, 5, -5, -5, 5, 5, -5, 1, -6, 6, 6, 3, -6, -12, -2, 6, 9, -6, 1, 7, -7, -7, -7, 7, 14, 7, 7, -7, -21, -7, 7, 14, -7, 1, -8, 8, 8, 8, 4, -8, -16, -16, -8, -8, 8, 24, 12, 24, 2, -8, -32, -16, 8, 20, -8, 1, 9, -9, -9, -9, -9, 9, 18, 18, 9, 9, 18, 3, -9, -27, -27, -27, -27, -9, 9, 36, 18, 54, 9, -9, -45, -30, 9, 27, -9, 1
Offset: 1

Views

Author

Wolfdieter Lang, Jan 13 2006

Keywords

Comments

Examples

			First few rows of triangle T(n,k) are as follows (see the link for rows 1..10):
   1;
  -2,  1;
   3, -3,  1;
  -4,  4,  2, -4, 1;
   5, -5, -5,  5, 5, -5, 1;
  ...
n=4: N*t^{(N)}_4 = -4*(sigma_4)^1 + 4*(sigma_1)*(sigma_3) + 2*(sigma_2)^2 -4*(sigma_1)^2*(sigma_2) + 1*(sigma_1)^4.
  (For 2 <= N < 4, one puts sigma_{N+1} = 0 = ... = sigma_4 = 0.) This becomes Sum_{k = 1..N} (x_k)^4 if the sigma functions are written in terms of the variables x_1, x_2, ..., x_N. E.g., for N=2: 0 + 0 + 2*(x_1*x_2)^2 -4*(x_1 + x_2)^2*(x_1*x_2) + 1*(x_1 + x_2)^4 = (x_1)^4 + (x_2)^4.
		

References

  • P. A. MacMahon, Combinatory Analysis, 2 vols., Chelsea, NY, 1960, see p. 5 (with a_k -> sigma_k).

Crossrefs

Cf. A210258 (in another ordering of partitions), A132460 (N=2), A325477 (N=3),
A324602 (N=4).

Formula

T(n,k) = (n/m(n,k))*A111786(n,k) for the k-th partition of n with m(n,k) parts in the Abramowitz-Stegun order for n >= 1 and k = 1..p(n), where p(n) := A000041(n).
Explicitly: T(n,k) = (-1)^(n + m(n,k)) * n * (m(n,k) - 1)!/(Product_{j = 1..n} e(k,j)!), where m(n,k):= Sum_{j = 1..n} e(k,j), with [1^e(k, 1), 2^e(k,2), ..., n^e(k,n)] being the k-th partition of n in the mentioned order. For m(n,k), see A036043.

Extensions

Various sections edited by Petros Hadjicostas, Dec 14 2019

A058694 Partial products p(0)*p(1)*...*p(n) of partition numbers A000041.

Original entry on oeis.org

1, 1, 2, 6, 30, 210, 2310, 34650, 762300, 22869000, 960498000, 53787888000, 4141667376000, 418308404976000, 56471634671760000, 9939007702229760000, 2295910779215074560000, 681885501426877144320000, 262525918049347700563200000, 128637699844180373275968000000
Offset: 0

Views

Author

N. J. A. Sloane, Dec 30 2000

Keywords

Comments

a(n) gives the number of partitions P(V(n)) of V(n)=[1,2,3,...,n]. A partition P(V(n)) acts on the components of V(n), i.e., the components of V(n) are partitioned. Therefore a(n) results as the product of the number of partitions P(i) of the component v(i)=i with i=1,...,n. For example, a(3) = 6 because we have 6 list partitions for the list V(n=3)=[1,2,3]: [[1], [1, 1], [2, 1]], [[1], [1, 1], [1, 1, 1]], [[1], [1, 1], [3]], [[1], [2], [2, 1]], [[1], [2], [1, 1, 1]], [[1], [2], [3]]. - Thomas Wieder, Sep 29 2007
Equals the eigensequence of triangle A174712; i.e., Triangle A174712 * A058694 preceded by a 1 shifts left. - Gary W. Adamson, Mar 27 2010

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember;
           combinat[numbpart](n)*`if`(n>0, a(n-1), 1)
        end:
    seq(a(n), n=0..40);  # Alois P. Heinz, Apr 21 2012
    #
    # The constant S in the Maple notation
    evalf(Zeta(0, -1/2, 23/24)*sqrt(2/3)*Pi - Zeta(0, 1/2, 23/24)*sqrt(3/2)/Pi+3*(D(GAMMA))(23/24)/(4*Pi^2*GAMMA(23/24)) - (Sum(Zeta(0, j/2, 23/24)*(sqrt(3/2)/Pi)^j/j, j=3..infinity)), 60); # Vaclav Kotesovec, Jun 24 2015
  • Mathematica
    Table[Product[PartitionsP[k], {k, 1, n}], {n, 1, 33}] (* Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
  • PARI
    a(n)=prod(k=2,n, numbpart(k)) \\ Charles R Greathouse IV, Jan 14 2017

Formula

a(n) ~ C * Product_{k=1..n} (exp(Pi*sqrt(2/3*(k-1/24))) / (4*sqrt(3)*(k-1/24)) * (1 - sqrt(3/(2*(k-1/24)))/Pi)), where C = 0.9110167313322499518... is the partition factorial constant A259314. - Vaclav Kotesovec, Jun 24 2015
a(n) ~ C * Gamma(23/24) / (n^(n + 11/24 + 3/(4*Pi^2)) * 2^(2*n) * 3^(n/2) * sqrt(2*Pi)) * exp(Pi*(2*n/3)^(3/2) + n + (11*Pi/(12*sqrt(6)) - sqrt(6)/Pi)*sqrt(n) + S), where C = A259314 and S = Zeta(-1/2, 23/24)*sqrt(2/3)*Pi - Zeta(1/2, 23/24)*sqrt(3/2)/Pi + 3*Gamma'(23/24)/(4*Pi^2*Gamma(23/24)) - Sum_{j>=3} Zeta(j/2, 23/24)*(sqrt(3/2)/Pi)^j/j = -0.02541933397793652709903012019225640813047573968579474..., Zeta is the Hurwitz Zeta Function, in Maple notation Zeta(0,z,v), in Mathematica notation Zeta[z,v], equivalently HurwitzZeta[z,v]. - Vaclav Kotesovec, Jun 24 2015

A083711 a(n) = A083710(n) - A000041(n-1).

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 5, 3, 7, 1, 14, 1, 13, 8, 20, 1, 33, 1, 40, 14, 44, 1, 85, 6, 79, 25, 117, 1, 181, 1, 196, 45, 233, 17, 389, 1, 387, 80, 545, 1, 750, 1, 839, 165, 1004, 1, 1516, 12, 1612, 234, 2040, 1, 2766, 48, 3142, 388, 3720, 1, 5295, 1, 5606, 663, 7038, 83, 9194, 1, 10379, 1005
Offset: 1

Views

Author

N. J. A. Sloane, Jun 16 2003

Keywords

Comments

Number of integer partitions of n with no 1's with a part dividing all the others. If n > 0, we can assume such a part is the smallest. - Gus Wiseman, Apr 18 2021

Examples

			From _Gus Wiseman_, Apr 18 2021: (Start)
The a(6) = 4 through a(12) = 13 partitions:
  (6)      (7)  (8)        (9)      (10)         (11)  (12)
  (3,3)         (4,4)      (6,3)    (5,5)              (6,6)
  (4,2)         (6,2)      (3,3,3)  (8,2)              (8,4)
  (2,2,2)       (4,2,2)             (4,4,2)            (9,3)
                (2,2,2,2)           (6,2,2)            (10,2)
                                    (4,2,2,2)          (4,4,4)
                                    (2,2,2,2,2)        (6,3,3)
                                                       (6,4,2)
                                                       (8,2,2)
                                                       (3,3,3,3)
                                                       (4,4,2,2)
                                                       (6,2,2,2)
                                                       (4,2,2,2,2)
                                                       (2,2,2,2,2,2)
(End)
		

References

  • L. M. Chawla, M. O. Levan and J. E. Maxfield, On a restricted partition function and its tables, J. Natur. Sci. and Math., 12 (1972), 95-101.

Crossrefs

Allowing 1's gives A083710.
The strict case is A098965.
The complement (except also without 1's) is counted by A338470.
The dual version is A339619.
A000005 counts divisors.
A000041 counts partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.
A339564 counts factorizations with a selected factor.

Programs

  • Maple
    with(combinat): with(numtheory): a := proc(n) c := 0: l := sort(convert(divisors(n), list)): for i from 1 to nops(l)-1 do c := c+numbpart(l[i]-1) od: RETURN(c): end: for j from 2 to 100 do printf(`%d,`,a(j)) od: # James Sellers, Jun 21 2003
    # second Maple program:
    a:= n-> max(1, add(combinat[numbpart](d-1), d=numtheory[divisors](n) minus {n})):
    seq(a(n), n=1..69);  # Alois P. Heinz, Feb 15 2023
  • Mathematica
    a[n_] := If[n==1, 1, Sum[PartitionsP[d-1], {d, Most@Divisors[n]}]];
    Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Feb 15 2023 *)

Formula

a(n) = Sum_{ d|n, dA000041(d-1).

Extensions

More terms from James Sellers, Jun 21 2003

A174713 Triangle read by rows, A173305 (A000009 shifted down twice) * A174712 (diagonalized variant of A000041).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 3, 2, 2, 4, 2, 2, 3, 5, 3, 4, 3, 6, 4, 4, 3, 5, 8, 5, 6, 6, 5, 10, 6, 8, 6, 5, 7, 12, 8, 10, 9, 10, 7, 15, 10, 12, 12, 10, 7, 11, 18, 12, 16, 15, 15, 14, 11, 22, 15, 20, 18, 20, 14, 11, 15
Offset: 0

Views

Author

Gary W. Adamson, Mar 27 2010

Keywords

Comments

Row sums = A000041, the partition numbers.
The current triangle is the 2nd in an infinite set, followed by A174714 (k=3), and A174715, (k=4); in which row sums of each triangle = A000041.
k-th triangle in the infinite set can be defined as having the sequence:
"Euler transform of ones: (1,1,1,...) interleaved with (k-1) zeros"; shifted down k times (except column 0) in successive columns, then multiplied * triangle A174712, the diagonalized variant of A000041, A174713 begins with A000009 shifted down twice (triangle A173305); where A000009 = the Euler transform of period 2 sequence: [1,0,1,0,...].
Similarly, triangle A174714 begins with A000716 shifted down thrice; where A000716 = the Euler transform of period 3 series: [1,1,0,1,1,0,...]. Then multiply the latter as an infinite lower triangular matrix * A174712, the diagonalized variant of A000041, obtaining triangle A174714 with row sums = A000041.
Case k=4 = triangle A174715 which begins with the Euler transform of period 4 series: [1,1,1,0,1,1,1,0,...], shifted down 4 times in successive columns then multiplied * A174712, the diagonalized variant of A000041.
All triangles in the infinite set have row sums = A000041.
The sequences: "Euler transform of ones interleaved with (k-1) zeros" have the following properties, beginning with k=2:
...
k=2, A000009: = Euler transform of [1,0,1,0,1,0,...] and satisfies
.....A000009. = p(x)/p(x^2), where p(x) = polcoeff A000041; and A000041 =
.....A000009(x) = r(x), then p(x) = r(x) * r(x^2) * r(x^4) * r(x^8) * ...
...
k=3, A000726: = Euler transform of [1,1,0,1,1,0,...] and satisfies
.....A000726(x): = p(x)/p(x^3), and given s(x) = polcoeff A000726, we get
.....A000041(x) = p(x) = s(x) * s(x^3) * s(x^9) * s(x^27) * ...
...
k=4, A001935: = Euler transform of [1,1,1,0,1,1,1,0,...] and satisfies
.....A001935(x) = p(x)/p(x^4) and given t(x) = polcoeff A001935, we get
.....A000041(x) = p(x) = t(x) * t(x^4) * t(x^16) * t(x^64) * ...
...
Also the number of integer partitions of n whose even parts sum to k, for k an even number from zero to n. The version including odd k is A113686. - Gus Wiseman, Oct 23 2023

Examples

			First few rows of the triangle =
1;
1;
1, 1;
2, 1;
2, 1, 2;
3, 2, 2;
4, 2, 2, 3;
5, 3, 4, 3;
6, 4, 4, 3, 5;
8, 5, 6, 6, 5;
10, 6, 8, 6, 5, 7;
12, 8, 10, 9, 10, 7;
15, 10, 12, 12, 10, 7, 11;
18, 12, 16, 15, 15, 14, 11;
22, 15, 20, 18, 20, 14, 11, 15;
...
From _Gus Wiseman_, Oct 23 2023: (Start)
Row n = 9 counts the following partitions:
  (9)          (72)        (54)       (63)      (81)
  (711)        (5211)      (522)      (6111)    (621)
  (531)        (3321)      (4311)     (432)     (441)
  (51111)      (321111)    (411111)   (42111)   (4221)
  (333)        (21111111)  (32211)    (3222)    (22221)
  (33111)                  (2211111)  (222111)
  (3111111)
  (111111111)
(End)
		

Crossrefs

Row sums are A000041.
The odd version is A365067.
The corresponding rank statistic is A366531, odd version A366528.
A000009 counts partitions into odd parts, ranks A066208.
A113685 counts partitions by sum of odd parts, even version A113686.
A239261 counts partitions with (sum of odd parts) = (sum of even parts).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Total[Select[#,EvenQ]]==k&]],{n,0,15},{k,0,n,2}] (* Gus Wiseman, Oct 23 2023 *)

Formula

As infinite lower triangular matrices, A173305 * A174712.
T(n,k) = A000009(n-2k) * A000041(k). - Gus Wiseman, Oct 23 2023
Previous Showing 31-40 of 3902 results. Next