|
%I #11 Sep 06 2020 06:41:19
%S 1,2,3,4,6,8,11,14,19,25,33,41,53,66,83,102,128,156,193,233,285,343,
%T 416,495,597,710,849,1003,1194,1404,1662,1946,2291,2675,3137,3646,
%U 4260,4939,5744,6637,7697,8868,10250,11778,13570,15558,17877,20437,23423,26727,30550,34781,39669,45068,51287,58157
%N Expansion of (1/(1 - x))*Product_{k>=1} (1 - x^prime(k))/(1 - x^k).
%C Partial sums of A002095.
%C Number of partitions of n into nonprime parts if there are two kinds of 1's.
%H David A. Corneth, <a href="/A303663/b303663.txt">Table of n, a(n) for n = 0..9999</a>
%H <a href="/index/Par#part">Index entries for sequences related to partitions</a>
%p b:= proc(n, i) option remember; `if`(n=0 or i=1, n+1,
%p b(n, i-1)+`if`(isprime(i), 0, b(n-i, min(n-i, i))))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..60); # _Alois P. Heinz_, Apr 28 2018
%t nmax = 55; CoefficientList[Series[1/(1 - x) Product[(1 - x^Prime[k])/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x]
%Y Cf. A000070, A002095, A018252, A023895, A034891.
%K nonn,easy
%O 0,2
%A _Ilya Gutkovskiy_, Apr 28 2018
|
