cp's OEIS Frontend

A022161 First row of spectral array W(e-1).

%I A022161 #12 Oct 24 2014 18:05:49
%S A022161 1,2,3,4,6,9,15,21,36,50,85,119,204,284,487,679,1166,1624,2790,3884,
%T A022161 6673,9291,15964,22226,38190,53169,91359,127191,218549,304267,522816,
%U A022161 727870,1250685,1741218,2991903,4165361,7157264,9964423,17121686
%N A022161 First row of spectral array W(e-1).
%H A022161 A. Fraenkel and C. Kimberling, <a href="http://dx.doi.org/10.1016/0012-365X(94)90259-3">Generalized Wythoff arrays, shuffles and interspersions</a>, Discrete Mathematics 126 (1994) 137-149.
%o A022161 (PARI)
%o A022161 \\ The first row of the generalized Wythoff array W(h),
%o A022161 \\   where h is an irrational number between 1 and 2.
%o A022161 row1(h, m) = {
%o A022161   my(
%o A022161     a=vector(m, n, floor(n*h)),
%o A022161     b=setminus(vector(m, n, n), a),
%o A022161     w=[a[1]^2, b[a[1]]],
%o A022161     j=3
%o A022161   );
%o A022161   while(1,
%o A022161     if(j%2==1,
%o A022161       if(w[j-1]<=#a, w=concat(w, a[w[j-1]]), return(w))
%o A022161     ,
%o A022161       if(w[j-2]<=#b, w=concat(w, b[w[j-2]]), return(w))
%o A022161     );
%o A022161     j++
%o A022161   );
%o A022161   w
%o A022161 }
%o A022161 allocatemem(10^9)
%o A022161 row1(exp(1)-1, 10^7) \\ _Colin Barker_, Oct 24 2014
%K A022161 nonn
%O A022161 0,2
%A A022161 _Clark Kimberling_
%E A022161 More terms from _Colin Barker_, Oct 24 2014