A022415 - OEIS

%I #15 Oct 15 2023 12:19:55

%S 3,10,13,15,18,21,23,26,28,31,34,36,39,42,44,47,49,52,55,57,60,62,65,

%T 68,70,73,76,78,81,83,86,89,91,94,97,99,102,104,107,110,112,115,117,

%U 120,123,125,128,131,133,136,138,141,144,146,149,151,154,157,159,162,165,167,170,172,175,178

%N Kim-sums: "Kimberling sums" K_n + K_4.

%D Posting to math-fun mailing list Jan 10 1997.

%H J. H. Conway, Allan Wechsler, Marc LeBrun, Dan Hoey, N. J. A. Sloane, <a href="/A269725/a269725.txt">On Kimberling Sums and Para-Fibonacci Sequences</a>, Correspondence and Postings to Math-Fun Mailing List, Nov 1996 to Jan 1997

%p Ki := proc(n,i)

%p     option remember;

%p     local phi ;

%p     phi := (1+sqrt(5))/2 ;

%p     if i= 0 then

%p         n;

%p     elif i=1 then

%p         floor((n+1)*phi) ;

%p     else

%p         procname(n,i-1)+procname(n,i-2) ;

%p     end if;

%p end proc:

%p Kisum := proc(n,m)

%p     local ks,a,i;

%p     ks := [seq( Ki(n,i)+Ki(m,i),i=0..5)] ;

%p     for i from 0 to 2 do

%p         for a from 0 do

%p             if Ki(a,0) = ks[i+1] and Ki(a,1) = ks[i+2] then

%p                 return a;

%p             end if;

%p             if Ki(a,0) > ks[i+1] then

%p                 break;

%p             end if;

%p         end do:

%p     end do:

%p end proc:

%p A022415 := proc(n)

%p     if n = 0 then

%p         3;

%p     else

%p         Kisum(n-1,3) ;

%p     end if;

%p end proc:

%p seq(A022415(n),n=0..80) ; # _R. J. Mathar_, Sep 03 2016

%t Ki[n_, i_] := Ki[n, i] = Which[i == 0, n, i == 1, Floor[(n + 1)* GoldenRatio], True, Ki[n, i - 1] + Ki[n, i - 2]];

%t Kisum[n_, m_] := Module[{ks, a, i}, ks = Table[Ki[n, i] + Ki[m, i], {i, 0, 5}]; For[i = 0, i <= 2, i++, For[a = 0, True, a++, If[Ki[a, 0] == ks[[i + 1]] && Ki[a, 1] == ks[[i + 2]], Return@a]; If[Ki[a, 0] > ks[[i + 1]], Break[]]]]];

%t a[n_] := If[n == 0, 3, Kisum[n - 1, 3]];

%t Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, Oct 15 2023, after _R. J. Mathar_ *)

%Y The "Kim-sums" K_n + K_i for i = 2 through 12 are given in A022413, A022414, A022415, A022416, ..., A022423.

%K nonn

%O 0,1

%A _Marc LeBrun_