OFFSET
0,1
REFERENCES
Posting to math-fun mailing list Jan 10 1997.
MAPLE
Ki := proc(n, i)
option remember;
local phi ;
phi := (1+sqrt(5))/2 ;
if i= 0 then
n;
elif i=1 then
floor((n+1)*phi) ;
else
procname(n, i-1)+procname(n, i-2) ;
end if;
end proc:
Kisum := proc(n, m)
local ks, a, i;
ks := [seq( Ki(n, i)+Ki(m, i), i=0..5)] ;
for i from 0 to 2 do
for a from 0 do
if Ki(a, 0) = ks[i+1] and Ki(a, 1) = ks[i+2] then
return a;
end if;
if Ki(a, 0) > ks[i+1] then
break;
end if;
end do:
end do:
end proc:
A022418 := proc(n)
if n = 0 then
6;
else
Kisum(n-1, 6) ;
end if;
end proc:
seq(A022418(n), n=0..80) ; # R. J. Mathar, Sep 03 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved