|
%I #12 Nov 16 2015 23:27:53
%S 1,1,0,1,2,1,0,2,2,1,2,2,3,4,3,2,6,6,1,5,9,6,5,9,9,9,11,8,13,17,11,12,
%T 22,19,13,23,25,22,26,28,30,37,34,31,47,45,35,50,61,52,56,69,68,76,74,
%U 72,95,100,82,100,130,112,113,139,144,149,154,156,183
%N Number of partitions of n into distinct positive Loeschian numbers (cf. A003136).
%H Reinhard Zumkeller, <a href="/A198727/b198727.txt">Table of n, a(n) for n = 0..1000</a>
%e a(20) = #{19+1, 16+4, 16+3+1, 13+7, 13+4+3, 12+7+1, 12+4+3+1, 9+7+4, 9+7+3+1} = 9;
%e a(21) = #{21, 16+4+1, 13+7+1, 13+4+3+1, 12+9, 9+7+4+1} = 6;
%e a(22) = #{21+1, 19+3, 13+9, 12+9+1, 12+7+3} = 5.
%o (Haskell)
%o import Data.MemoCombinators (memo2, list, integral)
%o a198727 n = a198727_list !! n
%o a198727_list = f 0 [] $ tail a003136_list where
%o f u vs ws'@(w:ws) | u < w = (p' vs u) : f (u + 1) vs ws'
%o | otherwise = f u (vs ++ [w]) ws
%o p' = memo2 (list integral) integral p
%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_, Nov 16 2015, Oct 30 2011
%Y Cf. A003136, A198726.
%K nonn
%O 0,5
%A _Reinhard Zumkeller_, Oct 30 2011
|
