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

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A193298 Gica-Panaitopol recursion: a(1) = 1; a(n+1) = 2*a(n) if a(n) <= n; otherwise a(n+1) = a(n) - 1. 5

%I #14 Mar 30 2012 19:00:09

%S 1,2,4,3,6,5,10,9,8,16,15,14,13,26,25,24,23,22,21,20,40,39,38,37,36,

%T 35,34,33,32,31,30,60,59,58,57,56,55,54,53,52,51,50,49,48,47,46,92,91,

%U 90,89,88,87,86,85,84,83,82,81,80,79,78,77,76,75,74,73,72,71,70,69

%N Gica-Panaitopol recursion: a(1) = 1; a(n+1) = 2*a(n) if a(n) <= n; otherwise a(n+1) = a(n) - 1.

%C Using the Prime Number Theorem, Gica and Panaitopol show that the sequence contains infinitely many primes.

%D A. Gica and L. Panaitopol, An application of the prime element theorem, Gazeta Matematica 21(100), No. 2 (2003), 113-115.

%H T. D. Noe, <a href="/A193298/b193298.txt">Table of n, a(n) for n = 1..10000</a>

%e The terms occur in disjoint blocks of decreasing consecutive numbers: 1; 2; 4, 3; 6, 5; 10, 9, 8; 16, 15, 14, 13; 26, 25, 24, 23, 22, 21, 20; . . .

%t a[1] = 1; a[n_] := a[n] = If[a[n-1] <= n-1, 2*a[n-1], a[n-1]-1]; Table[a[n], {n, 100}]

%Y Cf. A193299 (sorted sequence), A193300 (subset of primes), A193301 (complement of sorted sequence).

%K nonn

%O 1,2

%A _Jonathan Sondow_, Jul 21 2011

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 18 04:31 EDT 2024. Contains 371767 sequences. (Running on oeis4.)