login
Shifts 7 places left under Euler transform with a(0)=0 and a(n)=1 for n<7.
3

%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