OFFSET
0,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
a(n) ~ exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*n). - Vaclav Kotesovec, May 24 2018
MAPLE
with(numtheory):
b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
b(n, i-1) +`if`(i>n, 0, b(n-i, i)))
end:
a:= n-> b(n, phi(n)):
seq(a(n), n=0..100); # Alois P. Heinz, May 11 2015
MATHEMATICA
(* Requires version 6.0+ *) Table[Length[IntegerPartitions[n, n, Range[EulerPhi[n]]]], {n, 0, 47}] (* Ivan Neretin, May 11 2015 *)
intPartLen[n_, i_] := intPartLen[n, i] = If[n == 0 || i == 1, 1, intPartLen[n, i - 1] + If[i > n, 0, intPartLen[n - i, i]]]; intPartLenPhi[n_] := intPartLen[n, EulerPhi[n]]; Table[intPartLenPhi[n], {n, 0, 99}] (* Jean-François Alcover, Nov 11 2015, after Alois P. Heinz *)
PROG
(Haskell)
a227296 n = p [1 .. a000010 n] n where
p _ 0 = 1
p [] _ = 0
p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 05 2013
STATUS
approved