login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A133580 a(0)=a(1)=1; for n>1, a(n) = 2*a(n-1) + 1 if a(n-1) and n are coprime, otherwise a(n) = a(n-1)/gcd(a(n-1),n). 4
1, 1, 3, 1, 3, 7, 15, 31, 63, 7, 15, 31, 63, 127, 255, 17, 35, 71, 143, 287, 575, 1151, 2303, 4607, 9215, 1843, 3687, 1229, 2459, 4919, 9839, 19679, 39359, 78719, 157439, 314879, 629759, 1259519, 2519039, 5038079, 10076159, 20152319, 40304639 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

The initial value a(0)=1 is somehow artificial; using a(0)=0 would yield the same subsequent terms using the recurrence formula already for n=1. - M. F. Hasler, Feb 15 2015

LINKS

Harvey P. Dale, Table of n, a(n) for n = 0..1000

EXAMPLE

Write the GCD of a(n-1) and n under a(n-1):

n = : 0 1 2 3 4 5  6  7  8 9 ...

a(n): 1 1 3 1 3 7 15 31 63 7 ...

gcd : 1 1 3 1 1 1  1  1  9 1 ...

MATHEMATICA

a = {1, 1}; Do[If[GCD[n, a[[ -1]]] == 1, b = 2*a[[ -1]] + 1, b = a[[ -1]]/GCD[a[[ -1]], n]]; AppendTo[a, b], {n, 2, 50}]; a (* Stefan Steinerberger, Dec 31 2007 *)

nxt[{a_, b_}]:={a+1, If[CoprimeQ[b, a+1], 2b+1, b/GCD[b, a+1]]}; Join[{1}, Transpose[ NestList[nxt, {1, 1}, 50]][[2]]] (* Harvey P. Dale, Sep 16 2012 *)

PROG

(PARI) A=vector(1000, i, 1); for(n=2, #A, A[n]=if(gcd(A[n-1], n)>1, A[n-1]/gcd(A[n-1], n), A[n-1]*2+1)) \\ M. F. Hasler, Feb 15 2015

(PARI) a=0; #A133580=vector(1000, n, a=if(gcd(a, n)>1, a/gcd(a, n), a*2+1)) \\ M. F. Hasler, Feb 15 2015

CROSSREFS

Cf. A133058, A133579.

Sequence in context: A095732 A001644 A139123 * A019603 A171843 A132476

Adjacent sequences:  A133577 A133578 A133579 * A133581 A133582 A133583

KEYWORD

nonn

AUTHOR

Ctibor O. Zizka, Dec 26 2007

EXTENSIONS

Corrected and extended by Stefan Steinerberger, Dec 31 2007

Offset changed to 0 by N. J. A. Sloane, Feb 13 2015

The Mathematica programs are correct; b-file corrected by Harvey P. Dale, Feb 14 2015

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 April 16 05:26 EDT 2021. Contains 343030 sequences. (Running on oeis4.)