|
|
A282971
|
|
Number of compositions (ordered partitions) of n into primes of form x^2 + y^2 (A002313).
|
|
1
|
|
|
1, 0, 1, 0, 1, 1, 1, 2, 1, 3, 2, 4, 4, 6, 7, 9, 11, 15, 18, 24, 29, 37, 48, 58, 78, 92, 124, 149, 195, 243, 308, 393, 490, 629, 786, 1004, 1263, 1603, 2024, 2564, 3239, 4106, 5184, 6571, 8301, 10508, 13298, 16807, 21296, 26895, 34082, 43060, 54528, 68952, 87245, 110392, 139622, 176696, 223484, 282798, 357731
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,8
|
|
COMMENTS
|
Number of compositions (ordered partitions) of n into primes congruent to 1 or 2 mod 4.
Conjecture: every number > 16 is the sum of at most 4 primes of form x^2 + y^2.
|
|
LINKS
|
|
|
FORMULA
|
G.f.: 1/(1 - Sum_{k>=1} x^A002313(k)).
|
|
EXAMPLE
|
a(12) = 4 because we have [5, 5, 2], [5, 2, 5], [2, 2, 5] and [2, 2, 2, 2, 2, 2].
|
|
MATHEMATICA
|
nmax = 60; CoefficientList[Series[1/(1 - Sum[Boole[SquaresR[2, k] != 0 && PrimeQ[k]] x^k, {k, 1, nmax}]), {x, 0, nmax}], x]
|
|
PROG
|
(PARI) Vec(1/(1 - sum(k=1, 60, (isprime(k) && k%4<3)*x^k)) + O(x^61)) \\ Indranil Ghosh, Mar 15 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|