|
|
A075864
|
|
G.f. satisfies A(x) = 1 + Sum_{n>=0} (x*A(x))^(2^n).
|
|
3
|
|
|
1, 1, 2, 4, 10, 26, 72, 204, 594, 1762, 5318, 16270, 50360, 157392, 496016, 1574432, 5028962, 16152194, 52133154, 169004450, 550036778, 1796512970, 5886709502, 19346204982, 63751851400, 210605429496, 697337388556, 2313871053172, 7692939444640, 25623793107344
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Number of Dyck n-paths with all ascent lengths being a power of 2. [David Scambler, May 07 2012]
|
|
LINKS
|
|
|
FORMULA
|
G.f. A(x) satisfies x*A(x) = series_reversion( x / ( 1 + Sum_{k>=0} x^(2^k) ) ). - Joerg Arndt, Apr 01 2019
|
|
MAPLE
|
b:= proc(x, y, t) option remember; `if`(x<0 or y>x, 0,
`if`(x=0, 1, b(x-1, y+1, true)+`if`(t, add(
b(x-2^j, y-2^j, false), j=0..ilog2(y)), 0)))
end:
a:= n-> b(2*n, 0, true):
|
|
MATHEMATICA
|
seq = {};
f[x_, y_, d_] :=
f[x, y, d] =
If[x < 0 || y < x , 0,
If[x == 0 && y == 0, 1,
f[x - 1, y, 0] + f[x, y - If[d == 0, 1, d], If[d == 0, 1, 2*d]]]];
For[n = 0, n <= 27, n++, seq = Append[seq, f[n, n, 0]]]; seq
A[_] = 0; m = 32;
Do[A[x_] = 1+Sum[(x A[x])^(2^n)+O[x]^m, {n, 0, Log[2, m]//Ceiling}], {m}];
|
|
PROG
|
(PARI) N=66; K=ceil(log(N)/log(2))+1; x='x+O('x^N); Vec(serreverse(x/(1 + sum(k=0, K, x^(2^k) ) ) ) ) - Joerg Arndt, Apr 01 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|