login
Numbers which occur only once in A004001.
28

%I #29 Oct 12 2018 21:59:02

%S 3,5,6,9,10,11,13,17,18,19,20,22,23,25,28,33,34,35,36,37,39,40,41,43,

%T 44,46,49,50,52,55,59,65,66,67,68,69,70,72,73,74,75,77,78,79,81,82,84,

%U 87,88,89,91,92,94,97,98,100,103,107,108,110,113,117,122,129,130,131,132

%N Numbers which occur only once in A004001.

%C Out of the first one million terms (a(10^6) = 510403), 258661 occur only once.

%C Complement of A087686; A051135(a(n)) = 1. - _Reinhard Zumkeller_, Jun 03 2011

%C From _Antti Karttunen_, Jan 18 2016: (Start)

%C In general, out of the first 2^(n+1) terms of A004001, 2^(n-1) - 1 terms (a quarter) occur only once. See also illustration in A265332.

%C One more than the positions of ones in A093879.

%C (End)

%H Reinhard Zumkeller, <a href="/A088359/b088359.txt">Table of n, a(n) for n = 1..10000</a>

%F From _Antti Karttunen_, Jan 18 2016: (Start)

%F Other identities.

%F For all n >= 0, a(A000079(n)) = A000051(n+1), that is, a(2^n) = 2^(n+1) + 1.

%F For all n >= 1:

%F a(n) = A004001(A266399(n)).

%F (End)

%t a[1] = 1; a[2] = 1; a[n_] := a[n] = a[ a[n - 1]] + a[n - a[n - 1]]; hc = Table[ a[n], {n, 1, 261}]; RunLengthEncodeOne[x_List] := Length[ # ] == 1 & /@ Split[x]; r = RunLengthEncodeOne[hc]; Select[ Range[ Length[r]], r[[ # ]] == True &]

%o (Haskell)

%o import Data.List (elemIndices)

%o a088359 n = a088359_list !! (n-1)

%o a088359_list = map succ $ elemIndices 1 a051135_list

%o -- _Reinhard Zumkeller_, Jun 03 2011

%o (Scheme, with Antti Karttunen's IntSeq-library)

%o (define A088359 (ZERO-POS 1 1 (COMPOSE -1+ A051135)))

%o ;; _Antti Karttunen_, Jan 18 2016

%Y Positions of ones in A051135.

%Y Cf. A188163 (same sequence with prepended 1).

%Y Cf. A087686 (complement).

%Y Cf. A004001, A093879, A265332, A266399.

%Y Cf. also A267110, A267111, A267112.

%K nonn

%O 1,1

%A _Robert G. Wilson v_, Sep 26 2003