login
This site is supported by donations to The OEIS Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A117120 a(1)=1. a(n) is smallest positive integer not occurring earlier in the sequence where a(n) is congruent to -1 (mod a(n-1)). 11
1, 2, 3, 5, 4, 7, 6, 11, 10, 9, 8, 15, 14, 13, 12, 23, 22, 21, 20, 19, 18, 17, 16, 31, 30, 29, 28, 27, 26, 25, 24, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 63, 62, 61, 60, 59, 58, 57, 56, 55, 54, 53, 52, 51, 50, 49, 48, 95, 94, 93, 92, 91, 90, 89, 88, 87 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Sequence is a permutation of the positive integers.

The permutation is self-inverse. Except for fixed points 1, 2, 3 it consists completely of 2-cycles: (4,5), (6,7), (8,11), (9,10), (12,15), (13,14), (16,23), (17,22), ..., (24,31), ..., (32,47), ... . - Klaus Brockhaus

The permutation transforms enumeration system of positive irreducible fractions A071766/A229742 (HCS) into enumeration system A245325/A245326, and viceversa. Yosu Yurramendi, Jun 09 2015

A092569(a(n)) = a(A092569(n)), n > 0.

A258746(a(n)) = a(A258746(n)), n > 0.

A258996(a(n)) = a(A258996(n)), n > 0.

A054429(a(n)) = a(A054429(n)), n > 0.

a(n) = A054429(A063946(n)) = A063946(A054429(n)), n > 0.- Yosu Yurramendi, Mar 23 2017

LINKS

Robert Israel, Table of n, a(n) for n = 1..10000

FORMULA

For n >= 2: If a(n-1) = 2^m, m=positive integer, then a(n)= 2^(m+1)-1. If a(n-1) = 3*2^m, m= nonnegative integer, then a(n) = 3*2^(m+1)-1. Otherwise, a(n) = a(n-1) -1.

For n >= 2: a(2*n) = 2*a(n)+1, a(2*n+1) = 2*a(n). - Yosu Yurramendi, Jun 08 2015

MAPLE

A[1]:= 1: A[2]:= 2: B[1]:= 0: B[2]:= 0:

for n from 3 to 100 do

  for m from A[n-1]-1 by A[n-1] while assigned(B[m]) do od:

  A[n]:= m;

  B[m]:= 0;

od:

seq(A[n], n=1..100); # Robert Israel, Jun 09 2015

MATHEMATICA

f[n_] := Block[{a = {1}, i, k}, Do[k = 1; While[Or[Mod[k, a[[i - 1]]] != a[[i - 1]] - 1, MemberQ[a, k]], k++]; AppendTo[a, k], {i, 2, n}]; a]; f@ 120 (* Michael De Vlieger, Jun 11 2015 *)

A[n_]:= If[n<4, n, If[EvenQ[n], 2A[n/2] + 1, 2A[(n - 1)/2]]]; Table[A[n], {n, 100}] (* Indranil Ghosh, Mar 21 2017 *)

f[lst_List] := Block[{k = 2, m = lst[[-1]]}, While[ MemberQ[lst, k] || 1 + Mod[k, m] != m, k++]; Append[lst, k]]; Nest[f, {1}, 70] (* Robert G. Wilson v, Jan 22 2018 *)

PROG

(R)

a <- 1:3 # If it were c(1, 3, 2), it would be A054429

maxn <- 50 # by choice

#

for(n in 2:maxn){

  a[2*n  ] <- 2*a[n]+1

  a[2*n+1] <- 2*a[n]

}

#

a

# Yosu Yurramendi, Jun 08 2015

(PARI) A(n) = if(n<4, n, if(n%2, 2*A(n\2), 2*A(n/2)+1));

for(n=1, 50, print1(A(n), ", ")) \\ Indranil Ghosh, Mar 21 2017

CROSSREFS

Cf. A071766, A229742, A245325, A245326.

Sequence in context: A085790 A257339 A210882 * A181095 A276345 A257455

Adjacent sequences:  A117117 A117118 A117119 * A117121 A117122 A117123

KEYWORD

easy,nonn

AUTHOR

Leroy Quet, Apr 19 2006

EXTENSIONS

More terms from Klaus Brockhaus

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 16 15:53 EST 2019. Contains 319195 sequences. (Running on oeis4.)