The On-Line Encyclopedia of Integer Sequences (OEIS)

The OEIS is supported by the many generous donors to the OEIS Foundation.

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give

A116533 revision #6

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.

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

Table of n, a(n) for n=1..57.

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

editing