|
|
MAPLE
| f := proc(n, k) option remember; if k > n then RETURN(0); fi; if k= 0 then if n=0 then RETURN(1) else RETURN(0); fi; fi; if k = 1 then RETURN(1); fi; if n mod 2 = 1 then RETURN(f(n-1, k)); fi; f(n-1, k)+f(n/2, k-1); end; # present sequence is f(2m, 6)
GFF := k->x^(2^(k-2))/((1-x)*mul((1-x^(2^j)), j=0..k-2)); # present g.f. is GFF(6)/x^16
a:= proc(n) local m, r; m:= iquo (n, 16, 'r'); r:= r+1; [1, 2, 4, 6, 10, 14, 20, 26, 36, 46, 60, 74, 94, 114, 140, 166][r] +(((((128/5*m +8*(15+r))*m +(228 +[0, 32, 68, 104, 144, 184, 228, 272, 320, 368, 420, 472, 528, 584, 644, 704][r]))*m +(172 +[0, 43, 98, 153, 223, 293, 378, 463, 566, 669, 790, 911, 1053, 1195, 1358, 1521][r]))*m +(247/5 +[0, 22, 55, 88, 138, 188, 255, 322, 415, 508, 627, 746, 900, 1054, 1243, 1432][r]))*m)/3 end: seq (a(n), n=0..60); # Alois P. Heinz, Apr 17 2009
|
