login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A328325 Expansion of Product_{k>=0} 1/(1 - x^(k^k)). 2

%I #35 Sep 14 2022 07:26:24

%S 1,2,3,4,6,8,10,12,15,18,21,24,28,32,36,40,45,50,55,60,66,72,78,84,91,

%T 98,105,113,122,131,140,150,161,172,183,195,208,221,234,248,263,278,

%U 293,309,326,343,360,378,397,416,435,455,476,497,519,542,566,590,615,641,668,695

%N Expansion of Product_{k>=0} 1/(1 - x^(k^k)).

%C Partial sums of A328301.

%H Alois P. Heinz, <a href="/A328325/b328325.txt">Table of n, a(n) for n = 0..10000</a>

%F a(n) = Sum_{k=0..n} A328301(k).

%F G.f.: 1/(1-x) + Sum_{n>0} x^(n^n) / Product_{k=0..n} (1 - x^(k^k)).

%e G.f.: 1/(1-x) + x/(1-x)^2 + x^4/((1-x)^2*(1-x^4)) + x^27/((1-x)^2*(1-x^4)*(1-x^27)) + ... .

%p b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,

%p b(n, i-1)+(p-> `if`(p>n, 0, b(n-p, i)))(i^i))

%p end:

%p a:= proc(n) option remember; `if`(n<2, n+1, a(n-1)+

%p b(n, floor((t-> t/LambertW(t))(log(n)))))

%p end:

%p seq(a(n), n=0..100); # _Alois P. Heinz_, Oct 12 2019

%t b[n_, i_] := b[n, i] = If[n == 0 || i == 1, 1, b[n, i-1] + With[{p = i^i}, If[p > n, 0, b[n-p, i]]]];

%t a[n_] := a[n] = If[n < 2, n+1, a[n-1] + b[n, Floor[PowerExpand[Log[n]/ ProductLog[Log[n]]]]]];

%t Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Sep 14 2022, after _Alois P. Heinz_ *)

%o (PARI) N=99; x='x+O('x^N); m=1; while(N>=m^m, m++); Vec(1/prod(k=0, m-1, 1-x^k^k))

%Y Cf. A000312, A064985, A086858, A328301.

%K nonn

%O 0,2

%A _Seiichi Manyama_, Oct 12 2019

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)