OFFSET
0,10
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..15625
FORMULA
a(5*n) = a(n),
a(5*n+1) = a(n+1) - a(n),
a(5*n+2) = a(n+2) - 2*a(n+1) + a(n),
a(5*n+3) = a(n+3) - 3*a(n+2) + 3*a(n+1) - a(n),
a(5*n+4) = a(n+4) - 4*a(n+3) + 6*a(n+2) - 4*a(n+1) + a(n).
MAPLE
a:= proc(n) option remember; local m, q;
m:= irem(n, 5, 'q'); `if`(n<5, `if`(n=4, 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, 5]; If[n < 5, If[n == 4, 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, from Maple *)
CROSSREFS
KEYWORD
AUTHOR
Alois P. Heinz, Sep 27 2013
STATUS
approved