|
%I #12 Sep 16 2016 02:40:09
%S 0,1,1,1,2,1,4,2,6,4,9,12,13,48,17,192,21,768,26,3840,32,23040,39,
%T 161280,46,1128960,54,9031680,62,72253440,70,578027520,78,4624220160,
%U 87,41617981440,97,416179814400,108,4577977958400,120,54935735500800,132
%N List of pairs (a(n),b(n)): f(n) = A004001(n); a(n) = f(n) + a(n-1); b(n) = f(n)*b(n-1).
%H G. C. Greubel, <a href="/A154280/b154280.txt">Table of n, a(n) for n = 0..1001</a>
%t (*A004001*) f[0] = 0; f[1] = 1; f[2] = 1; f[n_] := f[n] = f[f[n - 1]] + f[n - f[n - 1]];
%t a[0] = 0; a[n_] := a[n] = f[n] + a[n - 1];
%t b[0] = 1; b[1] = 1; b[n_] := b[n] = (f[n])*b[n - 1];
%t Flatten[Table[{a[n], b[n]}, {n, 0, 30}]]
%Y Cf. A004001.
%K nonn,tabf
%O 0,5
%A _Roger L. Bagula_, Jan 06 2009
|
