login
a(n) = floor(n*phi), where phi = (1 + sqrt(5))/2.
27

%I #56 Aug 02 2024 05:03:45

%S 0,1,3,4,6,8,9,11,12,14,16,17,19,21,22,24,25,27,29,30,32,33,35,37,38,

%T 40,42,43,45,46,48,50,51,53,55,56,58,59,61,63,64,66,67,69,71,72,74,76,

%U 77,79,80,82,84,85,87,88,90,92,93,95,97,98,100,101,103,105,106

%N a(n) = floor(n*phi), where phi = (1 + sqrt(5))/2.

%C a(n) is the smallest number different from a(i) and a(i)+i for i < n.

%C The losing positions in the game of Wythoff-Nim are precisely the pairs (a(n), a(n)+n).

%H G. C. Greubel, <a href="/A066096/b066096.txt">Table of n, a(n) for n = 0..10000</a>

%F a(n) = A000201(n).

%F Duplicate values in A060143.

%F a(n) = 1 + A022342(n) = A000201(n).

%F a(n) = floor(n*phi), where phi = (1 + sqrt(5))/2. - _Peter Munn_, Jan 12 2018

%F a(n) = A026351(n) - 1. - _Philippe Deléham_, Jan 15 2023

%t Floor[GoldenRatio*Range[0, 80]] (* _G. C. Greubel_, Sep 12 2023 *)

%o (PARI) a(n) = (n+sqrtint(5*n^2))\2;

%o [a(n)|n<-[0..100]] \\ _Simon Strandgaard_, Jun 28 2022

%o (Magma) [Floor((1+Sqrt(5))*n/2): n in [0..80]]; // _G. C. Greubel_, Sep 12 2023

%o (SageMath) [floor(golden_ratio*n) for n in range(81)] # _G. C. Greubel_, Sep 12 2023

%Y Essentially the partial sums of A001468.

%Y Cf. A000201, A001622, A003622, A022342, A026351, A035336, A060143.

%Y Cf. A004919, A004920, A004921, A004922, A004923, A004924, A004925.

%Y Cf. A004926, A004927, A004928, A004929, A004930, A004931, A004932.

%Y Cf. A004933, A004934, A004935, A004976, A090909.

%K nonn,easy

%O 0,3

%A Michele Dondi (bik.mido(AT)tiscalenet.it), Dec 30 2001

%E Name corrected by _Peter Munn_, Dec 06 2017

%E New name using a formula from _Peter Munn_ by _Peter Luschny_, Jan 18 2023