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.

A090867 Number of partitions of n such that the set of even parts has only one element.

Original entry on oeis.org

0, 0, 1, 1, 3, 4, 6, 9, 13, 18, 23, 32, 42, 55, 69, 89, 112, 141, 175, 217, 266, 326, 396, 480, 581, 697, 834, 996, 1183, 1402, 1660, 1954, 2297, 2694, 3150, 3674, 4280, 4970, 5762, 6669, 7701, 8876, 10219, 11737, 13460, 15418, 17628, 20125, 22951, 26128, 29709
Offset: 0

Views

Author

Vladeta Jovovic, Feb 12 2004

Keywords

Comments

Conjecture: a(n) is also the difference between the number of parts in the odd partitions of n and the number of parts in the distinct partitions of n (offset 0). For example, if n = 5, there are 9 parts in the odd partitions of 5 (5, 311, 11111) and 5 parts in the distinct partitions of 5 (5, 41, 32), with difference 4. - George Beck, Apr 22 2017
George E. Andrews has kindly informed me that he has proved this conjecture and the result will be included in his article "Euler's Partition Identity and Two Problems of George Beck" which will appear in The Mathematics Student, 86, Nos. 1-2, January - June (2017). - George Beck, Apr 23 2017
a(n) is the number of partitions of n with exactly one repeated part. - Andrew Howroyd, Feb 14 2021

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n=0, t, `if`(i<1, 0,
          b(n, i-1, t)+`if`(i>n or t=1 and i::even, 0,
          add(b(n-i*j, i-1, `if`(i::even, 1, t)), j=1..n/i))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..70);  # Alois P. Heinz, Jun 17 2016
    A090867 := proc(n)
        add(numtheory[tau](k)*A000009(n-2*k),k=1..n/2) ;
    end proc: # R. J. Mathar, Jun 18 2016
  • Mathematica
    f[n_] := Count[ Plus @@@ Mod[ Union /@ IntegerPartitions[n] + 1, 2], 1]; Table[ f[n], {n, 0, 50}] (* Robert G. Wilson v, Feb 16 2004 *)
    a[n_] := Sum[DivisorSigma[0, k] PartitionsQ[n-2k], {k, 1, n/2}];
    a /@ Range[0, 70] (* Jean-François Alcover, May 24 2021, after R. J. Mathar *)
  • PARI
    seq(n)={Vec(sum(k=1, n\2, x^(2*k)/(1-x^(2*k)) + O(x*x^n))/prod(k=1, n\2, 1-x^(2*k-1) + O(x*x^n)), -(n+1))} \\ Andrew Howroyd, Feb 13 2021

Formula

G.f.: Sum_{m>0} x^(2*m)/(1-x^(2*m))/Product_{m>0} (1-x^(2*m-1)).
a(n) ~ 3^(1/4) * (2*gamma + log(3*n/Pi^2)) * exp(Pi*sqrt(n/3)) / (8*Pi*n^(1/4)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 25 2018
a(n) = A341494(n) + A341495(n) = A341496(n) + A341497(n). - Andrew Howroyd, Feb 14 2021

Extensions

More terms from Robert G. Wilson v, Feb 16 2004