login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A281545 Expansion of Sum_{k>=2} x^prime(k)/(1 + x^prime(k)) * Product_{k>=2} (1 + x^prime(k)). 1
0, 0, 1, 0, 1, 0, 1, 2, 0, 2, 1, 2, 1, 2, 3, 4, 1, 4, 4, 4, 6, 4, 7, 6, 6, 8, 9, 8, 10, 6, 13, 12, 12, 14, 15, 16, 16, 18, 23, 22, 19, 24, 24, 30, 28, 30, 33, 34, 34, 40, 44, 46, 44, 46, 58, 56, 60, 64, 65, 68, 70, 80, 86, 88, 87, 94, 101, 112, 114, 116, 125, 130, 132, 148, 159, 162, 163, 168, 190, 196 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,8
COMMENTS
Total number of parts in all partitions of n into distinct odd primes.
LINKS
FORMULA
G.f.: Sum_{k>=2} x^prime(k)/(1 + x^prime(k)) * Product_{k>=2} (1 + x^prime(k)).
EXAMPLE
a(23) = 7 because we have [23], [13, 7, 3], [11, 7, 5] and 1 + 3 + 3 = 7.
MATHEMATICA
nmax = 80; Rest[CoefficientList[Series[Sum[x^Prime[k]/(1 + x^Prime[k]), {k, 2, nmax}] Product[1 + x^Prime[k], {k, 2, nmax}], {x, 0, nmax}], x]]
PROG
(PARI)
sumparts(n, pred)={sum(k=1, n, 1 - 1/(1+pred(k)*x^k) + O(x*x^n))*prod(k=1, n, 1+pred(k)*x^k + O(x*x^n))}
{my(n=60); Vec(sumparts(n, v->v>2 && isprime(v)), -n)} \\ Andrew Howroyd, Dec 28 2017
CROSSREFS
Sequence in context: A329697 A029367 A192541 * A332033 A331981 A099302
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 23 2017
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 10:31 EDT 2024. Contains 371240 sequences. (Running on oeis4.)