cp's OEIS Frontend

A022159 First row of spectral array W(sqrt(3)).

%I A022159 #15 Oct 24 2014 18:05:34
%S A022159 1,2,3,4,6,9,15,21,36,49,84,115,199,272,471,643,1113,1521,2634,3598,
%T A022159 6231,8512,14743,20139,34881,47649,82530,112738,195267,266740,462007,
%U A022159 631113,1093119,1493229,2586348,3533017,6119364,8359207,14478571
%N A022159 First row of spectral array W(sqrt(3)).
%H A022159 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 A022159 (PARI)
%o A022159 \\ The first row of the generalized Wythoff array W(h),
%o A022159 \\   where h is an irrational number between 1 and 2.
%o A022159 row1(h, m) = {
%o A022159   my(
%o A022159     a=vector(m, n, floor(n*h)),
%o A022159     b=setminus(vector(m, n, n), a),
%o A022159     w=[a[1]^2, b[a[1]]],
%o A022159     j=3
%o A022159   );
%o A022159   while(1,
%o A022159     if(j%2==1,
%o A022159       if(w[j-1]<=#a, w=concat(w, a[w[j-1]]), return(w))
%o A022159     ,
%o A022159       if(w[j-2]<=#b, w=concat(w, b[w[j-2]]), return(w))
%o A022159     );
%o A022159     j++
%o A022159   );
%o A022159   w
%o A022159 }
%o A022159 allocatemem(10^9)
%o A022159 row1(sqrt(3), 10^7) \\ _Colin Barker_, Oct 24 2014
%K A022159 nonn
%O A022159 0,2
%A A022159 _Clark Kimberling_
%E A022159 More terms from _Colin Barker_, Oct 24 2014