OFFSET
0,6
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..19683
FORMULA
a(3*n) = a(n),
a(3*n+1) = a(n+1) - a(n),
a(3*n+2) = a(n+2) - 2*a(n+1) + a(n).
MAPLE
a:= proc(n) option remember; local m, q;
m:= irem(n, 3, 'q'); `if`(n<3, `if`(n=2, 1, 0),
add(a(q+m-j)*(-1)^j*binomial(m, j), j=0..m))
end:
seq(a(n), n=0..100);
MATHEMATICA
a[n_] := a[n] = Module[{m, q}, {q, m} = QuotientRemainder[n, 3]; If[n < 3, If[n == 2, 1, 0], Sum[a[q + m - j]*(-1)^j*Binomial[m, j], {j, 0, m}]]];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Jun 09 2018, translated from Maple *)
CROSSREFS
KEYWORD
AUTHOR
Alois P. Heinz, Sep 27 2013
STATUS
approved