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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A116533 a(1)=1, a(2)=2, for n > 2 if a(n-1) is prime, then a(n) = 2*a(n-1), otherwise a(n) = a(n-1) - 1. 13
 1, 2, 4, 3, 6, 5, 10, 9, 8, 7, 14, 13, 26, 25, 24, 23, 46, 45, 44, 43, 86, 85, 84, 83, 166, 165, 164, 163, 326, 325, 324, 323, 322, 321, 320, 319, 318, 317, 634, 633, 632, 631, 1262, 1261, 1260, 1259, 2518, 2517, 2516, 2515, 2514, 2513, 2512, 2511, 2510, 2509, 2508 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 COMMENTS For n >= 3, using Wilson's theorem, a(n) = a(n-1) + (-1)^r*gcd(a(n-1), W), where W = A038507(a(n-1) - 1), and r=1 if gcd(a(n-1), W) = 1 and r=0 otherwise. - Vladimir Shevelev, Aug 07 2009 LINKS MAPLE a[1]:=1: a[2]:=2: for n from 3 to 60 do if isprime(a[n-1])=true then a[n]:=2*a[n-1] else a[n]:=a[n-1]-1 fi od: seq(a[n], n=1..60); # Emeric Deutsch, Apr 02 2006 CROSSREFS Cf. A006992, A055496, A080359, A104272, A106108, A132199. - Vladimir Shevelev, Aug 07 2009 Sequence in context: A180625 A132340 A132666 * A087559 A193298 A168007 Adjacent sequences:  A116530 A116531 A116532 * A116534 A116535 A116536 KEYWORD nonn AUTHOR Rodolfo Kurchan, Mar 26 2006 EXTENSIONS More terms from Emeric Deutsch, Apr 02 2006 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.

Last modified September 26 05:00 EDT 2020. Contains 337346 sequences. (Running on oeis4.)