%I M3299 #57 Aug 26 2022 10:26:36
%S 1,4,7,9,12,15,17,20,22,25,28,30,33,36,38,41,43,46,49,51,54,56,59,62,
%T 64,67,70,72,75,77,80,83,85,88,91,93,96,98,101,104,106,109,111,114,
%U 117,119,122,125,127,130,132,135,138,140,143,145,148,151,153,156,159,161,164,166
%N a(n) = 1 + ceiling((n-1)*phi^2), phi = (1+sqrt(5))/2.
%C First column of dual Wythoff array.
%C Positions of 0's in A189479.
%C Skala (2016) asks if this sequence also gives the positions of the 0's in A283310. - _N. J. A. Sloane_, Mar 06 2017
%C Upper Wythoff sequence plus 2, when shifted by 1. - _Michel Dekking_, Aug 26 2019
%D Clark Kimberling, "Stolarsky interspersions," Ars Combinatoria 39 (1995) 129-138.
%D D. R. Morrison, "A Stolarsky array of Wythoff pairs," in A Collection of Manuscripts Related to the Fibonacci Sequence. Fibonacci Assoc., Santa Clara, CA, 1980, pp. 134-136.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Reinhard Zumkeller, <a href="/A007066/b007066.txt">Table of n, a(n) for n = 1..10000</a>
%H Clark Kimberling, <a href="http://faculty.evansville.edu/ck6/integer/intersp.html">Interspersions</a>
%H Matthew Skala, <a href="https://arxiv.org/abs/1604.04072">Graph Nimors</a>, arXiv preprint arXiv:1604.04072 [math.CO], 2016.
%H N. J. A. Sloane, <a href="/classic.html#WYTH">Classic Sequences</a>
%F a(n) = floor(1+phi*floor(phi*(n-1)+1)), phi=(1+sqrt(5))/2, n >= 2.
%F a(1)=1; for n>1, a(n)=a(n-1)+2 if n is already in the sequence, a(n)=a(n-1)+3 otherwise. - _Benoit Cloitre_, Mar 06 2003
%F a(n+1) = floor(n*phi^2) + 2, n>=1. - _Michel Dekking_, Aug 26 2019
%p Digits := 100: t := (1+sqrt(5))/2; A007066 := proc(n) if n <= 1 then 1 else floor(1+t*floor(t*(n-1)+1)); fi; end;
%t t = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {1, 0, 1}}] &, {0}, 6] (*A189479*)
%t Flatten[Position[t, 0]] (*A007066*)
%t Flatten[Position[t, 1]] (*A099267*)
%t With[{grs=GoldenRatio^2},Table[1+Ceiling[grs(n-1)],{n,70}]] (* _Harvey P. Dale_, Jun 24 2011 *)
%o (Haskell)
%o a007066 n = a007066_list !! (n-1)
%o a007066_list = 1 : f 2 [1] where
%o f x zs@(z:_) = y : f (x + 1) (y : zs) where
%o y = if x `elem` zs then z + 2 else z + 3
%o -- _Reinhard Zumkeller_, Sep 26 2014, Sep 18 2011
%o (Python)
%o from math import isqrt
%o def A007066(n): return (n+1+isqrt(5*(n-1)**2)>>1)+n if n > 1 else 1 # _Chai Wah Wu_, Aug 25 2022
%Y Cf. A064437. Apart from initial terms, same as A026356. Complement is (essentially) A026355. Equals 1 + A004957, also n + A004956.
%Y First differences give A076662.
%Y Complement of A099267. [_Gerald Hillier_, Dec 19 2008]
%Y Cf. A193214 (primes). Except for the first term equal to A001950 + 2.
%K nonn,easy,nice
%O 1,2
%A _N. J. A. Sloane_, _Mira Bernstein_
|