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

A001935 Number of partitions with no even part repeated; partitions of n in which no parts are multiples of 4.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 16, 22, 29, 38, 50, 64, 82, 105, 132, 166, 208, 258, 320, 395, 484, 592, 722, 876, 1060, 1280, 1539, 1846, 2210, 2636, 3138, 3728, 4416, 5222, 6163, 7256, 8528, 10006, 11716, 13696, 15986, 18624, 21666, 25169, 29190, 33808, 39104, 45164
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n where no part appears more than three times.
a(n) satisfies Euler's pentagonal number (A001318) theorem, unless n is in A062717 (see Fink et al.).
Also number of partitions of n in which the least part and the differences between consecutive parts is at most 3. Example: a(5)=6 because we have [4,1], [3,2], [3,1,1], [2,2,1], [2,1,1,1] and [1,1,1,1,1]. - Emeric Deutsch, Apr 19 2006
Equals A000009 convolved with its aerated variant, = polcoeff A000009 * A000041 * A010054 (with alternate signs). - Gary W. Adamson, Mar 16 2010
Equals left border of triangle A174715. - Gary W. Adamson, Mar 27 2010
The Cayley reference is actually to A083365. - Michael Somos, Feb 24 2011
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A000009 and A035457. - Vaclav Kotesovec, Aug 23 2015
Convolution inverse is A082303. - Michael Somos, Sep 30 2017
The g.f. in the form Sum_{n >= 0} x^(n*(n+1)/2) * Product_{k = 1..n} (1+x^k)/(1-x^k) = Sum_{n >= 0} x^(n*(n+1)/2) * Product_{k = 1..n} (1+x^k)/(1+x^k-2*x^k) == Sum_{n >= 0} x^(n*(n+1)/2) (mod 2). It follows that a(n) is odd iff n = k*(k + 1)/2 for some nonnegative integer k. Cf. A333374. - Peter Bala, Jan 08 2025

Examples

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 6*x^5 + 9*x^6 + 12*x^7 + 16*x^8 + 22*x^9 + ...
G.f. = q + q^9 + 2*q^17 + 3*q^25 + 4*q^33 + 6*q^41 + 9*q^49 + 12*q^57 + 16*q^65 + 22*q^73 + ...
a(5)=6 because we have [5], [4,1], [3,2], [3,1,1], [2,1,1,1] and [1,1,1,1,1].
		

References

  • A. Cayley, A memoir on the transformation of elliptic functions, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 9, p. 128.
  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (2.5.2).
  • M. D. Hirschhorn, The Power of q, Springer, 2017. See ped page 303ff.
  • R. Honsberger, Mathematical Gems III, M.A.A., 1985, p. 241.
  • 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, A010054. - Gary W. Adamson, Mar 16 2010
Cf. A174715. - Gary W. Adamson, Mar 27 2010
Cf. A082303.
Number of r-regular partitions for r = 2 through 12: A000009, A000726, A001935, A035959, A219601, A035985, A261775, A104502, A261776, A328545, A328546.

Programs

  • Haskell
    a001935 = p a042968_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, Sep 02 2012
  • Maple
    g:=product((1+x^j)*(1+x^(2*j)),j=1..50): gser:=series(g,x=0,55): seq(coeff(gser,x,n),n=0..48); # Emeric Deutsch, Apr 19 2006
    # second Maple program:
    with(numtheory):
    a:= proc(n) option remember; `if`(n=0, 1, add(a(n-j)*add(
         `if`(irem(d, 4)=0, 0, d), d=divisors(j)), j=1..n)/n)
        end:
    seq(a(n), n=0..50);  # Alois P. Heinz, Nov 24 2015
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, q] / EllipticTheta[ 2, Pi/4, q^(1/2)] / (16 q)^(1/8), {q, 0, n}]; (* Michael Somos, Jul 11 2011 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 - x^k, {k, 4, n, 4}] / Product[ 1 - x^k, {k, n}], {x, 0, n}]; (* Michael Somos, Jul 08 2011 *)
    CoefficientList[Series[Product[1+x^j+x^(2j)+x^(3j), {j,1,48}], {x,0,48}],x] (* Jean-François Alcover, May 26 2011, after Jon Perry *)
    QP = QPochhammer; CoefficientList[QP[q^4]/QP[q] + O[q]^50, q] (* Jean-François Alcover, Nov 24 2015 *)
    a[0] = 1; a[n_] := a[n] = Sum[a[n-j] DivisorSum[j, If[Divisible[#, 4], 0, #]&], {j, 1, n}]/n; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Feb 19 2016, after Alois P. Heinz *)
    Table[Count[IntegerPartitions@n, x_ /; ! MemberQ [Mod[x, 4], 0, 2] ], {n, 0, 49}] (* Robert Price, Jul 28 2020 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( eta(x^4 + x * O(x^n)) / eta(x + x * O(x^n)), n))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint( 8*n + 1) - 1)\2, prod(i=1, k, (1 + x^i) / (x^-i - 1), 1 + x * O(x^n))), n))}; /* Michael Somos, Jun 01 2004 */
    
  • PARI
    {a(n)=polcoeff(exp(sum(m=1, n+1, x^m/(1+(-x)^m+x*O(x^n))/m)),n)} \\ Paul D. Hanna, Jul 24 2013
    

Formula

Euler transform of period 4 sequence [ 1, 1, 1, 0, ...].
Expansion of q^(-1/8) * eta(q^4) / eta(q) in powers of q. - Michael Somos, Mar 19 2004
Expansion of psi(-x) / phi(-x) = psi(x) / phi(-x^2) = psi(x^2) / psi(-x) = chi(x) / chi(-x^2)^2 = 1 / (chi(x) * chi(-x)^2) = 1 / (chi(-x) * chi(-x^2)) = f(-x^4) / f(-x) in powers of x where phi(), psi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Jul 08 2011
G.f.: Product(j>=1, 1 + x^j + x^(2*j) + x^(3*j)). - Jon Perry, Mar 30 2004
G.f.: Product_{k>=1} (1+x^k)^(2-k%2). - Jon Perry, May 05 2005
G.f.: Product_{k>0} (1 + x^(2*k)) / (1 - x^(2*k-1)) = 1 + Sum_{k>0}(Product_{i=1..k} (x^i + 1) / (x^-i - 1)).
G.f.: Sum_{n>=0} ( x^(n*(n+1)/2) * Product_{k=1..n} (1+x^k)/(1-x^k) ). - Joerg Arndt, Apr 07 2011
G.f.: P(x^4)/P(x) where P(x) = Product_{k>=1} 1-x^k. - Joerg Arndt, Jun 21 2011
A083365(n) = (-1)^n a(n). Convolution square is A001936. a(n) = A098491(n) + A098492(n). a(2*n) = A081055(n). a(2*n + 1) = A081056(n).
G.f.: (1+ 1/G(0))/2, where G(k) = 1 - x^(2*k+1) - x^(2*k+1)/(1 + x^(2*k+2) + x^(2*k+2)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Jul 03 2013
G.f.: exp( Sum_{n>=1} (x^n/n) / (1 + (-x)^n) ). - Paul D. Hanna, Jul 24 2013
a(n) ~ Pi * BesselI(1, sqrt(8*n + 1)*Pi/4) / (2*sqrt(8*n + 1)) ~ exp(Pi*sqrt(n/2)) / (4 * (2*n)^(3/4)) * (1 + (Pi/(16*sqrt(2)) - 3/(4*Pi*sqrt(2))) / sqrt(n) + (Pi^2/1024 - 15/(64*Pi^2) - 15/128) / n). - Vaclav Kotesovec, Aug 23 2015, extended Jan 14 2017
a(n) = (1/n)*Sum_{k=1..n} A046897(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 25 2017
G.f. is a period 1 Fourier series which satisfies f(-1 / (256 t)) = 1/2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A082303. - Michael Somos, Sep 30 2017

Extensions

More terms from James Sellers

A042968 Numbers not divisible by 4.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98, 99, 101, 102
Offset: 1

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Equivalently, numbers whose square part is odd. Cf. A028982. - Peter Munn, Jul 14 2020
More generally the sequence of numbers not divisible by some fixed integer m >= 2 is given by a(n,m) = 1 + n + floor(n/(m-1)). - Benoit Cloitre, Jul 11 2009
Also a(n,m) = floor((m*n-1)/(m-1)) [with offset 1]. - Gary Detlefs, May 14 2011
Numbers not having more even than odd divisors: A048272(a(n)) >= 0. - Reinhard Zumkeller, Jan 21 2012
Extending the comments of Benoit Cloitre (Jul 11 2009) and Gary Detlefs (May 14 2011), the g.f. is A(m,x) = (1-x^m) / ((1-x^(m-1))*(1-x)^2) where m >= 2 is fixed. - Werner Schulte, Apr 26 2018

Examples

			G.f. = 1 + 2*x + 3*x^2 + 5*x^3 + 6*x^4 + 7*x^5 + 9*x^6 + 10*x^7 + 11*x^8 + ... - _Michael Somos_, Jun 17 2018
		

Crossrefs

Cf. A071619 (partial sums); A008586 (complement).
Numbers that are congruent to {k0,k1,k2} mod 4: A004772, A004773, A042965, a(n).

Programs

  • Haskell
    a042968 = (`div` 3) . (subtract 1) . (* 4)
    a042968_list = filter ((/= 0) . (`mod` 4)) [1..]
    -- Reinhard Zumkeller, Sep 02 2012
    
  • Magma
    [n+1+Floor(n/3): n in [0..80]]; // Vincenzo Librandi, Aug 03 2015
    
  • Maple
    seq(n+floor((n-1)/3), n=1..80); # Muniru A Asiru, Feb 17 2019
  • Mathematica
    Select[Table[n,{n,200}], Mod[#,4] != 0&] (* Vladimir Joseph Stephan Orlovsky, Feb 18 2011 *)
    LinearRecurrence[{1,0,1,-1},{1,2,3,5},80]  (* or *) Drop[Range[110],{4,-1,4}] (* Harvey P. Dale, Jan 07 2023 *)
  • PARI
    {a(n) = 1 + n + n\3};
    
  • Python
    def A042968(n): return n+(n-1)//3 # Chai Wah Wu, Apr 13 2025
  • Sage
    [1+n+floor(n/3) for n in (0..80)] # G. C. Greubel, Feb 17 2019
    

Formula

a(n) = a(n-1) + a(n-3) - a(n-4).
a(n) = a(n-3) + 4, with a(1) = 1.
G.f.: x * (1+x) * (1+x^2) / ( (1+x+x^2)*(1-x)^2 ). - Michael Somos, Jan 12 2000
A064680(A064680(a(n))) = a(n). - Reinhard Zumkeller, Oct 19 2001
Nearest integer to (Sum_{k>n} 1/k^4)/(Sum_{k>n} 1/k^5). - Benoit Cloitre, Jun 12 2003
a(n) = n + 1 + floor(n/3). - Benoit Cloitre, Jul 11 2009
a(n) = floor((4*n+3)/3). - Gary Detlefs, May 14 2011
A214546(a(n)) >= 0 for n > 0. - Reinhard Zumkeller, Jul 20 2012
a(n) = 2*n - ceiling(2*n/3) + 1. - Arkadiusz Wesolowski, Sep 21 2012
Sum_{k=0..n} a(n) = A071619(n+1). - L. Edson Jeffery, Jul 30 2014
The g.f. A(x) satisfies x*A(x)^2 = (B(x)/x)^2 + (B(x)/x), where B(x) is the o.g.f. of A042965. - Peter Bala, Apr 12 2017
a(n) = (12*n + 6 + 3*cos(2*n*Pi/3) + sqrt(3)*sin(2*n*Pi/3))/9. - Wesley Ivan Hurt, Sep 30 2017
Euler transform of length 4 sequence [2, 0, 1, -1]. - Michael Somos, Jun 17 2018
a(n) = -a(-1-n) for all n in Z. - Michael Somos, Jun 17 2018
E.g.f.: (2/3)*exp(x)*(1 + 2*x) + (1/9)*exp(-x/2)*(3*cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2)). - Stefano Spezia, Nov 16 2019
a(n) = (12*n + 6 + w^(2*n)*(w + 2) - w^n*(w - 1))/9 where w = (-1 + sqrt(-3))/2. - Guenther Schrack, Jun 07 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = (2*sqrt(2)-1)*Pi/8. - Amiram Eldar, Dec 05 2021

Extensions

Edited by Peter Munn, Nov 16 2019
I restored my original (1999) definition and offset, which in the intervening 21 years had been lost. - N. J. A. Sloane, Jun 12 2021

A003105 Schur's 1926 partition theorem: number of partitions of n into parts 6n+1 or 6n-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 3, 3, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20, 23, 26, 30, 34, 38, 42, 47, 53, 60, 67, 74, 82, 91, 102, 114, 126, 139, 153, 169, 187, 207, 228, 250, 274, 301, 331, 364, 399, 436, 476, 520, 569, 622, 679, 739, 804, 875, 953, 1038, 1128, 1224, 1327
Offset: 0

Views

Author

Keywords

Comments

There are many (at least 8) equivalent definitions of this sequence (besides the comments below, see also Schur, Alladi, Andrews). - N. J. A. Sloane, Jun 17 2011
Coefficients of replicable function number 72e. - N. J. A. Sloane, Jun 10 2015
Also number of partitions of n into odd parts in which no part appears more than twice, cf. A070048 and A096938. - Vladeta Jovovic, Jan 18 2005
Also number of partitions of n into distinct parts congruent to 1 or 2 modulo 3. (Follows from second g.f.) - N. Sato, Jul 20 2005
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A262928 and A261612. - Vaclav Kotesovec, Jan 13 2017
Convolution of A109702 and A109701. - Vaclav Kotesovec, Jan 21 2017

Examples

			G.f: A(x) = 1 + x + x^2 + x^3 + x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 3*x^8 + 3*x^9 + 4*x^10 + ...
T72e = 1/q + q^11 + q^23 + q^35 + q^47 + 2*q^59 + 2*q^71 + 3*q^83 + ...
The logarithm of the g.f. begins:
log(A(x)) = x + x^2/2 + x^3/3 + x^4/4 + 6*x^5/5 + x^6/6 + 8*x^7/7 + x^8/8 + x^9/9 + 6*x^10/10 + 12*x^11/11 + x^12/12 + ... + A186099(n)*x^n/n + ... . - _Paul D. Hanna_, Feb 17 2013
		

References

  • K. Alladi, Refinements of Rogers-Ramanujan type identities. In Special Functions, q-Series and Related Topics (Toronto, ON, 1995), 1-35, Fields Inst. Commun., 14, Amer. Math. Soc., Providence, RI, 1997.
  • G. E. Andrews, Schur's theorem, partitions with odd parts and the Al-Salam-Carlitz polynomials. In q-Series From a Contemporary Perspective (South Hadley, MA, 1998), 45-56, Contemp. Math., 254, Amer. Math. Soc., Providence, RI, 2000.
  • H. P. Robinson, Letter to N. J. A. Sloane, Jan 04 1974.
  • I. Schur, Zur Additiven Zahlentheorie, Ges. Abh., Vol. 2, Springer, pp. 43-50.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a003105 n = p 1 n where
       p k m | m == 0 = 1 | m < k = 0 | otherwise = q k (m-k) + p (k+2) m
       q k m | m == 0 = 1 | m < k = 0 | otherwise = p (k+2) (m-k) + p (k+2) m
    -- Reinhard Zumkeller, Nov 12 2011
  • Maple
    with(combinat);
    A:=proc(n) local i, j, t3, t2, t1;
        t2:=0;
        t1:=firstpart(n);
        for j from 1 to numbpart(n)+2 do
            t3:=1;
            for i from 1 to nops(t1) do
                if (t1[i] mod 6) <> 1 and (t1[i] mod 6) <> 5 then t3:=0; fi;
            od;
            if t3=1 then t2:=t2+1; fi;
            if nops(t1) = 1 then RETURN(t2); fi;
            t1:=nextpart(t1);
        od;
    end;
    # brute-force Maple program from N. J. A. Sloane, Jun 17 2011
  • Mathematica
    max = 63; f[x_] := 1/Product[1 - x^k + x^(2k), {k, 0, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 01 2011, after Vladeta Jovovic *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] / QPochhammer[ -x^3, x^3], {x, 0, n}]; (* Michael Somos, Jul 05 2014 *)
    nmax = 100; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[If[Mod[k, 3] != 0, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}];], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 13 2017 *)
    nmax = 63; kmax = nmax/6;
    s = Flatten[{Range[0, kmax]*6 + 1}~Join~{Range[kmax]*6 - 1}];
    Table[Count[IntegerPartitions@n, x_ /; SubsetQ[s, x]], {n, 0, nmax}] (* Robert Price, Jul 31 2020 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A) / (eta(x + A) * eta(x^6 + A)), n))}; /* Michael Somos, Jan 09 2005 */
    
  • PARI
    {S(n,x)=sumdiv(n,d,d*(1-x^d)^(n/d))}
    {a(n)=polcoeff(exp(sum(k=1,n,S(k,x)*x^k/k)+x*O(x^n)),n)}
    for(n=0,60,print1(a(n),", "))
    /* Paul D. Hanna, Feb 17 2013 */
    

Formula

G.f.: 1/Product_{k>=0} (1-x^(6*k+1))*(1-x^(6*k+5)) = Product_{k>=0} (1+x^(3*k+1))*(1+x^(3*k+2)) = 1/Product_{k>=0} (1-x^k+x^(2*k)). - Vladeta Jovovic, Jun 08 2003
Expansion of chi(-x^3) / chi(-x) in powers of x where chi() is a Ramanujan theta function. - Michael Somos, Mar 04 2012
Expansion of f(x, x^2) / f(-x^3) = f(-x^6) / f(-x, -x^5) in powers of x where f() is Ramanujan theta function. - Michael Somos, Jul 05 2014
Expansion of q^(1/12) * eta(q^2) * eta(q^3) / (eta(q) * eta(q^6)) in powers of q. - Michael Somos, Jan 09 2005
Euler transform of period 6 sequence [1, 0, 0, 0, 1, 0, ...]. - Michael Somos, Jan 09 2005
Given g.f. A(x), then B(q) = (A(q^12) / q)^4 satisfies 0 = f(B(q), B(q^2)) where f(u, v) = u*v^4 + (1 - u^3) * v^3 + 6*u^2*v^2 + (u^4 - u)*v + u^3. - Michael Somos, Jan 09 2005
The logarithmic derivative equals A186099. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} A186099(n) * x^n/n ) where A186099(n) = sum of divisors of n congruent to 1 or 5 mod 6. - Paul D. Hanna, Feb 17 2013
G.f.: exp( Sum_{n>=1} S(n,x) * x^n/n ) where S(n,x) = Sum_{d|n} d*(1-x^d)^(n/d). - Paul D. Hanna, Feb 17 2013
a(n) ~ Pi*sqrt(2) / sqrt(3*(12*n-1)) * BesselI(1, Pi*sqrt(12*n-1) / (3*sqrt(6))) ~ exp(Pi*sqrt(2*n)/3) / (2^(5/4) * sqrt(3) * n^(3/4)) * (1 - (9/(8*Pi) + Pi/36)/sqrt(2*n) + (5 - 135/(4*Pi^2) + Pi^2/81)/(64*n)). - Vaclav Kotesovec, Aug 23 2015, extended Jan 09 2017
a(n) = (1/n)*Sum_{k=1..n} A186099(k)*a(n-k), a(0) = 1. - Seiichi Manyama, Mar 21 2017

Extensions

More terms from Vladeta Jovovic, Jun 08 2003

A096938 McKay-Thompson series of class 60F for the Monster group.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 3, 4, 4, 6, 7, 8, 10, 12, 14, 16, 19, 22, 26, 30, 35, 41, 47, 54, 62, 70, 80, 92, 104, 118, 135, 152, 171, 194, 218, 244, 275, 308, 344, 386, 432, 481, 537, 598, 664, 738, 819, 908, 1006, 1114, 1232, 1362, 1503, 1658, 1828, 2012, 2214, 2436, 2676
Offset: 0

Views

Author

Noureddine Chair, Aug 18 2004

Keywords

Comments

The inverted graded parafermionic partition function.
Also number of partitions of n into parts congruent to {1,3,7,9} mod 10. Also number of partitions of n into odd parts parts in which no part appears more than 4 times.
Number of partitions of n into distinct parts in which no part is a multiple of 5.
This generating function is a generalization of the sequences A003105 and A006950. It arose in my recent work on partial supersymmetry in writing the graded parafermionic partition function in which I obtained a more general formula.

Examples

			a(8)=4, the number of partitions into distinct parts that exclude the number 5 because we can write 8=7+1=6+2=4+3+1.
T60F = 1/q + q^5 + q^11 + 2*q^17 + 2*q^23 + 2*q^29 + 3*q^35 + 4*q^41 +...
		

References

  • T. M. Apostol, An Introduction to Analytic Number Theory, Springer-Verlag, NY, 1976

Crossrefs

Cf. A133563.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A261770 (m=6), A097793 (m=7), A261771 (m=8), A112193 (m=9), A261772 (m=10).

Programs

  • Maple
    series(product(1/(1-x^k+x^(2*k)-x^(3*k)+x^(4*k)),k+1..150),x=0,100);
  • Mathematica
    CoefficientList[ Series[ Product[1/(1 - x^k + x^(2k) - x^(3k) + x^(4k)), {k, 70}], {x, 0, 60}], x] (* Robert G. Wilson v, Aug 19 2004 *)
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(5*k)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 31 2015 *)
    QP = QPochhammer; s = QP[q^2]*(QP[q^5]/(QP[q]*QP[q^10])) + O[q]^60; CoefficientList[s, q] (* Jean-François Alcover, Nov 12 2015 *)
  • PARI
    {a(n)=local(A); if(n<0,0, A=x*O(x^n); polcoeff( eta(x^2+A)*eta(x^5+A)/eta(x+A)/eta(x^10+A), n))} /* Michael Somos, Jan 18 2005 */

Formula

Euler transform of period 10 sequence [1, 0, 1, 0, 0, 0, 1, 0, 1, 0, ...]. - Vladeta Jovovic, Aug 19 2004
Expansion of q^(1/6)eta(q^2)eta(q^5)/(eta(q)eta(q^10)) in powers of q.
Given g.f. A(x), then B(x)=(A(x^6)/x)^2 satisfies 0=f(B(x), B(x^2)) where f(u, v)=(u^3+v^3)(1+uv)-uv(1-uv)^2. - Michael Somos, Jan 18 2005
G.f.: 1/product_{k>=1} (1-x^k+x^(2*k)-x^(3*k)+x^(4*k)) = 1/Product_{k>0} P10(x^k) where P10 is the 10th cyclotomic polynomial.
a(n) ~ exp(2*Pi*sqrt(n/15)) / (2 * 15^(1/4) * n^(3/4)) * (1 - (3*sqrt(15)/(16*Pi) + Pi/(6*sqrt(15))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017

Extensions

Definition corrected by Vladeta Jovovic, Aug 19 2004
More terms from Robert G. Wilson v, Aug 19 2004

A097793 McKay-Thompson series of class 56B for the Monster group.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 4, 5, 7, 8, 10, 12, 14, 17, 21, 24, 28, 34, 39, 46, 53, 61, 71, 82, 94, 108, 124, 142, 162, 185, 210, 238, 271, 306, 345, 390, 439, 494, 556, 623, 698, 783, 875, 977, 1092, 1216, 1354, 1508, 1674, 1859, 2064, 2286, 2532, 2803, 3098, 3424
Offset: 0

Views

Author

Michael Somos, Aug 24 2004

Keywords

Comments

Number of partitions of n into distinct parts not divisible by 7.
Also McKay-Thompson series of class 56C for Monster. - Michel Marcus, Feb 19 2014

Examples

			1 + x + x^2 + 2*x^3 + 2*x^4 + 3*x^5 + 4*x^6 + 4*x^7 + 5*x^8 + 7*x^9 + 8*x^10 +...
T56B = 1/q + q^3 + q^7 + 2q^11 + 2q^15 + 3q^19 + 4q^23 + 4q^27 +...
		

Crossrefs

Cf. A113297.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A096938 (m=5), A261770 (m=6), A261771 (m=8), A112193 (m=9), A261772 (m=10).

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(7*k)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 31 2015 *)
    QP = QPochhammer; s = QP[q^2]*(QP[q^7]/(QP[q]*QP[q^14])) + O[q]^60; CoefficientList[s, q] (* Jean-François Alcover, Nov 12 2015 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( prod( k=1, n, 1 + x^k, 1 + A) / prod( k=1, n\7, 1 + x^(7*k), 1 + A), n))}
    
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^7 + A) / (eta(x + A) * eta(x^14 + A)), n))}

Formula

Euler transform of period 14 sequence [ 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, ...].
Expansion of q^(1/4) * eta(q^2) * eta(q^7) / (eta(q) * eta(q^14)) in powers of q.
G.f.: Product_{k>0} (1 + x^k) / (1 + x^(7*k)).
a(n) ~ exp(Pi*sqrt(2*n/7)) / (2 * 14^(1/4) * n^(3/4)) * (1 - (3*sqrt(7)/ (8*Pi*sqrt(2)) + Pi/(4*sqrt(14))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017

A112193 Coefficients of replicable function number "54b".

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 5, 6, 7, 9, 11, 13, 16, 19, 23, 27, 32, 38, 44, 52, 61, 71, 82, 95, 110, 127, 145, 167, 191, 218, 249, 283, 322, 365, 414, 469, 529, 597, 673, 757, 851, 955, 1071, 1199, 1341, 1499, 1673, 1865, 2078, 2313, 2572, 2857, 3171, 3517, 3897
Offset: 0

Views

Author

Michael Somos, Aug 28 2005

Keywords

Comments

a(n) is the number of partitions of n into distinct parts where no part is a multiple of 9. - Joerg Arndt, Aug 31 2015
In general, if m > 1 and g.f. = Product_{k>=1} (1 + x^k)/(1 + x^(m*k)), then a(n) ~ exp(Pi*sqrt((m-1)*n/(3*m))) * (m-1)^(1/4) / (2^(3/2) * 3^(1/4) * m^(1/4) * n^(3/4)). - Vaclav Kotesovec, Aug 31 2015

Examples

			G.f. = 1 + x + x^2 + 2*x^3 + 2*x^4 + 3*x^5 + 4*x^6 + 5*x^7 + 6*x^8 + ... _Michael Somos_, Oct 06 2019
G.f. = q^-1 + q^2 + q^5 + 2*q^8 + 2*q^11 + 3*q^14 + 4*q^17 + 5*q^20 + ...
		

Crossrefs

Cf. A261733.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A096938 (m=5), A261770 (m=6), A097793 (m=7), A261771 (m=8), A261772 (m=10).

Programs

  • Maple
    b:= proc(n, i) option remember;  local r;
          `if`(2*n>i*(i+1)-(j-> 9*j*(j+1))(iquo(i, 9, 'r')), 0,
          `if`(n=0, 1, b(n, i-1)+`if`(i>n or r=0, 0, b(n-i, i-1))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, Aug 31 2015
  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(9*k)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 31 2015 *)
    b[n_, i_] := b[n, i] = Module[{q, r}, {q, r} = QuotientRemainder[i, 9]; If[2*n > i*(i+1) - 9*q*(q+1), 0, If[n == 0, 1, b[n, i-1] + If[i>n || r == 0, 0, b[n-i, i-1]]]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Oct 07 2016, after Alois P. Heinz *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ q^2] QPochhammer[ q^9] / (QPochhammer[ q] QPochhammer[ q^18]), {q, 0, n}]; (* Michael Somos, Oct 06 2019 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^9 + A) / (eta(x + A) * eta(x^18 + A)), n))}; /* Michael Somos, Oct 06 2019 */

Formula

a(n) ~ exp(2*Pi*sqrt(2*n/3)/3) / (6^(3/4) * n^(3/4)) * (1 - (9*sqrt(3)/ (16*Pi*sqrt(2)) + sqrt(2)*Pi/(9*sqrt(3))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017
From Michael Somos, Oct 06 2019: (Start)
Expansion of q^(1/3) * eta(q^2) * eta(q^9) / (eta(q) * eta(q^18)) in powers of q.
Euler transform of period 18 sequence [1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...].
G.f. is a period 1 Fourier Series which satisifies f(-1 / (18 t)) = f(t) where q = exp(2 Pi i t).
Given g.f. A(x), then B(q) = A(q^3) / q satisfies 0 = f(B(q), B(q^2)) where f(u, v) = (1 + u*v) * (u^3 + v^3) - u*v * (1 + u^2*v^2).
Given g.f. A(x), then B(q) = A(q^3) / q satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = (w^2 - v) * (u^2 - v) - 2*u*v*w.
Convolution inverse of A261733.
(End)

A261770 Expansion of Product_{k>=1} (1 + x^k) / (1 + x^(6*k)).

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 8, 9, 11, 13, 16, 19, 22, 26, 30, 35, 41, 47, 55, 63, 73, 84, 96, 110, 125, 143, 162, 184, 208, 235, 266, 300, 338, 380, 427, 479, 536, 600, 670, 748, 834, 929, 1034, 1149, 1277, 1417, 1571, 1740, 1925, 2129, 2351, 2596, 2863
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 31 2015

Keywords

Comments

a(n) is the number of partitions of n into distinct parts where no part is a multiple of 6. - Joerg Arndt, Aug 31 2015

Crossrefs

Cf. A261736.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A096938 (m=5), A097793 (m=7), A261771 (m=8), A112193 (m=9), A261772 (m=10).
Column k=6 of A290307.

Programs

  • Maple
    b:= proc(n, i) option remember;  local r;
          `if`(2*n>i*(i+1)-(j-> 6*j*(j+1))(iquo(i, 6, 'r')), 0,
          `if`(n=0, 1, b(n, i-1)+`if`(i>n or r=0, 0, b(n-i, i-1))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, Aug 31 2015
  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(6*k)), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(Pi*sqrt(5*n/2)/3) * 5^(1/4) / (2^(7/4) * sqrt(3) * n^(3/4)) * (1 - (9/(4*Pi*sqrt(10)) + 5*Pi*sqrt(5/2)/144) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017
G.f.: Product_{k>=1} (1 - x^(12*k-6))/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 07 2017

A261771 Expansion of Product_{k>=1} (1 + x^k) / (1 + x^(8*k)).

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 5, 5, 7, 9, 10, 13, 15, 18, 22, 26, 30, 36, 42, 49, 58, 67, 77, 89, 103, 118, 136, 156, 177, 203, 231, 263, 299, 338, 383, 433, 489, 550, 620, 696, 781, 877, 981, 1097, 1227, 1369, 1526, 1702, 1893, 2104, 2339, 2595, 2877, 3189, 3530
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 31 2015

Keywords

Comments

a(n) is the number of partitions of n into distinct parts where no part is a multiple of 8.

Crossrefs

Cf. A261735.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A096938 (m=5), A261770 (m=6), A097793 (m=7), A112193 (m=9), A261772 (m=10).
Column k=8 of A290307.

Programs

  • Maple
    with(numtheory):
    a:= proc(n) option remember; `if`(n=0, 1, add(add(d*
           [0, 1, 0, 1, 0, 1, 0, 1, -1, 1, 0, 1, 0, 1, 0, 1]
           [1+irem(d, 16)], d=divisors(j))*a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..80);  # Alois P. Heinz, Aug 31 2015
  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(8*k)), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(Pi*sqrt(7*n/6)/2) * 7^(1/4) / (4 * 6^(1/4) * n^(3/4)) * (1 - (3*sqrt(3)/ (2*Pi*sqrt(14)) + 7*Pi*sqrt(7)/(96*sqrt(6))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017
G.f.: Product_{k>=1} (1 - x^(16*k-8))/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 07 2017

A261772 Expansion of Product_{k>=1} (1 + x^k) / (1 + x^(10*k)).

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 9, 11, 14, 16, 20, 24, 28, 33, 40, 46, 54, 64, 74, 86, 100, 115, 133, 154, 176, 202, 231, 263, 300, 342, 388, 440, 499, 563, 636, 718, 808, 909, 1022, 1146, 1284, 1439, 1608, 1797, 2006, 2236, 2490, 2772, 3081, 3422, 3800, 4212
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 31 2015

Keywords

Comments

a(n) is the number of partitions of n into distinct parts where no part is a multiple of 10.
In general, if m > 1 and g.f. = Product_{k>=1} (1 + x^k)/(1 + x^(m*k)), then a(n) ~ exp(Pi*sqrt((m-1)*n/(3*m))) * (m-1)^(1/4) / (2^(3/2) * 3^(1/4) * m^(1/4) * n^(3/4)) * (1 - (3*sqrt(3*m)/(8*Pi*sqrt(m-1)) + (m-1)^(3/2)*Pi/(48*sqrt(3*m))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017

Crossrefs

Cf. A145707.
Cf. A000700 (m=2), A003105 (m=3), A070048 (m=4), A096938 (m=5), A261770 (m=6), A097793 (m=7), A261771 (m=8), A112193 (m=9).
Column k=10 of A290307.

Programs

  • Maple
    b:= proc(n, i) option remember;  local r;
          `if`(2*n>i*(i+1)-(j-> 10*j*(j+1))(iquo(i, 10, 'r')), 0,
          `if`(n=0, 1, b(n, i-1)+`if`(i>n or r=0, 0, b(n-i, i-1))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, Aug 31 2015
  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) / (1 + x^(10*k)), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(Pi*sqrt(3*n/10)) * 3^(1/4) / (2^(7/4) * 5^(1/4) * n^(3/4)) * (1 - (sqrt(15)/(4*Pi*sqrt(2)) + 3*Pi*sqrt(3)/(16*sqrt(10))) / sqrt(n)). - Vaclav Kotesovec, Aug 31 2015, extended Jan 21 2017
G.f.: Product_{k>=1} (1 - x^(20*k-10))/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 07 2017

A108961 Number of partitions that are "2-close" to being self-conjugate.

Original entry on oeis.org

1, 1, 2, 3, 3, 5, 7, 9, 12, 16, 20, 26, 33, 41, 51, 64, 79, 97, 119, 144, 175, 212, 254, 305, 365, 434, 516, 612, 722, 851, 1002, 1174, 1375, 1607, 1872, 2179, 2531, 2933, 3395, 3923, 4524, 5211, 5994, 6881, 7891, 9038, 10334, 11804, 13467, 15341, 17460, 19849
Offset: 0

Views

Author

John McKay (mckay(AT)cs.concordia.ca), Jul 22 2005

Keywords

Comments

Let (a1,a2,a3,...ad:b1,b2,b3,...bd) be the Frobenius symbol for a partition pi of n. Then pi is m-close to being self-conjugate if for all k, |ak-bk| <= m.
Convolution of A070048 and A035457. - Vaclav Kotesovec, Nov 13 2016

Examples

			1 + x + 2*x^2 + 3*x^3 + 3*x^4 + 5*x^5 + 7*x^6 + 9*x^7 + 12*x^8 + 16*x^9 + 20*x^10 + ...
1/q + q^23 + 2*q^47 + 3*q^71 + 3*q^95 + 5*q^119 + 7*q^143 + 9*q^167 + 12*q^191 + ...
		

References

  • D. M. Bressoud, Analytic and combinatorial generalizations of the Rogers-Ramanujan identities, Mem. Amer. Math. Soc. 24 (1980), no. 227, 54 pp.

Crossrefs

Cf. A000700 for m=0 (self-conjugate), A070047 for m=1, A108962 for m=3, A271661 for m=4, A280937 for m=5, A280938 for m=6.

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[(1 + x^k) * (1 + x^(2*k)) / (1 + x^(4*k)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 13 2016 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^4 + A)^2 / (eta(x + A) * eta(x^8 + A)), n))} /* Michael Somos, Oct 17 2006 */

Formula

Define the Dedekind eta function = q^1/24. Product(1-q^k), k >=1. Then the number of m-close partitions is q^(1/24).(m+2)^2/(1.(2m+4)) (where m denotes eta(q^m)).
Expansion of q^(1/24) * eta(q^4)^2 / (eta(q) * eta(q^8)) in powers of q. - Michael Somos, Oct 17 2006
Expansion of chi(x^2) * chi(x) / chi(-x^2) in powers of x where chi() is a Ramanujan theta function (see A000700). - Michael Somos, Oct 17 2006 [corrected by Peter Bala, Oct 09 2023]
Euler transform of period 8 sequence [ 1, 1, 1, -1, 1, 1, 1, 0, ...]. - Michael Somos, Oct 17 2006
G.f.: Product_{k>0} (1 + x^k) * (1 + x^(2*k)) / (1 + x^(4*k)). - Michael Somos, Oct 17 2006
a(n) ~ Pi * BesselI(1, Pi * sqrt(5*(24*n-1)/2)/12) / (2*sqrt((24*n-1)/5)) ~ 5^(1/4) * exp(sqrt(5*n/3)*Pi/2) / (2^(5/2) * 3^(1/4) * n^(3/4)) * (1 -(3*sqrt(3) / (4*Pi*sqrt(5)) + Pi*sqrt(5)/(96*sqrt(3)))/sqrt(n) + (5*Pi^2/55296 - 9/(32*Pi^2) + 5/128)/n). - Vaclav Kotesovec, Nov 13 2016, extended Jan 11 2017
Showing 1-10 of 16 results. Next