cp's OEIS Frontend

A145519 a(n) = Sum_{k=1..n} A145518(n,k).

%I A145519 #39 Jan 26 2025 14:29:28
%S A145519 1,2,7,19,54,134,354,838,2057,4794,11232,25412,58075,128670,286152,
%T A145519 625829,1365653,2941088,6331146,13474533,28642325,60404681,127082128,
%U A145519 265712673,554608226,1151374963,2385950536,4924685252,10145267212,20831428273,42708248451
%N A145519 a(n) = Sum_{k=1..n} A145518(n,k).
%C A145519 Row sums of A145518.
%C A145519 Also row sums of A129129, A215366.
%C A145519 a(n) = sum of the Heinz numbers of the partitions of n. The Heinz number of a partition p = [p_1, p_2, ..., p_r] is defined as Product(p_j-th prime, j=1...r) (concept used by _Alois P. Heinz_ in A215366 as an "encoding" of a partition). For example, for the 3 partitions of 3, namely [3], [1,2], and [1,1,1] we get 5, 2*3=6, and 2*2*2=8, respectively; their sum is a(3) = 19. - _Emeric Deutsch_, Jun 09 2015
%H A145519 Alois P. Heinz, <a href="/A145519/b145519.txt">Table of n, a(n) for n = 0..1000</a>
%H A145519 More terms in <a href="http://www.tilman-neumann.de/docs/A145518.txt">A145518 and A145519</a>
%F A145519 G.f.: 1/Product_{i>=1}(1-prime(i)*x^i). - _Vladeta Jovovic_, Nov 09 2008
%F A145519 a(n) ~ c * 2^n, where c = Product_{k>=2} 1/(1 - prime(k)/2^k) = 50.412394245500690832088704444961002125578414895935257436317... . - _Vaclav Kotesovec_, Sep 10 2014, updated Apr 11 2020
%p A145519 b:= proc(n, i) option remember; `if`(n=0 or i<2, 2^n,
%p A145519       add(b(n-i*j, i-1)*ithprime(i)^j, j=0..iquo(n, i)))
%p A145519     end:
%p A145519 a:= n-> b(n, n):
%p A145519 seq(a(n), n=0..35);  # _Alois P. Heinz_, Feb 19 2013
%t A145519 b[n_, i_] := b[n, i] = If[n == 0 || i < 2, 2^n, Sum[b[n-i*j, i-1]*Prime[i]^j, {j, 0, Quotient[n, i]}]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 35}] (* _Jean-François Alcover_, Feb 24 2015, after _Alois P. Heinz_ *)
%Y A145519 Cf. A129129, A145518, A215366, A334434.
%K A145519 nonn
%O A145519 0,2
%A A145519 _Tilman Neumann_, Oct 12 2008
%E A145519 a(0) inserted by _Alois P. Heinz_, Feb 19 2013