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-9 of 9 results.

A067589 Numbers k such that A067588(k) is an odd number.

Original entry on oeis.org

1, 5, 7, 15, 35, 51, 57, 77, 117, 145, 155, 187, 247, 287, 301, 345, 425, 477, 495, 551, 651, 715, 737, 805, 925, 1001, 1027, 1107, 1247, 1335, 1365, 1457, 1617, 1717, 1751, 1855, 2035, 2147, 2185, 2301, 2501, 2625, 2667, 2795, 3015, 3151, 3197, 3337
Offset: 1

Views

Author

Naohiro Nomoto, Jan 31 2002

Keywords

Comments

The terms are exactly the odd pentagonal numbers; that is, they are all the odd numbers of the form k*(3*k-1)/2 where k is an integer. - James Sellers, Jun 09 2007
Apparently groups of two odd pentagonal numbers (A000326, A014632) followed by two odd 2nd pentagonal numbers (A005449), which leads to the conjectured generating function x*(x^2+4*x+1)*(x^4-2*x^3+4*x^2-2*x+1)/((x^2+1)^2*(1-x)^3). - R. J. Mathar, Jul 26 2009
Odd generalized pentagonal numbers. - Omar E. Pol, Aug 19 2011
From Peter Bala, Jan 10 2025: (Start)
The sequence terms are the exponents in the expansion of Sum_{n >= 0} x^(2*n+1)/(Product_{k = 1..2*n+1} 1 + x^(2*k+1)) = x + x^5 - x^7 - x^15 + x^35 + x^51 - x^57 - x^77 + + - - ... (follows from Berndt et al., Theorem 3.3). Cf. A193828.
For positive integer m, define b_m(n) = Sum_{k = 1..n} k^(2*m+1)*A000009(k)*A000009(n-k). We conjecture that
i) for odd n, b(n)/ n is an integer
ii) b(2*n)/n is an integer, which is odd iff n is a member of this sequence.
Cf. A067567. (End)

Crossrefs

Programs

  • Mathematica
    With[{nn=50},Sort[Select[Table[(n(3n-1))/2,{n,-nn,nn}],OddQ]]] (* Harvey P. Dale, Feb 16 2014 *)

Formula

Sum_{n>=1} 1/a(n) = Pi/2. - Amiram Eldar, Aug 18 2022

Extensions

Corrected by T. D. Noe, Oct 25 2006

A038348 Expansion of (1/(1-x^2))*Product_{m>=0} 1/(1-x^(2m+1)).

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 8, 11, 14, 19, 24, 31, 39, 49, 61, 76, 93, 114, 139, 168, 203, 244, 292, 348, 414, 490, 579, 682, 801, 938, 1097, 1278, 1487, 1726, 1999, 2311, 2667, 3071, 3531, 4053, 4644, 5313, 6070, 6923, 7886, 8971, 10190, 11561
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n+2 with exactly one even part. - Vladeta Jovovic, Sep 10 2003
Also, number of partitions of n with at most one even part. - Vladeta Jovovic, Sep 10 2003
Also total number of parts, counted without multiplicity, in all partitions of n into odd parts, offset 1. - Vladeta Jovovic, Mar 27 2005
a(n) = Sum_{k>=1} k*A116674(n+1,k). - Emeric Deutsch, Feb 22 2006
Equals row sums of triangle A173305. - Gary W. Adamson, Feb 15 2010
Equals partial sums of A025147 (observed by Jonathan Vos Post, proved by several correspondents).
Conjecture: The n-th derivative of Gamma(x+1) at x = 0 has a(n+1) terms. For example, d^4/dx^4_(x = 0) Gamma(x+1) = 8*eulergamma*zeta(3) + eulergamma^4 + eulergamma^2*Pi^2 + 3*Pi^4/20 which has a(5) = 4 terms. - David Ulgenes, Dec 05 2023

Examples

			From _Gus Wiseman_, Sep 23 2019: (Start)
Also the number of integer partitions of n that are strict except possibly for any number of 1's. For example, the a(1) = 1 through a(7) = 11 partitions are:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (31)    (32)     (42)      (43)
             (111)  (211)   (41)     (51)      (52)
                    (1111)  (311)    (321)     (61)
                            (2111)   (411)     (421)
                            (11111)  (3111)    (511)
                                     (21111)   (3211)
                                     (111111)  (4111)
                                               (31111)
                                               (211111)
                                               (1111111)
(End)
		

Crossrefs

Programs

  • Maple
    f:=1/(1-x^2)/product(1-x^(2*j-1),j=1..32): fser:=series(f,x=0,62): seq(coeff(fser,x,n),n=0..58); # Emeric Deutsch, Feb 22 2006
  • Mathematica
    mmax = 47; CoefficientList[ Series[ (1/(1-x^2))*Product[1/(1-x^(2m+1)), {m, 0, mmax}], {x, 0, mmax}], x] (* Jean-François Alcover, Jun 21 2011 *)
  • SageMath
    # uses[EulerTransform from A166861]
    def g(n): return n % 2 if n > 2 else 1
    a = EulerTransform(g)
    print([a(n) for n in range(48)]) # Peter Luschny, Dec 04 2020

Formula

a(n) = A036469(n) - a(n-1) = Sum_{k=0..n} (-1)^k*A036469(n-k). - Vladeta Jovovic, Sep 10 2003
a(n) = A000009(n) + a(n-2). - Vladeta Jovovic, Feb 10 2004
G.f.: 1/((1-x^2)*Product_{j>=1} (1 - x^(2*j-1))). - Emeric Deutsch, Feb 22 2006
From Vaclav Kotesovec, Aug 16 2015: (Start)
a(n) ~ (1/2) * A036469(n).
a(n) ~ 3^(1/4) * exp(Pi*sqrt(n/3)) / (4*Pi*n^(1/4)). (End)
Euler transform of the sequence [1, 1, period(1, 0)] (A266591). - Georg Fischer, Dec 04 2020

A092268 Total number of smallest parts in all partitions of n into odd parts.

Original entry on oeis.org

1, 2, 4, 5, 8, 12, 15, 20, 29, 36, 46, 61, 74, 95, 122, 145, 180, 224, 268, 328, 399, 474, 567, 682, 807, 955, 1136, 1330, 1564, 1842, 2140, 2499, 2914, 3375, 3917, 4533, 5220, 6014, 6929, 7942, 9102, 10430, 11898, 13582, 15489, 17600, 19999, 22706, 25719
Offset: 1

Views

Author

Vladeta Jovovic, Feb 16 2004

Keywords

Examples

			Partitions of 6 into odd parts are: [1,1,1,1,1,1], [1,1,1,3], [3,3], [1,5]; thus a(6)=6+3+2+1=12.
		

Crossrefs

Programs

  • Mathematica
    nmax = 50; Rest[CoefficientList[Series[Sum[(x^(2*n - 1)/(1 - x^(2*n - 1))) / Product[(1 - x^(2*k - 1)), {k, n, nmax}], {n, 1, nmax}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, Jul 06 2019 *)

Formula

G.f.: Sum((x^(2*n-1)/(1-x^(2*n-1)))/Product((1-x^(2*k-1)), k=n..infinity), n=1..infinity).
a(n) ~ 3^(1/4) * exp(Pi*sqrt(n/3)) / (2*Pi*n^(1/4)). - Vaclav Kotesovec, Jul 07 2019

Extensions

More terms from Pab Ter (pabrlos(AT)yahoo.com), May 25 2004

A116676 Number of odd parts in all partitions of n into distinct parts.

Original entry on oeis.org

0, 1, 0, 2, 2, 3, 4, 5, 8, 10, 14, 16, 22, 26, 34, 43, 54, 64, 80, 96, 116, 142, 170, 202, 242, 288, 340, 404, 474, 556, 652, 762, 886, 1034, 1198, 1389, 1606, 1852, 2132, 2454, 2814, 3224, 3690, 4214, 4804, 5478, 6228, 7072, 8028, 9094, 10290, 11635, 13134
Offset: 0

Views

Author

Emeric Deutsch, Feb 22 2006

Keywords

Comments

a(n) = Sum(k*A116675(n,k), k>=0).

Examples

			a(9) = 10 because in the partitions of 9 into distinct parts, namely, [9], [81], [72], [6,3], [6,2,1], [5,4], [5,3,1] and [4,3,2], we have a total of 10 odd parts.
		

Crossrefs

Programs

  • Maple
    f:=product(1+x^j,j=1..64)*sum(x^(2*j-1)/(1+x^(2*j-1)),j=1..35): fser:=series(f,x=0,60): seq(coeff(fser,x,n),n=0..56);
    # second Maple program:
    b:= proc(n, i) option remember; local f, g;
          if n=0 then [1, 0] elif i<1 then [0, 0]
        else f:=b(n, i-1); g:=`if`(i>n, [0, 0], b(n-i, min(n-i, i-1)));
             [f[1]+g[1], f[2]+g[2] +irem(i, 2)*g[1]]
          fi
        end:
    a:= n-> b(n, n)[2]:
    seq(a(n), n=0..60);  # Alois P. Heinz, Nov 21 2012
  • Mathematica
    b[n_, i_] := b[n, i] = Module[{f, g}, Which [n == 0, {1, 0}, i<1 , {0, 0}, True, f = b[n, i-1]; g = If[i>n, {0, 0}, b[n-i, Min[n-i, i-1]]]; {f[[1]] + g[[1]],       f[[2]] + g[[2]] + Mod[i, 2]*g[[1]]}]]; a[n_] := b[n, n][[2]]; Table [a[n], {n, 0, 60}] (* Jean-François Alcover, May 22 2015, after Alois P. Heinz *)

Formula

G.f.: product(1+x^j, j=1..infinity)*sum(x^(2j-1)/(1+x^(2j-1)), j=1..infinity).
For n > 0, a(n) = A015723(n) - A116680(n). - Vaclav Kotesovec, May 26 2018
a(n) ~ 3^(1/4) * log(2) * exp(Pi*sqrt(n/3)) / (4*Pi*n^(1/4)). - Vaclav Kotesovec, May 26 2018

A116680 Number of even parts in all partitions of n into distinct parts.

Original entry on oeis.org

0, 0, 1, 1, 1, 2, 4, 5, 5, 8, 11, 14, 18, 23, 29, 37, 44, 55, 69, 83, 102, 124, 148, 178, 213, 253, 300, 356, 421, 494, 582, 680, 793, 926, 1074, 1246, 1446, 1668, 1922, 2215, 2545, 2918, 3345, 3823, 4366, 4982, 5668, 6445, 7321, 8300, 9401, 10639, 12021, 13566
Offset: 0

Views

Author

Emeric Deutsch, Feb 22 2006

Keywords

Examples

			a(9)=8 because in the partitions of 9 into distinct parts, namely, [9], [8,1], [7,2], [6,3], [6,2,1], [5,4], [5,3,1], and [4,3,2], we have a total of 8 even parts. [edited by _Rishi Advani_, Jun 07 2019]
		

Crossrefs

Programs

  • Magma
    m:=25; R:=PowerSeriesRing(Integers(), 3*m); [0,0] cat Coefficients(R!( (&*[1+x^j: j in [1..4*m]])*(&+[x^(2*k)/(1+x^(2*k)): k in [1..2*m]]) )); // G. C. Greubel, Jun 07 2019
    
  • Maple
    f:=product(1+x^j,j=1..70)*sum(x^(2*j)/(1+x^(2*j)),j=1..40): fser:=series(f,x=0,65): seq(coeff(fser,x,n),n=0..60);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(i*(i+1)/2 p+`if`(i::odd, 0, [0, p[1]]))(b(n-i, min(n-i, i-1)))))
        end:
    a:= n-> b(n$2)[2]:
    seq(a(n), n=0..60);  # Alois P. Heinz, May 24 2022
  • Mathematica
    With[{m = 25}, CoefficientList[Series[Product[1+x^j, {j,1,4*m}]* Sum[x^(2*k)/(1+x^(2*k)), {k,1,2*m}], {x,0,3*m}], x]] (* G. C. Greubel, Jun 07 2019 *)
  • PARI
    my(m=25); my(x='x+O('x^(3*m))); concat([0, 0], Vec( prod(j=1, 4*m, 1+x^j)*sum(k=1, 2*m, x^(2*k)/(1+x^(2*k))) )) \\ G. C. Greubel, Jun 07 2019
    
  • Sage
    m = 25
    R = PowerSeriesRing(ZZ, 'x')
    x = R.gen().O(3*m)
    s = product(1+x^j for j in (1..4*m))*sum(x^(2*k)/(1+x^(2*k)) for k in (1..2*m))
    [0, 0] + s.coefficients() # G. C. Greubel, Jun 07 2019

Formula

a(n) = Sum_{k >= 0} k*A116679(n,k).
G.f.: (Product_{j >= 1} (1+x^j)) * (Sum_{k >= 1} x^(2*k)/(1+x^(2*k))).
For n > 0, a(n) = A015723(n) - A116676(n). - Vaclav Kotesovec, May 26 2018
a(n) ~ 3^(1/4) * log(2) * exp(Pi*sqrt(n/3)) / (4*Pi*n^(1/4)). - Vaclav Kotesovec, May 26 2018

A305082 G.f.: Sum_{k>=1} x^k/(1-x^k) * Product_{k>=1} (1+x^k).

Original entry on oeis.org

0, 1, 3, 5, 9, 13, 20, 28, 39, 54, 71, 94, 124, 159, 201, 258, 322, 401, 499, 613, 750, 918, 1110, 1340, 1617, 1935, 2308, 2752, 3261, 3854, 4554, 5350, 6273, 7348, 8572, 9983, 11612, 13460, 15578, 18007, 20761, 23894, 27473, 31511, 36090, 41296, 47152, 53767
Offset: 0

Views

Author

Vaclav Kotesovec, May 25 2018

Keywords

Comments

Convolution of A000005 and A000009.
Apart from initial zero this is the convolution of A341062 and A036469. - Omar E. Pol, Feb 16 2021

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Sum[x^k/(1-x^k), {k, 1, nmax}]*Product[1+x^k, {k, 1, nmax}], {x, 0, nmax}], x]
    nmax = 50; CoefficientList[Series[((Log[1-x] + QPolyGamma[0, 1, x]) * QPochhammer[-1, x]) / (2*Log[x]), {x, 0, nmax}], x]

Formula

a(n) ~ 3^(1/4)*(2*gamma + log(12*n/Pi^2)) * exp(Pi*sqrt(n/3)) / (4*Pi*n^(1/4)), where gamma is the Euler-Mascheroni constant A001620.

A341497 Number of partitions of n with exactly one repeated part and that part is odd.

Original entry on oeis.org

0, 0, 1, 1, 2, 3, 5, 7, 9, 13, 17, 23, 30, 39, 49, 63, 78, 98, 122, 150, 184, 225, 272, 329, 397, 475, 567, 676, 802, 948, 1121, 1317, 1545, 1810, 2112, 2460, 2863, 3319, 3842, 4442, 5123, 5897, 6782, 7780, 8913, 10200, 11648, 13285, 15136, 17214, 19555, 22191, 25143
Offset: 0

Views

Author

Andrew Howroyd, Feb 13 2021

Keywords

Examples

			The a(2) = 1 partition is: 1+1.
The a(3) = 1 partition is: 1+1+1.
The a(4) = 2 partitions are: 1+1+2, 1+1+1+1.
The a(5) = 3 partitions are: 1+1+3, 1+1+1+2, 1+1+1+1+1.
		

Crossrefs

Programs

  • PARI
    seq(n)={Vec(sum(k=1, (n+2)\4, x^(4*k-2)/(1 - x^(4*k-2)) + O(x*x^n)) * prod(k=1, n, 1 + x^k + O(x*x^n)), -(n+1))}

Formula

G.f.: (Sum_{k>=1} x^(4*k-2)/(1 - x^(4*k-2))) * Product_{k>=1} (1 + x^k).
a(n) = A090867(n) - A341496(n).
a(n) = A116680(n) + A341496(n).
a(n) = A341495(n) for even n; a(n) = A341494(n) for odd n.
a(n) = (A067588(n) - A116676(n))/2. - Peter Bala, Jan 13 2025

A305105 G.f.: Sum_{k>=1} x^(2*k-1)/(1-x^(2*k-1)) * Product_{k>=1} (1+x^k)/(1-x^k).

Original entry on oeis.org

0, 1, 3, 8, 17, 34, 64, 114, 195, 325, 526, 832, 1292, 1970, 2958, 4384, 6413, 9276, 13283, 18836, 26478, 36924, 51096, 70210, 95844, 130019, 175350, 235192, 313802, 416618, 550540, 724250, 948719, 1237732, 1608508, 2082600, 2686857, 3454590, 4427144, 5655652
Offset: 0

Views

Author

Vaclav Kotesovec, May 25 2018

Keywords

Comments

Convolution of A066897 and A000009.
Convolution of A067588 and A000041.
Let A(x) = Sum_{k >= 1} x^(2*k-1)/(1 - x^(2*k-1)) * Product_{k >= 1} (1 + x^k)/(1 - x^k). Then A(x) = Sum_{k >= 1} x^(2*k-1)/(1 - x^(2*k-1)) * Product_{k >= 1} (1 + x^k)/(1 + x^k - 2*x^k) == Sum_{k >= 1} x^(2*k-1)/(1 - x^(2*k-1)) (mod 2). It follows from the comment in A001227 by Juri-Stepan Gerasimov, dated Jul 17 2016, that a(n) is odd iff n is a square or twice a square. - Peter Bala, Jan 10 2025

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Sum[x^(2*k-1)/(1-x^(2*k-1)), {k, 1, nmax}] * Product[(1+x^k)/(1-x^k), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ (2*gamma + log(16*n/Pi^2)) * exp(Pi*sqrt(n)) / (16*Pi*sqrt(n)), where gamma is the Euler-Mascheroni constant A001620.

A306925 Sum over all partitions of n into distinct parts of the bitwise XOR of the parts.

Original entry on oeis.org

0, 1, 2, 6, 6, 11, 16, 35, 36, 46, 50, 84, 94, 130, 158, 285, 338, 424, 460, 616, 672, 810, 816, 1162, 1346, 1680, 1754, 2308, 2562, 3164, 3288, 4486, 5306, 6838, 7522, 9627, 11006, 13496, 14200, 17462, 19682, 24036, 25650, 30842, 33884, 40302, 41644, 48896
Offset: 0

Views

Author

Alois P. Heinz, Mar 16 2019

Keywords

Crossrefs

Programs

  • Maple
    b:= proc(n, i, r) option remember; `if`(i*(i+1)/2 `if`(i b(n$2, 0):
    seq(a(n), n=0..51);

Formula

a(n) is odd <=> n in { A067589 }.
a(n) is odd <=> A067588(n) is odd.
Showing 1-9 of 9 results.