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.

A112021 Number of partitions of n into Chen primes.

This page as a plain text file.
%I A112021 #16 Dec 28 2017 21:30:34
%S A112021 0,1,1,1,2,2,3,3,4,5,6,7,9,10,12,14,17,19,23,26,30,35,40,46,52,60,67,
%T A112021 77,87,98,111,124,140,157,175,197,219,244,272,302,336,372,412,456,503,
%U A112021 556,613,675,742,816,896,983,1078,1180,1291,1411,1542,1683,1836,2001,2178
%N A112021 Number of partitions of n into Chen primes.
%C A112021 a(n) = A000607(n) for n <= 42.
%H A112021 Andrew Howroyd, <a href="/A112021/b112021.txt">Table of n, a(n) for n = 1..1000</a>
%F A112021 G.f.: Product_{k>=1} 1/(1 - x^A109611(k)). - _Andrew Howroyd_, Dec 28 2017
%t A112021 fQ[n_] := PrimeQ@n && (PrimeQ[n + 2] || 2 == Plus @@ Last /@ FactorInteger[n + 2]); f[n_] := Block[{c = k = 0, l = PartitionsP@n, p = Union /@ IntegerPartitions@n}, While[k++; k < l, If[Union[fQ /@ p[[k]]] == {True}, c++ ]]; c]; lst = {}; Do[ AppendTo[lst, f[n]], {n, 61}]; lst (* or *)
%t A112021 Rest@ CoefficientList[ Series[1/Times @@ (1 - x^Select[ Range@100, fQ@# &]), {x, 0, 61}], x] (* _Robert G. Wilson v_, Jun 16 2006 *)
%o A112021 (PARI)
%o A112021 ok(n)={isprime(n) && bigomega(n+2)<3}
%o A112021 {my(n=80); Vec(prod(k=1, n, if(ok(k), 1/(1-x^k) + O(x*x^n), 1))-1,-n)} \\ _Andrew Howroyd_, Dec 28 2017
%Y A112021 Cf. A112022, A101048, A109611.
%K A112021 nonn
%O A112021 1,5
%A A112021 _Reinhard Zumkeller_, Aug 26 2005