%I #23 May 20 2023 12:11:30
%S 1,1,0,1,0,2,0,4,1,5,1,11,0,17,4,13,13,37,2,53,13,51,35,103,10,135,78,
%T 167,89,255,4,339,253,378,306,542,121,759,558,872,498,1259,121,1609,
%U 1180,1677,1665,2589,808,3250,1969,3844,3325,5119,1850,6268,4758,7546,7070
%N Number of ways to partition n into distinct positive integers <= phi(n), where phi is Euler's totient function (A000010).
%D Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education, Vol. 31, No. 1, pp. 24-28, Winter 1997. MathEduc Database (Zentralblatt MATH, 1997c.01891).
%H Alois P. Heinz, <a href="/A079124/b079124.txt">Table of n, a(n) for n = 0..10000</a>
%H Mohammad K. Azarian, <a href="https://web.archive.org/web/20100727063142/http://www.math-cs.ucmo.edu/~mjms/2004.1/azar6.pdf">A Generalization of the Climbing Stairs Problem II</a>, Missouri Journal of Mathematical Sciences, Vol. 16, No. 1, Winter 2004, pp. 12-17. Zentralblatt MATH, Zbl 1071.05501.
%F a(n) = b(0, n), b(m, n) = 1 + sum(b(i, j): m<i<j<phi(n) & i+j=n).
%p with(numtheory):
%p b:= proc(n, i) option remember; `if`(n=0, 1,
%p `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i-1))))
%p end:
%p a:= n-> b(n, phi(n)):
%p seq(a(n), n=0..100); # _Alois P. Heinz_, May 11 2015
%t b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i-1]]]]; a[n_] := b[n, EulerPhi[n]]; Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Jun 30 2015, after _Alois P. Heinz_ *)
%o (Haskell)
%o a079124 n = p [1 .. a000010 n] n where
%o p _ 0 = 1
%o p [] _ = 0
%o p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
%o -- _Reinhard Zumkeller_, Jul 05 2013
%Y Cf. A079126, A000009, A079122, A079125, A067953.
%Y Cf. A227296, A036998.
%K nonn
%O 0,6
%A _Reinhard Zumkeller_, Dec 27 2002
%E a(0)=1 prepended by _Alois P. Heinz_, May 11 2015
|