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

A001156 Number of partitions of n into squares.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 4, 4, 4, 5, 6, 6, 6, 8, 9, 10, 10, 12, 13, 14, 14, 16, 19, 20, 21, 23, 26, 27, 28, 31, 34, 37, 38, 43, 46, 49, 50, 55, 60, 63, 66, 71, 78, 81, 84, 90, 98, 104, 107, 116, 124, 132, 135, 144, 154, 163, 169, 178, 192, 201, 209, 220, 235, 247, 256
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n such that number of parts equal to k is multiple of k for all k. - Vladeta Jovovic, Aug 01 2004
Of course p_{4*square}(n)>0. In fact p_{4*square}(32n+28)=3 times p_{4*square}(8n+7) and p_{4*square}(72n+69) is even. These seem to be the only arithmetic properties the function p_{4*square(n)} possesses. Similar results hold for partitions into positive squares, distinct squares and distinct positive squares. - Michael David Hirschhorn, May 05 2005
The Heinz numbers of these partitions are given by A324588. - Gus Wiseman, Mar 09 2019

Examples

			p_{4*square}(23)=1 because 23 = 3^2 + 3^2 + 2^2 + 1^2 and there is no other partition of 23 into squares.
G.f.: A(x) = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 2*x^6 + 2*x^7 +...
such that the g.f. A(x) satisfies the identity [_Paul D. Hanna_]:
A(x) = 1/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)*(1-x^25)*...)
A(x) = 1 + x/(1-x) + x^4/((1-x)*(1-x^4)) + x^9/((1-x)*(1-x^4)*(1-x^9)) + x^16/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)) + ...
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(14) = 6 integer partitions into squares are:
  (941)
  (911111)
  (44411)
  (44111111)
  (41111111111)
  (11111111111111)
while the a(14) = 6 integer partitions in which the multiplicity of k is a multiple of k for all k are:
  (333221)
  (33311111)
  (22222211)
  (2222111111)
  (221111111111)
  (11111111111111)
(End)
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000041, A000161 (partitions into 2 squares), A000290, A033461, A131799, A218494, A285218, A304046.
Cf. A078134 (first differences).
Row sums of A243148.
Euler trans. of A010052 (see also A308297).

Programs

  • Haskell
    a001156 = p (tail a000290_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, Aug 14 2011
    
  • Magma
    m:=70; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[1/(1-x^(k^2)): k in [1..(m+2)]]) )); // G. C. Greubel, Nov 11 2018
  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+ `if`(i^2>n, 0, b(n-i^2, i))))
        end:
    a:= n-> b(n, isqrt(n)):
    seq(a(n), n=0..120);  # Alois P. Heinz, May 30 2014
  • Mathematica
    CoefficientList[ Series[Product[1/(1 - x^(m^2)), {m, 70}], {x, 0, 68}], x] (* Or *)
    Join[{1}, Table[Length@PowersRepresentations[n, n, 2], {n, 68}]] (* Robert G. Wilson v, Apr 12 2005, revised Sep 27 2011 *)
    f[n_] := Length@ IntegerPartitions[n, All, Range@ Sqrt@ n^2]; Array[f, 67] (* Robert G. Wilson v, Apr 14 2013 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i^2>n, 0, b[n-i^2, i]]]]; a[n_] := b[n, Sqrt[n]//Floor]; Table[a[n], {n, 0, 120}] (* Jean-François Alcover, Nov 02 2015, after Alois P. Heinz *)
  • PARI
    {a(n)=polcoeff(1/prod(k=1, sqrtint(n+1), 1-x^(k^2)+x*O(x^n)), n)} \\ Paul D. Hanna, Mar 09 2012
    
  • PARI
    {a(n)=polcoeff(1+sum(m=1, sqrtint(n+1), x^(m^2)/prod(k=1, m, 1-x^(k^2)+x*O(x^n))), n)} \\ Paul D. Hanna, Mar 09 2012
    

Formula

G.f.: Product_{m>=1} 1/(1-x^(m^2)).
G.f.: Sum_{n>=0} x^(n^2) / Product_{k=1..n} (1 - x^(k^2)). - Paul D. Hanna, Mar 09 2012
a(n) = (1/n)*Sum_{k=1..n} A035316(k)*a(n-k). - Vladeta Jovovic, Nov 20 2002
a(n) = f(n,1,3) with f(x,y,z) = if xReinhard Zumkeller, Nov 08 2009
Conjecture (Jan Bohman, Carl-Erik Fröberg, Hans Riesel, 1979): a(n) ~ c * n^(-alfa) * exp(beta*n^(1/3)), where c = 1/18.79656, beta = 3.30716, alfa = 1.16022. - Vaclav Kotesovec, Aug 19 2015
From Vaclav Kotesovec, Dec 29 2016: (Start)
Correct values of these constants are:
1/c = sqrt(3) * (4*Pi)^(7/6) / Zeta(3/2)^(2/3) = 17.49638865935104978665...
alfa = 7/6 = 1.16666666666666666...
beta = 3/2 * (Pi/2)^(1/3) * Zeta(3/2)^(2/3) = 3.307411783596651987...
a(n) ~ 3^(-1/2) * (4*Pi*n)^(-7/6) * Zeta(3/2)^(2/3) * exp(2^(-4/3) * 3 * Pi^(1/3) * Zeta(3/2)^(2/3) * n^(1/3)). [Hardy & Ramanujan, 1917]
(End)

Extensions

More terms from Eric W. Weisstein
More terms from Gh. Niculescu (ghniculescu(AT)yahoo.com), Oct 08 2006

A000048 Number of n-bead necklaces with beads of 2 colors and primitive period n, when turning over is not allowed but the two colors can be interchanged.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 5, 9, 16, 28, 51, 93, 170, 315, 585, 1091, 2048, 3855, 7280, 13797, 26214, 49929, 95325, 182361, 349520, 671088, 1290555, 2485504, 4793490, 9256395, 17895679, 34636833, 67108864, 130150493, 252645135, 490853403, 954437120, 1857283155
Offset: 0

Views

Author

Keywords

Comments

Also, for any m which is a multiple of n, the number of 2m-bead balanced binary necklaces of fundamental period 2n that are equivalent to their complements. [Clarified by Aaron Meyerowitz, Jun 01 2024]
Also binary Lyndon words of length n with an odd number of 1's (for n>=1).
Also number of binary irreducible polynomials of degree n having trace 1.
Also number of binary irreducible polynomials of degree n having linear coefficient 1 (this is the same as the trace-1 condition, as the reciprocal of an irreducible polynomial is again irreducible).
Also number of binary irreducible self-reciprocal polynomials of degree 2*n; there is no such polynomial for odd degree except for x+1.
Also number of binary vectors (x_1,...x_n) satisfying Sum_{i=1..n} i*x_i = 1 (mod n+1) = size of Varshamov-Tenengolts code VT_1(n).
Also the number of dynamical cycles of period 2n of a threshold Boolean automata network which is a quasi-minimal negative circuit of size nq where q is odd and which is updated in parallel. - Mathilde Noual (mathilde.noual(AT)ens-lyon.fr), Mar 03 2009
Also the number of 3-elements orbits of the symmetric group S3 action on irreducible polynomials of degree 2n, n>1, over GF(2). - Jean Francis Michon, Philippe Ravache (philippe.ravache(AT)univ-rouen.fr), Oct 04 2009
Conjecture: Also the number of caliber-n cycles of Zagier-reduced indefinite binary quadratic forms with sum invariant equal to s, where (s-1)/n is an odd integer. - Barry R. Smith, Dec 14 2014
The Metropolis, Stein, Stein (1973) reference on page 31 Table II lists a(k) for k = 2 to 15 and is actually for sequence A056303 since there a(k) = 0 for k<2. - Michael Somos, Dec 20 2014

Examples

			a(5) = 3 corresponding to the necklaces 00001, 00111, 01011.
a(6) = 5 from 000001, 000011, 000101, 000111, 001011.
		

References

  • B. D. Ginsburg, On a number theory function applicable in coding theory, Problemy Kibernetiki, No. 19 (1967), pp. 249-252.
  • H. Kawakami, Table of rotation sequences of x_{n+1} = x_n^2 - lambda, pp. 73-92 of G. Ikegami, Editor, Dynamical Systems and Nonlinear Oscillations, Vol. 1, World Scientific, 1986.
  • Robert M. May, "Simple mathematical models with very complicated dynamics." Nature, Vol. 261, June 10, 1976, pp. 459-467; reprinted in The Theory of Chaotic Attractors, pp. 85-93. Springer, New York, NY, 2004. The sequences listed in Table 2 are A000079, A027375, A000031, A001037, A000048, A051841. - N. J. A. Sloane, Mar 17 2019
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Like A000013, but primitive necklaces. Half of A064355.
Equals A042981 + A042982.
Cf. also A001037, A056303.
Very close to A006788 [Fisher, 1989].
bisection (odd terms) is A131203

Programs

  • Maple
    with(numtheory); A000048 := proc(n) local d,t1; if n = 0 then RETURN(1) else t1 := 0; for d from 1 to n do if n mod d = 0 and d mod 2 = 1 then t1 := t1+mobius(d)*2^(n/d)/(2*n); fi; od; RETURN(t1); fi; end;
  • Mathematica
    a[n_] := Total[ MoebiusMu[#]*2^(n/#)& /@ Select[ Divisors[n], OddQ]]/(2n); a[0] = 1; Table[a[n], {n,0,35}] (* Jean-François Alcover, Jul 21 2011 *)
    a[ n_] := If[ n < 1, Boole[n == 0], DivisorSum[ n, MoebiusMu[#] 2^(n/#) &, OddQ] / (2 n)]; (* Michael Somos, Dec 20 2014 *)
  • PARI
    A000048(n) = sumdiv(n,d,(d%2)*(moebius(d)*2^(n/d)))/(2*n) \\ Michael B. Porter, Nov 09 2009
    
  • PARI
    L(n, k) = sumdiv(gcd(n,k), d, moebius(d) * binomial(n/d, k/d) );
    a(n) = sum(k=0, n, if( (n+k)%2==1, L(n, k), 0 ) ) / n;
    vector(55,n,a(n)) \\ Joerg Arndt, Jun 28 2012
    
  • Python
    from sympy import divisors, mobius
    def a(n): return 1 if n<1 else sum(mobius(d)*2**(n//d) for d in divisors(n) if d%2)//(2*n) # Indranil Ghosh, Apr 28 2017

Formula

a(n) = (1/(2*n)) * Sum_{odd d divides n} mu(d)*2^(n/d), where mu is the Mobius function A008683.
a(n) = A056303(n) for all integer n>=2. - Michael Somos, Dec 20 2014
Sum_{k dividing m for which m/k is odd} k*a(k) = 2^(m-1). (This explains the observation that the sequence is very close to A006788. Unless m has some nontrivial odd divisors that are small relative to m, the term m*a(m) will dominate the sum. Thus, we see for instance that a(n) = A006788(n) when n has one of the forms 2^m or 2^m*p where p is an odd prime with a(2^m) < p.) - Barry R. Smith, Oct 24 2015
A000013(n) = Sum_{d|n} a(d). - Robert A. Russell, Jun 09 2019
G.f.: 1 + Sum_{k>=1} mu(2*k)*log(1 - 2*x^k)/(2*k). - Ilya Gutkovskiy, Nov 11 2019

Extensions

Additional comments from Frank Ruskey, Dec 13 1999

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

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

A259793 Number of partitions of n^4 into fourth powers.

Original entry on oeis.org

1, 1, 2, 7, 36, 253, 1886, 14800, 118238, 955639, 7750456, 62777522, 506272363, 4056634991, 32252971687, 254209569990, 1985108901344, 15352968310930, 117579612410477, 891596419221856, 6694250497509934, 49768995849050468, 366423320400440927, 2671969175372760210
Offset: 0

Views

Author

N. J. A. Sloane, Jul 06 2015

Keywords

Crossrefs

A row of the array in A259799.

Programs

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

Formula

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

Extensions

More terms from Alois P. Heinz, Jul 10 2015

A259799 Array read by antidiagonals upwards: T(n,k) = number of partitions of k^n into n-th powers (n>=1, k>=0).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 4, 5, 1, 1, 2, 5, 8, 7, 1, 1, 2, 7, 17, 19, 11, 1, 1, 2, 9, 36, 62, 43, 15, 1, 1, 2, 13, 88, 253, 258, 98, 22, 1, 1, 2, 19, 218, 1104, 1886, 1050, 220, 30, 1, 1, 2, 27, 550, 5082, 15772, 14800, 4365, 504, 42, 1, 1, 2, 40, 1413, 24119, 140549, 241582, 118238, 18012, 1116, 56
Offset: 1

Views

Author

N. J. A. Sloane, Jul 06 2015

Keywords

Examples

			The array begins:
  1, 1, 2, 3, 5, 7, 11, 15, 22, 30, ...
  1, 1, 2, 4, 8, 19, 43, 98, 220, 504, ...
  1, 1, 2, 5, 17, 62, 258, 1050, 4365, 18012, ...
  1, 1, 2, 7, 36, 253, 1886, 14800, 118238, ...
  1, 1, 2, 9, 88, 1104, 15772, 241582, ...
  ...
		

Crossrefs

T(n,n) gives A331402.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0 or i=1, 1,
          `if`(i=2, 1+iquo(n, i^k), b(n, i-1, k)+
          `if`(i^k>n, 0, b(n-i^k, i, k))))
        end:
    T:= (n, k)-> b(k^n, k, n):
    seq(seq(T(d-k, k), k=0..d-1), d=1..12);  # Alois P. Heinz, Jul 10 2015
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n==0 || i==1, 1, If[i==2, 1+Quotient[n, i^k], b[n, i-1, k] + If[i^k>n, 0, b[n-i^k, i, k]]]]; T[n_, k_] := b[k^n, k, n]; Table[ Table[ T[d-k, k], {k, 0, d-1}], {d, 1, 12}] // Flatten (* Jean-François Alcover, Jul 15 2015, after Alois P. Heinz *)

Extensions

More terms from Alois P. Heinz, Jul 10 2015

A006788 a(n) = floor(2^(n-1)/n).

Original entry on oeis.org

1, 1, 1, 2, 3, 5, 9, 16, 28, 51, 93, 170, 315, 585, 1092, 2048, 3855, 7281, 13797, 26214, 49932, 95325, 182361, 349525, 671088, 1290555, 2485513, 4793490, 9256395, 17895697, 34636833, 67108864, 130150524, 252645135, 490853405, 954437176, 1857283155, 3616814565
Offset: 1

Views

Author

Keywords

Comments

Very close to A000048. [Fisher, 1989]
This is the number of nested polygons needed to produce a graph that is always concave, see the MathWorld article. - Jon Perry, Sep 15 2002

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

A027600 Number of partitions of 4^n-1 into n-th powers.

Original entry on oeis.org

3, 6, 15, 34, 86, 216, 548, 1411, 3677, 9620, 25295, 66735, 176535, 467996, 1242482, 3302288, 8784052, 23380492, 62260290, 165852146, 441922281, 1177759338, 3139292092, 8368647137, 22310795823, 59484268376
Offset: 1

Views

Author

Keywords

Crossrefs

Cf. A027601, = A027609(n) + 1.

A027609 Least non-partition into positive n-th powers.

Original entry on oeis.org

4, 7, 16, 35, 87, 217, 549, 1412, 3678, 9621, 25296, 66736, 176536, 467997, 1242483, 3302289, 8784053, 23380493, 62260291, 165852147, 441922282, 1177759339, 3139292093, 8368647138, 22310795824, 59484268377, 158602342445, 422894828107, 1127630055692
Offset: 1

Views

Author

Keywords

Crossrefs

Formula

a(n) = (partition of 4^n into n-th powers) - 1

Extensions

a(27)-a(29) from Alois P. Heinz, Jul 10 2015

A259794 Number of partitions of n^5 into fifth powers.

Original entry on oeis.org

1, 1, 2, 9, 88, 1104, 15772, 241582, 3869852, 63689650, 1065023018, 17948615155, 303219868652, 5116273886322, 86004191773864, 1437703756689091, 23869446608034827, 393225674878151704, 6423761195925513669, 104014146020398166139, 1668870762057827073994
Offset: 0

Views

Author

N. J. A. Sloane, Jul 06 2015

Keywords

Crossrefs

A row of the array in A259799.

Formula

a(n) = [x^(n^5)] Product_{k>=1} 1/(1 - x^(k^5)). - Ilya Gutkovskiy, Jan 29 2018

Extensions

More terms from Alois P. Heinz, Jul 10 2015
Showing 1-10 of 14 results. Next