|
%I #14 May 08 2019 16:47:47
%S 0,1,1,1,1,1,1,1,1,2,3,5,7,11,15,22,31,45,64,94,136,200,294,435,643,
%T 956,1420,2117,3157,4721,7064,10597,15909,23933,36038,54356,82059,
%U 124056,187707,284351,431114,654288,993780,1510785,2298471,3499653,5332313,8130576
%N Shifts 7 places left under Euler transform with a(0)=0 and a(n)=1 for n<7.
%H Alois P. Heinz, <a href="/A218024/b218024.txt">Table of n, a(n) for n = 0..1000</a>
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%F G.f.: x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 / Product_{n>=1} (1 - x^n)^a(n). - _Ilya Gutkovskiy_, May 08 2019
%p with(numtheory):
%p b:= proc(n) option remember; `if`(n=0, 1,
%p (add(add(d*a(d), d=divisors(j))*b(n-j), j=1..n))/n)
%p end:
%p a:= n-> `if`(n<7, signum(n), b(n-7)):
%p seq(a(n), n=0..50);
%Y Column k=7 of A144018.
%Y Cf. A316079.
%K nonn,eigen
%O 0,10
%A _Alois P. Heinz_, Oct 18 2012
|
