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 11-20 of 117 results. Next

A003108 Number of partitions of n into cubes.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 11, 11, 11, 12, 12, 13, 13, 13, 14, 14, 14, 15, 15, 17, 17, 17, 18, 18, 18, 19, 19, 21, 21, 21, 22, 22, 22, 23, 23, 25, 26, 26, 27, 27, 27, 28
Offset: 0

Views

Author

Keywords

Comments

The g.f. 1/(z+1)/(z**2+1)/(z**4+1)/(z-1)**2 conjectured by Simon Plouffe in his 1992 dissertation is wrong.

Examples

			a(16) = 3 because we have [8,8], [8,1,1,1,1,1,1,1,1] and [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1].
G.f.: A(x) = 1 + x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + 2*x^8 +...
such that the g.f. A(x) satisfies the identity [Paul D. Hanna]:
A(x) = 1/((1-x)*(1-x^8)*(1-x^27)*(1-x^64)*(1-x^125)*...)
A(x) = 1 + x/(1-x) + x^8/((1-x)*(1-x^8)) + x^27/((1-x)*(1-x^8)*(1-x^27)) + x^64/((1-x)*(1-x^8)*(1-x^27)*(1-x^64)) +...
		

References

  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • F. Smarandache, Sequences of Numbers Involved in Unsolved Problems, Hexis, Phoenix, 2006.

Crossrefs

Programs

  • Haskell
    a003108 = p $ tail a000578_list where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Oct 31 2012
    
  • Magma
    [#RestrictedPartitions(n,{d^3:d in [1..n]}): n in [0..150]]; // Marius A. Burtea, Jan 02 2019
    
  • Maple
    g:=1/product(1-x^(j^3),j=1..30): gser:=series(g,x=0,70): seq(coeff(gser,x,n),n=0..65); # Emeric Deutsch, Mar 30 2006
  • Mathematica
    nmax = 100; CoefficientList[Series[Product[1/(1 - x^(k^3)), {k, 1, nmax^(1/3)}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 19 2015 *)
    nmax = 60; cmax = nmax^(1/3);
    s = Table[n^3, {n, cmax}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Jul 31 2020 *)
  • PARI
    {a(n)=polcoeff(1/prod(k=1, ceil(n^(1/3)), 1-x^(k^3)+x*O(x^n)), n)} /* Paul D. Hanna, Mar 09 2012 */
    
  • PARI
    {a(n)=polcoeff(1+sum(m=1, ceil(n^(1/3)), x^(m^3)/prod(k=1, m, 1-x^(k^3)+x*O(x^n))), n)} /* Paul D. Hanna, Mar 09 2012 */
    
  • Python
    from functools import lru_cache
    from sympy import integer_nthroot, divisors
    @lru_cache(maxsize=None)
    def A003108(n):
        @lru_cache(maxsize=None)
        def a(n): return integer_nthroot(n,3)[1]
        @lru_cache(maxsize=None)
        def c(n): return sum(d for d in divisors(n,generator=True) if a(d))
        return (c(n)+sum(c(k)*A003108(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

Formula

G.f.: 1/Product_{j>=1} (1-x^(j^3)). - Emeric Deutsch, Mar 30 2006
G.f.: Sum_{n>=0} x^(n^3) / Product_{k=1..n} (1 - x^(k^3)). - Paul D. Hanna, Mar 09 2012
a(n) ~ exp(4 * (Gamma(1/3)*Zeta(4/3))^(3/4) * n^(1/4) / 3^(3/2)) * (Gamma(1/3)*Zeta(4/3))^(3/4) / (24*Pi^2*n^(5/4)) [Hardy & Ramanujan, 1917]. - Vaclav Kotesovec, Dec 29 2016

A109298 Primal codes of finite idempotent functions on positive integers.

Original entry on oeis.org

1, 2, 9, 18, 125, 250, 1125, 2250, 2401, 4802, 21609, 43218, 161051, 300125, 322102, 600250, 1449459, 2701125, 2898918, 4826809, 5402250, 9653618, 20131375, 40262750, 43441281, 86882562, 181182375, 362364750, 386683451, 410338673, 603351125, 773366902, 820677346
Offset: 1

Views

Author

Jon Awbrey, Jul 06 2005

Keywords

Comments

Finite idempotent functions are identity maps on finite subsets, counting the empty function as the idempotent on the empty set.
From Gus Wiseman, Mar 09 2019: (Start)
Also numbers whose ordered prime signature is equal to the distinct prime indices in increasing order. A prime index of n is a number m such that prime(m) divides n. The ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base. The case where the prime indices are taken in decreasing order is A324571.
Also numbers divisible by prime(k) exactly k times for each prime index k. These are a kind of self-describing numbers (cf. A001462, A304679).
Also Heinz numbers of integer partitions where the multiplicity of m is m for all m in the support (counted by A033461). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of distinct elements of A062457. For example, 43218 = prime(1)^1 * prime(2)^2 * prime(4)^4.
(End)

Examples

			Writing (prime(i))^j as i:j, we have the following table of examples:
Primal Codes of Finite Idempotent Functions on Positive Integers
` ` ` 1 = { }
` ` ` 2 = 1:1
` ` ` 9 = ` ` 2:2
` ` `18 = 1:1 2:2
` ` 125 = ` ` ` ` 3:3
` ` 250 = 1:1 ` ` 3:3
` `1125 = ` ` 2:2 3:3
` `2250 = 1:1 2:2 3:3
` `2401 = ` ` ` ` ` ` 4:4
` `4802 = 1:1 ` ` ` ` 4:4
` 21609 = ` ` 2:2 ` ` 4:4
` 43218 = 1:1 2:2 ` ` 4:4
`161051 = ` ` ` ` ` ` ` ` 5:5
`300125 = ` ` ` ` 3:3 4:4
`322102 = 1:1 ` ` ` ` ` ` 5:5
`600250 = 1:1 ` ` 3:3 4:4
From _Gus Wiseman_, Mar 09 2019: (Start)
The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2) has prime signature {1,2} and the distinct prime indices are also {1,2}.
       1: {}
       2: {1}
       9: {2,2}
      18: {1,2,2}
     125: {3,3,3}
     250: {1,3,3,3}
    1125: {2,2,3,3,3}
    2250: {1,2,2,3,3,3}
    2401: {4,4,4,4}
    4802: {1,4,4,4,4}
   21609: {2,2,4,4,4,4}
   43218: {1,2,2,4,4,4,4}
  161051: {5,5,5,5,5}
  300125: {3,3,3,4,4,4,4}
  322102: {1,5,5,5,5,5}
  600250: {1,3,3,3,4,4,4,4}
(End)
		

Crossrefs

Programs

  • Mathematica
    Select[Range[10000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]==k]&]
  • PARI
    is(n) = my(f = factor(n)); for(i = 1, #f~, if(prime(f[i, 2]) != f[i, 1], return(0))); 1 \\ David A. Corneth, Mar 09 2019

Formula

Sum_{n>=1} 1/a(n) = Product_{n>=1} (1 + 1/prime(n)^n) = 1.6807104966... - Amiram Eldar, Jan 03 2021

Extensions

Offset set to 1, missing terms inserted and more terms added by Alois P. Heinz, Mar 08 2019

A037444 Number of partitions of n^2 into squares.

Original entry on oeis.org

1, 1, 2, 4, 8, 19, 43, 98, 220, 504, 1116, 2468, 5368, 11592, 24694, 52170, 108963, 225644, 462865, 941528, 1899244, 3801227, 7550473, 14889455, 29159061, 56722410, 109637563, 210605770, 402165159, 763549779, 1441686280, 2707535748, 5058654069, 9404116777
Offset: 0

Views

Author

Keywords

Comments

Is lim_{n->inf} a(n)^(1/n) > 1? - Paul D. Hanna, Aug 20 2002
The limit above is equal to 1 (see formula by Hardy & Ramanujan for A001156). - Vaclav Kotesovec, Dec 29 2016

Crossrefs

Entries with square index in A001156.
A row or column of the array in A259799.

Programs

  • Haskell
    a037444 n = p (map (^ 2) [1..]) (n^2) where
       p _      0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 14 2011
  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
          b(n, i-1) +`if`(i^2>n, 0, b(n-i^2, i)))
        end:
    a:= n-> b(n^2, n):
    seq(a(n), n=0..40);  # Alois P. Heinz, Apr 15 2013
  • Mathematica
    max=33; se = Series[ Product[1/(1-x^(k^2)), {k, 1, max}], {x, 0, max^2}]; a[n_] := Coefficient[se, x^(n^2)]; a[0] = 1; Table[a[n], {n, 0, max}] (* Jean-François Alcover, Oct 18 2011 *)

Formula

a(n) = A001156(n^2) = coefficient of x^(n^2) in the series expansion of Prod_{k>=1} 1/(1 - x^(k^2)).
a(n) ~ 3^(-1/2) * (4*Pi)^(-7/6) * Zeta(3/2)^(2/3) * n^(-7/3) * exp(2^(-4/3) * 3 * Pi^(1/3) * Zeta(3/2)^(2/3) * n^(2/3)) [Hardy & Ramanujan, 1917, modified from A001156]. - Vaclav Kotesovec, Dec 29 2016

A078134 Number of ways to write n as sum of squares > 1.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 2, 1, 1, 0, 2, 1, 1, 0, 2, 3, 1, 1, 2, 3, 1, 1, 3, 3, 3, 1, 5, 3, 3, 1, 5, 5, 3, 3, 5, 7, 3, 3, 6, 8, 6, 3, 9, 8, 8, 3, 9, 10, 9, 6, 9, 14, 9, 8, 11, 15, 12, 9, 15, 15, 16, 9, 18, 18, 18, 13, 19, 23, 18, 17, 21, 28, 22, 19, 26, 30, 28, 19, 31, 34, 34
Offset: 1

Views

Author

Reinhard Zumkeller, Nov 19 2002

Keywords

Comments

a(A078135(n))=0; a(A078136(n))=1; a(A078137(n))>0;
Conjecture (lower bound): for all k exists b(k) such that a(n)>k for n>b(k); see b(0)=A078135(12)=23 and b(1)=A078136(15)=39. This is true - see comments by Hieronymus Fischer.
Also first difference of A001156 (number of partitions of n into squares). - Wouter Meeussen, Oct 22 2005
Comments from Hieronymus Fischer, Nov 11 2007 (Start): First statement of monotony: a(n+k^2)>=a(n) for all k>1. Proof: we restrict ourselves on a(n)>0 (the case a(n)=0 is trivial). Let T(i), 1<=i<=a(n), be the a(n) different sum expressions of squares >1 representing n. Then, adding k^2 to those expressions, we get a(n) sums of squares T(i)+k^2, obviously representing n+k^2, thus a(n+k^2) cannot be less than a(n).
Second statement of monotony: a(n+m)>=max(a(n),a(m)) for all m with a(m)>1. Proof: let T(i), 1<=i<=a(n), be the a(n) different sum expressions of squares >1 representing n; let S(i), 1<=i<=a(m), be the a(m) different sum expressions of squares >1 representing m. Then, adding those expressions, we get a(n) sums of squares T(i)+S(1), representing n+m, further we get a(m) sums T(1)+S(i), also representing n+m, thus a(n+m) cannot be less than the maximum of a(n) and a(m).
The author's conjecture holds true. Proof by induction: b(0) exists; if b(k) exists, then a(j)>k for all j>b(k). Setting m:=b(k)+1, we find that there are k+1 sums B(0,i) of squares >1, 1<=i<=k+1, with m=B(0,i). Further there are k+1 such sum expressions B(1,i), B(2,i) and B(3,i), 1<=i<=k+1, representing m+1, m+2 and m+3, respectively. For n>b(k) we have n=m+4*floor((n-m)/4)+(n-m) mod 4.
Thus n=m+r+s*2^2, where r=0,1,2 or 3. Hence n can be written B(r,i)+s*2^2 and there are k+1 such representations. Let q be the maximal number (to be squared) occurring as a term within those sum expressions B(r,i), 0<=r<=3,1<=i<=k+1. We select a number p>q and we set c:=b(k)+p^2. For n>c, we have the k+1 representations B(r(n),i)+s(n)*2^2.
Additionally, for n-p^2 (which is >b(k)) there are also k+1 representations B(r_p,i)+s_p*2^2, where r_p:=r(n-p^2), s_p:=s(n-p^2). Thus n can be written B(r(n),i)+s(n)*2^2, 1<=i<=k+1 and B(r_p,i)+s_p*2^2+p^2, 1<=i<=k+1.
By choice of p all these sum representations of n are different, which implies, that there are 2k+2 such representations. It follows a(n)>2k+2>k+1 for all n>c, which implies, that b(k+1) exists.
A more precise formulation of the author's conjecture is "b(k):=min( n | a(j)>k for all j>n) exists for all k>=0". (End)
A033183(n) <= a(n). [From Reinhard Zumkeller, Nov 07 2009]

Examples

			a(42)=3: 2*3^2+6*2^2 = 4^2+2*3^2+2*2^2 = 5^2+3^2+2*2^2.
		

Crossrefs

See A134754 for the sequence representing b(k).

Programs

  • Haskell
    a078134 = p $ drop 2 a000290_list where
       p _          0 = 1
       p ks'@(k:ks) x = if x < k then 0 else p ks' (x - k) + p ks x
    -- Reinhard Zumkeller, May 04 2013
  • Mathematica
    Join[{1}, Table[Length[PowersRepresentations[n, n, 2]], {n, 1, 90}]] // Differences
    (* or *)
    m = 91; CoefficientList[Product[1/(1 - x^(k^2)), {k, 1, m}] + O[x]^m, x] // Differences (* Jean-François Alcover, Mar 01 2019 *)

Formula

a(n) = 1/n*Sum_{k=1..n} (A035316(k)-1)*a(n-k), a(0) = 1. - Vladeta Jovovic, Nov 20 2002
G.f. g(x)=product{k>1, 1/(1-x^(k^2))}-1 = 1/((1-x^4)*(1-x^9)*(1-x^16)*(1-x^25)*(1-x^36)*...)-1. - Hieronymus Fischer, Nov 19 2007
a(n) ~ exp(3*Pi^(1/3) * Zeta(3/2)^(2/3) * n^(1/3) / 2^(4/3)) * Zeta(3/2)^(4/3) / (2^(11/3) * sqrt(3) * Pi^(5/6) * n^(11/6)). - Vaclav Kotesovec, Jan 05 2017

A087153 Number of partitions of n into nonsquares.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 3, 3, 5, 5, 8, 9, 13, 15, 20, 24, 30, 37, 47, 55, 71, 83, 103, 123, 151, 178, 218, 257, 310, 366, 440, 515, 617, 722, 857, 1003, 1184, 1380, 1625, 1889, 2214, 2570, 3000, 3472, 4042, 4669, 5414, 6244, 7221, 8303, 9583, 10998, 12655, 14502
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 21 2003

Keywords

Comments

Also, number of partitions of n where there are fewer than k parts equal to k for all k. - Jon Perry and Vladeta Jovovic, Aug 04 2004. E.g. a(8)=5 because we have 8=6+2=5+3=4+4=3+3+2.
Convolution of A276516 and A000041. - Vaclav Kotesovec, Dec 30 2016
From Gus Wiseman, Apr 02 2019: (Start)
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The Heinz numbers of the integer partitions described in Perry and Jovovic's comment are given by A325128, while the Heinz numbers of the integer partitions described in the name are given by A325129. In the former case, the first 10 terms count the following integer partitions:
() (2) (3) (4) (5) (6) (7) (8) (9)
(32) (33) (43) (44) (54)
(42) (52) (53) (63)
(62) (72)
(332) (432)
while in the latter case they count the following:
() (2) (3) (22) (5) (6) (7) (8) (63)
(32) (33) (52) (53) (72)
(222) (322) (62) (333)
(332) (522)
(2222) (3222)
(End)

Examples

			n=7: 2+5 = 2+2+3 = 7: a(7)=3;
n=8: 2+6 = 2+2+2+2 = 2+3+3 = 3+5 = 8: a(8)=5;
n=9: 2+7 = 2+2+5 = 2+2+2+3 = 3+3+3 = 3+6: a(9)=5.
		

References

  • G. E. Andrews, K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004. See page 48.

Crossrefs

Programs

  • Haskell
    a087153 = p a000037_list where
       p _          0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Apr 25 2013
    
  • Maple
    g:=product((1-x^(i^2))/(1-x^i),i=1..70):gser:=series(g,x=0,60):seq(coeff(gser,x^n),n=1..53); # Emeric Deutsch, Feb 09 2006
  • Mathematica
    nn=54; CoefficientList[ Series[ Product[ Sum[x^(i*j), {j, 0, i - 1}], {i, 1, nn}], {x, 0, nn}], x] (* Robert G. Wilson v, Aug 05 2004 *)
    nmax = 100; CoefficientList[Series[Product[(1 - x^(k^2))/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Dec 29 2016 *)
  • PARI
    first(n)=my(x='x+O('x^(n+1))); Vec(prod(m=1,sqrtint(n), (1-x^m^2)/(1-x^m))*prod(m=sqrtint(n)+1,n,1/(1-x^m))) \\ Charles R Greathouse IV, Aug 28 2016

Formula

G.f.: Product_{m>0} (1-x^(m^2))/(1-x^m). - Vladeta Jovovic, Aug 21 2003
a(n) = (1/n)*Sum_{k=1..n} (A000203(k)-A035316(k))*a(n-k), a(0)=1. - Vladeta Jovovic, Aug 21 2003
G.f.: Product_{i>=1} (Sum_{j=0..i-1} x^(i*j)). - Jon Perry, Jul 26 2004
a(n) ~ exp(Pi*sqrt(2*n/3) - 3^(1/4) * Zeta(3/2) * n^(1/4) / 2^(3/4) - 3*Zeta(3/2)^2/(32*Pi)) * sqrt(Pi) / (2^(3/4) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, Dec 30 2016

Extensions

Zeroth term added by Franklin T. Adams-Watters, Jan 25 2010

A243148 Triangle read by rows: T(n,k) = number of partitions of n into k nonzero squares; n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1
Offset: 0

Views

Author

Alois P. Heinz, May 30 2014

Keywords

Examples

			T(20,5) = 2 = #{ (16,1,1,1,1), (4,4,4,4,4) } since 20 = 4^2 + 4 * 1^2 = 5 * 2^2.
Triangle T(n,k) begins:
  1;
  0, 1;
  0, 0, 1;
  0, 0, 0, 1;
  0, 1, 0, 0, 1;
  0, 0, 1, 0, 0, 1;
  0, 0, 0, 1, 0, 0, 1;
  0, 0, 0, 0, 1, 0, 0, 1;
  0, 0, 1, 0, 0, 1, 0, 0, 1;
  0, 1, 0, 1, 0, 0, 1, 0, 0, 1;
  0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1;
  0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1;
  0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1;
  (...)
		

Crossrefs

Columns k = 0..10 give: A000007, A010052 (for n>0), A025426, A025427, A025428, A025429, A025430, A025431, A025432, A025433, A025434.
Row sums give A001156.
T(2n,n) gives A111178.
T(n^2,n) gives A319435.
T(n,k) = 1 for n in A025284, A025321, A025357, A294675, A295670, A295797 (for k = 2..7, respectively).

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n=0, `if`(t=0, 1, 0),
          `if`(i<1 or t<1, 0, b(n, i-1, t)+
          `if`(i^2>n, 0, b(n-i^2, i, t-1))))
        end:
    T:= (n, k)-> b(n, isqrt(n), k):
    seq(seq(T(n, k), k=0..n), n=0..14);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+(s-> `if`(s>n, 0, expand(x*b(n-s, i))))(i^2)))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, isqrt(n))):
    seq(T(n), n=0..14);  # Alois P. Heinz, Oct 30 2021
  • Mathematica
    b[n_, i_, k_, t_] := b[n, i, k, t] = If[n == 0, If[t == 0, 1, 0], If[i < 1 || t < 1, 0, b[n, i-1, k, t] + If[i^2 > n, 0, b[n-i^2, i, k, t-1]]]]; T[n_, k_] := b[n, Sqrt[n] // Floor, k, k]; Table[Table[T[n, k], {k, 0, n}], {n, 0, 14}] // Flatten (* Jean-François Alcover, Jun 06 2014, after Alois P. Heinz *)
    T[n_, k_] := Count[PowersRepresentations[n, k, 2], r_ /; FreeQ[r, 0]]; T[0, 0] = 1; Table[T[n, k], {n, 0, 14}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 19 2016 *)
  • PARI
    T(n,k,L=n)=if(n>k*L^2, 0, k>n-3, k==n, k<2, issquare(n,&n) && n<=L*k, k>n-6, n-k==3, L=min(L,sqrtint(n-k+1)); sum(r=0,min(n\L^2,k-1),T(n-r*L^2,k-r,L-1), n==k*L^2)) \\ M. F. Hasler, Aug 03 2020

Formula

T(n,k) = [x^n y^k] 1/Product_{j>=1} (1-y*x^A000290(j)).
Sum_{k=1..n} k * T(n,k) = A281541(n).
Sum_{k=1..n} n * T(n,k) = A276559(n).
Sum_{k=0..n} (-1)^k * T(n,k) = A292520(n).

A259792 Number of partitions of n^3 into cubes.

Original entry on oeis.org

1, 1, 2, 5, 17, 62, 258, 1050, 4365, 18012, 73945, 301073, 1214876, 4852899, 19187598, 75070201, 290659230, 1113785613, 4224773811, 15866483556, 59011553910, 217410395916, 793635925091, 2871246090593, 10297627606547, 36620869115355, 129166280330900
Offset: 0

Views

Author

N. J. A. Sloane, Jul 06 2015

Keywords

Crossrefs

A row of the array in A259799.
Cf. A279329.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
          b(n, i-1) +`if`(i^3>n, 0, b(n-i^3, i)))
        end:
    a:= n-> b(n^3, n):
    seq(a(n), n=0..26);  # Alois P. Heinz, Jul 10 2015
  • Mathematica
    $RecursionLimit = 1000; b[n_, i_] := b[n, i] = If[n==0 || i==1, 1, b[n, i-1] + If[ i^3>n, 0, b[n-i^3, i]]]; a[n_] := b[n^3, n]; Table[a[n], {n, 0, 26}] (* Jean-François Alcover, Jul 15 2015, after Alois P. Heinz *)

Formula

a(n) = [x^(n^3)] Product_{j>=1} 1/(1-x^(j^3)). - Alois P. Heinz, Jul 10 2015
a(n) = A003108(n^3). - Vaclav Kotesovec, Aug 19 2015
a(n) ~ exp(4 * (Gamma(1/3)*Zeta(4/3))^(3/4) * n^(3/4) / 3^(3/2)) * (Gamma(1/3)*Zeta(4/3))^(3/4) / (24*Pi^2*n^(15/4)) [after Hardy & Ramanujan, 1917]. - Vaclav Kotesovec, Dec 29 2016

Extensions

More term from Alois P. Heinz, Jul 10 2015

A276516 Expansion of Product_{k>=1} (1-x^(k^2)).

Original entry on oeis.org

1, -1, 0, 0, -1, 1, 0, 0, 0, -1, 1, 0, 0, 1, -1, 0, -1, 1, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, -1, 1, -1, 1, 1, 0, -1, 0, 0, 0, 0, 0, 0, -2, 1, 1, 1, 0, 0, -1, -1, 1, 1, -1, 0, 0, -1, 1, -1, 2, -1, 0, 1, -2, 0, 1, 0, 1, 0, -1, 0, -2, 2, 0
Offset: 0

Views

Author

Vaclav Kotesovec, Dec 12 2016

Keywords

Comments

The difference between the number of partitions of n into an even number of distinct squares and the number of partitions of n into an odd number of distinct squares. - Ilya Gutkovskiy, Jan 15 2018

Crossrefs

Programs

  • Mathematica
    nn = 15; CoefficientList[Series[Product[(1-x^(k^2)), {k, nn}], {x, 0, nn^2}], x]
    nmax = 200; nn = Floor[Sqrt[nmax]]+1; poly = ConstantArray[0, nn^2 + 1]; poly[[1]] = 1; poly[[2]] = -1; poly[[3]] = 0; Do[Do[poly[[j + 1]] -= poly[[j - k^2 + 1]], {j, nn^2, k^2, -1}];, {k, 2, nn}]; Take[poly, nmax+1]

Formula

a(n) = Sum_{k>=0} (-1)^k * A341040(n,k). - Alois P. Heinz, Feb 03 2021
a(n) = A033461(n) - 2*A339367(n). - R. J. Mathar, Jul 29 2025

A324572 Number of integer partitions of n whose multiplicities (where if x < y the multiplicity of x is counted prior to the multiplicity of y) are equal to the distinct parts in decreasing order.

Original entry on oeis.org

1, 1, 0, 0, 2, 0, 1, 0, 1, 1, 2, 0, 3, 0, 2, 0, 4, 1, 2, 1, 4, 1, 3, 1, 5, 3, 5, 1, 6, 2, 6, 1, 7, 2, 7, 2, 11, 4, 8, 3, 11, 5, 10, 4, 13, 5, 11, 5, 16, 8, 14, 5, 19, 8, 18, 6, 22, 8, 22, 7, 26, 10, 25, 8, 33, 12, 29, 11, 36, 13, 34, 12, 40, 16, 41, 14, 47, 17, 45, 16, 55
Offset: 0

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing partitions (cf. A001462, A304679).
The Heinz numbers of these partitions are given by A324571.
The case where the distinct parts are taken in increasing order is counted by A033461, with Heinz numbers given by A109298.

Examples

			The first 19 terms count the following integer partitions:
   1: (1)
   4: (22)
   4: (211)
   6: (3111)
   8: (41111)
   9: (333)
  10: (511111)
  10: (322111)
  12: (6111111)
  12: (4221111)
  12: (33222)
  14: (71111111)
  14: (52211111)
  16: (811111111)
  16: (622111111)
  16: (4444)
  16: (442222)
  17: (43331111)
  18: (9111111111)
  18: (7221111111)
  19: (533311111)
		

Crossrefs

Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

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

Extensions

More terms from Alois P. Heinz, Mar 08 2019

A090677 Number of ways to partition n into sums of squares of primes.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Dec 19 2003

Keywords

Comments

From Hieronymus Fischer, Nov 11 2007: (Start)
First statement of monotony: a(n+p^2)>=a(n) for all primes p. Proof: we restrict ourselves on a(n)>0 (the case a(n)=0 is trivial). Let T(i), 1<=i<=a(n), be the a(n) different sums of squares of primes representing n. Then, adding p^2 to those expressions, we get a(n) sums of squares of primes T(i)+p^2, obviously representing n+p^2, thus a(n+p^2) cannot be less than a(n).
Second statement of monotony: a(n+m)>=max(a(n),a(m)) for all m with a(m)>1. Proof: let T(i), 1<=i<=a(n), be the a(n) different sums of squares of primes representing n; let S(i), 1<=i<=a(m), be the a(m) different sums of squares of primes representing m. Then, adding these expressions, we get a(n) sums of squares of primes T(i)+S(1), representing n+m, further we get a(m) sums T(1)+S(i), also representing n+m. Thus a(n+m) cannot be less than the maximum of a(n) and a(m).
The minimum b(k):=min( n | a(j)>k for all j>n) exists for all k>=0. See A134755 for that sequence representing b(k). (End)

Examples

			a(25)=2 because 25 = 5^2 = 4*(2^2)+3^2.
a(83)=8 because 83 = 3^2+5^2+7^2 = 4*(2^2)+2*(3^2)+7^2
                   = 2*(2^2)+3*(5^2) = 6*(2^2)+3^2+2*(5^2)
                   = 2^2+6*(3^2)+5^2 = 10*(2^2)+2*(3^2)+5^2
                   = 5*(2^2)+7*(3^2) = 14*(2^2)+3*(3^2).
		

References

  • R. F. Churchouse, Representation of integers as sums of squares of primes. Caribbean J. Math. 5 (1986), no. 2, 59-65.

Crossrefs

Programs

  • Mathematica
    CoefficientList[ Series[ Product[1/(1 - x^Prime[i]^2), {i, 111}], {x, 0, 101}], x] (* Robert G. Wilson v, Sep 20 2004 *)

Formula

G.f.: 1/((1-x^4)*(1-x^9)*(1-x^25)*(1-x^49)*(1-x^121)*(1-x^169)*(1-x^289)...).
G.f.: 1 + Sum_{i>=1} x^(prime(i)^2) / Product_{j=1..i} (1 - x^(prime(j)^2)). - Ilya Gutkovskiy, May 07 2017
Previous Showing 11-20 of 117 results. Next