OFFSET
0,1
REFERENCES
M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956, p. 11.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..300
M. S. Cheema and H. Gupta, Tables of Partitions of Gaussian Integers. National Institute of Sciences of India, Mathematical Tables, Vol. 1, New Delhi, 1956 (Annotated scanned pages from, plus a review)
FORMULA
From Vaclav Kotesovec, Jan 07 2017: (Start)
G.f.: (4 - x - 3*x^2 + x^3) / ((1-x)^3 * (1+x)) * Product_{k>=1} 1/(1-x^k).
a(n) ~ exp(Pi*sqrt(2*n/3)) * 3*sqrt(n)/(2*sqrt(2)*Pi^3).
(End)
MAPLE
with(numtheory):
b:= proc(n, k) option remember;
`if`(n>k, 0, 1) +`if`(isprime(n), 0,
add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
end:
a:= n-> b((45*2^n)$2):
seq(a(n), n=0..50); # Alois P. Heinz, May 26 2013
MATHEMATICA
b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, DeleteCases[Divisors[n], 1|n]}]]; a[n_] := b[45*2^n, 45*2^n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Mar 20 2014, after Alois P. Heinz *)
nmax = 100; CoefficientList[Series[(4 - x - 3*x^2 + x^3) / ((1 - x)^3 * (1 + x)) / Product[1 - x^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 07 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Extended beyond a(25) by Alois P. Heinz, May 26 2013
STATUS
approved