OFFSET
1,2
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
For n > 1: a(n+1) = a(n) + a(n - A249040(n)) + a(n - A249041(n)) by mutual recursion. - Reinhard Zumkeller, Nov 11 2014
MAPLE
M:=100;
v[1]:=1; v[2]:=2; w[1]:=0; w[2]:=1; x[1]:=1; x[2]:=1;
for n from 3 to M do
v[n]:=v[n-1]+v[n-1-w[n-1]]+v[n-1-x[n-1]];
if v[n] mod 2 = 0 then w[n]:=w[n-1]+1; x[n]:=x[n-1];
else w[n]:=w[n-1]; x[n]:=x[n-1]+1; fi;
od:
[seq(v[n], n=1..M)]; # A249039
[seq(w[n], n=1..M)]; # A249040
[seq(x[n], n=1..M)]; # A249041
PROG
(Haskell)
import Data.List (genericIndex)
a249039 n = genericIndex a249039_list (n - 1)
a249039_list = 1 : 2 : f 2 2 1 1 where
f x u v w = y : f (x + 1) y (v + 1 - mod y 2) (w + mod y 2)
where y = u + a249039 (x - v) + a249039 (x - w)
-- Reinhard Zumkeller, Nov 11 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Oct 26 2014
STATUS
approved