|
| |
|
|
A088359
|
|
Numbers which occur only once in A004001.
|
|
6
| |
|
|
3, 5, 6, 9, 10, 11, 13, 17, 18, 19, 20, 22, 23, 25, 28, 33, 34, 35, 36, 37, 39, 40, 41, 43, 44, 46, 49, 50, 52, 55, 59, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 77, 78, 79, 81, 82, 84, 87, 88, 89, 91, 92, 94, 97, 98, 100, 103, 107, 108, 110, 113, 117, 122, 129, 130, 131, 132
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,1
|
|
|
COMMENTS
| Out of the first one million terms (a(10^6)= 510403), 258661 occur only once.
Complement of A087686; A051135(a(n)) = 1. [Reinhard Zumkeller, Jun 03 2011]
|
|
|
LINKS
| Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
|
|
|
MATHEMATICA
| 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 &]
|
|
|
PROG
| (Haskell)
import Data.List (elemIndices)
a088359 n = a088359_list !! (n-1)
a088359_list = map succ $ elemIndices 1 a051135_list
-- Reinhard Zumkeller, Jun 03 2011
|
|
|
CROSSREFS
| Cf. A004001.
Sequence in context: A056875 A087757 A188163 * A184413 A187345 A179185
Adjacent sequences: A088356 A088357 A088358 * A088360 A088361 A088362
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Robert G. Wilson v (rgwv(AT)rgwv.com), Sep 26 2003
|
| |
|
|
